✍️ 田冬冬 • 📅 2025-04-24

# 绘制区域地形图

使用 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