如何用GMT绘制GTOPO30的DEM地形图.docx

《如何用GMT绘制GTOPO30的DEM地形图.docx》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、[转载]如何用GMT绘制GTOPO30的DEM地形图 (2012-05-1719:06:57)转载▼标签： 转载多谢分享！原文地址：如何用GMT绘制GTOPO30的DEM地形图作者：forever上次提到了GTOPO30的命名法则以及下载地址，弄完，当然就要开始绘图了。一、解压缩把下载好的*.tar.gz文件放到一个文件夹里，然后在终端进入，输入ls*.tar.gz

2、gawk'{print"tarzxvf",$1;}'

3、sh，就会在当前目录下把各个压缩包里的东东解压出来。一般来说，每个包都含有8个文件，后缀分别是DE

4、M,DMW,GIF,HDR,PRJ,SCH,SRC,STX，在使用GMT绘图时只需用到DEM文件，但建议其它格式文件也都保留，另外的一些绘图工具要用到（比如MATLAB好像要用到HDR文件）。二、定义文件环境之所以这样做是因为要用到grdraster，它从DEM文件中抽取需求范围内的数据并形成grd文件，而使用该工具的前提是要对DEM数据的文件代码、别称、及数据类型进行描述，这些描述必须存储在grdraster.info文件中。假设GMT的安装目录为/usr/local/gmt，那么grdraster.info会在/

5、usr/local/gmt/share/dbase里。进入到这个目录，并把解压缩好了的DEM数据拷贝进该目录，然后编辑grdraser.info文件，按要求加入DEM文件的描述。对于中国来说，E060N90,E100N90,E060N40,E100N40这4个DEM文件就足够了。我的是这么写的：29"GTOPO30E060N90"            "m"     -R60/100/40/90          -I0.5m          Pi1           0       -9999   E060N

6、90.DEMB30"GTOPO30E060N40"            "m"     -R60/100/-10/40         -I0.5m          Pi1           0       -9999   E060N40.DEMB32"GTOPO30E100N90"            "m"     -R100/140/40/90         -I0.5m          Pi1           0       -9999   E100N90.DEMB33"GTOPO30E100

7、N40"            "m"     -R100/140/-10/40        -I0.5m          Pi1           0       -9999   E100N40.DEMB分别是文件代码、文件别称、单位、范围等等，倒数第二列是文件名，最后一个字符“B”在我的系统上必须要加，不然出的图就像是电视没信号，都是雪花（似乎见过一些朋友没加但也用得挺好的）。写完后保存退出。三、绘图假设绘一个四川的地图，范围为97/110/25/35，这个范围涉及到两个DEM数据文件——E060N40,E

8、100N40，可先用grdraster分块取出然后用grdpaste合并之。步骤是grdraster30-R97/100/25/35-I0.5m-Gsc1.grdgrdraster33-R100/110/25/35-I0.5m-Gsc2.grd30和33是grdraster.info文件中对应的文件代码，-R中所指定的范围不要超过DEM文件的范围，且二者需要有共同的边界，-I指定采样间隔，-G输出文件grdpastesc1.grdsc2.grd-Gsc.grd嗯，将两个文件合并在一起，输出为sc.grd。然后制作调色

9、板（CPT）文件，如果已经有了grd文件，则可使用grd2cpt工具来制作grd2ctpsc.grd-Ctopo-S-200/6000/200-Z>sc.cpt-C指定生产的cpt文件类型，GMT里内置了很多的，可以查看帮助文档，-S指定Z（高程）的开始和结束的范围以及间隔，-Z生产连续颜色的CPT文件。然后绘图：grdimage-JM6i-R97/110/25/35-B2/2-Csc.cptsc.grd-K-P>sc.ps-J指定投影方式，-R指定绘图范围，-B地图旁的刻度，-C调色板文件然后加入色标，用pssca

10、le来加：psscale-D6.5i/2i/7.5c/0.75c-Csc.cpt-I-E-B400-K  -P-O>>sc.ps-Dx位置/y位置/长度/宽度最后生成的图是这个样子的：