sconstruct

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

from rsfproj import *

Flow('aliasp',None,'aliasp cycles=25')
Plot('aliasp','window n1=575 | grey title=Input crowd=0.85')

Flow('dip2','aliasp',
     'twodip2 order=3 nj1=4 nj2=4 p0=2 q0=-1 eps=4 gauss=n')
Flow('dip1','aliasp dip2',
     'dip idip=${SOURCES[1]} order=3 nj1=4 rect1=7 rect2=3')
Flow('dip','dip1',
     'transp | spline n1=96 o1=0 d1=0.25 | transp')

Flow('aliasp4 ones4','aliasp','lpad jump=4 mask=${TARGETS[1]}')

for prec in [0,1]:
    dealias = 'dealias%d' % prec
    
    Flow(dealias,'aliasp4 dip ones4',
         '''
         planemis2 dip=${SOURCES[1]} mask=${SOURCES[2]} order=3 verb=y prec=%d
         ''' % prec)
    Plot(dealias,'window n1=575 | grey title=Output crowd=0.85')

    Result('aliasp%d' % prec,['aliasp',dealias],'SideBySideAniso')

End()