sconstruct

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

from rsfproj import *

Fetch('dynamics.H','toldi')

Result('toldi','dynamics.H',
       'window n1=375 n2=1 f2=9 | dd form=native type=float | grey wanttitle=n allpos=y')

Flow('cmp','dynamics.H',
     'window f3=5 n3=1 | dd form=native type=float | add add=-128 | nmostretch v0=1.5 half=n inv=y')

near = [0,5,9]
far = map(lambda x: 27-x,near)

for f2 in near+far:
    win = 'win%d' % f2
    Flow(win,'cmp','window f2=%d n2=1' % f2)

for case in range(2):
    plot = 'case%d' % case
    Plot(plot,'win%d  win%d win%d' % map(tuple,[near,far])[case],
         '''
         cat axis=2 ${SOURCES[1:3]} |
         wiggle wanttitle=n plotcol=%d zplot=.4 wantaxis2=n xll=0 xur=13 yll=.5 yur=10.
         ''' % (6-case))

Result('reciptrace','case0 case1','Overlay')

for s in range(2):
    slice = 'slice%d' % s
    Plot(slice,'dynamics.H',
         'window min1=%g n1=1 | dd form=native type=float | grey wanttitle=n allpos=y' % (1,2.5)[s])

Result('recipslice','slice0 slice1','OverUnderAniso')

End()