学习GMT也有几天了 ，GMT作为地学领域的一门优秀绘图语言，一直深深的吸引着我，联系了几天也绘制了一副中国区域的地形图，分辨率为5m。在此感谢李建涛（涛）的文章http://blog.sciencenet.cn/home.php?mod=space&uid=3391834，还有GMT中文社区https://gmt-china.org/#docs还有周老师的一些例子http://blog.sciencenet.cn/home.php?mod=space&uid=350278，  非常感谢刘珠妹老师，九天学者郭志道博士的热心讲解，这次我是在linux平台gmt5.4.5下绘制 代码如下，在此不做注释（虽说没有注释就是没有灵魂的代码）

#!/bin/bash

R=70/135/15/55

J=M110/16c

PS=china.ps

D=earth_relief_05m.grd

gmt set MAP_FRAME_TYPE=fancy

#gmt set MAP_GRID_PEN_PRIMARY 0.25p,black,2_2:1

#gmt set FROMAT_GEO_MAP ddd.xxF

gmt set MAP_FRAME_WIDTH 2.5p

gmt set FONT_ANNOT_PRIMARY 8p

gmt set MAP_FONT_OFFSET_PRIMARY 0.5c

gmt set FONT_TITLE 15P

gmt set FONT_TITLE_OFFSET_FRAME 2c

gmt set MAP_TICK_LENGTH_PRIMARY -2p

gmt grdcut $D -R65/140/12/65 -Gchina2.grd

gmt grdgradient china2.grd -A0 -Nt -Gint.grd

gmt psbasemap -R$R -J$J -B10g10 -BESwn+t"relief Map of China" -U"Made by Hongbin Li" -P -V -K > $PS

gmt grdimage -R$R -J$J china2.grd -Iint.grd  -Cetopo1.cpt -K -P -O >> $PS

gmt pscoast -R$R -J$J -B10g10   -Ia/0.15p,177/178/183 -I1/0.5p,61/99/172 -C81/174/254 -Lg102/12+c30+w1000+u+f -K -O >> $PS

gmt psxy CN-border-La.dat   -R$R -J$J -K -O >> $PS

gmt psxy CN-capitals.dat -J$J -R$R -Sc0.08c -Gblack -K -O >> $PS

gmt psxy CN-capitals.dat -J$J -R$R -Sc0.15c -W0.2p,black -K -O >> $PS

gmt pstext CN-capitals.dat -J$J -R$R -F+f7p,35+j -Dj0.15c/0.15c -K -O >> $PS

绘制的地形图如下：