使用Madagascar绘制剖面图

1.查看要绘制的segy文件的总道数，以确定绘制范围

sfsegyread < 文件名.segy endian=1 | sfin

n1 : 每一道的采样点数
n2 : 总道数

2.快速预览从第a道起的b道剖面（注意:a+b<n2）

sfsegyread < 文件名.segy endian=1 | sfput d1=0.01 unit1=km label1=Depth title= 标题 |sfpow powl=0 | sfwindow f2=a n2=b | sfgrey pclip=97 | sfpen &

如果想加色标，可以这么写:

sfsegyread < 文件名.segy endian=1 | sfput d1=0.01 unit1=km label1=Depth title=标题 |sfpow powl=0 | sfwindow f2=a n2=b | sfgrey pclip=97 scalebar=y | sfpen &

3.将剖面图转为vpl文件

sfsegyread < 文件名.segy endian=1 | sfput d1=0.01 unit1=km label1=Depth title=标题 |sfpow powl=1 | sfwindow f2=a n2=b | sfgrey pclip=97 > 文件名.vpl

4.把vpl文件转换为pdf文件

vpconvert 文件名.vpl 文件名.pdf bgcolor=white

如果想查看pdf，可以输入命令 evince 文件名.pdf 进行查看

5.把pdf文件转换为png文件

convert 文件名.pdf 图片名.png

6.显示png图片

display 图片名.png &