葡萄 发表于 2013-4-13 22:28:31

GMT范例一

首先到http://topex.ucsd.edu/cgi-bin/get_data.cgi下载要绘图范围的1-minute grids in ASCII XYZ-format地形数据，输入经纬度范围E117-125，N19-28，文件保存为tw.dxbat文件代码：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效果：

页: [1]
查看完整版本: GMT范例一