利用GMT的點(X,Y)如果有Z值,如何畫成contour/Use coordinate to plot the contour map by coordinate point and z value./Pakai detik koordinat melukis suatu peta contour sama data nilai Z.

繪製contour,利用blockmean, blockmode等,triangulate, nearneighbor以及 surface 等功能. Drawing contour, by using blockmean, blockmode and triangulate, nearneighbor.
Melukis Peta contour, sama dengan fungsi blockmean, blockmode, triangulate, nearneigbor dan lain-lain

 

分別由X,Y,Z的資料,Z資料可以為深度,高度等等 

XYZ data represented the coordinate (X,Y,Z)

Data XYZ bermakna detik koordinat (X,Y,Z)

 沒有加入色調的等深線圖
no color bathymetry map
Peta dalam dasar laut  tidak berwarna

加入等深線色調的圖
with color pallett
Peta yang ber warna

 

Triangulate

不連續色表
 Uncontinuous color pallett
Warna tidak berterusan

程式碼/script/skript

blockmean petroleum_game.csv -I0.1 -R0/20/0/30 -V > nf_filt.dat
#triangulate petroleum_game.csv -Gpetroleum_game1.grd -I0.2 -R0/20/0/30 -V
#nearneighbor petroleum_game.csv -Gpetroleum_game1.grd -I0.2 -N2 -R0/20/0/30 -S50 -V
surface nf_filt.dat -Gpetroleum_game1.grd -I0.1 -R0/20/0/30 -S25 -T0.2 -V
makecpt -Crainbow -T-6900/-4500/100 -Z -V > nf.cpt
psbasemap -JX4/6i -R0/20/0/30 -Y-3i -B10f1g1nwSE -V -K > Zuolan_formation.ps
grdimage petroleum_game1.grd.grd -JX4i/6i -R -Cnf.cpt -Inf_shade.grd -V -K -O >> Zuolan_formation.ps
psxy petroleum_game.csv -Sc0.05 -Gblack -R -J -V -K -O >> Zuolan_formation.ps
grdcontour petroleum_game1.grd -JX4/6i -B10f1g1nwSE -R -C50 -A200 -V -K -O >> Zuolan_formation.ps
pause

也可嘗試利用triangulate 或nearneighbor, blockmean 也可改成 blockmode, blockmedian等.
You can try another function like triangulate, nearneighbor, blockmean, blockmode, blockmedian
Anda boleh menggunakan fungsi yang lain bagi memuaskan peta yang lukis. Seperti triangulate, nearneighbor, blockmean, blockmode, blockmedian

地形圖鋸齒狀的修正和精化/smooth the topographic by low resolusion DEM/ Menghaluskan peta topographic dari resolusi rendah DEM

地形圖原始資料的解析度不好,或者是有鋸齒狀,修改的的方法則是在grdimage 中加入-Evalue,例如說-E10, -E50, -E100
Low resolution of DEM raw data, we can use grdimage -Evalue, eg. -E10, -E50, -E100 to smooth the rugged boundary.
Data DEM resolusi rendah boleh menggunakan  fungsi grdimage, tambah nilai -Evalue, misalnya. -E10, -E50, -E100 bisa menghaluskan garisan yang bengkang-bengkok.

圖1/Figure1/Gambar1

圖2/Figure2/ Gambar 2

圖3/Figure3/Gambar3

圖1 為原始讀入的DEM資料,圖二和圖三則是在grdimage中分別加入-E50和-E100的地形圖
Figure 1 represent the DEM raw data, Figure2 and Figure 3 is using the function grdimage add the -E50 and -E100 value respectively.
Gambar 1 menunjuk Data DEM asli, Gambar 2 dan Gambar 3 pakai fungsi grdimage, tambah nilai -E50 dan -E100.

程式碼/ script/ skript

makecpt -Cglobe -T-199/200/50 -Z -V > milun.cpt
psbasemap  -R121.58/121.65/23.95/24.06 -JM3i -B0.05f0.05NWes -G200/200/225 -V -K > magnifymilun.ps
grdimage "/Numerical process/GMT_DTM_data/twdtm_small.grd" -R121.58/121.65/23.95/24.06 -JM3i -E100 -Cmilun.cpt -V -K -O >> magnifymilun.ps
psxy hualian01_1.txt -R121.58/121.65/23.95/24.06 -JM3i -Wwhite -G0/200/200 -I1 -V -O -K >> magnifymilun.ps
grdcontour "/Numerical process/GMT_DTM_data/twdtm_small.grd" -R121.58/121.65/23.95/24.06 -JM3i -A20+w60 -C10 -G20c -L1/100 -Wared -V -O -K >> magnifymilun.ps
psxy meilun_river.asc -R121.58/121.65/23.95/24.06 -JM3i -Y0.1 -W10/blue -V -K -O >> magnifymilun.ps
psxy moea_fault.gmt -R121.57/121.65/23.95/24.06 -JM3i -MX -W5/red -Sf0.5/0.1irt -X0.1 -Y0.4  -Gred -V -K -O >> magnifymilun.ps
psscale -Cmilun.cpt -D5.1i/1.35i/2.88i/0.4i -Ac -B200:Topographic:/:m: -V -O >> magnifymilun.ps

台灣東部北花蓮地形/Peta permukaan darat Hualien, timur Taiwan/Contour map nothern Hualien, eastern Taiwan

修改色盤的方法-利用word pad 開啟.cpt的色盤後修改.
menukarakan warna color pallett yang di pakai, buka kan file .cpt dengan word pad.
Change the color pallett by word pad,

 根據以上的的設定,例如1公尺到5公尺範圍,設定RGB為相同的值,
mengikut contoh setting di atas, 1meter sampai 5meter ketinggian, pakai nilai RGB sama,
 According the example, the RGB value are setting from 1 meter to 5 meter high with same value.

GTM 程式碼/Skript GMT/ Script GMT

makecpt -Cglobe -T-490/500/50 -Z -V > milun.cpt  %移除-Z也可以變成非連續的色標, nilai -Z boleh dibuangkan kalau tidak mahu warna berterusan yang pakai untuk peta%
psbasemap -R121.5/121.65/23.84/24.1 -JM3i -B0.05f0.05NWes -G200/200/225 -V -K > shuilien.ps
grdimage "/Numerical process/GMT_DTM_data/twdtm_small.grd" -R121.5/121.65/23.84/24.1 -JM3i -I"/Numerical process/GMT_DTM_data/twdtm_small_temp_shade.grd" -Cmilun4.cpt -V -K -O >> shuilien.ps
#grdimage "/Numerical process/GMT_DTM_data/twdtm_small.grd" -R121.5/121.65/23.84/24.1 -JM3i -Csea1.cpt -V -K -O >> shuilien.ps
psxy hualian01_1.txt  R121.5/121.65/23.84/24.1 -JM3i -Wwhite -G0/200/200 -I1 -V -O -K >> shuilien.ps
grdcontour "/Numerical process/GMT_DTM_data/twdtm_small.grd" -R121.5/121.65/23.84/24.1 -JM3i -A20+w60 -C10 -G20c -L1/100 -V -O -K >> shuilien.ps
psxy moea_fault.gmt -R121.5/121.65/23.84/24.1 -JM3i -MX -W5/red -X0.1 -Y-0.1 -Sf0.5/0.1irt -Gred -V -O >> shuilien.ps
#pscoast -R121.5/121.65/23.85/24.06 -JM3i -Df -W2/0 -Na -V -O >> shuilien.ps
pause

%移除-Z也可以變成非連續的色表, nilai -Z boleh dibuangkan kalau tidak mahu warna berterusan yang pakai untuk peta%
makecpt -CCrainbow -T-5/55/5  -V > milun.cpt
psbasemap  -R121.58/121.65/23.95/24.06 -JM3i -B0.05f0.05NWes -G200/200/225 -V -K > magnifymilun.ps
grdimage "/Numerical process/GMT_DTM_data/twdtm_small.grd" -R121.58/121.65/23.95/24.06 -JM3i -I"/Numerical process/GMT_DTM_data/twdtm_small_temp_shade.grd" -C milun.cpt -V -K -O >> magnifymilun.ps
#grdimage "/Numerical process/GMT_DTM_data/twdtm_small.grd" -R121.58/121.65/23.95/24.06 -JM3i -Csea1.cpt -V -K -O >> magnifymilun.ps
psxy hualian01_1.txt -R121.58/121.65/23.95/24.06 -JM3i -Wwhite -G0/200/200 -I1 -V -O -K >> magnifymilun.ps
grdcontour "/Numerical process/GMT_DTM_data/twdtm_small.grd" -R121.58/121.65/23.95/24.06 -JM3i -A20+w60 -C10 -G20c -L1/100 -Wared -V -O -K >> magnifymilun.ps
psxy meilun_river.asc -R121.58/121.65/23.95/24.06 -JM3i -Y0.1 -W10/blue -V -K -O >> magnifymilun.ps
psxy moea_fault.gmt -R121.57/121.65/23.95/24.06 -JM3i -MX -W5/red -Sf0.5/0.1irt -X0.1 -Y0.4  -Gred -V -O >> magnifymilun.ps
pause