sconstruct

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

from rsfproj import *

def plot(title):
    return 'grey titlesz=15 title="%s" wantaxis=n' % title

Flow('make',None,'make n2=20 n=5 p=11 t1=3 t2=3')
Plot('make',plot('original'))

Flow('win','hole','window j2=2 f2=1')
Flow('strip','win','math output=0 | interleave axis=2 $SOURCE')
Plot('strip',plot('gapped'))

Flow(['pef','lag'],'strip',
     'mspef maskin=$SOURCE a=5,3 ns=2 jump=1,2 lag=${TARGETS[1]}')
Flow('gap',['strip','pef'],'msmiss filt=${SOURCES[1]}')
Plot('gap',plot('restored'))

Result('mshole','make strip gap','SideBySideAniso',vppen='txscale=2')

Flow('hole','make','hole')

outs = []
for jump in ([1],[1,4]):
    jumps = string.join(map(str,jump),',')
    pef = 'pef' + jumps
    lag = 'lag' + jumps
    Flow([pef,lag],'hole',
         'mspef maskin=$SOURCE a=5,3 ns=%d jump=%s lag=${TARGETS[1]}' %
         (len(jump),jumps))
    for niter in [1,30]:
        out = 'out%s-%d' % (jumps,niter)
        outs.append(out)
        Flow(out,['hole',pef],'msmiss filt=${SOURCES[1]} niter=%d' % niter)
        Plot(out,plot('%s, iter=%d' %
                      (('singlescale','multiscale')[len(jump)-1],niter))) 

Result('msiter',outs,'TwoRows')

###########################################################################
End()