一种利用光谱遥感图像识别海上人工设施的方法与流程

本发明涉及一种利用光谱遥感图像识别海上人工设施的方法，属于遥感图像智能识别的技术领域。

背景技术：

管理、维护与利用海域中的岛礁是国家主权的重要内容之一，探明海域内的海上人工设施，对于保障海上正常航行、开发利用与维护国家海洋权益都有重要意义。由于历史的种种原因，我国南海海上人工设施的了解还远没有满足实际的需求。随着东海、南海等海域内主权争端和开发利用需求的升温，探明我海域海上人工设施情况的需求也越来越迫切。

探明我海域海上人工设施情况有两大类方法：一类是基于船只抵近获取的传统探测方法，由于该方法必须利用船只抵近勘探，难以在船只无法到达的海域开展探测(例如在存在主权归属争议的岛屿附近)。另一类是基于遥感的探测方法，该方法无需抵近勘探，适用于任何海域。

目前，基于遥感探测人工设施的方法主要还是依赖人工判图，无法大规模大范围判明人工设施，特别是小型人工设施受遥感空间分辨率影响，可能会存在几个像素甚至亚像素的情况，这时人工判图基本无效。

技术实现要素：

本发明的目的就在于克服现有技术的不足，利用现有的遥感成像手段，提供一种利用光谱遥感图像识别海上人工设施的方法。

为了实现上述目的，本发明提供了一种利用光谱遥感图像识别海上人工设施的方法，包括如下步骤：

(1)背景分析：获取所要检测人工设施的相关信息，包括大小、形状和材质；

(2)异常检测：根据步骤(1)中获取的所要检测人工设施的大小与遥感图像空间分辨率大小的对比，分为纯像素异常检测和亚像素异常检测，将检测出的异常结果根据人工设施的大小进行初步筛选，滤除过大和过小的结果，剩余的检测结果即为疑似目标；

(3)光谱特征提取：由于海水的光谱曲线较为平缓，而水下礁滩在近红外部分几乎不反射任何光线，因此水下礁滩的光谱曲线在近红外部分会有明显下降梯度,所以采用近红外波段与红波段的波段比值用于鉴别水下礁滩与海水；由于海上人工设施与水上礁滩的光谱没有固定的差异，且不同材质的礁滩光谱曲线变化复杂，因此采用光谱导数挖掘海上人工设施与水上礁滩的光谱鉴别特征；最后，采用特征选择方法对波段比值与光谱导数特征进行筛选，得到有效的光谱特征；

(4)图像特征提取：海上人工设施具有明显的几何特征，由于红波段图像受大气散射的影响较小，地物影像清晰，对人文地物判读有利，不同地质构造的边界在图像上有显示的反映，因此只在红波段图像对海上人工设施进行几何特征的提取；

(5)由于海上人工设施必须同时具备露出海面以及一定形状的特点，因此采用步骤(3)得到的光谱特征和步骤(4)得到的图像特征与步骤(2)得到的疑似目标进行光谱匹配和图像匹配，当同时满足光谱匹配和图像匹配条件的疑似目标即为人工设施。

所述背景分析前对高光谱数据进行数据预处理，主要采用包括去除水吸收波段的方法去除噪声过大的波段。

步骤(3)中，特征选择方法主要是选择海上人工设施内部差异小和与非人工设施间差异大的特征，非人工设施主要包括海水、水下礁滩、水上礁滩，其中选择与非人工设施间差异和海上人工设施内部差异比值较大的一组特征参与后续的光谱匹配。

步骤(3)中的一阶光谱导数定义为：

其中δλ＝λj-λi，s(λj)是指波段j处的光谱值，s(λi)是指波段i处的光谱值；

二阶光谱导数定义为：

其中δλ＝λk-λj-λi；s(λk)是指波段k处的光谱值

其他高阶的光谱导数也可以用类似的方法获取，定义如下：

其中：s(λi+n)是指波段i+n处的光谱值。

步骤(4)中海上人工设施进行几何特征提取的方法：

主要利用霍夫变换提取线特征和圆特征，霍夫变换线特征提取主要利用极坐标系直线用参数极径和极角(r，θ)表示，由于对于每一个给定点，所有通过它的直线将在极坐标系得到一条正弦曲线，如果比较多的正弦曲线相交于一点，则认为这些点构成了一条直线；

霍夫变换圆特征提取主要利用极坐标系直线用参数极径和圆心点(r，x,y)表示，由于对于每一个给定点，所有通过它的直线将在极坐标系得到一条正弦曲线，如果比较多的正弦曲线相交于一点，则认为这些点构成了一条直线。

步骤(5)中：光谱匹配主要是根据步骤(3)得到的已知海上人工设施的光谱特征与疑似目标进行欧式距离计算匹配，如果距离较近，在已知海上人工设施内部动态范围内，则在光谱方面具有高度相似性；

图像匹配主要是根据步骤(4)得到的已知海上人工设施的几何特征，如果疑似目标的几何特征在已知海上人工设施内部动态范围内，则在图像方面具有高度相似性；

如果疑似目标在光谱和图像方面均与已知海上人工设施具有高度相似性，则认为疑似目标为人工设施，否则为非人工设施。

与现有技术相比，本发明的有益效果：本发明所提供的方法填补了利用遥感图像自动识别海上人工设施的空白。适用于高光谱传感器对海上人工设施的检测识别，在多光谱图像上也可适用，具有良好的应用前景。

附图说明

为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。

图1为本发明实施例的流程图。

图2是本发明实施例中的极坐标示意图。

具体实施方式

下面结合附图对发明进一步说明，但不用来限制本发明的范围。

实施例

如图1所示，本发明提供的一种利用光谱遥感图像识别海上人工设施的方法，包括如下步骤：

(1)背景分析：获取所要检测人工设施的相关信息，包括大小、形状和材质等，以为后续步骤提供输入信息，其大小用于步骤(2)中确定异常检测的结果是纯像素还是亚像素以及大概在图像中有多少像素，形状主要用于步骤(5)中的图像匹配，材质主要用于步骤(5)的光谱匹配。

所述背景分析前对高光谱数据进行数据预处理，主要采用包括去除水吸收波段的方法去除噪声过大的波段。

(2)异常检测：考虑到在全图进行特征提取和分析耗时比较长，为了减轻后续处理压力，首先进行异常检测；根据步骤(1)中获取的所要检测人工设施的大小与遥感图像空间分辨率大小的对比，分为纯像素异常检测和亚像素异常检测，将检测出的异常结果根据人工设施的大小进行初步筛选，滤除过大和过小的结果，剩余的检测结果即为疑似目标；

如果光谱图像单个像素的实际大小比人工设施大，则为亚像素，表示人工设施只占有不到一个像素。其中纯像素异常检测和亚像素异常检测可以采用常见的rxd等异常检测算法；rxd算法是高光谱异常检测中的一个基准算法，是一种基于广义似然比检验的恒虚警率检测算法，检测过程中不需要知道待检测目标的光谱，只需要估计背景的协方差矩阵，在背景波动范围外的标记为异常目标。

(3)光谱特征提取：由于海水的光谱曲线较为平缓，而水下礁滩在近红外部分几乎不反射任何光线，因此水下礁滩的光谱曲线在近红外部分会有明显下降梯度,所以采用近红外波段与红波段的波段比值用于鉴别水下礁滩与海水；由于海上人工设施与水上礁滩的光谱没有固定的差异，且不同材质的礁滩光谱曲线变化复杂，因此采用光谱导数挖掘海上人工设施与水上礁滩的光谱鉴别特征；最后，采用特征选择方法对波段比值与光谱导数特征进行筛选，得到有效的光谱特征；

(4)图像特征提取：海上人工设施具有明显的几何特征，其几何特征包括目标的形状、长度、宽度、长宽比、面积等，由于红波段图像受大气散射的影响较小，地物影像清晰，对人文地物判读有利，不同地质构造的边界在图像上有显示的反映，因此只在红波段图像对海上人工设施进行几何特征的提取；

(5)由于海上人工设施必须同时具备露出海面以及一定形状的特点，因此采用步骤(3)得到的光谱特征和步骤(4)得到的图像特征与步骤(2)得到的疑似目标进行光谱匹配和图像匹配，当同时满足光谱匹配和图像匹配条件的疑似目标即为人工设施。

步骤(3)中，特征选择方法主要是选择海上人工设施内部差异小和与非人工设施间差异大的特征，非人工设施主要包括海水、水下礁滩、水上礁滩等，其中选择与非人工设施间差异和海上人工设施内部差异比值较大的一组特征参与后续的光谱匹配。

步骤(3)中的一阶光谱导数定义为：

其中δλ＝λj-λi，s(λj)是指波段j处的光谱值，s(λi)是指波段i处的光谱值；

二阶光谱导数定义为：

其中δλ＝λk-λj-λi；s(λk)是指波段k处的光谱值

其他高阶的光谱导数也可以用类似的方法获取，定义如下：

其中：s(λi+n)是指波段i+n处的光谱值。

步骤(4)中海上人工设施进行几何特征提取的方法：

如图2所示，主要利用霍夫变换提取线特征和圆特征，霍夫变换线特征提取主要利用极坐标系直线用参数极径和极角(r，θ)表示，由于对于每一个给定点，所有通过它的直线将在极坐标系得到一条正弦曲线，如果比较多的正弦曲线相交于一点，则认为这些点构成了一条直线；

霍夫变换圆特征提取主要利用极坐标系直线用参数极径和圆心点(r，x,y)表示，由于对于每一个给定点，所有通过它的直线将在极坐标系得到一条正弦曲线，如果比较多的正弦曲线相交于一点，则认为这些点构成了一条直线。

步骤(5)中：光谱匹配主要是根据步骤(3)得到的已知海上人工设施的光谱特征与疑似目标进行欧式距离计算匹配，如果距离较近，在已知海上人工设施内部动态范围内，则在光谱方面具有高度相似性；

图像匹配主要是根据步骤(4)得到的已知海上人工设施的几何特征，如果疑似目标的几何特征在已知海上人工设施内部动态范围内，则在图像方面具有高度相似性；

如果疑似目标在光谱和图像方面均与已知海上人工设施具有高度相似性，则认为疑似目标为人工设施，否则为非人工设施。

本发明公开的一种利用光谱遥感图像识别海上人工设施的方法。对高光谱数据预处理后进行异常检测，获取疑似目标。然后在光谱域进行提取波段比值和光谱导数等光谱特征，在图像域利用红波段的几何纹理清晰特点，提取目标图像直线和圆特征，最后通过融合分析识别感兴趣目标的方法。本发明所提供的方法填补了利用遥感图像自动识别海上人工设施的空白。适用于高光谱传感器对海上人工设施的检测识别，在多光谱图像上也可适用，具有良好的应用前景。

以上显示和描述了本发明的基本原理、主要特征和本发明的优点。本行业的技术人员应该了解，本发明不受上述实施例的限制，上述实施例和说明书中描述的只是说明本发明的原理，在不脱离本发明精神和范围的前提下本发明还会有各种变化和改进，这些变化和改进都落入要求保护的本发明范围内。本发明要求保护范围由所附的权利要求书及其等同物界定。