GMT制图
合集下载
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
GMT绘图
报告人：王浩璠 学号：201421043 专业：地球探测与信息技术
常用的画图工具： Excel Origin CorelDRAW ......
GMT绘图 GMT主要是由Paul Wessel和Walter H.F.Smith以及全球 的一群志愿者组织开发并且定时更新的绘图软件。它是 国际地学界普遍使用的绘图软件。 优点：数据量大，以数万计； 分辨率高，高于Excel、Origin等软件； 功能强大，可以绘制二维、三维图件； 图件规范且标准； 缺点：不支持中文
常用图件：
1. 二维线性图 2. 二维散点图 3. 频率分布直方图 4. 平面图 5. 地震剖面、地震道显示 6. 测井曲线图 7. 三维立体图
1. 二维线性图
GMT编写： psxy data1.txt -R0/50/0/0.15 -JX6i/4i -B10:"@a":/0.1:"Rpp":enWS -W0.8p,black -P -K > data.ps psxy data2.txt -R -JX -W0.8p,red -P -O >> data.ps
x,y
7. 三维立体图
GMT编写： ①创建cpt文件： makecpt -Crainbow -T0/11/3 -Z > E_f.cpt ②创建网格： gmtconvert data.asc -F0,1,3 | surface -R7/24/0/100 -I0.01/0.01 -GE_f.grd
x,y,z
What is more?
坐标轴设置： gmtset PAPER_MEDIA A4 （纸张大小） gmtset ANOT_FONT_SIZE 12 (坐标轴刻度字体) gmtset LABEL_OFFSET 10p (刻度与轴名的距离)
......
GMT是基于linux系统上完成的。
sh 89well.sh 运行GMT程序
7. 三维立体图
③成图： grdview E_f.grd -R7/24/0/100/0/11 -JX5i/5i -JZ3i -CE_f.cpt -Wthin -P -K -Ba2g2:"Pore Pressure,MPa":/a5g5:"CO@-2@Saturation":/a2g2:"E_f":wSEnZ+ -K > E_f.ps ④色棒： psscale -CE_f.cpt -D7.5i/2i/3i/0.2i -P -O -B/:"E_f": >> E_f.ps
x,y,z
-N: 显示小于0的值
6. 测井曲线图
GMT编写： gmtconvert g89.txt -F1,0 | psxy -R7.0/26.0/-400/-300 -JX1i/7i B20:"CAL,cm":/a100f100:"Dep th":eNWs -W0.2p -K -X1i Y10i -P > g89_welllog.ps gmtconvert g89.txt -F2,0 | psxy -R0/113/-400/-300 -JX1i/7i B50:"GR,API":/f100::eNws -W0.2p -X1i -O -P -K >> g89_welllog.ps .........
x
Байду номын сангаас 4. 平面图
GMT编写： ①创建cpt文件： makecpt -Crainbow -T0/10000/5000 -I5000 >amp.cpt ②创建网格： gmtconvert data.asc -F0,1,6 | xyz2grd -Gg89_t7.grd -I25/25 -R580575/585575/4111975/4116950
x,y
2. 二维散点图
GMT编写：
x, y, 井号
① 投点：psxy data.asc R581000/585000/4.113e+06/4.1157e +06 -Jx4i/6i -B1000:"X (m)":/500:"Y (m)":enWS -W0.5p,blue -Gblue -Sc0.08i -P -K V >data.ps
x,y,z
4. 平面图
③成图： grdimage -P g89_t7.grd -JX4i/4i R580575/585575/4111975/4116950 B2000:"X":/2000:"Y"::.AMP_t3: WNse -Camp.cpt -K > AMP_t3.ps ④色棒： psscale -D2i/-0.2i/3i/0.2ih -Camp.cpt -B5000/:"AMP": -O >> AMP_t3.ps
5. 地震剖面、地震道显示
GMT编写： gmtconvert traces.txt -F0,4,5 |pswiggle -R0/34/-162/0/-0.055/0.075 -JX1i/5i -B10:"Angle":/20:"":esWN -Z0.9 -W2ta -G0 -P > traces.ps
......
2. 二维散点图
x, y, 图例标注
②图例： psxy xy_data.asc -R -Jx -Sc0.08i Gblue -O -K -V >> data.ps awk '{printf "%g %s 9 0 1 5 %15s\n", $1+80.0, $2, $3}' xy_data.asc | pstext -R -Jx Gblack -O -K -V -N >>data.ps ......
Thank you !
2. 二维散点图
③ 井号： awk '{printf "%10g %s 9 0 1 5 %15s\n", $1-20, $2-100, $3}' data.asc | pstext -R -Jx -O -Gblack -K -N >>data.ps
......
3. 频率分布直方图
GMT编写： gmtconvert data.asc -F2 | pshistogram Ba400f100:"Max_Amp":/a10f5:"F requency"::,%::."orig_1992_cap": WSne -R0/2500/0/50 -JX4.8i/2.4i -G200 -X1.5i -Y28i L0.5p -W200 -P >Max_Amp_histogram.ps
报告人：王浩璠 学号：201421043 专业：地球探测与信息技术
常用的画图工具： Excel Origin CorelDRAW ......
GMT绘图 GMT主要是由Paul Wessel和Walter H.F.Smith以及全球 的一群志愿者组织开发并且定时更新的绘图软件。它是 国际地学界普遍使用的绘图软件。 优点：数据量大，以数万计； 分辨率高，高于Excel、Origin等软件； 功能强大，可以绘制二维、三维图件； 图件规范且标准； 缺点：不支持中文
常用图件：
1. 二维线性图 2. 二维散点图 3. 频率分布直方图 4. 平面图 5. 地震剖面、地震道显示 6. 测井曲线图 7. 三维立体图
1. 二维线性图
GMT编写： psxy data1.txt -R0/50/0/0.15 -JX6i/4i -B10:"@a":/0.1:"Rpp":enWS -W0.8p,black -P -K > data.ps psxy data2.txt -R -JX -W0.8p,red -P -O >> data.ps
x,y
7. 三维立体图
GMT编写： ①创建cpt文件： makecpt -Crainbow -T0/11/3 -Z > E_f.cpt ②创建网格： gmtconvert data.asc -F0,1,3 | surface -R7/24/0/100 -I0.01/0.01 -GE_f.grd
x,y,z
What is more?
坐标轴设置： gmtset PAPER_MEDIA A4 （纸张大小） gmtset ANOT_FONT_SIZE 12 (坐标轴刻度字体) gmtset LABEL_OFFSET 10p (刻度与轴名的距离)
......
GMT是基于linux系统上完成的。
sh 89well.sh 运行GMT程序
7. 三维立体图
③成图： grdview E_f.grd -R7/24/0/100/0/11 -JX5i/5i -JZ3i -CE_f.cpt -Wthin -P -K -Ba2g2:"Pore Pressure,MPa":/a5g5:"CO@-2@Saturation":/a2g2:"E_f":wSEnZ+ -K > E_f.ps ④色棒： psscale -CE_f.cpt -D7.5i/2i/3i/0.2i -P -O -B/:"E_f": >> E_f.ps
x,y,z
-N: 显示小于0的值
6. 测井曲线图
GMT编写： gmtconvert g89.txt -F1,0 | psxy -R7.0/26.0/-400/-300 -JX1i/7i B20:"CAL,cm":/a100f100:"Dep th":eNWs -W0.2p -K -X1i Y10i -P > g89_welllog.ps gmtconvert g89.txt -F2,0 | psxy -R0/113/-400/-300 -JX1i/7i B50:"GR,API":/f100::eNws -W0.2p -X1i -O -P -K >> g89_welllog.ps .........
x
Байду номын сангаас 4. 平面图
GMT编写： ①创建cpt文件： makecpt -Crainbow -T0/10000/5000 -I5000 >amp.cpt ②创建网格： gmtconvert data.asc -F0,1,6 | xyz2grd -Gg89_t7.grd -I25/25 -R580575/585575/4111975/4116950
x,y
2. 二维散点图
GMT编写：
x, y, 井号
① 投点：psxy data.asc R581000/585000/4.113e+06/4.1157e +06 -Jx4i/6i -B1000:"X (m)":/500:"Y (m)":enWS -W0.5p,blue -Gblue -Sc0.08i -P -K V >data.ps
x,y,z
4. 平面图
③成图： grdimage -P g89_t7.grd -JX4i/4i R580575/585575/4111975/4116950 B2000:"X":/2000:"Y"::.AMP_t3: WNse -Camp.cpt -K > AMP_t3.ps ④色棒： psscale -D2i/-0.2i/3i/0.2ih -Camp.cpt -B5000/:"AMP": -O >> AMP_t3.ps
5. 地震剖面、地震道显示
GMT编写： gmtconvert traces.txt -F0,4,5 |pswiggle -R0/34/-162/0/-0.055/0.075 -JX1i/5i -B10:"Angle":/20:"":esWN -Z0.9 -W2ta -G0 -P > traces.ps
......
2. 二维散点图
x, y, 图例标注
②图例： psxy xy_data.asc -R -Jx -Sc0.08i Gblue -O -K -V >> data.ps awk '{printf "%g %s 9 0 1 5 %15s\n", $1+80.0, $2, $3}' xy_data.asc | pstext -R -Jx Gblack -O -K -V -N >>data.ps ......
Thank you !
2. 二维散点图
③ 井号： awk '{printf "%10g %s 9 0 1 5 %15s\n", $1-20, $2-100, $3}' data.asc | pstext -R -Jx -O -Gblack -K -N >>data.ps
......
3. 频率分布直方图
GMT编写： gmtconvert data.asc -F2 | pshistogram Ba400f100:"Max_Amp":/a10f5:"F requency"::,%::."orig_1992_cap": WSne -R0/2500/0/50 -JX4.8i/2.4i -G200 -X1.5i -Y28i L0.5p -W200 -P >Max_Amp_histogram.ps