2011年7月3日 星期日

不同投影法的台灣,加部分河川流域/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

沒有留言:

張貼留言