sconstruct

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

from rsfproj import *

data = 'wz.24.H'
Fetch(data,'wz')

Flow('data',data,
     '''
     dd form=native | window n1=400 n2=24 j1=2 | tpow tpow=2 | reverse which=2
     ''')
Flow('nmo','data','nmostretch v0=1.44 half=n')

Plot('data','grey wanttitle=n labelsz=12')
Plot('nmo','grey wanttitle=n labelsz=12')
Result('stretch','data nmo','SideBySideAniso')

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