xyz2grd的功能以及GMT在windows的安裝/ Fungsi xyz2grd dan intall GMT dalam windows/Function xyz2grd and install GMT for windows

DEM具備有經緯度以及高層三個資料，如果我們從別處獲得資料後要將不同的資料格式給GMT讀必須要做轉檔的方法，例如 .XYZ 的資料格式

DEM ialah salah satu file yang memiliki 3 data, iaitu data longtitude dan latitude dan ketinggian. Kalau kita dapat data peta bukan GMT boleh dibaca, jadi kita mesti buat konversi file itu kepada format yang boleh dibaca, contohnya data format .XYZ.
 


利用Global mapper開啟的XYZ檔的DEM
然後選擇overlay control center進入conntrol center之後則選擇metadata
可看到以下的視窗

pakai software Global mapper, buka file format XYZ,file ini ialah DEM
selepas itu pilih overlay control center, masuk ke dalam control center dan pilih metadata,
sterusnya kita dapat melihat windows dibawah.
 

將Data 的Upper left ，Upper right，Lower left ，Lower right PIXEL WIDTH, PIXEL HEIGHT
抄下 最後將資訊寫入GMT程式碼內

Sila salinkan data Upper left, Upper right, Lower left dan Lower right PIXEL WIDTH, PIXEL HEIGHT, selepas itu maklumat yang salin itu tulis ke dalam skript GMT seperti dibawah.

GMT程式碼/Skript GMT
xyz2grd Taiwan_topo.xyz -GTaiwan_topo01.grd -R117/125/18/27 -I/0.0166666666666667/0.0166666666666667 -V %-I為PIXEL的大小/ -I bermakna saiz PIXEL%
pause

這樣就轉成了GMT可讀的grd檔了
Jadi sudah berjaya pernukaran format kepada file format grd, grd boleh dibaca oleh GMT.

另外則是轉貼/
Ini pautan yang macam mana install GMT dalam Windows.
GMT在window 下的安裝

地質圖和圖例/Peta geologi dan petunjuk

一般上地質圖繪製主要是數化圖層的部份以及圖層上下疊伏的關係，將圖層以及圖例資料先在別的軟體進行數化。最後再利用GMT畫圖。

Biasanya lukis sesuatu peta geologi dalam komputer berkaitan dengan buat digitize lapisan(layer) dan hubungan atas dan bawah antara lapisan-lapisan. Sebelum pakai GMT baca data peta, kita harus guna software lain buat digitize, selepas itu pakai GMT skript lukis peta geologi dan petunjuk.

GMT程式碼如下/Skript GMT

地質圖/Peta geologi
set R= -R120.25/21.8/122.8/24r
set J=-Joc121.1/23.5/21/108/2.2i
makecpt -Cglobe -T-4000/4000/100 -Z -V > usa.cpt
makecpt -Cglobe -T-8000/8000/100 -Z -V > sea.cpt
psscale -D14c/1c/3c/0.5c -Cusa1.cpt -I0.5 -B1000:m: -V -K > studyarea4.ps
psbasemap %R% %J% -B0.5f0.5nwES --ANNOT_FONT_SIZE=10 --ANNOT_FONT=4 -G241/252/255 -V -K -O >> studyarea4.ps
pscoast %R% %J% -Na -T122.25/23.7/1 --HEADER_FONT=33 --HEADER_FONT_SIZE=15 -Lf121.5/21.55/21.55/40k+u -V -O -K >> studyarea4.ps
pscoast %R% %J% -Df -W2 -G255/246/143 -Na -V -O -K >> studyarea4.ps
psxy pingtungredsoil.txt %R% %J% -W1 -Mx -G139/71/38 -V -O -K >> studyarea4.ps
psxy argillite02.txt %R% %J% -X-0.02 -W1 -Mx -G50/205/50 -V -O -K >> studyarea4.ps
psxy argilite.txt %R% %J% -W1 -Mx -G46/139/87 -V -O -K >> studyarea4.ps
psxy ms.txt %R% %J% -W1 -Mx -G154/205/54 -V -O -K >> studyarea4.ps
psxy redsoil.txt %R% %J% -W1 -Mx -G139/71/38 -V -O -K >> studyarea4.ps
psxy B.txt %R% %J% -W1 -Mx -X0.12 -Y0.0075 -G255/215/0 -V -O -K >> studyarea4.ps
psxy P3.txt %R% %J% -W1 -Mx -G255/215/0 -V -O -K >> studyarea4.ps
psxy D1.txt %R% %J% -W1 -Mx -G0/191/255 -V -O -K >> studyarea4.ps
psxy G.txt %R% %J% -W1 -Mx -G210/105/30 -V -O -K >> studyarea4.ps
psxy D3.txt %R% %J% -W1 -Mx -G0/191/255 -V -O -K >> studyarea4.ps
psxy slate.txt %R% %J% -G24/116/205 -X-0.12 -Y-0.0075 -W1 -Mx -V -O -K >> studyarea4.ps
psxy greenschist.txt %R% %J% -Mx -W1 -G200/50/155 -V -O -K >> studyarea4.ps
psxy yulibelt.txt %R% %J% -W1 -Mx -G255/131/250 -V -O -K >> studyarea4.ps
psxy marble.txt %R% %J% -G122/197/205 -W1/176/226/255 -Mx -V -O -K >> studyarea4.ps
psxy pinanshan_Wuheconglomerate.gmt %R% %J% -W1 -Mx -Gp400/62:B238/216/174 -V -O -K >> studyarea4.ps
psxy gneiss.txt %R% %J% -G255/211/155 -W1/255/69/0 -Mx -V -O -K >> studyarea4.ps
psxy lichi.txt %R% %J% -W1 -Mx -G105/139/34 -V -O -K >> studyarea4.ps
psxy northfanshuliao.txt %R% %J% -G207 -W2/122/139/139 -Mx -V -O -K >> studyarea4.ps
psxy paliwan.txt %R% %J% -W1 -Mx -Gp400/34:B147/112/216 -V -O -K >> studyarea4.ps
psxy Tuluanshan.txt %R% %J% -W1/122/139/139 -Mx -Gp450/28:Bred -V -O -K >> studyarea4.ps
psxy pinanshan_Wuheconglomerate.gmt %R% %J% -W1 -Mx -Gp400/62:B238/216/174 -V -O -K >> studyarea4.ps
psxy shuilienconglomerate.txt %R% %J% -W1/122/139/139 -Mx -Gp400/62:B238/118/0 -V -O -K >> studyarea4.ps
psxy middle_southfanshuliao.txt %R% %J% -G207 -W1/122/139/139 -Mx -V -O -K >> studyarea4.ps
psxy lichi02.txt %R% %J% -W1 -Mx -G105/139/34 -V -O -K >> studyarea4.ps
psxy exoticrock.txt %R% %J% -W1 -Mx -Gp300/62:B0/50/0 -V -O -K >> studyarea4.ps
psxy shuilien.asc %R% %J% -W5/blue -V -O -K >> studyarea4.ps
psxy xiukuluan.asc %R% %J% -W5/blue -V -O -K >> studyarea4.ps
psxy madagida.asc %R% %J% -W5/blue -V -O -K >> studyarea4.ps
psxy mukua.asc %R% %J% -W5 -V -O -K >> studyarea4.ps
psxy luye.asc %R% %J% -W5 -V -O -K >> studyarea4.ps
psxy xuhai.asc -X-0.01 %R% %J% -W5 -V -O -K >> studyarea4.ps
psxy liwu.asc %R% %J% -W5 -V -O -K >> studyarea4.ps
psxy jinlun.asc %R% %J% -W5 -V -O -K >> studyarea4.ps
psxy lele.asc %R% %J% -W5 -V -O -K >> studyarea4.ps
psxy 01.asz %R% %J% -G0 -Mx -V -O -K >> studyarea4.ps
psxy Lanshu.asz %R% %J% -W1/122/139/139 -Mx -Gp450/28:Bred -X0.08 -Y0.01 -V -O -K >> studyarea4.ps
psxy moea_fault.gmt %R% %J% -MX -X-0.05 -W2/red -V -O >> studyarea4.ps
pause

圖例/Petunjuk
set R= -R121.59/121.7/23.701/23.85
set J=-JX4/7
makecpt -Cglobe -T-3999/4000/100 -Z -V > usa2.cpt
psbasemap %R% %J% -Ba0.3f0.3nwes --ANNOT_FONT_SIZE=9 -V -K > studyarea5.ps
psxy 1.asc %R% %J% -W2 -G255/246/143 -V -O -K >> studyarea5.ps
psxy 2.asc %R% %J% -W2 -G139/71/38 -V -O -K >> studyarea5.ps
psxy 3.asc %R% %J% -W2 -Gp100/62:B238/216/174 -V -O -K >> studyarea5.ps
psxy 4.asc %R% %J% -W2 -Gp100/62:B238/118/0 -V -O -K >> studyarea5.ps
psxy 5.asc %R% %J% -W2 -Gp150/34:B147/112/216 -V -O -K >> studyarea5.ps
psxy 6asc.csv %R% %J% -W2 -G207 -V -O -K >> studyarea5.ps
psxy 7.csv %R% %J% -W2 -Gp200/28:Bred -V -O -K >> studyarea5.ps
psxy 8.csv %R% %J% -W2 -G105/139/34-V -O -K >> studyarea5.ps
psxy 9.csv %R% %J% -G154/205/54 -W2 -V -O -K >> studyarea5.ps
psxy 10.csv %R% %J% -W2 -G50/205/50 -V -O -K >> studyarea5.ps
psxy 11.csv %R% %J% -W2 -G46/139/87 -V -O -K >> studyarea5.ps
psxy 12.csv %R% %J% -W2 -G24/116/205 -V -O -K >> studyarea5.ps
psxy 13.csv %R% %J% -W2 -G255/131/250 -V -O -K >> studyarea5.ps
psxy 14.csv %R% %J% -W2 -G200/50/155 -V -O -K >> studyarea5.ps
psxy 15.csv %R% %J% -W2 -G122/197/205 -V -O -K >> studyarea5.ps
psxy 16.csv %R% %J% -W2/255/69/0 -G255/211/155 -V -O -K >> studyarea5.ps
psxy 17.csv %R% %J% -W2/255/69/0 -G255/211/155 -V -O -K >> studyarea5.ps
pstext textlegend1.asc %R% %J% -V -O >> studyarea5.ps
pause

利用NOAA(美國國家海洋及大氣觀測局)的資料畫SST(海水表面溫度)圖/Peta suhu lapisan air laut teratas pakai GMT,data dari NOAA AMERIKA SYARIKAT

網路上有很多開放的性的資料庫，我們可以善加利用，如果要畫一張海水表面溫度圖，可以選擇美國國家海洋大氣觀測局的資料庫 http://ghrsst.nodc.noaa.gov/acessdata.html 的資料，但是在讓GMT畫圖的時候，我們必須對資料進行轉換.

Dalam internet, kita boleh dapat banyak data-data dari katalog atau FTP terbuka, Kita boleh guna katalog terbuka buat peta. Kalau nak buat sesuatu tajuk peta, seperti  suhu permukaan air laut dunia.
Di sini ialah satu contoh NOAA(National Oceanographic Data Center Amerika syarikat) katalog pautan http://ghrsst.nodc.noaa.gov/acessdata.html,Sebelum pakai data yang didapat dari katalog terbuka ini, kita harus buat konversi file format.

1.下載的介面


http://ghrsst.nodc.noaa.gov/acessdata.html > HTTP: http://data.nodc.noaa.gov/ghrsst/
>L4 > Glob > Navo > K10_SST > 可以選擇年度有資料天數的SST溫度資料，點選天數之後下載檔案格式為 ******SST.nc.bz2的資料

2.將資料下載完畢後解壓縮。然後開啟Global Mapper，選擇所解壓縮的.nc檔案
會出現一個視窗，選擇ANALYSED_SST


開啟之後選擇 >Tools > control centre 就會出現以下的畫面

點選圖層，滑鼠左鍵按兩下會出現 elevation options 的視窗 選擇 alter elevation values > Valid elvation range的Minimum data將無資料的點轉換成0。

最後則將.nc檔案轉成XYZ檔
流程為File> export Raster and elevation data > export XYZ grid > options > coordinate seperator > space, 就可以轉成XYZ檔了。

3.其餘的部份則是利用GMT來畫圖
GMT的程式碼如下

minmax analysised20110414.xyz %找出XYZ的最大和最小值%
gawk "{print $1, $2 , $3*0.1}" analysised20110414.xyz > sst01.xyz %Z溫度值為給予*0.1%
makecpt -Crainbow -T0/35/5 -V -Z > temp.cpt%製作調色盤%
xyz2grd sst01.xyz -Gsst.grd -I0.1 -R-180/180/-90/90 %將xyz值轉成grd檔案格式%
grdimage sst.grd -R-180/180/-80/80 -Ctemp.cpt -Jj150/0.025i -V -K > sst01.ps %將轉好的grd檔案讀入%
pscoast -R-180/180/-80/80 -Jj150/0.025i -B30 -G0 -O -V >> sst01.ps %畫海岸線，海岸以內-G填為黑色%
pause

不同投影法的台灣,加部分河川流域/Jenis projection yang lain dan tambah basin sungai

程式碼/ Skript GMT
set R= -R120.45/22.23/122.2/24.2r
set J=-Joc121.1/23.5/21/108/2.5i
makecpt -Cglobe -T-4000/4000/100 -Z -V > usa.cpt
makecpt -Cglobe -T-8000/8000/100 -Z -V > sea.cpt
psscale -D9c/5c/5c/0.5c -Cusa1.cpt -I0.5 -B1000 -V -K > qixing_coast1.ps
psbasemap %R% %J% -B1f1nwES -G241/252/255 -V -K -O >> qixing_coast1.ps
pscoast %R% %J% -Na --HEADER_FONT=33 -T121.8/23.8/1 -Lf121.15/22.2/22.2/40k -V -O -K >> qixing_coast1.ps
pscoast %R% %J% -Df -Gc -V -K -O >> qixing_coast1.ps
grdimage "/Numerical process/GMT_DTM_data/twdtm_small.grd" %R% %J% -I"/Numerical process/GMT_DTM_data/twdtm_small_temp_shade.grd" -Cusa.cpt -V -K -O >>
qixing_coast1.ps
pscoast %R% %J% -Q -V -O -K >> qixing_coast1.ps
pscoast %R% %J% -Df -W2/red -Na -V -O -K >> qixing_coast1.ps
psxy liwu.asc %R% %J% -W5/red -V -O -K >> qixing_coast1.ps
psxy lele.asc %R% %J% -W5/red -V -O -K >> qixing_coast1.ps
psxy jinlun.asc %R% %J% -W5/red -V -O -K >> qixing_coast1.ps
psxy luye.asc %R% %J% -W5/red -V -O -K >> qixing_coast1.ps
psxy mukua.asc %R% %J% -W5/red -V -O -K >> qixing_coast1.ps
psxy xuhai.asc %R% %J% -W5/red -X-0.02 -Y-0.02 -V -O -K >> qixing_coast1.ps
psxy shuilien.asc %R% %J% -W5/blue -V -O -K >> qixing_coast1.ps
psxy xiukuluan.asc %R% %J% -W5/blue -V -O -K >> qixing_coast1.ps
psxy madagida.asc %R% %J% -W5/blue -V -O -K >> qixing_coast1.ps
pause

3D台灣加構造線/ dimensi 3 dan garisan struktur Taiwan

GMT程式碼/Skript GMT/Script GMT
makecpt -Cglobe -T-11000/8000/100 -Z -V > sea.cpt
psscale -D20c/12.75c/8c/0.5c --ANNOT_FONT_SIZE=8 -B5/:meter: -E -Csea1.cpt -I0.5 -B2000 -V -K > qixing_coast2.ps
#grdsample "/Numerical process/GMT_DTM_data/Taiwan_topo01.grd" -L -I0.02 -T -R119/124/21/26 -Ql -fg -GTaiwan_topo01_1by1.grd -V
grdgradient Taiwan_topo01_1by1.grd -G"/Numerical process/GMT_DTM_data/Taiwan_topo01_1by1_temp_shade.grd" -A60 -Ne0.9 -V
grdview Taiwan_topo01_1by1.grd -R119/124/21/26/-20000/4000 -Jm1 -Jz0.000075 -N-20000/gray -I"/Numerical process/GMT_DTM_data/Taiwan_topo01_1by1_temp_shade.grd" -E140/40 -B1/1/10000:"meter":SEwnZ -Qi -Csea1.cpt -V -K -O >> qixing_coast2.ps
pscoast -R -Jm1 -Jz0.000075 -W3 -E140/40 -Df -Tx123/25/5 -Z10 -O -V -K >>qixing_coast2.ps
psxyz deformationfront01.xyz -R -Jm1 -Jz0.000075 -W4/red -Sf0.5/0.1ilt -E140/40 -Gred -Z10 -O -V -K >> qixing_coast2.ps
psxyz deformationfront02.xyz -R -Jm1 -Jz0.000075 -W4,red,- -Sf0.5/0.1ilt -E140/40 -Gred -Z10 -O -V -K >> qixing_coast2.ps
psxyz deformationfront03.xyz -R -Jm1 -Jz0.000075 -W4/red -Sf0.5/0.1ilt -E140/40 -Gred -Z10 -O -V -K >> qixing_coast2.ps
psxyz lishanfault02.xyz -R -Jm1 -Jz0.000075 -W4/red -Sf0.5/0.1ilt -E140/40 -Gred -Z10 -O -V -K >> qixing_coast2.ps
psxyz Luzontrench.xyz -R -Jm1 -Jz0.000055 -W4,red,- -E140/40 -Z6500 -O -V -K >> qixing_coast2.ps
psxyz Ryukyutrench01.xyz -R -Jm1 -Jz0.000095 -W4,red,- -Sf0.6/0.1irt -E140/40 -Gred -Z-2750 -O -V -K >> qixing_coast2.ps
psxyz Ryukyutrench02.xyz -R -Jm1 -Jz0.000065 -W4,red -Sf0.6/0.1irt -E140/40 -Gred -Z-5600 -O -V -K >> qixing_coast2.ps
psxyz chuchifault.xyz -R -Jm1 -Jz0.000070 -W4/red -Sf0.5/0.1ilt -E140/40 -Gred -O -V -Z2100 -K >> qixing_coast2.ps
psxyz lishanfault.xyz -R -Jm1 -Jz0.000068 -W4/red -Sf0.5/0.1irt -E140/40 -Gred -Z3500 -O -V >> qixing_coast2.ps
pause

利用地震資料來畫地震分布,深度以及規模/Pakai data gempa bumi dari katalog terbuka di internet

主要是利用地震資料(可以利用美國地震資料中心 http://earthquake.usgs.gov/earthquakes/eqarchives/epic/epic_rect.php，台灣地震科學中心TEC http://tecdc.earth.sinica.edu.tw/ 下載資料),選取所需要的經緯度，所下載的資料。資料將會有不同的屬性，其中會有經緯度，規模以及深度。再將這些資料利用GMT畫出來。

Pakai data Gempa bumi dari katalog terbuka(Pautan USGS amerika syarikat http://earthquake.usgs.gov/earthquakes/eqarchives/epic/epic_rect.php，Taiwan TEC http://tecdc.earth.sinica.edu.tw/ ) ,pilih keluasan longtitude dan latitude, download data. Dalam data yang dapat itu ada properti longtitude, latititude, kekuatan gempa bumi(magnitude), dalam(depth).Pakai GMT skript boleh baca file ini.

以下是GMT程式碼/ Skript GMT

makecpt -Crainbow -T0/300/5 -V -Z > depth.cpt %先建立調色盤/ tetapkan color pallett%
psbasemap -JM6i -R114/128/3/19 -B2f2NWes -V -K >philipines.ps
%設定圖框大小，經緯度範圍，投影方法等參數/ setting keluasan peta dan jenis projection%
psxy -JM6i -R114/128/3/19 eq_nu_phi2.txt -Sc -Cdepth.cpt -: -V -K -O >> philipines.ps
%讀入所下載的地震資料檔案，-Sc 為simbol的意思，c為circle/ baca data download dari katalog, -Sc bermakna simbol, c bererti bulat%
pscoast -JM6i -R114/128/3/19 -Dh -W3 -Na -Ia -T115.5/10.25/1.25 -L115.5/8.5/8.5/100k -V -K -O >> philipines.ps
%畫出海岸線, -Dh為線條的解析度高 ，另外還有Di，Dc，Dh，Df可選擇，-W海岸線的線條粗細；後面加/─類似-W3/100為灰階色，也可以改成彩色的RGB的色階可寫成-W3/0/0/255；-T為指北，前兩個參數為經緯度，第三個為箭頭大小；-L則為比例尺前面三個為經緯度，第三個為所要設定的間隔為多少/Fungsi ini melukis garisan pantai, -Dh, Di, Dc atau Df boleh dipilih, ini ialah resolusi garisan yang hendak ditunjuk, -W bermakna garisan pantai tebal atau nipis, contoh nya -W3, belakangnya boleh tambah warna dalam jenis kelabu gred(gray grade), atau jenis RGB. -T bermakna penunjuk arah utara, nilai di petak depan dan kedua bermakna longtitude dan latitude, ketiga iaitu saiz tanda penunjuk; -L ialah skela yang ingin dipakai, petak depan dan kedua longtitude dan latitude juga, ketiga selang skela.%
psscale -D1.25c/14c/5c/0.25 -Cdepth.cpt -B50 -V -O >> philipines.ps
%表示深度的scale，將所建好的色盤讀入/skela dalam gempa bumi, baca color pallett yang setting didepan.%
pause