0%

GMT绘制中国各省DEM地图

发表于 更新于 阅读次数： Valine：

今天分享如何使用GMT绘制DEM，并使用省界进行裁剪。

DEM的获取

DEM可以从国内的中科大镜像网站下载：
https://mirrors.ustc.edu.cn/gmt/data/

为了保持更高的细节分辨率，建议下载15秒数据（约2.7G）。

目前还有更高分辨率的DEM数据，如需要可以从网络上检索到30m分辨率（1秒）的SRTM内插数据。

Code

本例子主要用到：

  • pscoast 导出省界。省界代码如-ECN.41（河南）可从GMT目录下dcw-states.txt找到。
  • psclip 使用省界裁切地图
  • grdimage 添加DEM绘图，并添加方位光照效果
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
#!/usr/bin/env bash
ps=hnan.ps

gmt gmtset FORMAT_GEO_MAP = dddF MAP_FRAME_WIDTH=2p
gmt gmtset FONT_ANNOT_PRIMARY 7p,Helvetica,black FONT_LABEL 7p,Helvetica,black 

gmt pscoast -ECN.41 -M -Df >sd.txt

gmt psclip sd.txt `gmt gmtinfo sd.txt -I0.000001` -JM2i  -K > $ps # Begin of the clip

gmt makecpt -Cdem2 -T0/2000/1  >mss.cpt
gmt grdimage -J -R C:/Users/yangleir/Documents/dtu/*.grd -Cmss.cpt   -O -K  -I+a0 -Q >>$ps

gmt psclip -C -O  -R -J >> $ps # End of the clip

gmt psconvert $ps  -A -P -Tg

部分结果