bitsmithy.net

So it turns out it absolutely is possible to approximate a cosine squared taper in Seismic Unix. The trick is to do a sine or cosine taper, square the whole thing, and then divide by the original.

Interestingly, I’d already tried this – it’s just that the particular view of the image I was using had a big anomaly in about the same region as the taper that I hadn’t seen using my default colour scales. I thought there were problems with it: I was wrong.

Very quick and dirty example:

< model1.su suwind key=tracl min=1 max=400 > model1.start.su
< model1.su suwind key=tracl min=401 max=500 > model1.end.su
< model2.su suwind key=tracl min=401 max=500 > model2.start.su
< model2.su suwind key=tracl min=501 max=901 > model2.end.su

< model1.end.su suflip flip=2 | sutxtaper taper=2 tbeg=0 tend=0 tr1=101 tr2=0 | suflip flip=2 | suop op=ssqr > model1.end.taper.su
< model2.start.su sutxtaper taper=2 tbeg=0 tend=0 tr1=101 tr2=0 | suop op=ssqr > model2.start.taper.su

suop2 model2.start.taper.su model2.start.su op=quo > inter1.su
suop2 model1.end.taper.su model1.end.su op=quo > inter2.su
susum inter1.su inter2.su > intermediate.su

cat model1.start.su intermediate.su model2.end.su | suximage