GMT語句命令總結(jié)

1、1. psbasemap：畫底圖基礎命令語法構(gòu)成：psbasemap -R10/70/-3/8 -JX4i/3i -B10/1 >(輸出或者生成新文件)plot.ps（顏色標注）-B10/1是-Ba10/a1的縮寫，a為刻度標記，這里表示X軸的刻度標記為10，Y軸的刻度標記為1-Bf，刻度；-Bg網(wǎng)格線2.psxy：畫點和線（這一條是描點）例如：圖1 圖1語法構(gòu)成：psbasemap -R -B -J （繪制底圖）>······psxy *.text -J（不加任何數(shù)字，默認和上行一致）-R（同理） -Sc -G -W &

2、gt;>······【注意：第二行不加-B，為了防止坐標疊加若第一行加-Xc -Yc，那么第二行不加-Xc -Yc】其中：-S選項：-Sc表示圓圈，-St表示三角形，-S表示五角星，后面的數(shù)字表示大小，如-Sc0.2-G選項：（圓圈的填充顏色）黑 0/0/0白 255/255/255藍 0/0/255紅 255/0/0或者-Gred表示紅色new.text為經(jīng)緯度文件-W表示圓圈邊緣的粗細以及圓圈邊緣的顏色2. pscoast：繪制海岸線，湖泊和國界線還是如圖1所示3. gawk管道文件：用于合并命令一般是gawk“print&

3、#183;··” ···.text（其中列與列用，隔開）4. psmeca：用于畫震源機制解的球語法形式：psmecad+文件+-J+-R+-B+-Sa或者-Sc其中：-S選項有兩個-Sa和-Sc-Sa（10列）：X（經(jīng)度），Y（緯度），depth（深度），strike（走向），dip（傾角），rake（滑動角），mag（震級），New X（指定的假的位置），New Y（指定的假的位置），title（標題）在這里使用-Sa，由于沒有假使位置，故使用0,0代替。-Sc：X（經(jīng)度），Y（緯度），depth（深度），strike1（走向），dip1（

4、傾角），rake1（滑動角），strike2（走向），dip2（傾角），rake2（滑動角），mantissa（地震矩尾數(shù)），exponent（地震矩指數(shù)部分），New X（指定的假的位置），New Y（指定的假的位置），title（標題）如果已知scalar moment=9.56e+26，則地震矩尾數(shù)為9.56，指數(shù)為265. psvelo：繪制GPS速度場（畫箭頭畫線） 第一二行進行缺省設置，第一行設置外框粗細，第二行設置字體大小。第三行進行psbasemap底圖繪制第四行進行海岸線和國界的繪制第五行進行斷層的繪制第六行進行psvelo：Psvelo基本選項及原則（用于畫P軸和T軸，所以

5、只有方向沒有大?。┗具x項： -Sx+輸入（五列：經(jīng)度、緯度、0.001（無窮?。?、傾角、方位角）注意區(qū)別，畫P軸是使用的是3.6（$17+90）所以后面的-Sx0.04（較小）；畫T軸是使用的是3.6/（$17+90）所以后面的-Sx5（較大）（3.6*···，后面-Sx數(shù)值較??；3.6/···，后面-Sx數(shù)值較大）6. psxy畫線注意兩幅圖先畫斷層后畫地震和先畫地震后畫斷層的區(qū)別6. Grdcut：剪切地形數(shù)據(jù)grdcut+.grd（輸入文件）+-R+-G（輸出文件）7. Grdgradient：梯度，畫光照文件grdgradi

6、ent+.grd（輸入文件）+-A310（其中A為azimuth，310為光照的方位，范圍為0-360）+-Ne0.8（一般為定制參數(shù)，在0.5-1之間。也可以是-Nt。其中，-Ne為拉普拉斯變換，-Nt為反正切變換）+-G（輸出文件，注意其文件后綴為.）8.Grdimage：繪制彩色影像圖像：grdgradient+.grd（輸入文件）+-C（將彩色文件拷貝到此文件夾下，如-Crainbow.cpt。這里是將gltopo.cpt.文件放入了）+-R+-JM+-I（輸入文件為grdgradient中的最后的輸出文件.）>(一個大于)生成第一個新圖（剪切圖）.c

7、pt為調(diào)色板文件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列經(jīng)緯度，第三列最大主張應力值、第四列最大主壓應力值，第五列壓應力軸方位這里-Sx后跟的數(shù)值0.1為方向的大小粗細-W這里為設置線寬和顏色11. surface：生成.grd 文件（梯度文件）Surface的用法：1、 為什么要用surface？GMT只提供一個.grd文件，所以提供其他的.grd文件，需要自己用生成一個.grd文件。（.grd文件為進行梯度運算的文件）2、 Surface的語法規(guī)則Surface+.txt+-I（0.1c表示0.1弧秒；0.1m表示0.1弧分）+-C（固定參數(shù)）+-R+-J+-T(固定參數(shù)0-1之間)+-Ll（lower）+-Lu（upper）3、 函數(shù)的查詢GMT/doc/html/man4、 由于這個給的文件stress_.txt表示Z軸方向的拉壓應力，所以選取最大最小值打開文件，第三