GMT常用绘图模板及技巧
合集下载
相关主题
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
echo G 0.05i >legend_china.txt echo S 0.15i c 0.3 blue 0.25p 0.3i 20151009-20151015(14)>> legend_china.txt echo S 0.15i c 0.3 red 0.25p 0.3i @;red;20151015-20151022(9)@;;>> legend_china.txt gmt pslegend legend_china.txt -Dx0/0/3i/0.6i/BL -F+gwhite+pblack -R J -K -O >> %PS%
房立华，2013
6
尝试读一个样例脚本：
数据来源：中国地震局地震研究 所，GNSS中心
7
set PS=JP.ps set R=-R123/149/20/47 set J=-JM4i gmt gmtset FORMAT_GEO_MAP dddF gmt gmtset FONT_ANNOT_PRIMARY 7 gmt pscoast %R% %J% -BWSEN -B5g5 -N2 -W1 -A1000 -Glightbrown Slightblue -Lf128/45/35/400+lkm+jt -K > %PS%
只画文件尾
gmt psxy %R% %J% -T -O >> a.p
17
绘图的基本要求
符合规范：字体大小，线条粗细 遵守约定俗成的规则：高速蓝色，低速红色；震中圆圈 信息突出、数据准确 清晰、清楚、美观 自己看着舒服，同行看着也舒服：颜色搭配
图: 黑白图和彩色图的分辨率不能低于600 dpi, 图中线条要清晰, 线条粗细约0.5-
— 图 层 BODY
— 文件尾 TRAILER
15
打“开”文件 忽略文件尾
psxy …-K > a.ps
psxy …-O –K >> a.ps
“O”verlay 忽略文件头
psxy…-O >> a.ps
16
只画文件头
gmt psxy %R% %J% -T -K > a.ps
gmt psxy …-O –K >> a.ps gmt psxy …-O –K >> a.ps
37
1、读懂脚本
http://icgem.gfzpotsdam.de/ICGEM/
26
布格异常图
本图主要涉及命令为 grdimage 、psxy、 pslegend; 布格异常数据来源： BGI（ http://bgi.omp.obsmip.fr/ ），World Gravity Map (参考文献： Spherical harmonic modelling to ultra-high degree of Bouguer and isostatic anomalies, J.Geod., Balmino et al., 2011.） 亦可根据DEM模型，对自 由空气异常进行地形改正 获得布格异常。
27
28
gawk "{print $5,$4,0.08*$3 }" 上周.txt |gmt psxy -R -J -Sc -W0.5,black, Gblue -K -O >> %PS% gawk "{print $5,$4,0.08*$3 }" 本周.txt |gmt psxy -R -J -Sc -W0.5,black, Gred -K -O >> %PS%
9
2.1 分析图层
文字 标注 站点 断层 行政 边界 陆地、 海洋
10
文字 标注 （n） 站点 断层 （2） 地形 陆地、 海洋
Shapiro,2005
11
2.2 准备数据
① 海岸线数据 （GMT自带）
② DEM数据 ③ 颜色表CPT （30″） （GMT自带）
④ 断层经纬度、名称（自行准备） ⑤ 站点经纬度、名称（自行准备）
12
2.3 找到对应模块（module） pstext psxy psxy grdgradient + grdimage
文字 标注 （n） 站点 断层 （2） 地形
pscoast
陆地、 海洋
13
2.4 组织顺序，查文档，编写脚本
图层叠加顺序 PostScript生成顺序
14
— 文件头 HEADER
重力管理部制图
34
gmt psvelo %D% -Sx0.01 %R% %J% -O -K -W1.5p,blue >> %PS% echo 90 24 40 40 45 |gmt psvelo -Sx0.01 %R% %J% -W1p,0/0/255 -O -K>>%PS% echo 90.5 24 40 nanostrain/yr|gmt pstext %R% %J% -F+f10+jLM -O >> %PS%
35
MAP_VECTOR_SHA PE [0,2)
gmt psvelo %D% -Se0.03/0.95/0 %R% %J% -G0/0/255 -h1 -K -O >>%PS%
36
psmeca -R -J -Sm0.7/0/0 -L0.5p,black -Ggreen -O -K >> %PS%
房立华，2013
18
常用投影
直角坐标系：线性投影 -JX 地理坐标系： 全球： 半球： 区域： -JN Robinson 投影 /-JY/-JR -JA –JS 美国 全国100万： 省区图： 一些特殊投影： -JP 极坐标 不知道用什么： 等角圆柱墨卡托（-JM）；流行、通用 等距圆柱投影（-JQ）：适合纬度跨度大
31
剪应变率场
32
33
重力变化点值时间序列
gmt psbasemap -R$x1\T/$x2\T/$y1/$y2 -JX16c/7c -Bsxa1Y -Bpxa3of1o – Bya10f2 -BSn -P -K > obser.ps gmt psxy obser.dat -R -JX -Sd0.25c -Ey -Gred -K -O >> obser.ps gmt psxy mean.dat -R -J -W1.0,green,- -K -O >> obser.ps
推荐阅读：GMT添加断层名 https:///plot-fault-names.html
30
重力变化图
gmt surface xyg -R -I5m G%GRD% gmt grdimage %GRD% Ccpt.cpt -J -R -K -O >> %outputps% gmt grdcontour %GRD% -R -J C1.5 -A3 -O -K >> %outputps%
3、常用绘图模板
4、常用绘图技巧 5、学习方法及学习资源
3
Window s? or Linux?
温馨提示：本报告中脚本均为windows平台下使用
4
GMT4? or GMT5?
5
模块（Module）+ 选项（Arguments）= 命令行
gmt psbasemap -R0/10/0/10 -JX15c -BWSNE+tTitle -B2g2 > test.ps
gawk "{print $3, $2}" Japan.llh | gmt psxy %R% %J% -Sc0.07 -Gdarkred -O >> %PS%
gmt psconvert -A+r %PS% -Tg -P
8
1、读懂脚本
提纲
CONTENT
2、编写自己的绘图脚本
3、常用绘图模板
4、常用绘图技巧 5、学习方法及学习资源
0.6 mm. 中文图中的汉字为7 pt的幼圆字体, 英文图中的文字和阿拉伯数字为8 pt 的Times New Roman字体. 图的宽度分两种: 单栏图宽4-8 cm, 通栏图宽12-17
cm. 图的长度一般不超过20 cm. 如果出现中国地图, 请注意中国地图的完整性(如国
界线和南海诸岛等). 南海诸岛作插图时, 务必完整表达附图中与主图相同区域内的 专题要素。
地壳形变资料分析与应用培训班 2016 · 09
/
GMT常用绘图模板及技巧
刘珠妹 中国地震局地震研究所
1
1、读懂脚本
提纲
CONTENT
2、编写自己的绘图脚本
3、常用绘图模板
4、常用绘图技巧 5、学习方法及学习资源
2
1、读懂脚本
提纲
CONTENT
2、编写自己的绘图脚本
19
-JB Albers等面积投影 –JA Lambert –JB Albers等面积投影
-JQ -JN -JR -JA -JY
20
2.5 调试
“pause” → CMD Notepad++
推荐阅读：Windows 下使用 GMT 的正确姿势 https:///how-to-use-gmt-under-windows.html
24
gmt psxy %R% %J% -T -K -X2i> Tibet.ps
gmt grdimage %R% %J% -Ba10g10 Tibet.grd -Iint.grad -CTibet.cpt -K -O >> Tibet.ps gmt pscoast %R% %J% -N1/0.5p,white -Ia/0.15p,177/178/183 I1/0.5p,61/99/172 -C81/174/254 -Lg86/18+c30+w1000+u+f -K -O >> Tibet.ps gmt psscale -D7i/2.7i/2.5i/0.15i -Ccpt1.cpt -Bx1000f200 -By+l"m" -K -O >> Tibet.ps gmt psscale -D7i/1.6/1.3i/0.15i -Ccpt2.cpt -L -O -K>> Tibet.ps gmt psxy %R% %J% -T -O>> Tibet.ps
29
REM 断层 gmt psxy %R% %J% D:\04数据\断层\Asia_faults.Lnormal -Sf-2/0.1+l+f W1.5p,brown -Gbrown -O -K >> %PS% echo 100.1 37 Qinghai Lake |gmt pstext %R% %J% F+a0+jCM+f11,4,blue -O -K >> %PS% echo 100.7 36.6 Nanshan Fault |gmt pstext %R% %J% -F+a26+jLM+f11,6 -O -K >> %PS%
21
2.6 输出
psconvert
gmt psconvert -Tj -C-sFONTPATH=C:\windows\fonts %PS% -A+r
-E ?
22
1、ห้องสมุดไป่ตู้懂脚本
提纲
CONTENT
2、编写自己的绘图脚本
3、常用绘图模板
4、常用绘图技巧 5、学习方法及学习资源
23
gmt psbasemap -R110/111.5/30.5/31.3 -JM6.5i -B0.2 -K > sx.ps gmt psxy 三峡库首区_县界.txt -R -J -K -O -W0.8p,gray,.>> sx.ps gmt psxy 三峡库首区_居民点.txt -R -J -Sc0.1 -K -O >> sx.ps gmt pstext 三峡库首区_居民点.txt -R -J -F+f7,35+jCM -O -D0/0.2 >>sx.ps
25
自由空气异常图
本图主要涉及命令 为grdimage 、psxy、 pslegend; 本图自由空气异常 优先采用，直接用img2grd读取 数据来源（1分格网）
/ma rine_grav/mar_grav.htm l其中，陆地重力异
常来自EGM2008； 海洋重力异常为卫 星测高数据； 亦可下载重力场模 型（球谐系数）
命令行 + 命令行 =
脚本文件
gmt psbasemap -R0/10/0/10 -JX15c -BWSNE+tTitle -B2g2 -K > test.ps gmt psxy drawDif.txt -R -J -Sc0.3 -W0.5p,blue -Gred -O -K >> test.ps echo 5 5 测试| gmt pstext -R -J -F+f16,35,red+jLM -O >> test.ps
房立华，2013
6
尝试读一个样例脚本：
数据来源：中国地震局地震研究 所，GNSS中心
7
set PS=JP.ps set R=-R123/149/20/47 set J=-JM4i gmt gmtset FORMAT_GEO_MAP dddF gmt gmtset FONT_ANNOT_PRIMARY 7 gmt pscoast %R% %J% -BWSEN -B5g5 -N2 -W1 -A1000 -Glightbrown Slightblue -Lf128/45/35/400+lkm+jt -K > %PS%
只画文件尾
gmt psxy %R% %J% -T -O >> a.p
17
绘图的基本要求
符合规范：字体大小，线条粗细 遵守约定俗成的规则：高速蓝色，低速红色；震中圆圈 信息突出、数据准确 清晰、清楚、美观 自己看着舒服，同行看着也舒服：颜色搭配
图: 黑白图和彩色图的分辨率不能低于600 dpi, 图中线条要清晰, 线条粗细约0.5-
— 图 层 BODY
— 文件尾 TRAILER
15
打“开”文件 忽略文件尾
psxy …-K > a.ps
psxy …-O –K >> a.ps
“O”verlay 忽略文件头
psxy…-O >> a.ps
16
只画文件头
gmt psxy %R% %J% -T -K > a.ps
gmt psxy …-O –K >> a.ps gmt psxy …-O –K >> a.ps
37
1、读懂脚本
http://icgem.gfzpotsdam.de/ICGEM/
26
布格异常图
本图主要涉及命令为 grdimage 、psxy、 pslegend; 布格异常数据来源： BGI（ http://bgi.omp.obsmip.fr/ ），World Gravity Map (参考文献： Spherical harmonic modelling to ultra-high degree of Bouguer and isostatic anomalies, J.Geod., Balmino et al., 2011.） 亦可根据DEM模型，对自 由空气异常进行地形改正 获得布格异常。
27
28
gawk "{print $5,$4,0.08*$3 }" 上周.txt |gmt psxy -R -J -Sc -W0.5,black, Gblue -K -O >> %PS% gawk "{print $5,$4,0.08*$3 }" 本周.txt |gmt psxy -R -J -Sc -W0.5,black, Gred -K -O >> %PS%
9
2.1 分析图层
文字 标注 站点 断层 行政 边界 陆地、 海洋
10
文字 标注 （n） 站点 断层 （2） 地形 陆地、 海洋
Shapiro,2005
11
2.2 准备数据
① 海岸线数据 （GMT自带）
② DEM数据 ③ 颜色表CPT （30″） （GMT自带）
④ 断层经纬度、名称（自行准备） ⑤ 站点经纬度、名称（自行准备）
12
2.3 找到对应模块（module） pstext psxy psxy grdgradient + grdimage
文字 标注 （n） 站点 断层 （2） 地形
pscoast
陆地、 海洋
13
2.4 组织顺序，查文档，编写脚本
图层叠加顺序 PostScript生成顺序
14
— 文件头 HEADER
重力管理部制图
34
gmt psvelo %D% -Sx0.01 %R% %J% -O -K -W1.5p,blue >> %PS% echo 90 24 40 40 45 |gmt psvelo -Sx0.01 %R% %J% -W1p,0/0/255 -O -K>>%PS% echo 90.5 24 40 nanostrain/yr|gmt pstext %R% %J% -F+f10+jLM -O >> %PS%
35
MAP_VECTOR_SHA PE [0,2)
gmt psvelo %D% -Se0.03/0.95/0 %R% %J% -G0/0/255 -h1 -K -O >>%PS%
36
psmeca -R -J -Sm0.7/0/0 -L0.5p,black -Ggreen -O -K >> %PS%
房立华，2013
18
常用投影
直角坐标系：线性投影 -JX 地理坐标系： 全球： 半球： 区域： -JN Robinson 投影 /-JY/-JR -JA –JS 美国 全国100万： 省区图： 一些特殊投影： -JP 极坐标 不知道用什么： 等角圆柱墨卡托（-JM）；流行、通用 等距圆柱投影（-JQ）：适合纬度跨度大
31
剪应变率场
32
33
重力变化点值时间序列
gmt psbasemap -R$x1\T/$x2\T/$y1/$y2 -JX16c/7c -Bsxa1Y -Bpxa3of1o – Bya10f2 -BSn -P -K > obser.ps gmt psxy obser.dat -R -JX -Sd0.25c -Ey -Gred -K -O >> obser.ps gmt psxy mean.dat -R -J -W1.0,green,- -K -O >> obser.ps
推荐阅读：GMT添加断层名 https:///plot-fault-names.html
30
重力变化图
gmt surface xyg -R -I5m G%GRD% gmt grdimage %GRD% Ccpt.cpt -J -R -K -O >> %outputps% gmt grdcontour %GRD% -R -J C1.5 -A3 -O -K >> %outputps%
3、常用绘图模板
4、常用绘图技巧 5、学习方法及学习资源
3
Window s? or Linux?
温馨提示：本报告中脚本均为windows平台下使用
4
GMT4? or GMT5?
5
模块（Module）+ 选项（Arguments）= 命令行
gmt psbasemap -R0/10/0/10 -JX15c -BWSNE+tTitle -B2g2 > test.ps
gawk "{print $3, $2}" Japan.llh | gmt psxy %R% %J% -Sc0.07 -Gdarkred -O >> %PS%
gmt psconvert -A+r %PS% -Tg -P
8
1、读懂脚本
提纲
CONTENT
2、编写自己的绘图脚本
3、常用绘图模板
4、常用绘图技巧 5、学习方法及学习资源
0.6 mm. 中文图中的汉字为7 pt的幼圆字体, 英文图中的文字和阿拉伯数字为8 pt 的Times New Roman字体. 图的宽度分两种: 单栏图宽4-8 cm, 通栏图宽12-17
cm. 图的长度一般不超过20 cm. 如果出现中国地图, 请注意中国地图的完整性(如国
界线和南海诸岛等). 南海诸岛作插图时, 务必完整表达附图中与主图相同区域内的 专题要素。
地壳形变资料分析与应用培训班 2016 · 09
/
GMT常用绘图模板及技巧
刘珠妹 中国地震局地震研究所
1
1、读懂脚本
提纲
CONTENT
2、编写自己的绘图脚本
3、常用绘图模板
4、常用绘图技巧 5、学习方法及学习资源
2
1、读懂脚本
提纲
CONTENT
2、编写自己的绘图脚本
19
-JB Albers等面积投影 –JA Lambert –JB Albers等面积投影
-JQ -JN -JR -JA -JY
20
2.5 调试
“pause” → CMD Notepad++
推荐阅读：Windows 下使用 GMT 的正确姿势 https:///how-to-use-gmt-under-windows.html
24
gmt psxy %R% %J% -T -K -X2i> Tibet.ps
gmt grdimage %R% %J% -Ba10g10 Tibet.grd -Iint.grad -CTibet.cpt -K -O >> Tibet.ps gmt pscoast %R% %J% -N1/0.5p,white -Ia/0.15p,177/178/183 I1/0.5p,61/99/172 -C81/174/254 -Lg86/18+c30+w1000+u+f -K -O >> Tibet.ps gmt psscale -D7i/2.7i/2.5i/0.15i -Ccpt1.cpt -Bx1000f200 -By+l"m" -K -O >> Tibet.ps gmt psscale -D7i/1.6/1.3i/0.15i -Ccpt2.cpt -L -O -K>> Tibet.ps gmt psxy %R% %J% -T -O>> Tibet.ps
29
REM 断层 gmt psxy %R% %J% D:\04数据\断层\Asia_faults.Lnormal -Sf-2/0.1+l+f W1.5p,brown -Gbrown -O -K >> %PS% echo 100.1 37 Qinghai Lake |gmt pstext %R% %J% F+a0+jCM+f11,4,blue -O -K >> %PS% echo 100.7 36.6 Nanshan Fault |gmt pstext %R% %J% -F+a26+jLM+f11,6 -O -K >> %PS%
21
2.6 输出
psconvert
gmt psconvert -Tj -C-sFONTPATH=C:\windows\fonts %PS% -A+r
-E ?
22
1、ห้องสมุดไป่ตู้懂脚本
提纲
CONTENT
2、编写自己的绘图脚本
3、常用绘图模板
4、常用绘图技巧 5、学习方法及学习资源
23
gmt psbasemap -R110/111.5/30.5/31.3 -JM6.5i -B0.2 -K > sx.ps gmt psxy 三峡库首区_县界.txt -R -J -K -O -W0.8p,gray,.>> sx.ps gmt psxy 三峡库首区_居民点.txt -R -J -Sc0.1 -K -O >> sx.ps gmt pstext 三峡库首区_居民点.txt -R -J -F+f7,35+jCM -O -D0/0.2 >>sx.ps
25
自由空气异常图
本图主要涉及命令 为grdimage 、psxy、 pslegend; 本图自由空气异常 优先采用，直接用img2grd读取 数据来源（1分格网）
/ma rine_grav/mar_grav.htm l其中，陆地重力异
常来自EGM2008； 海洋重力异常为卫 星测高数据； 亦可下载重力场模 型（球谐系数）
命令行 + 命令行 =
脚本文件
gmt psbasemap -R0/10/0/10 -JX15c -BWSNE+tTitle -B2g2 -K > test.ps gmt psxy drawDif.txt -R -J -Sc0.3 -W0.5p,blue -Gred -O -K >> test.ps echo 5 5 测试| gmt pstext -R -J -F+f16,35,red+jLM -O >> test.ps