使用Madagascar实现三维地震数据体的动态展示

#!/bin/sh
plot()
{   segyFile=3D.segy
    segyDir=.
    imgDir=.
    title=3D
    ny= number Of Inline
    nx= number Of Crossline
    crl0= the first Crossline to show
    inl0= the first Inline to show
    gain=98
    inlshow= the default inline to show
    movinl= how many inline to show
   ((movfrom=$movinl-$inl0))
   ((windowf3=$inlshow-$inl0))
    sfsegyread verb=1 endian=y tape=$segyFile | sfwindow f1=0 f2=0 f3=0 n1=trace length\
    | sfput n2=$ny n3=$nx label1=Sample label2=Crossline label3=Inline title=""\
     d1=1 d2=1 d3=1 unit1="" unit2="" unit3=""  o1=0 o2=${crl0} o3=${inl0} \
    | sfpow pow1=1 >$title.rsf 
    sfbyte <$title.rsf gainpanel=a  verb=1 >$title-byte.rsf
    sfbar <$title.rsf verb=1 >bar.rsf                                                                                    sfgrey3 <$title-byte.rsf  flat=n point1=0.8 point2=0.7 frame1=300 frame2=300 frame3=${movfrom} dframe=20 color=i movie=3 mean=y pclip=$gain scalebar=y bar=bar.rsf labelfat=2 verb=1 bgcolor=white screentatio=0.75 screanht=8  >$title-mov.vpl
    gdpen type=gif $title-mov.vpl bgcolor=white >$title-mov.gif  
    display $title-mov.gif &
}
plot 
rm *.rsf
exit