绘制区域地形图

示例贡献者:

田冬冬（作者）

---

使用 GMT 内置的行政区划边界裁剪出区域并绘制地形图。

警告

GMT内置数据提供的中国国界数据不符合中国的领土主张，在正式刊物中发表使用此类数据的图件时 都可能存在问题。

#!/usr/bin/env bash
#
# Earth relief in Anhui, China
#

gmt begin ex013
    # Clip the Anhui province
    gmt coast -ECN.AH+c -RCN.AH -JM10c
    # Plot the earth relief
    gmt grdimage @earth_relief_15s -I+d
    # Add rivers and lakes
    gmt coast -Ia/1p,lightblue -Slightblue
    # Undo the clipping
    gmt clip -C
gmt end show