本博文的示例采用 SConstruct 脚本来书写,关于如何采用 SConstruct 写 Seismic Unix 脚本详见另一篇博文 http://nickwhyy.top/su-scons-together/

二维数据

灰度图

from rsf.proj import *
import rsf.suproj as sp

sp.Flow('velocity.eps','velocity.dat',
    '''
    psimage
    n1=300 d1=12.192 d2=12.192
    d1num=1000 f1num=0 n1tic=1
    d2num=2000 f2num=0 n2tic=1
    width=3.33 height=1.6
    label2='Distance (m)' label1='Depth (m)' labelsize=8
    legend=1 lstyle=vertright units='Velocity (m/s)'
    lwidth=0.1 lheight=1.6
    ''')

End()

捕获.PNG

彩色图

from rsf.proj import *
import rsf.suproj as sp

sp.Flow('vel-color.eps','velocity.dat',
    '''
    psimage
    n1=300 d1=12.192 d2=12.192
    d1num=1000 f1num=0 n1tic=1
    d2num=2000 f2num=0 n2tic=1
    width=3.33 height=1.6
    label2='Distance (m)' label1='Depth (m)' labelsize=8
    legend=1 lstyle=vertright units='Velocity (m/s)'
    lwidth=0.1 lheight=1.6
    threecolor=1 wrgb=1.0,0,0 grgb=1.0,1.0,1.0 brgb=0,0,1.0
    ''')

End()

捕获1.PNG

波形图

sp.Flow('dn.eps','dn.bin',
    '''
    pswigb
    n1=501 d1=0.002 d2=1
    wbox=3.33 hbox=4
    label2='Trace #' label1='Time (s)' labelsize=8
    ''')

image.png

三维数据

灰度图

sp.Flow('salt3d_dat_txy.eps','salt3d_dat_txy.TRC',
    '''
    pscube n1=501 n2=250 n3=201
    d1=0.008 f1=0
    d2=0.04 f2=0
    d3=0.04 f3=0
    f3num=0.5
    d3num=2
    n3tic=1
    size1=4 size2=3.33 size3=1.5
    label1='Time (s)'
    label2='x (km)'
    label3='y (km)'
    labelsize=8
    perc=98
    bps=24
    legend=1
    lstyle=vertright
    units='Amplitude'
    lwidth=0.1
    lheight=4.0
    ''')

image.png

彩色图

sp.Flow('salt3d_vel_zxy.eps','salt3d_vel_zxy.TRC',
    '''
    pscube n1=201 n2=250 n3=250
    d1=0.04 f1=0
    d2=0.04 f2=0
    d3=0.04 f3=0
    f3num=0.5
    d3num=2
    n3tic=1
    size1=4 size2=3.33 size3=1.5
    label1='Depth (km)'
    label2='x (km)'
    label3='y (km)'
    labelsize=8
    perc=98
    bps=24
    legend=1
    lstyle=vertright
    units='Amplitude'
    lwidth=0.1
    lheight=4.0
    wrgb=0.0,0.0,1.0 grgb=1.0,0.0,0.0 brgb=1.0,1.0,1.0
    bps=24
    ''')

捕获4.PNG

参考资料

http://blog.sciencenet.cn/blog-2834901-1084725.html