5-1-在GMT中使用高精度地形数据SRTM3

28 Jul 2018

在GMT中,最常用的地图就是地形图,本文要讲的就是高精度地形图如何在Windows中使用

etopo1数据成图效果展示

高精度地形图简介

  1. etopo1:精度为1弧分,约为1~2km,覆盖范围是全球(即包括南北极)
  2. etopo5:与etopo1类似,精度为5弧分,现在用的人比较少了。
  3. GEBCO,精度为30弧秒,覆盖范围是全球
  4. SRTM30和STRM15:精度分别为30弧秒和15弧秒,范围是纬度-81~81
  5. SRTM3:精度为3弧秒,大约80m,只覆盖纬度-60~60的区域,中国完全覆盖,所以研究国内区域可以使用该高精度地形数据
  6. ASTER GDEM:精度1弧秒,不太了解。可以参考:地形起伏数据 数据下载及介绍参见:GMT中文社区地形起伏数据

    SRTM3介绍

    美国航空航天局NASA在2000时利用奋进号航天飞机上的雷达测观测所得,是以前用得最多的高程数据,覆盖了全球南北纬60度以内的区域。有SRTM1和SRTM3两种,即分别是1角秒和3角秒精度的数据,对应精度为30米和90米。谷歌地球所使用高程数据即为SRTM3。公开出来的覆盖中国区域的只有90米精度。中科院镜像可以下载
    参考资料:DLR、ASTER GDEM、SRTM3、GMTED2010等5种全球高程数据对比
    问题
    SRTM数据把全球分成很多不同的块,其数据也相应打包成不同的压缩包,若要绘制一个较大区域的地形,则需要选择相应的数据包,解压和拼接。

SRTM-reCoordinate.png

拼接工具安装

max_logngitude表示感兴趣区域的最小经度值,后面三个参数类似,要用具体数值代替。outfilename要用具体文件名代替。。

  1. 本目录下可以找到outfilename.grd这个文件,这就是地形数据生成的网格文件,可以直接在GMT下绘制地形图。

更新:
发现GMT自带一种更简单下载的方式, 可以直接在命令行输入。

Some data sets are ubiquitous and used by nearly all GMT users. At the moment this set is limited to Earth relief grids. If you reference files called earth_relief_res.grd on a command line then that grid will automatically be downloaded from the GMT Data Site and placed in $GMT_USERDIR [~/.gmt]. The resolution res allows a choice among 13 command grid spacings: 60m, 30m, 20m, 15m, 10m, 6m, 5m, 4m, 3m, 2m, 1m, 30s, and 15s (with file sizes 111 kb, 376 kb, 782 kb, 1.3 Mb, 2.8 Mb, 7.5 Mb, 11 Mb, 16 Mb, 27 Mb, 58 Mb, 214 Mb, 778 Mb, and 2.6 Gb respectively). Once one of these have been downloaded any future reference will simply obtain the file from $GMT_USERDIR (except if explicitly removed by the user). Note: The four highest resolutions are the original data sets SRTM15+, SRTM30+, ETOPO1 and ETOPO2V2. Lower resolutions are spherically Gaussian-filtered versions of ETOPO1.

即可以直接用erath_relief_30s.grd作为输入文件,gmt会自动下载SRTM 3的数据,存放在~/.gmt中,~在windows下一般是c:\Users\yourname
考虑到在国内,估计网速……
另一种可行的方法,下载ASCII的SRTM3数据,就是数据比较占空间,但可以直接用xyz2grd网格化,在GMT下使用很简单,就不需要cutsrtm3-new.sh脚本了。