sconstruct

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

from rsfproj import *

dat = {'dat':'gomshot.h',
       'zof':'gomnearoffset.H'}

for case in dat.keys():
    Fetch(dat[case],'gom')
    Flow(case,dat[case],'dd form=native')

def showdat(title,max1,f2=1,flip='',wantaxis1=1):
    return '''
    window min1=1.8 max1=%g max2=10000 f2=%d | %s
    tpow tpow=1 | grey  crowd=0.85 pclip=99 
    label1=Time unit1=s label2=Offset unit2=ft
    title="%s" wantaxis1=%d
    ''' % (max1,f2,flip,title,wantaxis1)

Flow('mask','dat','math output=1 | mutter half=n v0=4800 tp=1.6')

for filt in range(2):
    pef = 'pef%d' % filt
    lag = 'lag%d' % filt
    decon = 'decon%d' % filt
    data = 'shot' + str(filt)
    Flow([pef,lag],'dat mask',
         'pef a=%s maskin=${SOURCES[1]} lag=${TARGETS[1]}' %
         ['30,1','20,5'][filt])
    Flow(decon,['dat',pef],'helicon filt=${SOURCES[1]}')
    Plot(decon,showdat('PEF*Data',2.8+filt,f2=1+3*filt,wantaxis1=0))
    Plot(data,'dat',showdat('Data',2.8+filt,flip='reverse which=2 |'))
    Result(decon,[data,decon],'SideBySideAniso')

def showspec(title):
    return '''
    spectra all=y |
    graph title="%s" crowd1=0.85
    label1=Frequency unit1=Hz label2=Amplitude unit2=
    wheretitle=bot wherexlabel=top
    ''' % title

Plot('specdat','dat',showspec('Spectrum Data'))
Plot('specdecon','decon0',showspec('Spectrum PEF*Data'))
Result('specdecon','specdat specdecon','OverUnderAniso')

def showdat(title):
    return '''
    window f2=10 f1=10 |
    grey crowd1=0.85  crowd2=0.8
    label1=Time unit1=s label2=X unit2=ft title="%s"
    ''' % title

Flow(['zofpef','zoflag'],'zof',
     'window max3=6000 | pef a=10,2 lag=${TARGETS[1]}')
Flow('zofdecon','zof zofpef','helicon filt=${SOURCES[1]}')
Plot('zof','zof',showdat('Data'))
Plot('zofdecon',showdat('10x2 PEF * Data'))
Result('zof','zof zofdecon','OverUnderAniso')

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