绘制海岸线¶

GMT中使用 coast 模块绘制海岸线。 coast 模块，顾名思义是海岸线、湖岸线的意思。该模块不仅可以绘制这些岸线，还可以绘制比例尺和行政边界等。

这一节我们将通过绘制一张美国地图来介绍 coast 模块的基本用法。

绘制海岸线和湖岸线¶

下面的脚本用最简单的命令绘制了一张美国地图。其中 -R 选项设置了绘图区域，-JM15c 设置投影方式和图片尺寸， -Baf 设置底图边框属性，这些在前一节都已经做了介绍。

coast 模块的选项 -W0.5p,black 表示绘制岸线，并设置岸线为 0.5p宽的黑色线条。

gmt begin coastline png,pdf
    gmt coast -R-130/-50/20/60 -JM15c -Baf -W0.5p,black
gmt end show

美国海岸线

上图看上去有些复杂。这是因为默认情况下 -W 选项会绘制所有水体与陆地的界限，包括海岸线、湖岸线、湖中岛的边界等。根据水体或陆地的不同，coast模块将岸线分为从1到4的四个级别，依次指海岸线、湖岸线、湖中岛，以及湖中岛内的湖边界。

因而，我们可以在使用 -W 选项时指定要绘制哪一个等级的岸线。下面的脚本中，我们使用 -W1/0.5p,black 表示用 0.5p 宽的黑色线条绘制1级海岸线。

gmt begin coastline png,pdf
    gmt coast -R-130/-50/20/60 -JM15c -Baf -W1/0.5p,black
gmt end show

绘制海岸线

可以在一个命令中多次使用 -W 选项，分别用不同的线条绘制不同等级的岸线。这里我们使用 -W1/0.5p,black 绘制0.5p黑色的海岸线，同时使用 -W2/1p,lightred 绘制1p宽的浅红色湖岸线。

gmt begin coastline png,pdf
    gmt coast -R-130/-50/20/60 -JM15c -Baf -W1/0.5p,black -W2/1p,lightred
gmt end show

绘制海岸线和湖岸线

上图看上去依然有些乱。主要是因为美国有大大小小的很多岛屿和湖泊，GMT默认会把这些岛屿和湖泊都画出来。如果能够只绘制比较大的岛屿或湖泊就好了，我们可以通过 -A 选项设置要绘制的对象的最小面积（单位为平方千米）。

gmt begin coastline png,pdf
    gmt coast -R-130/-50/20/60 -JM15c -Baf -W1/0.5p,black -W2/1p,lightred -A5000
gmt end show

设置要绘制的对象的最小面积

这里我们使用 -A5000 表示只绘制面积大于5000平方千米的湖泊或岛屿。这样子得到的图看上去更加干净简洁了。

上面介绍了如何使用 -W 绘制海岸线和湖岸线。我们还可以不绘制岸线，而是为陆地和水体设置不同的填充色。

gmt begin coastline png,pdf
    gmt coast -R-130/-50/20/60 -JM15c -Baf -A5000 -Gred -Slightblue -Clightred
gmt end show

填充颜色

其中，-G 设置了陆地区域的填充色，-S 设置水体的颜色， -C 则设置湖泊的颜色（若不指定 -C，则湖泊颜色由 -S 控制）。

当然，你也可以同时使用 -W 选项和 -G 选项，即绘制岸线并填充颜色。

使用 -N 选项可以绘制国界/州界等行政边界。-N1 表示绘制国界线， -N2 表示绘制州界/省界线（目前只有美洲各国以及澳大利亚的国界的数据）。

gmt begin coastline png,pdf
    gmt coast -R-130/-50/20/60 -JM15c -Baf -A5000 -Gred -Slightblue -Clightred -N1/1p -N2/0.25p
gmt end show

绘制国界与州界

除了 -N 选项之外，-E 选项也可以用于绘制州界、国界和省界，在 DCW: 世界数字图表一节会详细介绍。本入门教程中暂时略过。

最后，我们还需要为地图添加比例尺。为了绘制比例尺，我们需要提供如下参数：

在下面的例子中，我们使用了 -Lg-60/25+c25+w1000k+f+u 增加比例尺，其中：

gmt begin coastline png,pdf
gmt coast -R-130/-50/20/60 -JM15c -Baf -A5000 -Gred -Slightblue -Clightred -Lg-60/25+c25+w1000k+f+u
gmt end show

添加比例尺