GMT常用绘图模板及技巧.pptx
- 文档编号:18811081
- 上传时间:2023-11-29
- 格式:PPTX
- 页数:56
- 大小:13.01MB
GMT常用绘图模板及技巧.pptx
《GMT常用绘图模板及技巧.pptx》由会员分享,可在线阅读,更多相关《GMT常用绘图模板及技巧.pptx(56页珍藏版)》请在冰点文库上搜索。
,GMT常用绘图模板及技巧,刘珠妹中国地震局地震研究所,地壳形变资料分析与应用培训班201609,http:
/gmt-china.org/,1,CONTENT,1、读懂脚本2、编写自己的绘图脚本3、常用绘图模板4、常用绘图技巧5、学习方法及学习资源,提纲,2,CONTENT,1、读懂脚本2、编写自己的绘图脚本3、常用绘图模板4、常用绘图技巧5、学习方法及学习资源,提纲,3,Windows?
orLinux?
温馨提示:
本报告中脚本均为windows平台下使用,4,GMT4?
orGMT5?
5,模块(Module)+选项(Arguments)=命令行,gmtpsbasemap-R0/10/0/10-JX15c-BWSNE+tTitle-B2g2test.ps,gmtpsbasemap-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,房立华,2013,6,命令行+命令行=脚本文件,尝试读一个样例脚本:
数据来源:
中国地震局地震研究所,GNSS中心,7,setPS=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%gawkprint$3,$2Japan.llh|gmtpsxy%R%J%-Sc0.07-Gdarkred-O%PS%gmtpsconvert-A+r%PS%-Tg-P,8,CONTENT,1、读懂脚本2、编写自己的绘图脚本3、常用绘图模板4、常用绘图技巧5、学习方法及学习资源,提纲,9,2.1分析图层,文字标注站点断层行政边界陆地、海洋,10,Shapiro,2005,11,文字标注(n)站点断层
(2)地形陆地、海洋,12,2.2准备数据,海岸线数据(GMT自带)DEM数据(30)颜色表CPT(GMT自带)断层经纬度、名称(自行准备)站点经纬度、名称(自行准备),13,2.3找到对应模块(module),文字标注(n)站点断层
(2)地形陆地、海洋,pstextpsxypsxygrdgradient+grdimagepscoast,14,2.4组织顺序,查文档,编写脚本,图层叠加顺序PostScript生成顺序,文件头HEADER图层BODY文件尾TRAILER,15,psxy-Ka.ps,psxy-OKa.ps,psxy-Oa.ps,打“开”文件忽略文件尾,“O”verlay忽略文件头,16,gmtpsxy%R%J%-T-Ka.ps,gmtpsxy-OKa.ps,gmtpsxy%R%J%-T-Oa.ps,只画文件头,只画文件尾,gmtpsxy-OKa.ps,17,绘图的基本要求符合规范:
字体大小,线条粗细遵守约定俗成的规则:
高速蓝色,低速红色;震中圆圈信息突出、数据准确清晰、清楚、美观自己看着舒服,同行看着也舒服:
颜色搭配图:
黑白图和彩色图的分辨率不能低于600dpi,图中线条要清晰,线条粗细约0.5-0.6mm.中文图中的汉字为7pt的幼圆字体,英文图中的文字和阿拉伯数字为8pt的TimesNewRoman字体.图的宽度分两种:
单栏图宽4-8cm,通栏图宽12-17cm.图的长度一般不超过20cm.如果出现中国地图,请注意中国地图的完整性(如国界线和南海诸岛等).南海诸岛作插图时,务必完整表达附图中与主图相同区域内的专题要素。
房立华,2013,18,常用投影直角坐标系:
线性投影-JX地理坐标系:
全球:
-JNRobinson投影/-JY/-JR半球:
-JAJS区域:
美国-JBAlbers等面积投影全国100万:
JALambert省区图:
JBAlbers等面积投影一些特殊投影:
-JP极坐标不知道用什么:
等角圆柱墨卡托(-JM);流行、通用等距圆柱投影(-JQ):
适合纬度跨度大,19,-JQ-JN-JR-JA-JY,20,21,2.5调试,“pause”CMDNotepad+,推荐阅读:
Windows下使用GMT的正确姿势https:
/seisman.info/how-to-use-gmt-under-windows.html,22,2.6输出,gmtpsconvert-Tj-C-sFONTPATH=C:
windowsfonts%PS%-A+r,psconvert,-E?
CONTENT,1、读懂脚本2、编写自己的绘图脚本3、常用绘图模板4、常用绘图技巧5、学习方法及学习资源,提纲,23,24,gmtpsbasemap-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.ps,25,gmtpsxy%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分格网)http:
/topex.ucsd.edu/marine_grav/mar_grav.html其中,陆地重力异常来自EGM2008;海洋重力异常为卫星测高数据;亦可下载重力场模型(球谐系数)http:
/icgem.gfz-potsdam.de/ICGEM/,26,布格异常图,本图主要涉及命令为grdimage、psxy、pslegend;布格异常数据来源:
BGI(http:
/bgi.omp.obs-mip.fr/),WorldGravityMap(参考文献:
Sphericalharmonicmodellingtoultra-highdegreeofBouguerandisostaticanomalies,J.Geod.,Balminoetal.,2011.)亦可根据DEM模型,对自由空气异常进行地形改正获得布格异常。
27,28,29,gawkprint$5,$4,0.08*$3上周.txt|gmtpsxy-R-J-Sc-W0.5,black,-Gblue-K-O%PS%gawkprint$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%,30,REM断层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添加断层名https:
/seisman.info/plot-fault-names.html,重力变化图,31,gmtsurfacexyg-R-I5m-G%GRD%gmtgrdimage%GRD%-Ccpt.cpt-J-R-K-O%outputps%gmtgrdcontour%GRD%-R-J-C1.5-A3-O-K%outputps%,剪应变率场,32,33,重力变化点值时间序列,gmtpsbasemap-R$x1T/$x2T/$y1/$y2-JX16c/7c-Bsxa1Y-Bpxa3of1oBya10f2-BSn-P-Kobser.psgmtpsxyobser.dat-R-JX-Sd0.25c-Ey-Gred-K-Oobser.psgmtpsxymean.dat-R-J-W1.0,green,-K-Oobser.ps,34,重力管理部制图,35,gmtpsvelo%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_SHAPE,36,gmtpsvelo%D%-Se0.03/0.95/0%R%J%-G0/0/255-h1-K-O%PS%,37,psmeca-R-J-Sm0.7/0/0-L0.5p,black-Ggreen-O-K%PS%,CONTENT,1、读懂脚本2、编写自己的绘图脚本3、常用绘图模板4、常用绘图技巧5、学习方法及学习资源,提纲,38,安装GMT4,使用gmtpsxy安装GMT5,使用psxy不知道自己安装的是4还是5,一条命令用gmtpsxy一条用psxy,https:
/seisman.info/multiple-versions-of-gmt.html,gmtswitch?
39,X,GMT4和GMT5的兼容,删除GMT5对4的兼容,直接安装两个版本,同一文件同一风格。
实(TOU)用(LAN)派,C:
programsgmt5sharepostscriptlightPSL_custom_fonts.txt,STSong-Light-GB-EUC-H0.7001STFangsong-Light-GB-EUC-H0.7001STHeiti-Regular-GB-EUC-H0.7001STKaiti-Regular-GB-EUC-H0.7001,gmtpstext-L,40,中文支持,原(XIA)理(CAI)党,GMTpostscriptPS,gsgs9.07examplescjkgscjk_ag.ps,X,41,中文支持,深(CAICOU)究(BUFUZE)党,X,.gsgs9.07libcidfmap、cidfmap.ps,42,中文支持,GMTpostscriptPS,X,setR=-R73/135/17/54setJ=-JM6isetlocalenabledelayedexpansionfor/f%iin(dir/s/b*.DAT)do(setdata=%isetps=%i.pssurface!
R!
-I5m!
data!
-Gtmp.grdgrd2cpttmp.grd-Chaxbytmp.cptgrdimagetmp.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)*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-Gcrap_img.ps,cpt,G);gmt(destroy);,44,二次开发与自己的工作对接,45,二次开发与自己的工作对接,更多资源:
http:
/soliton.vm.bytemark.co.uk/pub/cpt-city/,46,颜色表的制作和修改,两种方式设定颜色表:
grd2cpt:
根据数据范围产生;(gmtgrd2cptmydataCrainbowmycpt)makecpt:
人为指定范围。
(gmtmakecpt-T-100/100/50Crainbowmycpt),47,我以为我会画成这样,实际我画成了这样.,48,49,gmtmakecpt-T-100/100/50Crainbowmycpt,(YiShuang,2014),50,51,CONTENT,1、读懂脚本2、编写自己的绘图脚本3、常用绘图模板4、常用绘图技巧5、学习方法及学习资源,提纲,52,5、学习方法及学习资源,对照样例,修改参数。
遇到不明白的参数,查找系统自带帮助文档。
具体问题具体分析,不需要从头到尾读某本教程。
个人推荐参考网站:
http:
/gmt-china.org/GMT中文社区http:
/seisman.info/seismanhttp:
/科学网郑晖博客,脚本样例http:
/www.gl.ntu.edu.tw/geodesy/index.php?
option=com_content&task=view&id=30&Itemid=49台湾大学胡植庆教授GMT课件QQ群:
地学GMT学习218905582教程、帮助文档等:
官网百度google,53,数据资源,行政边界数据GSHHG行政边界(连续)国家测绘局基础地理信息中心提供的标准数据GMT5的DCW地形数据:
全球:
etopo5(5分)全国:
GTOPO30(30秒)/SRTM30_PLUS省市:
SRTM/ASTERGDEM海底:
GEBCO自己的专题数据:
MIF转gmt,https:
/seisman.info/introduction-to-dcw-gmt5.html,https:
/seisman.info/global-relief-models.html,54,本讲义引用了台湾大学地质科学系胡植庆教授的GMT系列课件中的图片,房立华博士常用绘图软件报告,GMT之路2.0(YiShuang)以及SeisMan博客(http:
/seisman.info/)、GMT中文社区(http:
/www.gmt-china.org/)中的内容,在此一并致谢。
谢谢!
55,参考sample1,画一个底色为黄色,线条为红色,y轴标注间隔为20,x轴格网密度为3个月的时序图。
参考sample5,将震中点改为三角形表示,颜色不变。
参考sample2,用自己的数据画一幅空间分布图。
56,
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- GMT 常用 绘图 模板 技巧