1987-2017年海南岛海岸线数据集

毕京鹏，张丽，宋茜茜，隋燕，温礼

1.海南省地球观测重点实验室，海南三亚 572000

2.中国科学院遥感与数字地球研究所，数字地球重点实验室，北京 100094

3.山东科技大学，测绘科学与工程学院，山东青岛 266590

4.江西理工大学，建筑与测绘工程学院，江西赣州 341000

5.中国土地勘测规划院，北京 100035

数据库（集）基本信息简介

数据库（集）名称 1987-2017年海南岛海岸线数据集数据作者 毕京鹏，张丽，宋茜茜，隋燕，温礼数据通信作者 张丽（zhangli@radi.ac.cn ）数据时间范围 1987-2017年地理区域海岸线数据集覆盖范围为海南岛，包括海南岛沿海12个市县（海口市、澄迈县、临高县、儋州市、昌江黎族自治县、东方市、乐东黎族自治县、三亚市、陵水黎族自治县、万宁市、琼海市、文昌市）。空间分辨率 30 m数据量 1.75 MB数据格式 *.shp数据服务系统网址 http://www.sciencedb.cn/dataSet/handle/669基金项目 海南省重大科技计划项目（ZDKJ2016021）；国家发展和改革委员会促进大数据发展重大工程（2016-999999-65-01-000696-01）。数据库（集）组成本数据集主要包括海南岛8期海岸线数据产品，这些数据保存为1个压缩文件（1987-2017年海南岛8期海岸线数据集.rar），每期海岸线数据单独存贮在以年份命名的文件夹里，总数据量为1.75 MB。海岸线的类型和长度信息在数据表中分别以“类型”“Length”值字段记录。

引 言

海南岛是“21世纪海上丝绸之路”的重要节点，海南自由贸易试验区建设等国家政策的实施逐渐加快了海南岛的建设与发展。海南岛独特的地理位置使得其拥有丰富的海岸线资源，海岸线的时序变迁监测对保护海岸带的生物多样性、合理开发岸线资源以及沿海生态环境的保护与可持续发展具有重要意义，同时能为政府部门加强海岸带的资源监管和治理工作提供决策支持与依据，对沿海地区的经济发展及海岸资源的管理和使用具有重要意义。本文以多时相遥感数据为基础数据，结合实地调查记录构建海岸线类型解译标志，利用遥感数字化方法得到1987-2017年共8期的海岸线数据产品，并提供数据共享服务。相关的研究成果作为海南省海岸带生态环境监测等成果的一部分，反映海南岛时序变迁历史下的生态环境发展状况。

1 数据采集和处理方法

1.1 数据来源

本海岸线产品数据集以美国地质勘探局（http://glovis.usgs.gov/）提供下载的海南岛1987年、1990年、1995年、2000年、2005年、2010年、2015年和2017年的30 m空间分辨率的多光谱影像Landsat TM/OLI为数据源，共获取无云和条带影响的40景影像。为了从无云层遮挡的遥感影像数据中准确提取海岸线信息产品，本文采用对应时相相邻一年的遥感数据作为补充和替换，所用到的遥感影像信息如表1所示。根据成像时间对Landsat影像进行拼接，得到完整覆盖海南岛的1987-2017年共8期的遥感影像图。

表1 海南岛遥感影像数据列表

序号 成像日期 卫星 传感器 轨道号16 2001-04-09 Landsat-5 TM 123-046 17 2001-04-09 Landsat-5 TM 123-047 18 2000时相2000-11-07 Landsat-5 TM 124-046 19 2000-03-28 Landsat-5 TM 124-047 20 2000-07-25 Landsat-5 TM 125-047 21 2005-04-20 Landsat-5 TM 123-046 22 2005-04-20 Landsat-5 TM 123-047 23 2005时相2005-05-13 Landsat-5 TM 124-046 24 2005-07-16 Landsat-5 TM 124-047 25 2005-08-24 Landsat-5 TM 125-047 26 2011-06-08 Landsat-5 TM 123-046 27 2011-07-26 Landsat-5 TM 123-047 28 2010时相2010-03-24 Landsat-5 TM 124-046 29 2011-02-07 Landsat-5 TM 124-047 30 2011-07-08 Landsat-5 TM 125-047 31 2015-04-16 Landsat-8 OLI 123-046 32 2015-04-16 Landsat-8 OLI 123-047 33 2015时相2015-11-17 Landsat-8 OLI 124-046 34 2015-10-16 Landsat-8 OLI 124-047 35 2015-11-24 Landsat-8 OLI 125-047 36 2017-04-21 Landsat-8 OLI 123-046 37 2017-04-21 Landsat-8 OLI 123-047 38 2017时相2016-12-05 Landsat-8 OLI 124-046 39 2017-03-11 Landsat-8 OLI 124-047 40 2017-03-02 Landsat-8 OLI 125-047

1.2 数据处理方法

不同类型的海岸线有不同的解译和判读方式，建立海岸线的界定和判读方法是研究海岸线的重要方面。本数据集生产过程中根据海岸线所处的地理环境、海岸物质组成以及海岸开发状况，同时考虑河口岸线的特殊性，确定海岸线的位置和类型，在遥感影像上依据地物的颜色、色调、纹理结构、阴影、形状、斑块形状、位置等构建不同解译判读标志[1-2]。

从海岸线类型上划分为：自然岸线和人工岸线。其中自然岸线划分为河口岸线、基岩岸线、砂质岸线、生物岸线、淤泥质岸线；人工岸线包括港口码头岸线和人工建筑岸线等人为参与的海岸线，例如按照海岸线的用途可以细分为养殖围堤岸线、盐田围堤岸线、农田围堤岸线、建设围堤岸线、港口码头岸线、交通围堤岸线、护岸和海堤岸线、丁坝岸线等[3-5]。

结合海岸线在遥感影像上的光学、纹理、下垫面、地理位置等特征形成的“色”“形”“位”解译标志，借助遥感增强方法（如图像融合、图像拉伸、主成分变换等），在ArcGIS 10.1软件中将Landsat TM/OLI影像按照最优的波段组合及图像增强方式显示，结合实地踏勘地点，建立海岸线的类型解译标志[6-7]（表 2）。

表2 海岸线类型解译标志

海岸线一级类 海岸线二级类 图像特点（TM 642假彩色波段组合） 空间分布与解译标志人工岸线 （图示为养殖围塘岸线）人工岸线是指改变原有的自然状态，完全由人工作用而成的海岸，分布于有港口、码头、鱼塘、盐田、防潮堤、防波堤、防潮闸、人工陆连岛等建筑物的地区。

海岸线可通过多种海水和陆地的交界线进行划分，如低潮线、高潮线（平均高潮线、多年大潮高潮线等）、平均海平面与陆地交界线、瞬时水边线、痕迹线或干湿分界线等[8-9]作为海岸线的指示边界线。在不同的国家和地区之间、不同的学科分支和不同的研究目之间又有各自的定义和使用方法[10]。我国相关规范和国内外学者研究多以多年平均大潮高潮线作为海岸线的界定标准[9,11-14]，同时受限于遥感影像资料以及潮汐洋流等影响，准确划定海岸线有一定的困难，因此往往采用各种技术手段与方法提取指示海岸线的基础上，通过海岸线类型构造解译标志、依数据源和地区进行潮位校正、海岸线特征点修正和重新位置定位以及结合多种海岸线精度指标和标准评估后得到海岸线的理论准确位置。本文对覆盖海南岛的6种海岸线类型以遥感影像解译的标志为基础，对其位置与类型判别标准进行了界定[6-7]（表3）。

表3 海南岛海岸线类型位置界定标准

海岸线一级类 海岸线二级类 海岸线位置与类型界定标准淤泥质岸线已开发的淤泥质海岸常用于养殖，在影像上近海一侧修筑防浪堤坝，海岸线位置较稳定。自然淤泥质海岸的岸滩面积较大，靠陆地一侧生长有耐盐植物，在TM 432（OLI 543）波段组合的假彩色影像上呈暗红色调，向陆一侧植被与向海一侧岸滩的分界线即为淤泥质岸线。人工岸线在TM 321（OLI 432）波段组合的真彩色影像上，港口码头内部有许多工厂和居民区，建筑物多为灰白色，且有清晰的条带状道路，码头外部形状规则，海水不能淹没其外缘线，因此建筑物的外边缘线即为人工岸线；防潮堤、防波堤等是修建于沿海滩涂、河流入海口处的堤坝，用以抵挡海水入侵和风暴潮，堤外是淤泥质滩涂，堤内可能是未开发的淤泥质滩涂或者吹填沙体，因此堤坝的外围即为人工岸线所在处；养殖区和盐田蒸发池呈深蓝色调，盐田结晶池呈亮白色调，养殖塘和盐田的外边缘线即为人工岸线；此外海南岛特殊的人工岛、围填海工程建设，在形成一定的规模后，陆连岛外边缘拓展了海岸线的曲折度和利用长度，因此取其外边缘作为人工岸线。

在获取海岸线的遥感影像解译标志和界定标准后，以2017年的Landsat OLI遥感影像数据为基础，采用自动化方法获取2017年的海岸线位置和类型信息，并通过精度验证、位置纠正、类型纠正等处理后，得到精确的2017年的海岸线数据（图1）。以2017年的海岸线数据为本底数据，叠加其余时相的遥感影像，通过目视解译与判读，采用数字矢量化方法修改获取对应时相（1987年、1990年、1995年、2000年、2005年、2010年、2015年）的海岸线的位置与类型信息。

其中2017时相的海岸线提取方法如下：基于遥感和GIS相结合的技术，首先运用改进的归一化差异水体指数MNDWI对预处理后的遥感影像进行海陆分割，设置全局阈值，实现图像的二值化；其次对二值化影像进行灰度值填充来去除较小的噪声点，在形态学处理的基础上，利用MATLAB语言实现Canny算子对二值图像的海岸线提取[6]；这样得到的海岸线实际只是水边线，最后根据海图、地形图和潮汐资料等数据对水边线进行校正和修正处理，获取海南岛2017时相的海岸线数据。

图1 海岸线产品精度验证方法

2 数据样本描述

2.1 数据组成

本数据集包括海南岛1987-2017年共计8期海岸线位置分布及类型数据，这些数据保存为1个压缩文件（“1987-2017年海南岛8期海岸线数据集.rar”），每期数据单独存贮在以年份命名的文件夹里，总数据量为1.75 MB。每个文件夹中单独存贮对应的shp矢量数据文件，矢量数据采用WGS 84-UTM坐标投影系统。每个以年份为名字的文件夹中包含有2份数据，其中一份是包含位置信息、长度信息、类型信息等全部信息属性，例如命名为“2010coastline_class.shp”的海岸线数据；另外一份是只包含位置信息、长度信息，例如命名为“2010coastline.shp”的海岸线数据，可供不同领域的研究需要。

2.2 数据样本

以2017年的海岸线的提取结果为基准，在ArcGIS软件中固定显示比例尺，通过1.2章节的海岸线类型解译标志和界定方法，获取海南岛1987年、1990年、1995年、2000年、2005年、2010年和2015年海岸线位置及类型信息。海南岛海岸线的分布特征如图2所示。

图2 1987-2017年海南岛海岸线分布

3 数据质量控制和评估

3.1 数据结果与精度验证技术路线

2017年海岸线数据产品采用自动化提取方法结合数字化解译提取得到。通过借助ArcGIS软件中的DSAS（Digital Shoreline Analysis System，数字岸线分析系统）扩展插件，将提取结果与目视解译提取结果进行定量比较分析，对海岸线位置出现错误偏差的地方进行修正。此外，在2018年开展海岸线野外环境调研（http://www.radi.cas.cn/dtxw/kjyw/201802/t20180201_4941129.html，中国科学院遥感与数字地球研究所，数字地球重点实验室），对2017时相海岸线数据进行海岸线类型验证（图1）。

其余时相的海岸线（1987年、1990年、1995年、2000年、2005年、2010年、2015年）均以2017年海岸线为本底数据，结合目标时相的遥感影像，对海岸线的位置与类型信息进行修改完善得到[14-15]（图 2）。

3.2 数据结果验证方法

海岸线的类型解译及验证方法采用实地考察验证，范围主要包括海南岛沿海岸线变迁重点区域，其中包括海口市（如东寨港红树林自然保护区）、儋州市（如海花岛人工岛建设）、三亚市（如大东海、小东海以及西玳瑁岛珊瑚礁自然保护区）、文昌市（如会文区近海养殖区沿线）4个主要地区，海岸线类型的确定和历史变迁状况采用实地勘探采样、与当地居民交谈问询、相关机构情报交流等方式获取。

测量坐标点位采用手持GPS仪（型号为Trimble JUNO SB）确定点位信息，海岸线记录表采用我国近海海洋综合调查与评价专项（简称908专项）的《岸线测量登记表》[16]（表4）。

表4 海岸线测量记录表

3.3 数据结果验证精度

选取2017时相文昌地区的自动化海岸线提取结果与真实海岸线的结果进行对比，通过DSAS基线法对海岸线位置误差统计分析。通过设置约 2 000条有序断面线，对比统计分析海岸线的位置偏差（图3）。验证结果表明，海岸线位置的绝对位置误差在像元级别的缓冲区内浮动，实现了95%的区域的海岸线自动化高精度提取。最后针对误差存在的复杂河口位置、淤泥质岸线处以及近海养殖区等，采用人工修正的方法，获取海岸线的准确位置。

图3 海岸线位置验证方法

海南岛环岛野外实地调查共获取82个点位的海岸线类型信息，共拍摄380张典型海岸线照片。通过野外调研，以沿岸红树林为代表的生物岸线、人工修筑的人工岸线以及基岩岸线、砂质岸线等6类海岸线类型均得到了实地的考察验证。

4 数据使用方法和建议

1987-2017年海南岛8期海岸线数据集均为shp格式，可利用ArcGIS等地理信息系统软件对本数据集进行编辑、查询、属性显示以及后续分析工作。在海南岛海岸线研究方面，大多数学者选择局部海岸线变化较大的地区进行典型性分析，而针对海南岛全岛的海岸线变迁研究较少，本文可作为海南岛海岸线生态环境监测的基础数据集，应用于海岸线时空分布特征、陆海格局演变等研究。