首先到http://topex.ucsd.edu/cgi-bin/get_data.cgi下载要绘图范围的1-minute grids in ASCII XYZ-format地形数据,输入经纬度范围E117-125,N19-28,文件保存为tw.dx bat文件代码: gmtdefaults -D > .gmtdefaults4
gmtset BASEMAP_TYPE Plain
set filename=gmt2011_01a.ps
surface tw.dx -R117.0083/125.0083/19.003/28.0051 -I1m -Gtw.grd
grdgradient tw.grd -A0/45 -Gtw.grd.int -Ne1 -M
grd2cpt tw.grd -Csealand -L4000/8000 -S-8000/4000/500 > topo_sealand.cpt
grdimage tw.grd -Jm3 -R118.5/123.5/20/26.5 -Ctopo_sealand.cpt -Itw.grd.int -P -K -V -U"xuxiaoqiang"> %filename%
pscoast -R -Jm -Di -Ba2f1 -W0.8p,0/0/0 -P -O -V >> %filename%
del .gmt*
del *.int
del *.cpt
del *.grd 效果:
|