海洋论坛▏面向对象的无人机遥感影像海岸线提取方法研究
海岸线是海岸带的重要标志之一，也是海洋开发战略的基准线和“起跑线”，它不仅标识了陆地与海洋的分界线，而且还蕴含了丰富的环境信息，并对沿海的滩涂利用、湿地生态系统的兴衰等具有重要的指示作用。在我国，海岸线系指多年大潮平均高潮位时海陆分界痕迹线，是重要的基础地理信息要素，也是一种特殊的自然资源，快速而又准确地获取海岸线的位置及其动态变化信息，对于国土空间规划、海岸带资源可持续利用与综合管理以及生态环境保护与修复等工作具有十分重要的意义。
对海岸线探测和获取的方法通常是通过实地测量或卫星遥感影像解译。由于部分区域难以到达，人工现场测量困难，越来越多的专家学者开始利用遥感手段快速提取海岸线，自动解译成为海岸线提取技术发展的主流，各种算法不断呈现。由于海岸带区域地形复杂和卫星遥感影像分辨率的不足，常规方法存在一定的局限性，随着无人机遥感系统的普及，成为低空遥感应用领域的一个崭新发展方向。面向对象思想也随着无人机遥感的广泛应用得到的进一步发展，面向对象的方法将形状、光谱、纹理、语义和空间结构等相同特征的像元划分为同一类对象，其分析过程更加接近人类的认知过程，不仅可以降低光谱特征变化对精度的影响，而且还能够充分地利用像元之间丰富的细节信息，实现遥感图像分类和目标地物提取，信息提取精度更有保证。
因此，本研究以无人机(UAV)遥感影像为基础数据，提出一种面向对象的海岸线提取方法，开展人工海岸线和砂质海岸线识别的应用实验，并结合现场实测数据，对实验结果进行精度分析与评定。
一、研究方法
⒈海岸线识别方法
根据《全国海岸线技术规程》(自然资办函[2019]1187号)，海岸线分为人工海岸线、自然海岸线和其他海岸线。人工海岸线识别通常以海岸工程等的外沿为其位置线。自然海岸线又包括砂质岸线、泥质岸线、基岩岸线和生物岸线。对砂质岸线来说，发育有滩脊的，海岸线一般确定在滩脊的顶部向海一侧，而具陡崖的海滩通常与基岩陡岸相接，其位置界定为崖下滩与崖的相交线。泥质岸线应根据海岸植被生长变化状况、大潮平均高潮位时的海水痕迹线以及植物碎屑、贝壳碎片、杂物垃圾分布的痕迹线等综合分析界定。基岩海岸线位置界定在陡崖的基部。生物海岸线界定为毗邻或穿越珊瑚礁、红树林和海草床的海岸线。
⒉海岸线提取方法
在分析国内外相关研究方法、成果和经验的基础上，提出了一种面向对象的海岸线提取算法，该方法只需要提供少量的标记信息，而无需对海陆进行参数化分析。首先，利用多尺度分割算法对无人机遥感影像进行初始分割；然后，进行人工标记，在此基础上利用直方图为特征度量进行对象相似度描述；最后，通过对象合并算法将水陆分离，提取目标海岸线。
⑴图像分割
借助eCognition软件的图像多尺度分割算法对无人机遥感影像边界控制进行初始分割，得到均匀同质的对象，形成影像对象层次网络。
⑵目标标记
对分割的图像，通过专家目视判断，人工标记海岸和海域，为对象相似性测度提供背景参考。如图1a和图2a所示。
养殖池为海域部分
图1 IM1的人工岸线的手动标记(a)
和半自动提取结果(b)
⑶对象相似度描述
鉴于光谱直方图能够描述蕴含丰富特征信息的对象全局特征，因而，通过构造光谱直方图来描述对象特征。首先将每个光谱向量量化为r个箱格，把量化后的光谱矢量转化为反映光谱矢量信息的标量索引I。然后计算每个对象的直方图HrQ：
其中V表示对象标号，r表示箱格号，L表示对象V中像元的位置，I1，I2为箱格r的上下界，I为像元的光谱索引值，δ为狄拉克函数。利用对象光谱直方图，选择巴氏系数表征区域对象C和D的相似度：
其中HrC和HrD表示对象C和D的光谱直方图，两个区域的光谱相似度越大，直方图越相似，巴氏系数就越高，B(C，D)能够充分描述对象之间的相似性。
⑷对象合并
利用具有全局属性的极值优化对象合并方法，基于测量两个区域对象边界上像素之间的差异性进行合并。
其中F(P，Q)是最小权边缘，P和Q代表区域分量；m，n是对应于图像元素的一组节点，W是连接相邻节点对的边缘集合，每个边缘(m，n)有一个对应的权重来度量由该边连接的两个节点的相异性U(m，n)。
通过循环迭代合并最相似的对象，并把未标识对象逐步合并到已标识对象中，将水陆分离，实现目标的提取，相应的合并判定依据满足下列条件：
其中，GXs表示对象X的任一相邻对象。如果Y是X的最相似对象，那么Y＝Y∪X；反之，X与Y对象合并取消。
二、结果与分析
⒈人工海岸线
人工海海岸线通常受潮汐作用影响较小，仅需提供目视判断信息，对此，只需要大致的标记海域和陆地的位置。本研究采用了两幅具有人工海岸线的无人机图像IM1，IM2，进行试验。IM1为养殖池区域人工海岸线，IM2为码头人工海海岸线，大小为1581像素×1588像素，两幅图像分辨率均为0.05m。
由于人工海岸线构成复杂，根据不同的应用需求，本研究提取了两种海岸线结果，一种是大范围的海岸线获取，另外一种为小范围局部海岸线提取，并以小范围局部海岸线提取结果精度进行精确评估。
图1是IM1养殖池作为海域的一部分，如图1a所示，陆地区域设置为红色目标，海洋设置为蓝色背景。图1b是本研究方法提取的海岸线结果；图2是养殖池边界作为水陆分界，同样的，如图2a所示，陆地设置为目标红色，海洋设置为背景蓝色，图2b是本研究方法提取的海岸线结果。两种不同情况下，得到的大范围区域海岸线同样具有相对平滑、连续的特点。
养殖池为陆域部分
图2 IM1的人工岸线的手动标记(a)
和半自动提取结果(b)
图3是T2人工建筑海岸线，以建筑上最高水位作为水路分界，水位线作为海岸线的人工海岸线提取结果。如图3a所示，陆地作为目标设置成蓝色，海洋作为背景设置成绿色。图3b是本研究方法提取的海岸线结果，由于水位线不够平滑，故提取的海岸线虽然位置较为准确、连续，但是缺少平滑的特点。
图3 IM2的人工海岸线标记(a)及提取结果(b)
在此基础上，对提取的海岸线进行了多项式拟合，得到最终提取的拟合的海岸线如图4a所示，图4b是专家目视解译的参考海岸线位置。从图4可以看出，最终得到的海岸线接近人工目视解译的结果，而且平滑连续。
图4 IM2的人工海岸线拟合结果(a)与参考海岸线(b)对比
⒉砂质海岸线
受潮汐作用的影响，在进行无人机航摄获取遥感影像的时刻通常是瞬时水位的痕迹影像，而在实际进行海岸线的提取往往得到的是瞬时水边线，因此应进行潮位改正。一般根据无人机航摄时刻潮位信息、多年平均大潮高潮位高度、海岸DEM等参数因子，计算瞬时水边线到平均大潮高潮线之间的距离，从而对瞬时水边线进行校正，最终获得海岸线的实际位置。
故，在进行海岸线提取时，有必要融合潮位、海岸地貌、DEM等多源数据信息进行辅助，以精确推算海岸线。对于砂质海岸线提取，本研究采用了一幅IM3无人机遥感影像进行试验，其分辨率为0.05m。
本方法可根据海岸地形地貌信息灵活设置辅助标志，以提高海岸线识别的速度和精度，如图5a所示，在潮滩区域进行绿色标识，植被以上的陆地区域划分蓝色标记。提取的海岸线结果如图5b所示，为了得到较为自然的平滑效果，对提取的海岸线进行二项式拟合，如图6a所示。图6b是专家根据丰富的经验知识，结合现场海岸线的实际位置，经过目视解译得到的参考海岸线，从图6综合对比可以看出，在潮滩区域存在的模糊边界和潮沟等并没有影响本方法的提取结果，两者获得海岸线走势基本一致。
图5 IM3砂质海岸线标记(a)及提取结果(b)
图6 IM3砂质海岸线拟合结果(a)与参考海岸线(b)对比
⒊精度评定
为验证本研究提出的海岸线提取方法的精度，确保实验结果的可靠性，我们利用现场实测的数据，结合目视解译，形成一条参考背景海岸线。通过概率边缘指数PRI和边缘定位误差BDE两个指标进行定量计算对比，以评价海岸线提取结果与参考背景值之间的符合性。
⑴PRI
概率边缘指数PRI是指实际分割结果和参考分割结果中具有相同标签的像素对数与总的像素对数的比值。取值范围为[0，1]，值越接近1，表明实际分割结果与人工参考分割结果越接近。
⑵BDE
边缘定位误差BDE度量实际分割边缘和参考边缘的边缘像素的平均定位误差，一个边界像素的误差定义为其与另一分割图像的相邻最近的边缘像素的距离。取值范围为[0，+∞)，值越小，表示分割结果与人工解译结果边缘越相似。
在本实验中，PRI用来检验实际分割结果与参考结果之间的整体一致性。BDE从局部检验边缘定位的准确度。一般来讲，PRI的值越大，BDE的值越小，表明海岸线的提取精度越好。
由表1可以看出，人工岸线的PRI指数数值大于0.97，即提取海岸线与参考海岸线的结果达到了97%的重合度；BDE小于4.5个像素，根据图像分辨率可知，也就是提取海岸线与参考海岸线的实际定位误差仅小于0.22m，不足0.5m。砂质岸线的PRI指数数值大于0.88，即提取的海岸线与参考海海岸线相比，整体精度达到了88%；BDE小于2.9像素，提取海岸线与参考海海岸线的实际定位误差仅小于0.14m，即海岸线定位误差不足0.5m。本方法提取人工海岸线和砂质海岸线与实测海海岸线结果整体上相似，仅在局部细微处存在微小差异，海岸线识别精度较高。综合两种指标，提取的海岸线结果可以满足海岸线位置精度的提取要求。
表1 海岸线提取结果精度评定
序号 | 岸线类型 | 评价指标 | 指标值 |
1 | 人工岸线 | PRI | 0.974 |
BDE | 4.334 | ||
2 | 砂质岸线 | PRI | 0.883 |
BDE | 2.842 |
三、结论
本研究将面向对象的思想方法融入到无人机遥感影像海岸线提取中，分别对人工海岸线和砂质海岸线进行了提取实验，并通过概率边缘指数PRI和边缘定位误差BDE进行精度评估，人工海岸线和砂质海岸线与实测海岸线结果整体结构上基本相似，仅在局部细微处存在微小差异，表明该方法海岸线识别速度快、精度较高，是一种有价值的无人机遥感影像海岸线提取算法，能够满足大比例尺成图要求，可在实际岸线资源管理中推广使用。
1
END
1
【作者简介】文/麻德明 刘焱雄 金永德 程宝权，分别来自自然资源部第一海洋研究所、中国海洋大学环境科学与工程学院和青建集团股份公司。第一作者麻德明，1982年出生，男，山东胶南人，工程师，在职博士生，研究方向为海洋测绘与GIS应用。本文为基金项目，海洋公益性科研专项(201405028-4)；文章来自《海洋科学》（2020年第10期），参考文献略，用于学习与交流，版权归作者及出版社共同拥有，转载也请备注由“溪流之海洋人生”微信公众号编辑与整理。
相关阅读推荐
论文专区▏海岸带遥感影像与海图融合显示技术研究海测论坛▏LiDAR技术用于我国海岸线测量的优势对比分析
论文专区▏江苏滩涂机载LiDAR测量数据处理技术
海测论坛▏基于遥感影像的海岸线提取方法综述
论文专区▏基于潮汐分带校正的海岸线遥感推算研究
论文专区▏基于水边线的遥感影像防波堤提取方法
公众号
溪流之海洋人生
微信号▏xiliu92899
用专业精神创造价值
用人文关怀引发共鸣
您的关注就是我们前行的动力
投稿邮箱▏12163440@qq.com