GMT常用绘图模板及技巧

GMT常用绘图模板及技巧刘珠妹中国地震局地震研究所地壳形变资料分析与应用培训班2016·09、读懂脚本2、编写自己的绘图脚本3、常用绘图模板4、常用绘图技巧5、学习方法及学习资源提纲2CONTENT1、读懂脚本2、编写自己的绘图脚本3、常用绘图模板4、常用绘图技巧5、学习方法及学习资源提纲3Windows?orLinux?温馨提示：本报告中脚本均为windows平台下使用4GMT4?orGMT5?5模块（Module）+选项（Arguments）=命令行gmtpsbasemap-R0/10/0/10-JX15c-BWSNE+tTitle-B2g2test.psgmtpsbasemap-R0/10/0/10-JX15c-BWSNE+tTitle-B2g2-Ktest.psgmtpsxydrawDif.txt-R-J-Sc0.3-W0.5p,blue-Gred-O-Ktest.psecho55测试|gmtpstext-R-J-F+f16,35,red+jLM-Otest.ps房立华，20136命令行+命令行=脚本文件尝试读一个样例脚本：数据来源：中国地震局地震研究所，GNSS中心7setPS=JP.pssetR=-R123/149/20/47setJ=-JM4igmtgmtsetFORMAT_GEO_MAPdddFgmtgmtsetFONT_ANNOT_PRIMARY7gmtpscoast%R%%J%-BWSEN-B5g5-N2-W1-A1000-Glightbrown-Slightblue-Lf128/45/35/400+lkm+jt-K%PS%gawk{print$3,$2}Japan.llh|gmtpsxy%R%%J%-Sc0.07-Gdarkred-O%PS%gmtpsconvert-A+r%PS%-Tg-P8CONTENT1、读懂脚本2、编写自己的绘图脚本3、常用绘图模板4、常用绘图技巧5、学习方法及学习资源提纲92.1分析图层文字标注站点断层行政边界陆地、海洋10Shapiro,200511文字标注（n）站点断层（2）地形陆地、海洋122.2准备数据①海岸线数据（GMT自带）②DEM数据（30″）③颜色表CPT（GMT自带）④断层经纬度、名称（自行准备）⑤站点经纬度、名称（自行准备）132.3找到对应模块（module）文字标注（n）站点断层（2）地形陆地、海洋pstextpsxypsxygrdgradient+grdimagepscoast142.4组织顺序，查文档，编写脚本图层叠加顺序PostScript生成顺序—文件头HEADER—图层BODY—文件尾TRAILER15psxy…-Ka.pspsxy…-O–Ka.pspsxy…-Oa.ps打“开”文件忽略文件尾“O”verlay忽略文件头16gmtpsxy%R%%J%-T-Ka.psgmtpsxy…-O–Ka.psgmtpsxy%R%%J%-T-Oa.ps只画文件头只画文件尾gmtpsxy…-O–Ka.ps17绘图的基本要求符合规范：字体大小，线条粗细遵守约定俗成的规则：高速蓝色，低速红色；震中圆圈信息突出、数据准确清晰、清楚、美观自己看着舒服，同行看着也舒服：颜色搭配图:黑白图和彩色图的分辨率不能低于600dpi,图中线条要清晰,线条粗细约0.5-0.6mm.中文图中的汉字为7pt的幼圆字体,英文图中的文字和阿拉伯数字为8pt的TimesNewRoman字体.图的宽度分两种:单栏图宽4-8cm,通栏图宽12-17cm.图的长度一般不超过20cm.如果出现中国地图,请注意中国地图的完整性(如国界线和南海诸岛等).南海诸岛作插图时,务必完整表达附图中与主图相同区域内的专题要素。房立华，201318常用投影直角坐标系：线性投影-JX地理坐标系：全球：-JNRobinson投影/-JY/-JR半球：-JA–JS区域：美国-JBAlbers等面积投影全国100万：–JALambert省区图：–JBAlbers等面积投影一些特殊投影：-JP极坐标不知道用什么：等角圆柱墨卡托（-JM）；流行、通用等距圆柱投影（-JQ）：适合纬度跨度大19-JQ-JN-JR-JA-JY20212.5调试“pause”→CMDNotepad++推荐阅读：Windows下使用GMT的正确姿势输出gmtpsconvert-Tj-C-sFONTPATH=C:\windows\fonts%PS%-A+rpsconvert-E?CONTENT1、读懂脚本2、编写自己的绘图脚本3、常用绘图模板4、常用绘图技巧5、学习方法及学习资源提纲2324gmtpsbasemap-R110/111.5/30.5/31.3-JM6.5i-B0.2-Ksx.psgmtpsxy三峡库首区_县界.txt-R-J-K-O-W0.8p,gray,.sx.psgmtpsxy三峡库首区_居民点.txt-R-J-Sc0.1-K-Osx.psgmtpstext三峡库首区_居民点.txt-R-J-F+f7,35+jCM-O-D0/0.2sx.ps25gmtpsxy%R%%J%-T-K-X2iTibet.psgmtgrdimage%R%%J%-Ba10g10Tibet.grd-Iint.grad-CTibet.cpt-K-OTibet.psgmtpscoast%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-OTibet.psgmtpsscale-D7i/2.7i/2.5i/0.15i-Ccpt1.cpt-Bx1000f200-By+lm-K-OTibet.psgmtpsscale-D7i/1.6/1.3i/0.15i-Ccpt2.cpt-L-O-KTibet.psgmtpsxy%R%%J%-T-OTibet.ps自由空气异常图本图主要涉及命令为grdimage、psxy、pslegend;本图自由空气异常数据来源（1分格网）其中，陆地重力异常来自EGM2008；海洋重力异常为卫星测高数据；亦可下载重力场模型（球谐系数）优先采用，直接用img2grd读取26布格异常图本图主要涉及命令为grdimage、psxy、pslegend;布格异常数据来源：BGI（），WorldGravityMap(参考文献：Sphericalharmonicmodellingtoultra-highdegreeofBouguerandisostaticanomalies,J.Geod.,Balminoetal.,2011.）亦可根据DEM模型，对自由空气异常进行地形改正获得布格异常。272829gawk{print$5,$4,0.08*$3}上周.txt|gmtpsxy-R-J-Sc-W0.5,black,-Gblue-K-O%PS%gawk{print$5,$4,0.08*$3}本周.txt|gmtpsxy-R-J-Sc-W0.5,black,-Gred-K-O%PS%echoG0.05ilegend_china.txtechoS0.15ic0.3blue0.25p0.3i20151009-20151015(14)legend_china.txtechoS0.15ic0.3red0.25p0.3i@;red;20151015-20151022(9)@;;legend_china.txtgmtpslegendlegend_china.txt-Dx0/0/3i/0.6i/BL-F+gwhite+pblack-R-J-K-O%PS%30REM断层gmtpsxy%R%%J%D:\04数据\断层\Asia_faults.Lnormal-Sf-2/0.1+l+f-W1.5p,brown-Gbrown-O-K%PS%echo100.137QinghaiLake|gmtpstext%R%%J%-F+a0+jCM+f11,4,blue-O-K%PS%echo100.736.6NanshanFault|gmtpstext%R%%J%-F+a-26+jLM+f11,6-O-K%PS%推荐阅读：GMT添加断层名重力变化图31gmtsurfacexyg-R-I5m-G%GRD%gmtgrdimage%GRD%-Ccpt.cpt-J-R-K-O%outputps%gmtgrdcontour%GRD%-R-J-C1.5-A3-O-K%outputps%剪应变率场3233重力变化点值时间序列gmtpsbasemap-R$x1\T/$x2\T/$y1/$y2-JX16c/7c-Bsxa1Y-Bpxa3of1o–Bya10f2-BSn-P-Kobser.psgmtpsxyobser.dat-R-JX-Sd0.25c-Ey-Gred-K-Oobser.psgmtpsxymean.dat-R-J-W1.0,green,--K-Oobser.ps34重力管理部制图35gmtpsvelo%D%-Sx0.01%R%%J%-O-K-W1.5p,blue%PS%echo9024404045|gmtpsvelo-Sx0.01%R%%J%-W1p,0/0/255-O-K%PS%echo90.52440nanostrain/yr|gmtpstext%R%%J%-F+f10+jLM-O%PS%[0,2)MAP_VECTOR_SHAPE36gmtpsvelo%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%PS%CONTENT1、读懂脚本2、编写自己的绘图脚本3、常用绘图模板4、常用绘图技巧5、学习方法及学习资源提纲38安装GMT4，使用gmtpsxy…安装GMT5，使用psxy…不知道自己安装的是4还是5，一条命令用gmtpsxy…一条用psxy…？39GMT4和GMT5的兼容删除GMT5对4的兼容，直接安装两个版本，同一文件同一风格。实（TOU）用（LAN）派C:\programs\gmt5\share\postscriptlight\PSL_custom_fonts.txtSTSong-Light--GB-EUC-H0.7001STFangsong-Light--GB-EUC-H0.7001STHeiti-Regular--GB-EUC-H0.7001STKaiti-Regular--GB-EUC-H0.7001gmtpstext-L40中文支持原(XIA)理(CAI)党GMTpostscriptPS…\gs\gs9.07\examples\cjk\gscjk_ag.ps41中文支持深（CAICOU）究（BUFUZE）党...\gs\gs9.07\lib\cidfmap、cidfmap.ps42中文支持GMTpostscriptPSsetR=-R73/135/17/54setJ=-JM6isetlocalenabledelayedexpansionfor/f%%iin('dir/s/b*.DAT