陈桂廷的个人博客分享 http://blog.sciencenet.cn/u/jinshi2015 Reproducible Research Philosophy

博文

如何用Seismic Unix画出漂亮的地震剖面图

已有 10479 次阅读 2017-11-10 21:40 |个人分类:地球物理学软件|系统分类:科研笔记| 画图

如何画出SCI文章里漂亮的地震剖面, 不需要借助别的软件,强大的SU画图功能就能实现,看如下代码:


psimage<interp_dft.bin  n1=500 n2=500 f1=0 d1=12 f2=0 d2=12 title="DFT interpolation" perc=99  legend=1 lstyle=vertright    lbeg=-0.1 lend=0.1 ldnum=0.05  lgrid=dot height=6.0 width=6.0 label1="Depth(m)" label2="Lateral(m)" verbose=1 lwidth=0.3 lheight=6  n1tic=4 f1num=0 d1num=1000 n2tic=4 f2num=0 d2num=1000  axeswidth=2 labelsize=20 >interp_dft.eps

interp_dft.bin 为二进制文件,大小500*500

这里采用的psimage命令,输出EPS格式。

n1,n2分别是x和y方向的采样个数。

f1,f2两轴起点坐标。

d1,d2,采样间隔。

这里legend参数非常重要,详细请看:

有时候需要几幅图片的scalebar范围一致,用lbeg和lend非常方便,这里ldnum定义标尺的间隔。

lgrid=dot 标尺内部风格。

lwidth=0.3 lheight=6定义了标尺的长宽。

n1tic=4则定义了大刻度里面有4个小刻度,默认是没有刻度的。

f1num=0 d1num=1000 则定义了大刻度的起点和间隔。

axeswidth=2 坐标轴加粗,可以使图片更加美观。

更多参数请参看psimage 命令。


这样就画出了非常漂亮的剖面图,psimage简单易学,输入只需要二进制文件,输出直接eps,比madagascar的sfgrey更灵活漂亮,注意这里y轴的刻度数值是横着显示的,马达的为竖着的。



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

上一篇:kirchhoff叠前时间偏移源代码(SU修改版)
收藏 IP: 159.226.117.*| 热度|

0

该博文允许注册用户评论 请点击登录 评论 (1 个评论)

数据加载中...

Archiver|手机版|科学网 ( 京ICP备07017567号-12 )

GMT+8, 2024-11-29 19:36

Powered by ScienceNet.cn

Copyright © 2007- 中国科学报社

返回顶部