GMT常用繪圖模板及技巧

1、GMT繪圖模板及技巧劉珠地殼形變資料分析與應(yīng)用培訓(xùn)班 2016 091CONTENT1、讀懂腳本2、編寫自己的繪圖腳本3、常用繪圖模板4、常用繪圖技巧5、學(xué)習(xí)方法及學(xué)習(xí)資源提 綱2CONTENT1、讀懂腳本2、編寫自己的繪圖腳本3、常用繪圖模板4、常用繪圖技巧5、學(xué)習(xí)方法及學(xué)習(xí)資源提 綱3Windows?orLinux?溫馨提示：本報(bào)告中腳本均為windows平臺(tái)下使用4GMT4?orGMT5?5模塊（Module）+ 選項(xiàng)（Arguments）= 命令行g(shù)mt psbasemap -R0/10/0/10 -JX15c -BWSNE+tTitle -B2g2 test.psgmt psbase

2、map -R0/10/0/10 -JX15c -BWSNE+tTitle -B2g2 -K test.psgmt psxy drawDif.txt -R -J -Sc0.3 -W0.5p,blue -Gred -O -K test.psecho 5 5 測試| gmt pstext -R -J -F+f16,35,red+jLM -O test.ps房立華，20136命令行 + 命令行 = 腳本文件嘗試讀一個(gè)樣例腳本：數(shù)據(jù)來源：中國地震局地震研究所，GNSS中心7set PS=JP.psset R=-R123/149/20/47set J=-JM4igmt gmtset FORMAT_GEO_

3、MAP dddFgmt gmtset FONT_ANNOT_PRIMARY 7gmt pscoast %R% %J% -BWSEN -B5g5 -N2 -W1 -A1000 -Glightbrown -Slightblue -Lf128/45/35/400+lkm+jt -K %PS%gawk print $3, $2 Japan.llh | gmt psxy %R% %J% -Sc0.07 -Gdarkred -O %PS%gmt psconvert -A+r %PS% -Tg -P8CONTENT1、讀懂腳本2、編寫自己的繪圖腳本3、常用繪圖模板4、常用繪圖技巧5、學(xué)習(xí)方法及學(xué)習(xí)資源提 綱

4、92.1 分析圖層文字標(biāo)注站點(diǎn)斷層行政邊界陸地、海洋10Shapiro,200511文字標(biāo)注（n）站點(diǎn)斷層（2）地形陸地、海洋122.2 準(zhǔn)備數(shù)據(jù)海岸線數(shù)據(jù) （GMT自帶）DEM數(shù)據(jù) （30）顏色表CPT （GMT自帶）斷層經(jīng)緯度、名稱（自行準(zhǔn)備）站點(diǎn)經(jīng)緯度、名稱（自行準(zhǔn)備） 132.3 找到對(duì)應(yīng)模塊（module）文字標(biāo)注（n）站點(diǎn)斷層（2）地形陸地、海洋pstext psxypsxygrdgradient + grdimagepscoast142.4 組織順序，查文檔，編寫腳本圖層疊加順序PostScript生成順序 文件頭 HEADER 圖 層 BODY 文件尾 TRAILER15psx

5、y -K a.pspsxy -O K a.pspsxy-O a.ps打“開”文件忽略文件尾“O”verlay忽略文件頭16gmt psxy %R% %J% -T -K a.psgmt psxy -O K a.psgmt psxy %R% %J% -T -O a.ps只畫文件頭只畫文件尾gmt psxy -O K a.ps17繪圖的基本要求 符合規(guī)范：字體大小，線條粗細(xì) 遵守約定俗成的規(guī)則：高速藍(lán)色，低速紅色；震中圓圈 信息突出、數(shù)據(jù)準(zhǔn)確 清晰、清楚、美觀 自己看著舒服，同行看著也舒服：顏色搭配 圖: 黑白圖和彩色圖的分辨率不能低于600 dpi, 圖中線條要清晰, 線條粗細(xì)約0.5-0.6 m

6、m. 中文圖中的漢字為7 pt的幼圓字體, 英文圖中的文字和阿拉伯?dāng)?shù)字為8 pt的Times New Roman字體. 圖的寬度分兩種: 單欄圖寬4-8 cm, 通欄圖寬12-17 cm. 圖的長度一般不超過20 cm. 如果出現(xiàn)中國地圖, 請(qǐng)注意中國地圖的完整性(如國界線和南海諸島等). 南海諸島作插圖時(shí), 務(wù)必完整表達(dá)附圖中與主圖相同區(qū)域內(nèi)的專題要素。房立華，201318常用投影直角坐標(biāo)系：線性投影 -JX地理坐標(biāo)系： 全球：-JN Robinson 投影 /-JY/-JR 半球：-JA JS 區(qū)域：美國-JB Albers等面積投影 全國100萬： JA Lambert 省區(qū)圖： JB

7、Albers等面積投影一些特殊投影：-JP 極坐標(biāo)不知道用什么：等角圓柱墨卡托（-JM）；流行、通用 等距圓柱投影（-JQ）：適合緯度跨度大19-JQ -JN -JR-JA -JY20212.5 調(diào)試“pause” CMD Notepad+推薦閱讀：Windows 下使用 GMT 的正確姿勢(shì) https:/how-to-use-gmt-under-windows.html222.6 輸出gmt psconvert -Tj -C-sFONTPATH=C:windowsfonts %PS% -A+rpsconvert-E ?CONTENT1、讀懂腳本2、編寫自己的繪圖腳本3、常用繪圖模板4、常用繪

8、圖技巧5、學(xué)習(xí)方法及學(xué)習(xí)資源提 綱2324gmt psbasemap -R110/111.5/30.5/31.3 -JM6.5i -B0.2 -K sx.psgmt psxy 三峽庫首區(qū)_縣界.txt -R -J -K -O -W0.8p,gray,. sx.psgmt psxy 三峽庫首區(qū)_居民點(diǎn).txt -R -J -Sc0.1 -K -O sx.psgmt pstext 三峽庫首區(qū)_居民點(diǎn).txt -R -J -F+f7,35+jCM -O -D0/0.2 sx.ps25gmt psxy %R% %J% -T -K -X2i Tibet.psgmt grdimage %R% %J% -B

9、a10g10 Tibet.grd -Iint.grad -CTibet.cpt -K -O Tibet.psgmt pscoast %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 -O Tibet.psgmt psscale -D7i/2.7i/2.5i/0.15i -Ccpt1.cpt -Bx1000f200 -By+lm -K -O Tibet.psgmt psscale -D7i/1.6/1.3i/0.15i -Ccpt2.cpt

10、-L -O -K Tibet.psgmt psxy %R% %J% -T -O Tibet.ps自由空氣異常圖本圖主要涉及命令為grdimage 、psxy、pslegend;本圖自由空氣異常數(shù)據(jù)來源（1分格網(wǎng)）其中，陸地重力異常來自EGM2008；海洋重力異常為衛(wèi)星測高數(shù)據(jù)；亦可下載重力場模型（球諧系數(shù)）優(yōu)先采用，直接用img2grd讀取26布格異常圖本圖主要涉及命令為grdimage 、psxy、pslegend;布格異常數(shù)據(jù)來源：BGI（ ），World Gravity Map (參考文獻(xiàn)：Spherical harmonic modelling to ultra-high degree

11、 of Bouguer and isostatic anomalies, J.Geod., Balmino et al., 2011.）亦可根據(jù)DEM模型，對(duì)自由空氣異常進(jìn)行地形改正獲得布格異常。272829gawk print $5,$4,0.08*$3 上周.txt |gmt psxy -R -J -Sc -W0.5,black, -Gblue -K -O %PS%gawk print $5,$4,0.08*$3 本周.txt |gmt psxy -R -J -Sc -W0.5,black, -Gred -K -O %PS% echo G 0.05i legend_china.txtech

12、o S 0.15i c 0.3 blue 0.25p 0.3i 20155(14) legend_china.txtecho S 0.15i c 0.3 red 0.25p 0.3i ;red;20152(9); legend_china.txtgmt pslegend legend_china.txt -Dx0/0/3i/0.6i/BL -F+gwhite+pblack -R -J -K -O %PS% 30REM 斷層gmt psxy %R% %J% D:04數(shù)據(jù)斷層Asia_faults.Lnormal -Sf-2/0.1+l+f -W1.5p,brown -Gbrown -O -K %

13、PS%echo 100.1 37 Qinghai Lake |gmt pstext %R% %J% -F+a0+jCM+f11,4,blue -O -K %PS%echo 100.7 36.6 Nanshan Fault |gmt pstext %R% %J% -F+a-26+jLM+f11,6 -O -K %PS%推薦閱讀：GMT添加斷層名 https:/plot-fault-names.html重力變化圖31gmt surface xyg -R -I5m -G%GRD% gmt grdimage %GRD% -Ccpt.cpt -J -R -K -O %outputps%gmt grdco

14、ntour %GRD% -R -J -C1.5 -A3 -O -K %outputps%剪應(yīng)變率場3233重力變化點(diǎn)值時(shí)間序列g(shù)mt psbasemap -R$x1T/$x2T/$y1/$y2 -JX16c/7c -Bsxa1Y -Bpxa3of1o Bya10f2 -BSn -P -K obser.psgmt psxy obser.dat -R -JX -Sd0.25c -Ey -Gred -K -O obser.psgmt psxy mean.dat -R -J -W1.0,green,- -K -O obser.ps34重力管理部制圖35gmt psvelo %D% -Sx0.01 %R

15、% %J% -O -K -W1.5p,blue %PS%echo 90 244040 45 |gmt psvelo -Sx0.01 %R% %J% -W1p,0/0/255 -O -K%PS%echo 90.5 24 40 nanostrain/yr|gmt pstext %R% %J% -F+f10+jLM -O %PS%0,2)MAP_VECTOR_SHAPE36gmt psvelo %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

16、%PS%CONTENT1、讀懂腳本2、編寫自己的繪圖腳本3、常用繪圖模板4、常用繪圖技巧5、學(xué)習(xí)方法及學(xué)習(xí)資源提 綱38 安裝GMT4，使用 gmt psxy 安裝GMT5，使用psxy 不知道自己安裝的是4還是5，一條命令用gmt psxy 一條用psxygmtswitch ？39XGMT4和GMT5的兼容刪除GMT5對(duì)4的兼容，直接安裝兩個(gè)版本，同一文件同一風(fēng)格。實(shí)（TOU）用（LAN）派C:programsgmt5sharepostscriptlightPSL_custom_fonts.txtSTSong-Light-GB-EUC-H 0.7001STFangsong-Light-GB-

17、EUC-H 0.7001STHeiti-Regular-GB-EUC-H 0.7001STKaiti-Regular-GB-EUC-H 0.7001gmt pstext -L40中文支持原(XIA)理(CAI)黨GMT postscript PSgsgs9.07examplescjk gscjk_ag.psX41中文支持深（CAICOU）究（BUFUZE）黨X.gsgs9.07libcidfmap、cidfmap.ps42中文支持GMT postscript PSXset R=-R73/135/17/54set J=-JM6isetlocal enabledelayedexpansionfor

18、 /f %i in (dir /s/b *.DAT) do (set data=%iset ps=%i.pssurface !R! -I5m !data! -Gtmp.grd grd2cpt tmp.grd -Chaxby tmp.cptgrdimage tmp.grd !R! !J! -Ctmp.cpt -B10f5 -P !ps!ps2raster !ps! -Tj -Adel .gmt* tmp.*)變量延遲 ！能夠接收變量的動(dòng)態(tài)變化43腳本的批量運(yùn)行 GMT的bin目錄加入Matlab的搜索路徑。% 創(chuàng)建一個(gè)100*3矩陣，xyz值均為0150之間的隨機(jī)數(shù)t= rand(100,3)*

19、150% 利用GMT的surface命令對(duì)t進(jìn)行格網(wǎng)化，輸出為結(jié)構(gòu)體GG = gmt(surface -R0/150/0/150 -I1, t );% 利用Grd2cpt創(chuàng)建顏色表文件，輸出為顏色表結(jié)構(gòu)體cptcpt = gmt(grd2cpt -Cjet, G);% 利用Grdimage繪制格網(wǎng)化結(jié)果gmt(grdimage -JX8c -Ba -P -C -G crap_img.ps, cpt, G);gmt(destroy);44二次開發(fā)與自己的工作對(duì)接45二次開發(fā)與自己的工作對(duì)接更多資源：http:/soliton.vm.bytemark.co.uk/pub/cpt-city/46顏色表的制作和修改兩種方式設(shè)定顏色表:grd2cpt：根據(jù)數(shù)據(jù)范圍產(chǎn)生；（gmt grd2cpt mydata Crainbow mycpt）makecpt：人為指定范圍。（gmt makecpt -T-100/100/50 Crainbow mycpt ）47我以為我會(huì)畫成這樣實(shí)際我畫成了這樣.4849gmt makecpt -T-100/100/50 Crainbow mycpt (Yi Shuang ,2014)5051CONTENT1、讀懂腳本2、編寫自己的繪圖腳本3、常用繪圖模板4