资源描述：
《gmt绘制gps速度场的应用_占伟》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、第28卷第3期华北地震科学Vol.28,No.32010年9月NORTHCHINAEARTHQUAKESCIENCESSep.,2010文章编号:1003-1375(2010)03-0061-04GMT绘制GPS速度场的应用占伟,孟宪刚,刘志广(中国地震局第一监测中心,天津300180)摘要:GMT是一款功能强大、开源的绘图软件。介绍了应用GMT绘制GPS速度场的方法,以及如何实现读入MapInfo格式数据。实践证明,上述方法具有较强的应用价值,适用于多信息的图形表达,能绘制高质量的图形。关键词:GMT;GPS;速度场;数据格式中图分类号:P321.72文献标
2、识码:A装使用,用户可以到GMT主页(http://gmt.soest.0引言hawaii.edu/)下载安装包。GMT(GenericMappingTools)是在美国国家在网络的支持下,Ubuntu8.04系统可以通过科学基金会的资助下,于1987年由PaulWessel和aptgetinstallgmt命令实现GMT的自动安装,软WalterSmithH.F.Smith共同开发的一款免费的、件安装完毕后需进行系统配置,详细的方法在参考开源的绘图软件。该软件支持30多种地图投影参文献[3]有介绍,故本文不再重复。数,能在几乎所有的操作系统下安装使用。GMT1.2GPS速
3、度场的分步绘制可以由命令行、脚本和用户程序调用,具有较强的适GMT绘图命令较多,该软件提供了较为详尽[4]应性。GMT具有很强的灵活性,能够结合地形数的使用手册。参考使用手册,笔者总结得出应用据、配色方案,将复杂的信息以高质量的ps(PostGMT绘制GPS速度场的主要步骤:Script)格式输出,绘制出非常漂亮的2D和3D图(1)绘制国界、海岸线形。GMT已在地理、大气、海洋、地震等研究领域pscoastR75/135/18/55Jm0.4cB10g0WSen[1-2]得到较为广泛的应用。目前,GMT的最新版本W1N1DiVKX1.5Y1.5>$psfil
4、e为2010年1月15日发布的4.5.2版本。主要的参数有J,投影方法,用户可根据绘图需与目前常用的绘图软件MapInfo相比,GMT要自主选择;R,地图的经纬度范围,按西/东/南/操作更为简便,绘制的图形质量更高,在绘制GPS北赋给数值;B,设定图边界尺标的间隔,a表示文速度场的应用中具有较强的可选性。但由于GMT字批注的间隔,f表示边框标记的间隔,g表示隔线绘图命令、参数选项、步骤众多且使用文本形式输入的间隔,如果要针对X轴或Y轴做不同的设定,就命令,对于初学者而言,顺利地绘制一幅完整的在参数间加斜线/;K表示后面还有绘图命令;GPS速度场还是具有一定难度。因此,笔者
5、在$psfile代表输出图形文件名。Ubuntu8.04系统安装了GMT4.2版,用于GPS(2)绘制地形图速度场制图,并加入了MapInfo数据。本文结合实grdimage$grdfileRJmBC$cptfilePK例详细介绍GMT绘制GPS速度场的使用方法和X1.5Y1.5>$psfile处理流程。$grdfile代表地形数据文件。本例使用的地形数据为GAMIT/GLOBK软件包内的55的1GPS速度场的绘制全球地形数据epoto5.grd。读者可以到美国国家地1.1GMT安装和配置球物理数据中心(NGDC)网站(http://www.ngdc.GM
6、T软件可在Linux、Windows等系统下安noaa.gov/mgg/global/)获取更高分辨率的地形数收稿日期:20100407基金项目:中国地震局第一监测中心人才培养基金项目专项资助作者简介:占伟(1983-),男,江西金溪人,助理工程师,从事GPS应用技术研究.Email:zw000373@163.com62华北地震科学28卷据;C后接调色文件(cpt格式);O,意为把这一行指令绘制的图加到前面的图上。(3)绘制点位速度psvelo$velfileA0.04/0.12/0.06H0RJmD2W1/0/0/0G0/0/0Se
7、0.05/0.95/8KPOV>>$psfile$velfile代表GPS点位速度文件,文件格式由S选项决定。在经过GAMIT/GLOBK平差计算得到速度场结果后,可使用sh_exglk提取平差结果org文件得到GPS点位速度文件;A,设置箭头属性;Se,速度尺寸/置信区间/字体大小,如字体大小为零,则不显示点名。该选项要求GPS点位文件每行格式必须相同,且每行都分为8列,格式为!XYVxVySigxSigyCorxyName∀,即表示!经度纬度E方向速度N方向速度
