R: Convolve spikes with wavelets

sigconv {RSEIS} | R Documentation

## Convolve spikes with wavelets

### Description

Convolve spikes with wavelets

### Usage

    sigconv(wigmat, wavepulse)

### Arguments

`wigmat` | matrix, spikes
`wavepulse` | wavelet for convolution

### Details

Convolution is done in Frequency domain on each trace

### Value

Matrix, waveforms

### Author(s)

Jonathan M. Lees<jonathan.lees@unc.edu>

### See Also

wiggleimage, symshot1, genrick

### Examples

    S1 <- symshot1()

    ##############   S1$THEORY$treflex

    d <- dim(S1$smograms)
    G1 <- matrix( rep(0, length=d[1]*d[2]), ncol=d[2], nrow=d[1])

    ############   set up the spike set for reflexions
    for(i in 1:3){
    p <- round( S1$THEORY$treflex[i,]/S1$dt );

    G1[cbind(p , 1:d[2])  ] <- 1

    }

    #### plot the spikes
    wiggleimage(0.1*G1, dt = -S1$dt, dx = S1$x, col = "black")

    ###  make a ricker wavelet
    wavelet <- genrick(25,S1$dt,35)
    klem <- 11
    ###
    nwave <- RPMG::RESCALE(wavelet, 0, 1, wavelet[1], max(wavelet))

    ##############  convolve the wavelet with the set of spikes
    H1 <- sigconv(G1, nwave)

    ############  plot
    wiggleimage(0.1*H1, dt = -S1$dt, dx = S1$x, col = "black")

[Package RSEIS version 3.7-4 Index]