GMT利用gtopo30绘制陕西省地形图
一、准备工作(参照http://xpku.blog.163.com/blog/static/239650020102224449595/)1、解压缩把下载好的*.tar.gz文件放到一个文件夹里，用winrar解压开。一般来说，每个包都含有8个文件，后缀分别是DEM, DMW, GIF, HDR, PRJ, SCH, SRC, STX，在使用GMT绘图时只需用到DEM文件，但建议其它格式文件也都保留，另外的一些绘图工具要用到（比如MATLAB好像要用到HDR文件）。
2、定义文件环境
之所以这样做是因为要用到grdraster，它从DEM文件中抽取需求范围内的数据并形成grd文件，而使用该工具的前提是要对DEM数据的文件代码、 别称、及数据类型进行描述，这些描述必须存储在grdraster.info文件中。假设GMT的安装目录为/usr/local/gmt，那么 grdraster.info会在/usr/local/gmt/share/dbase里。进入到这个目录，并把解压缩好了的DEM数据拷贝进该目录， 然后编辑grdraser.info文件，按要求加入DEM文件的描述。对于中国来说，E060N90, E100N90, E060N40, E100N40这4个DEM文件就足够了。我的是这么写的：
29 "GTOPO30 E060N90" "m" -R60/100/40/90 -I0.5m P i 1 0 -9999 E060N90.DEM B
30 "GTOPO30 E060N40" "m" -R60/100/-10/40 -I0.5m P i 1 0 -9999 E060N40.DEM B
32 "GTOPO30 E100N90" "m" -R100/140/40/90 -I0.5m P i 1 0 -9999 E100N90.DEM B
33 "GTOPO30 E100N40" "m" -R100/140/-10/40 -I0.5m P i 1 0 -9999 E100N40.DEM B
分别是文件代码、文件别称、单位、范围等等，倒数第二列是文件名，最后一个字符“B”在我的系统上必须要加二、bat文件代码：set filegrd=shaaxi.grd
set range=-R105/112/31/40
REM 陕西的范围只涉及一个DEM文件（E100N40），grdraster-extract subregion from a binary raster and write a grid filegrdraster 33 %range% -I0.5m -G%filegrd%
grd2cpt %filegrd% -Ctopo -S50/3700/100 -Z > sx.cpt
grdgradient %filegrd% -A0 -Nt -Gtmp.grd
grdimage %filegrd% -JM10c %range% -Itmp.grd -B2/2 -Csx.cpt -K -P > sx.ps
psscale -D12c/5c/7c/0.6ch -Csx.cpt -I -E -B400 -K -P -O >> sx.ps
gawk '{print $1,$2}' city.dm | psxy -JM10c %range% -K-O -Sc0.15c -W2 -V >> sx.ps
del *.cpt
del *.gmt*
del *.grd三、enjoy成果图原图效果：
能不能把goto中国的高程数据发给我一下？谢啦 看起来不错
页:
[1]