sconstruct

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

from rsfproj import *

Fetch('mona.img','imgs')

Flow('mona','mona.img',
     '''
     echo n1=512 n2=513 in=$SOURCE data_format=native_uchar |
     dd type=float
     ''',stdin=0)

colors = {
    'i': 'Greyscale',
    'j': 'Jet',
    't': 'Traffic',
    'h': 'Hot',
    'p': 'Pink',
    'b': 'Bone',
    'f': 'Flag',
    'a': 'Rainbow',
    'c': 'Cool',
    'l': 'Linear',
    'e': 'Blue-white-red',
    'g': 'Black-white-red',
    'H': 'Hot Reversed',
    'iC': 'Greyscale Clipped'
    }

monas=[]
for col in colors.keys():
    if col != 'H' and col != 'iC':
        mona = 'mona-%s' % col
        monas.append(mona)
        Plot(mona,'mona',
             '''
             window n2=100 f2=130 |
             grey allpos=y pclip=100 scalebar=y crowd1=0.8
             color=%s title="%s (color=%s)"
             transp=n wantaxis=n titlesz=24
             ''' % (col,colors[col],col))
        
Plot('mona1',monas[:6],'OverUnderAniso')
Plot('mona2',monas[6:],'OverUnderAniso')
Result('mona','mona1 mona2','SideBySideAniso',vppen='txscale=1.5')

for col in ('h','H'):
    mona = 'mona+%s' % col
    Plot(mona,'mona',
         '''
         grey allpos=y pclip=100 
         color=%s title="%s (color=%s)"
         transp=n wantaxis=n
         ''' % (col,colors[col],col))
Result('monacap','mona+h mona+H','SideBySideIso',vppen='ysize=10 xsize=10')

for col in ('i','iC'):
    mona = 'mona+%s' % col
    Plot(mona,'mona',
         '''
         grey allpos=y pclip=98 screenratio=1
         color=%s title="%s (color=%s)"
         transp=n wantaxis=n
         ''' % (col,colors[col],col))
Result('monaclip','mona+i mona+iC','SideBySideIso',vppen='ysize=10 xsize=10')
    
End()