seismic.convolve {RGeostats}  R Documentation 
Perform convolution of the seismic attributes
seismic.convolve(db, type = 0, ntw = 5, option = 0, tindex = 0, fpeak = 1, period = 1, amplitude = 1, distort = 0, flag.operate = TRUE, flag.contrast = FALSE, wavelet = NA, val.before=0., val.middle=0., val.after=0., radix = "Convolve", modify.target = db.locmod())
db 
The 
type 
Type of convolution

ntw 
Number of samples in the halflength of the wavelet excluding the center. If the wavelet is provided (type<0), this parameter is calculated automatically. 
option 
Egde processing option:

tindex 
Time index to locate the spike 
fpeak 
Peak frequency of the Ricker and AKB wavelets (f) 
period 
Wavelet period for Ricker wavelets (λ) 
amplitude 
Wavelet amplitude for Ricker & AKB wavelets (α) 
distort 
Wavelet distortion factor for Ricker wavelet (β) 
flag.operate 
When TRUE, the convolution is performed (the output is the modified Db). When FALSE, the discretized wavelet is printed out (no additional output is provided) 
flag.contrast 
When TRUE, the initial value is first converted into a contrast before convolution. 
wavelet 
Array containing the discretized wavelet used when type=1. Its dimension must be 2*ntw+1. 
val.before 
Assigned to undefined values of each trace before the first defined sample 
val.middle 
Assigned to undefined values of each trace between two defined samples 
val.after 
Assigned to undefined values of each trace after the last defined sample 
radix 
Radix of the name given to the variables storing the results in the target Db. 
modify.target 
Decides whether or not the newly created variables will have their
locator defined or not. For more information, see 
The Ricker1 wavelet equation is:
w(t) = ( 1  2 * (π*f*dt*(tt0))^2 ) / exp[ (π*f*dt*(tt0))^2 ]
The Ricker2 wavelet equation is:
w(t) = α * (1  2 * ( 2.5 * dt*(tt0)/λ)^2 ) / exp[ (2.5 * dt*(tt0)/λ)^2 ] * ( 1  2 * β * dt*(tt0)/λ )
The Alford, Kelly and Boore wavelet equation is:
w(t) = α * (tt0) * exp[ 2 * (f*dt*(tt0))^2]
with:
w(t) : Wavelet equation
dt : Grid mesh along the vertical axis
t0 : Center of the wavelet
The output Db with as many output variables as there are seismic attributes in the input Db.
seismic.grid.define
seismic.convert
seismic.operate