sconstruct

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

from rsfproj import *

Fetch('ida.H','ida')

Flow('ida','ida.H','dd form=native | window n1=600')
Plot('ida','grey gpow=0.5 crowd=0.85 title="stack (gpow=0.5)" ')

Flow('npef nlag','ida','lopef a=5,1 w=81,11 k=12,30 lag=${TARGETS[1]}')
Flow('spef slag','ida',
     'lopef a=9,3 w=81,11 k=12,30 gap=0,1 lag=${TARGETS[1]}')

eps = 2

Flow('sneps','ida spef npef',
     '''
     losignoi sfilt=${SOURCES[1]} nfilt=${SOURCES[2]}
     niter=30 eps=%g
     ''' % eps)
Plot('sneps','grey gpow=0.5 crowd=0.85 title="epsilon=%g" ' % eps)

Result('sneps','ida sneps','SideBySideAniso')

End()