2012年7月18日 星期三

台灣東部北花蓮地形/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