[发明专利]一种基于遥感图像的沙波体范围检测方法

说明书

技术领域

本发明涉及遥感技术应用和浅海地形监测领域，特别是一种基于遥感图像的沙波体范围检测方法。

背景技术

海底地形测绘是探索和研究海洋的第一步，是开发和保护海洋的首要条件，也是现今海洋学研究的迫切需求。在回声测深仪发明之前，主要靠测深杆和测深锤来测量水深，测量精度较差。20世纪20年代回声测深仪出现之后，现代意义上海图的绘制才得以实现。但早期的测深仪为单波束发射，一次发射只能得到测量船正下方的水深，因而只能实现点、线测量，无法反映测线之间的地形地貌。1970年代中期出现的多波束测深技术实现了带状测量(最大宽度可达水深的7倍)，显著提高了测深效率。尽管如此，限于测量周期长、人力消耗大和资金需求高等方面的劣势，在地形高变化的浅海区域进行大范围的动态变化监测任务中，仅依赖多波束测深方法仍然很难满足需求。对于浅海区域广泛分布的浅海沙波地形动态变化监测，需要掌握沙波的位置和沙波体范围等信息的基础，通过多期比较分析，确定沙波的新生、移动、形态改变和消亡等动态变化信息。如何充分发挥遥感技术的大范围高动态观测能力，对浅海沙波区的沙波范围进行动态监测，对于掌握地形演化规律和特征都有十分重要的意义。

本发明针对广泛分布的浅海沙波分布区进行地形变化动态监测的需求，利用遥感技术方法的大范围覆盖和高频率重访的能力，基于遥感图像中所包含的因沙波地形对水动力调制而引起的灰度差异特征信息(亮暗条纹信息)，提取浅海沙波体的范围，用于浅海沙波的识别管理与沙波地形动态变化监测。从检索到的公开资料看，尚未有基于遥感方法进行浅海沙波体范围快速确定的方法。

发明内容

本发明所要解决的技术问题是提供一种新的基于遥感图像的沙波体范围检测方法。

本发明解决上述技术问题所采用的技术方案是：

一种基于遥感图像的沙波体范围检测方法，包括以下步骤：

(1)通过遥感器获取沙波纹理分布的遥感图像，提取沙波波脊线；

(2)以沙波波脊线的一个端点为起点，等间距地添加检测点；

(3)基于检测点添加检测线；

(4)将检测线叠加于遥感图像上，以检测点为起点，获得检测线位置上的遥感图像像元的灰度值，经过归一化处理得到检测线的遥感图像剖面线；

(5)对获得的每一条检测线的遥感图像归一化剖面线采用均值滤波法进行平滑处理，获得平滑灰度剖面线；

(6)对平滑灰度剖面线进行梯度计算，获得梯度剖面线；

(7)从沙波波脊线出发，依次读取梯度值进行判定，当梯度值小于判断阈值S时，则此梯度值所在数据点作为该检测线的范围分界点；

(8)将全部检测线的范围分界点连线，得到沙波体范围。

作为优选，所述步骤(1)中提取沙波波脊线，根据遥感图像上沙波纹理分布特征，亮暗条纹相接处为沙波波脊线的位置，采用人机交互方法提取相应的沙波波脊线。

作为优选，所述步骤(2)中等间距地添加检测点，检测点的分布间距为5到10个遥感图像像元点，检测点采用均匀分布的原则生成，具体步骤如下：

(21)根据沙波波脊线的长度、遥感图像的像元尺寸，计算确定检测点的间距r，间距r为遥感图像空间分辨率的5-10倍；

(22)以沙波波脊线的一个端点为起点，以距离r为步长，依次确定检测点的位置，添加检测点。

作为优选，所述步骤(3)中基于检测点添加检测线，首先基于每一个检测点分别向两侧添加垂直于沙波波脊线的检测线；其次在沙波波脊线两端，分别添加至少3条接近于或等于均匀间隔的检测线。

作为优选，所述步骤(4)中，将检测线叠加于遥感图像上，从遥感图像上获得与检测线相交的所有像元点，并读取其灰度值，经过归一化处理，得到检测线的遥感图像剖面线；所述归一化处理方法为：

(41)通过比较和计算获得该检测线上所有像元灰度值的平均值(DN_a)、最大值(DN_max)和最小值(DN_min)；

(42)计算任意像元点的归一化像元值，计算公式如下：