Wavelet Transform

#!/bin/sh

project=/hou/data/D_152_001/SU_Class/class_project

cd ${project}

cp int_datasets/template.dat int_datasets/hilb_wavelet_cwt.dat

bhpio delete=yes pathlist=int_datasets/hilb_wavelet_cwt.dat filename=hilb_wavelet_cwt

# make a wavelet to use in the transform program

bhpwavelet > seismic/wavelet.su

# make wavelet transform of every 200th line and trace

bhpread filename=hilb_seismic pathlist=int_datasets/hilb_seismic.dat \

 keys=ep,cdp keylist=7025-9263[200]:18858-21494[200] | \

bhpcwt first=0 inc=0.02 last=1.0 wavelet=seismic/wavelet.su | \

bhpsmooth npass=3 length=20 key=cdp nxmax=51 | \

bhpwritecube init=yes key2=cdp,18858,200,15 key1=ep,7025,200,13 key3=cdpt,1,1,51 \

 filename=hilb_wavelet_cwt pathlist=int_datasets/hilb_wavelet_cwt.dat

# read back on transform and display

bhpread filename=hilb_wavelet_cwt pathlist=int_datasets/hilb_wavelet_cwt.dat keys=ep,cdp \

 keylist=7025:19058 | suximage &