[SMS]获取gshhs岸线数据

大家好！第2期，和大家一起交流一下如何获取和生成岸线数据。

最近半年，组里陆陆续续来过几个短期交流的人，有从缅因州开车来的，也有专程从日本坐飞机来的，都是想来学习FVCOM的。结果，老师们很忙，所有事情就落在了我身上。起初，我以为是关于模拟结果方面的讨论，可是，当真正见面后我才知道，我需要向来宾从零开始讲解FVCOM的使用。结果就是，短短的一周过后，我们才刚刚完成FVCOM网格的构建。几次类似经历之后，我觉得有必要把这部分内容总结一下，希望能为刚开始接触FVCOM的初学者节约时间。

如果你想使用FVCOM模拟某海域，第一件事就是构建三角形网格。而构建三角形网格的第一步是生成岸线。就像FVCOM manual里写的一样，通常，我们使用SMS软件构建三角形网格。SMS是由AQUAVEO公司制作的一款水动力模型软件，功能比较多。我只使用过其中的构建三角形网格和插值，所以在这里不展开介绍了。下载网址是

https://www.aquaveo.com/。

接下来是准备岸线数据。岸线数据的获取方式有两种。

第一种是手动勾画，即将Google Earth图片或者其他岸线图片导入SMS，然后用鼠标人为的勾勒出岸线。第二种方法是导入岸线数据，即从网站上下载岸线数据，然后改写成SMS可读的格式（比如*.cst）并导入。多数情况下，我们会使用第二种方法，因为这种方法高效快速，同时可以避免出现投影方面的错误。然而，有时我们不得不选择手动勾画的方式。比如，我们获取的岸线数据与实际情况不符，或者，项目是为了模拟岸线有所改变之后的情况。

我最常用的一个数据库是GSHHG（A Global Self-consistent, Hierarchial, High-resolution Geography Database），该数据库是由夏威夷大学和NOAA共同开发和维护的全球地理信息数据集。数据下载网址是

http://www.soest.hawaii.edu/wessel/gshhg/。

其中，岸线数据分为以下五个分辨率（摘自www.52ocean.cn）：

对于以上五个分辨率，使用过m_map的朋友一定不会陌生。没错，m_map中的高分辨率岸线也是使用这个数据集。获取数据之后，可以将研究区域的岸线提取出来，然后写成cst格式。如下：

以渤海黄海为例，我们设定经纬度范围为117-127E、33-41N，使用gshhs_h岸线数据，便可得到下图岸线：

图中蓝线为最长、最主要的亚洲大陆岸线，红线表示各个小岛屿。根据项目目的，我们可以选择性的将部分岛屿删除。写成cst文件后，我们便可以使用SMS打开岸线文件，效果如下图：

至此，我们构建FVCOM三角形网格的第一步就完成了。