sconstruct

国外免费地震资料处理软件包

from rsfproj import *

Flow('plane',None,
     '''
     spike n1=64 n2=32 d2=1 o2=1 label2=Trace
     nsp=3 k1=8,20,32 k2=4 l2=28 p2=2,1,0 |
     noise var=0.03
     ''')

for plot in ('wiggle','grey'):
    Result(plot,'plane',plot + '''
    clip=1 transp=y yreverse=y poly=y
    wanttitle=n wheretitle=b wherexlabel=t
    ''')

End()