GMT常用绘图模板及技巧

1、GMT绘图模板及技巧刘珠地壳形变资料分析与应用培训班 2016 091CONTENT1、读懂脚本2、编写自己的绘图脚本3、常用绘图模板4、常用绘图技巧5、学习方法及学习资源提 纲2CONTENT1、读懂脚本2、编写自己的绘图脚本3、常用绘图模板4、常用绘图技巧5、学习方法及学习资源提 纲3Windows?orLinux?温馨提示：本报告中脚本均为windows平台下使用4GMT4?orGMT5?5模块（Module）+ 选项（Arguments）= 命令行gmt psbasemap -R0/10/0/10 -JX15c -BWSNE+tTitle -B2g2 test.psgmt psbase

2、map -R0/10/0/10 -JX15c -BWSNE+tTitle -B2g2 -K test.psgmt psxy drawDif.txt -R -J -Sc0.3 -W0.5p,blue -Gred -O -K test.psecho 5 5 测试| gmt pstext -R -J -F+f16,35,red+jLM -O test.ps房立华，20136命令行 + 命令行 = 脚本文件尝试读一个样例脚本：数据来源：中国地震局地震研究所，GNSS中心7set PS=JP.psset R=-R123/149/20/47set J=-JM4igmt gmtset FORMAT_GEO_

3、MAP dddFgmt gmtset FONT_ANNOT_PRIMARY 7gmt pscoast %R% %J% -BWSEN -B5g5 -N2 -W1 -A1000 -Glightbrown -Slightblue -Lf128/45/35/400+lkm+jt -K %PS%gawk print $3, $2 Japan.llh | gmt psxy %R% %J% -Sc0.07 -Gdarkred -O %PS%gmt psconvert -A+r %PS% -Tg -P8CONTENT1、读懂脚本2、编写自己的绘图脚本3、常用绘图模板4、常用绘图技巧5、学习方法及学习资源提 纲

4、92.1 分析图层文字标注站点断层行政边界陆地、海洋10Shapiro,200511文字标注（n）站点断层（2）地形陆地、海洋122.2 准备数据海岸线数据 （GMT自带）DEM数据 （30）颜色表CPT （GMT自带）断层经纬度、名称（自行准备）站点经纬度、名称（自行准备） 132.3 找到对应模块（module）文字标注（n）站点断层（2）地形陆地、海洋pstext psxypsxygrdgradient + grdimagepscoast142.4 组织顺序，查文档，编写脚本图层叠加顺序PostScript生成顺序 文件头 HEADER 图 层 BODY 文件尾 TRAILER15psx

5、y -K a.pspsxy -O K a.pspsxy-O a.ps打“开”文件忽略文件尾“O”verlay忽略文件头16gmt psxy %R% %J% -T -K a.psgmt psxy -O K a.psgmt psxy %R% %J% -T -O a.ps只画文件头只画文件尾gmt psxy -O K a.ps17绘图的基本要求 符合规范：字体大小，线条粗细 遵守约定俗成的规则：高速蓝色，低速红色；震中圆圈 信息突出、数据准确 清晰、清楚、美观 自己看着舒服，同行看着也舒服：颜色搭配 图: 黑白图和彩色图的分辨率不能低于600 dpi, 图中线条要清晰, 线条粗细约0.5-0.6 m

6、m. 中文图中的汉字为7 pt的幼圆字体, 英文图中的文字和阿拉伯数字为8 pt的Times New Roman字体. 图的宽度分两种: 单栏图宽4-8 cm, 通栏图宽12-17 cm. 图的长度一般不超过20 cm. 如果出现中国地图, 请注意中国地图的完整性(如国界线和南海诸岛等). 南海诸岛作插图时, 务必完整表达附图中与主图相同区域内的专题要素。房立华，201318常用投影直角坐标系：线性投影 -JX地理坐标系： 全球：-JN Robinson 投影 /-JY/-JR 半球：-JA JS 区域：美国-JB Albers等面积投影 全国100万： JA Lambert 省区图： JB

7、Albers等面积投影一些特殊投影：-JP 极坐标不知道用什么：等角圆柱墨卡托（-JM）；流行、通用 等距圆柱投影（-JQ）：适合纬度跨度大19-JQ -JN -JR-JA -JY20212.5 调试“pause” CMD Notepad+推荐阅读：Windows 下使用 GMT 的正确姿势 https:/how-to-use-gmt-under-windows.html222.6 输出gmt psconvert -Tj -C-sFONTPATH=C:windowsfonts %PS% -A+rpsconvert-E ?CONTENT1、读懂脚本2、编写自己的绘图脚本3、常用绘图模板4、常用绘

8、图技巧5、学习方法及学习资源提 纲2324gmt psbasemap -R110/111.5/30.5/31.3 -JM6.5i -B0.2 -K sx.psgmt psxy 三峡库首区_县界.txt -R -J -K -O -W0.8p,gray,. sx.psgmt psxy 三峡库首区_居民点.txt -R -J -Sc0.1 -K -O sx.psgmt pstext 三峡库首区_居民点.txt -R -J -F+f7,35+jCM -O -D0/0.2 sx.ps25gmt psxy %R% %J% -T -K -X2i Tibet.psgmt grdimage %R% %J% -B

9、a10g10 Tibet.grd -Iint.grad -CTibet.cpt -K -O Tibet.psgmt 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.psgmt psscale -D7i/2.7i/2.5i/0.15i -Ccpt1.cpt -Bx1000f200 -By+lm -K -O Tibet.psgmt psscale -D7i/1.6/1.3i/0.15i -Ccpt2.cpt

10、-L -O -K Tibet.psgmt psxy %R% %J% -T -O Tibet.ps自由空气异常图本图主要涉及命令为grdimage 、psxy、pslegend;本图自由空气异常数据来源（1分格网）其中，陆地重力异常来自EGM2008；海洋重力异常为卫星测高数据；亦可下载重力场模型（球谐系数）优先采用，直接用img2grd读取26布格异常图本图主要涉及命令为grdimage 、psxy、pslegend;布格异常数据来源：BGI（ ），World Gravity Map (参考文献：Spherical harmonic modelling to ultra-high degree

11、 of Bouguer and isostatic anomalies, J.Geod., Balmino et al., 2011.）亦可根据DEM模型，对自由空气异常进行地形改正获得布格异常。272829gawk 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% echo G 0.05i legend_china.txtech

12、o S 0.15i c 0.3 blue 0.25p 0.3i 20155(14) legend_china.txtecho S 0.15i c 0.3 red 0.25p 0.3i ;red;20152(9); legend_china.txtgmt pslegend legend_china.txt -Dx0/0/3i/0.6i/BL -F+gwhite+pblack -R -J -K -O %PS% 30REM 断层gmt psxy %R% %J% D:04数据断层Asia_faults.Lnormal -Sf-2/0.1+l+f -W1.5p,brown -Gbrown -O -K %

13、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+a-26+jLM+f11,6 -O -K %PS%推荐阅读：GMT添加断层名 https:/plot-fault-names.html重力变化图31gmt surface xyg -R -I5m -G%GRD% gmt grdimage %GRD% -Ccpt.cpt -J -R -K -O %outputps%gmt grdco

14、ntour %GRD% -R -J -C1.5 -A3 -O -K %outputps%剪应变率场3233重力变化点值时间序列gmt psbasemap -R$x1T/$x2T/$y1/$y2 -JX16c/7c -Bsxa1Y -Bpxa3of1o Bya10f2 -BSn -P -K obser.psgmt psxy obser.dat -R -JX -Sd0.25c -Ey -Gred -K -O obser.psgmt psxy mean.dat -R -J -W1.0,green,- -K -O obser.ps34重力管理部制图35gmt psvelo %D% -Sx0.01 %R

15、% %J% -O -K -W1.5p,blue %PS%echo 90 244040 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%0,2)MAP_VECTOR_SHAPE36gmt psvelo %D% -Se0.03/0.95/0 %R% %J% -G0/0/255 -h1 -K -O %PS%37psmeca -R -J -Sm0.7/0/0 -L0.5p,black -Ggreen -O -K

16、%PS%CONTENT1、读懂脚本2、编写自己的绘图脚本3、常用绘图模板4、常用绘图技巧5、学习方法及学习资源提 纲38 安装GMT4，使用 gmt psxy 安装GMT5，使用psxy 不知道自己安装的是4还是5，一条命令用gmt psxy 一条用psxygmtswitch ？39XGMT4和GMT5的兼容删除GMT5对4的兼容，直接安装两个版本，同一文件同一风格。实（TOU）用（LAN）派C:programsgmt5sharepostscriptlightPSL_custom_fonts.txtSTSong-Light-GB-EUC-H 0.7001STFangsong-Light-GB-

17、EUC-H 0.7001STHeiti-Regular-GB-EUC-H 0.7001STKaiti-Regular-GB-EUC-H 0.7001gmt pstext -L40中文支持原(XIA)理(CAI)党GMT postscript PSgsgs9.07examplescjk gscjk_ag.psX41中文支持深（CAICOU）究（BUFUZE）党X.gsgs9.07libcidfmap、cidfmap.ps42中文支持GMT postscript PSXset R=-R73/135/17/54set J=-JM6isetlocal enabledelayedexpansionfor

18、 /f %i in (dir /s/b *.DAT) do (set data=%iset ps=%i.pssurface !R! -I5m !data! -Gtmp.grd grd2cpt tmp.grd -Chaxby tmp.cptgrdimage tmp.grd !R! !J! -Ctmp.cpt -B10f5 -P !ps!ps2raster !ps! -Tj -Adel .gmt* tmp.*)变量延迟 ！能够接收变量的动态变化43脚本的批量运行 GMT的bin目录加入Matlab的搜索路径。% 创建一个100*3矩阵，xyz值均为0150之间的随机数t= rand(100,3)*

19、150% 利用GMT的surface命令对t进行格网化，输出为结构体GG = gmt(surface -R0/150/0/150 -I1, t );% 利用Grd2cpt创建颜色表文件，输出为颜色表结构体cptcpt = gmt(grd2cpt -Cjet, G);% 利用Grdimage绘制格网化结果gmt(grdimage -JX8c -Ba -P -C -G crap_img.ps, cpt, G);gmt(destroy);44二次开发与自己的工作对接45二次开发与自己的工作对接更多资源：http:/soliton.vm.bytemark.co.uk/pub/cpt-city/46颜色表的制作和修改两种方式设定颜色表:grd2cpt：根据数据范围产生；（gmt grd2cpt mydata Crainbow mycpt）makecpt：人为指定范围。（gmt makecpt -T-100/100/50 Crainbow mycpt ）47我以为我会画成这样实际我画成了这样.4849gmt makecpt -T-100/100/50 Crainbow mycpt (Yi Shuang ,2014)5051CONTENT1、读懂脚本2、编写自己的绘图脚本3、常用绘图模板4