2010年2月1日 星期一

利用GMT -pscoast clip的方法/Satu fungsi dalam GMT- pscoast clip


主要是利用clip 的方法 , 將兩張圖疊在一起 , 以上為例台灣的底圖較世界的底圖解析度較高
故利用此方法做疊合 , 如果遇到的問題是解析度較好的台灣底圖沒有海底的資料 , 這樣海洋的部份則以調色盤的-1 - 1 之間變化範圍內的-1為起始點 ,如果在相同色盤的情況下,所出現的情況是解析度較高的台灣島較低處和海底地形較低處的顏色相似 , 這似乎在顯示地形高低的情況是不可能發生的 , 故做兩個色表再加以套疊.

Fungsi ini ialah melakukan cara clip, terutamanya pakai dua lapisan yang berkaitan atas bawah, diatas itu resolusi tinggi(resolution),dan dibawah itu lebih rendah tetapi resolusi tinggi lapisan itu data yang meliputi itu kecil daripada resolusi yang rendah itu. Kita mungkin hanya dapat permukaan dasar laut pixel saiz besar.

以下是程式碼/Skript GMT

set R=-R118/127/17/26 %設定範圍/ masukan keluasan%
set J=-JM4i %設定投影方法/ masukan projection yang mahu di pakai%
makecpt -Cpanoply -T-5500/4000/100 -Z -V > usa.cpt %設定色盤/ pilih color pallett %
makecpt -Crainbow -T-8000/2000/100 -Z -V > sea.cpt %設定色盤/ pilih color pallett %
psbasemap %R% %J% -B2f2nwES -V -K > qixing_coast1.ps
grdimage "/Numerical process/GMT_DTM_data/etopo2.grd" %R% %J% -I"/Numerical process/GMT_DTM_data/etopo2_temp_shade.grd" -Csea.cpt -V -K -O >> qixing_coast1.ps
grdcontour "/Numerical process/GMT_DTM_data/etopo2.grd" %R% %J% -A1000+w60+s4+ap -G10c -L-11000/-1000 -Wared -V -K -O >> qixing_coast1.ps %此為加入等高線的指令/ Fungsi ini tambah garisan kontor(contour)%
pscoast %R% %J% -Df -Gc -V -K -O >> qixing_coast1.ps %加入-Gc 之後 , 就可以進行剪圖的動作/ Fungsi ini ialah fungsi Clip, tambah simbol bendera(flag) -Gc dapat melakukan fungsi klip %
grdimage "/Numerical process/GMT_DTM_data/etopo2.grd" %R% %J% -I"/Numerical process/GMT_DTM_data/etopo2_temp_shade.grd" -Cusa.cpt -V -K -O >> 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% -Df -W2 -Na -V -O >> qixing_coast1.ps
pause

沒有留言:

張貼留言