GMT语句命令总结

1. psbasemap：画底图基础命令语法构成：psbasemap -R10/70/-3/8 -JX4i/3i -B10/1 (输出或者生成新文件)plot.ps（颜色标注）-B10/1是-Ba10/a1的缩写，a为刻度标记，这里表示X轴的刻度标记为10，Y轴的刻度标记为1-Bf，刻度；-Bg网格线2.psxy：画点和线（这一条是描点）例如：图1 图1语法构成：psbasemap -R -B -J （绘制底图）psxy *.text -J（不加任何数字，默认和上行一致）-R（同理） -Sc -G -W 【注意：第二行不加-B，为了防止坐标叠加若第一行加-Xc -Yc，那么第二行不加-Xc -Yc】其中：-S选项：-Sc表示圆圈，-St表示三角形，-S表示五角星，后面的数字表示大小，如-Sc0.2-G选项：（圆圈的填充颜色）黑 0/0/0白 255/255/255蓝 0/0/255红 255/0/0或者-Gred表示红色new.text为经纬度文件-W表示圆圈边缘的粗细以及圆圈边缘的颜色2. pscoast：绘制海岸线，湖泊和国界线还是如图1所示3. gawk管道文件：用于合并命令一般是gawk“print” .text（其中列与列用，隔开）4. psmeca：用于画震源机制解的球语法形式：psmecad+文件+-J+-R+-B+-Sa或者-Sc其中：-S选项有两个-Sa和-Sc-Sa（10列）：X（经度），Y（纬度），depth（深度），strike（走向），dip（倾角），rake（滑动角），mag（震级），New X（指定的假的位置），New Y（指定的假的位置），title（标题）在这里使用-Sa，由于没有假使位置，故使用0,0代替。-Sc：X（经度），Y（纬度），depth（深度），strike1（走向），dip1（倾角），rake1（滑动角），strike2（走向），dip2（倾角），rake2（滑动角），mantissa（地震矩尾数），exponent（地震矩指数部分），New X（指定的假的位置），New Y（指定的假的位置），title（标题）如果已知scalar moment=9.56e+26，则地震矩尾数为9.56，指数为265. psvelo：绘制GPS速度场（画箭头画线） 第一二行进行缺省设置，第一行设置外框粗细，第二行设置字体大小。第三行进行psbasemap底图绘制第四行进行海岸线和国界的绘制第五行进行断层的绘制第六行进行psvelo：Psvelo基本选项及原则（用于画P轴和T轴，所以只有方向没有大小）基本选项： -Sx+输入（五列：经度、纬度、0.001（无穷小）、倾角、方位角）注意区别，画P轴是使用的是3.6（$17+90）所以后面的-Sx0.04（较小）；画T轴是使用的是3.6/（$17+90）所以后面的-Sx5（较大）（3.6*，后面-Sx数值较小；3.6/，后面-Sx数值较大）6. psxy画线注意两幅图先画断层后画地震和先画地震后画断层的区别6. Grdcut：剪切地形数据grdcut+.grd（输入文件）+-R+-G（输出文件）7. Grdgradient：梯度，画光照文件grdgradient+.grd（输入文件）+-A310（其中A为azimuth，310为光照的方位，范围为0-360）+-Ne0.8（一般为定制参数，在0.5-1之间。也可以是-Nt。其中，-Ne为拉普拉斯变换，-Nt为反正切变换）+-G（输出文件，注意其文件后缀为.）8.Grdimage：绘制彩色影像图像：grdgradient+.grd（输入文件）+-C（将彩色文件拷贝到此文件夹下，如-Crainbow.cpt。这里是将gltopo.cpt.文件放入了）+-R+-JM+-I（输入文件为grdgradient中的最后的输出文件.）(一个大于)生成第一个新图（剪切图）.cpt为调色板文件8. Grd2cpt：生成 .cpt文件Grd2cpt命令：grd2cpt+-SZ(zstart最低值（此次为高程）)/Zstop（最大值（此次为高程）/采样间隔+-S+-L（可以视为下面.cpt的最小最大值）见GMT第四天9. psscale：绘制色棒10.pstext：绘制图例10. psvelo：绘制P轴和T轴注意，此处的P轴和T轴无大小，只有方向-Sx对应文档有五列，分别为：1、2列经纬度，第三列最大主张应力值、第四列最大主压应力值，第五列压应力轴方位这里-Sx后跟的数值0.1为方向的大小粗细-W这里为设置线宽和颜色11. surface：生成.grd 文件（梯度文件）Surface的用法：1、 为什么要用surface？GMT只提供一个.grd文件，所以提供其他的.grd文件，需要自己用生成一个.grd文件。（.grd文件为进行梯度运算的文件）2、 Surface的语法规则Surface+.txt+-I（0.1c表示0.1弧秒；0.1m表示0.1弧分）+-C（固定参数）+-R+-J+-T(固定参数0-1之间)+-Ll（lower）+-Lu（upper）3、 函数的查询GMT/doc/html/man4、 由于这个给的文件stress_.txt表示Z轴方向的拉压应力，所以选取最大最小值打开文件，第三列进行选取（这个