CN107657272A - 一种遥感图像海上目标检测方法 - Google Patents

Abstract

本发明公开了一种遥感图像海上目标检测方法，涉及遥感领域中目标显著性检测、目标检测识别等应用。它首先利用SIFT特征提取获得计算显著度的基础；再根据特征点的位置计算目标所占外接矩形面积；最后保留每一维特征值，对其进行累加作为单点显著度，累加各个点的显著度并除以目标的矩形面积，获得目标的显著度度量。本发明方法具有效率高、准确率高、鲁棒性好等特点，是对现有技术的一种重要改进。

Description

一种遥感图像海上目标检测方法

技术领域

本发明属于目标检测技术领域，特别是指一种遥感图像海上目标检测方法。

背景技术

在遥感领域中目标检测的主要依据是目标的显著性，显著性的形式有如下几种，但它们在性能上均存在某些缺陷。

(1)简单的像素累加计算，虽然可以形成初步的显著性表示，但是缺乏鲁棒性和稳定性。

(2)Sliency模型使用了初级视觉特征如颜色、亮度和方位，并将每个特征图分成几个不同尺度的高斯金字塔图，用中央周边差分的方法在细尺度和粗尺度中模拟视觉接受域，这些特征通过不同尺度间的相减被合并成亮度、颜色和方位3个特征图，这些特征图被归一化后联合生成一个显著图，该显著图中利用生物学中赢者取全的机制得到图像中最显著的空间位置，用来向导注意位置的选取，最后采用返回抑制的方法来完成注意焦点的转移。但由于该模型没有利用任何的先验知识同时无法体现目标的属性，难以直接应用到遥感图像目标检测中。

发明内容

有鉴于此，本发明提供一种遥感图像海上目标检测方法，该方法具有高效率、高准确率的特点，并具有较好的鲁棒性，适合用于可见光遥感影像处理等领域。

为了实现上述目的，本发明所采用的技术方案是：

一种遥感图像海上目标检测方法，其包括以下步骤：

(1)在遥感图像上进行SIFT特征提取；

(2)利用SIFT特征的128维向量，以特征点的特征值密度作为目标区域的显著度；

(3)根据步骤(2)所获得的显著度以及预设阈值，检测海上目标。

可选的，所述步骤(1)包括：

(101)生成系列不同尺度空间的差分图像；

(102)将各差分图像的局部极值点确定为初始特征点p_i，其中i为特征点顺序号；

(103)去除不稳定的特征点；

(104)通过特征点描述获得SIFT特征点。

可选的，所述步骤(2)包括：

(201)在诸多SIFT特征点坐标位置中，选取其中横纵坐标的最大值和最小值，形成目标外接矩形的四个顶点，并求算目标外接矩形的面积；

(202)每一个SIFT特征点具有128维特征，保留每一维特征值，并对所有128个特征值进行累加，将累加的结果作为单点显著性的评价标准；

(203)将目标区域内所有特征点的单点显著性累加，并将累加结果与所述目标外接矩形的面积的比值作为目标区域的显著度。

与现有技术相比，本发明具有如下有益效果：

1、本发明克服了现有方法中缺乏鲁棒性和稳定性的缺点，改善了显著性检测的性能。

2、本发明具有较好的实际应用和海上环境检测的性能，能够满足海上环境下的显著目标检测应用要求。

附图说明

图1是本发明实施例中SIFT特征点的邻域及其划分示意图。

图2是本发明实施例中目标特征点集的外接矩形示意图。

图3a、图3b是本发明实施例中军舰目标的SIFT特征点提取示意图。

图4a、图4b是本发明实施例中潜艇目标的SIFT特征点提取示意图。

图5a、图5b是本发明实施例中民用船只目标的SIFT特征点提取示意图。

图6是本发明实施例中军舰目标的显著度计算的示意图。

图7是本发明实施例中潜艇目标的显著度计算的示意图。

图8是本发明实施例中民用船只目标的显著度计算的示意图。

图9是本发明实施例中有云遥感图像显著度计算的示意图。

图10是本发明实施例中有海浪的遥感图像显著度计算的示意图。

图11是本发明实施例的一种方法流程图。

具体实施方式

为使本发明的目的、技术方案和优点更加清楚明白，以下结合具体实施例，并参照附图，对本发明做进一步的详细说明。

如图11所示，一种遥感图像海上目标检测方法，其包括以下步骤：

(1)在遥感图像上进行SIFT特征提取；

(2)利用SIFT特征的128维向量，以特征点的特征值密度作为目标区域的显著度；

(3)根据步骤(2)所获得的显著度以及预设阈值，检测海上目标。

其中阈值的预设值要根据实际情况进行选取，该选取方式为现有技术，此处不再赘述。

可选的，所述步骤(1)包括：

(101)生成系列不同尺度空间的差分图像；

设原图像为I₁(x，y)，原图的增倍插值图像为I₀(x，y)、隔点采样图像为I_r(x，y)，r＝1，2，其中x、y为图像坐标。

首先通过高斯函数G(x，y，kσ)对图像I_r(x，y)生成不同尺度空间的图像：

L_rk(x，y)＝G(x，y，kσ)*I_r(x，y) k＝1，2，…，K，r＝0，1，2 (1)

将r为定值的相邻L_rk(x，y)相减生成差分图像

D_rk＝L_r(k+1)-L_rk k＝1，2，…，K-1，r＝0，1，2 (2)

(102)将各差分图像的局部极值点确定为初始特征点p_i，其中i为特征点顺序号；p_i应满足如下条件：

①p_i位于图像D_rk，k＝2，…K-2之中，r＝0，1，2，该点在图像D_rk中的取值为D_rk(p_i)；

②设p_il，l＝1，2，…8为p_i的八邻域点，则

D_rk(p_i)＝max{D_rk(p_i)，D_rk(p_il)，D_r(k-1)(p_i)，D_r(k-1)(p_il)，D_r(k+1)(p_i)，D_r(k+1)(p_il)} (3)

或

D_rk(p_i)＝min{D_rk(p_i)，D_rk(p_il)，D_r(k-1)(p_i)，D_r(k-1)(p_il)，D_r(k+1)(p_i)，D_r(k+1)(p_il)} (4)

(103)去除不稳定的特征点；具体包括：

①去掉弱极值点。对差分图像函数D_rk进行泰勒展开，将D_rk中初始特征点的极值点p_i带入展开式，得到其函数值D_rk(p_i)。若|D_rk(p_i)|＜0.03，则视D_rk中p_i为不稳定点，剔除之。

②去掉边缘特征点。位于边缘的特征点会在边缘方向有很大的主曲率而在垂直边缘的方向主曲率很小。借助对D_rk的计算得到2×2海森矩阵H，根据(Tr(H))²与Det(H)的比值进行边缘不稳定点的剔除，其中Tr()为求矩阵的迹，Det()为求矩阵的行列式。

(104)通过特征点描述获得SIFT特征点；具体方式为：

①遍历各D_rk中所有特征点，以特征点p_i为中心，圈定在L_rk上的17×17的邻域Ω_i；

②对Ω_i内的每一点(x，y)，按照下式求其变化幅值m和幅角θ：

③将Ω_i划为16个4×4的子区域(ψ_i，i＝1，2，…，16)，如图1所示。将每个子域的16个点的变化方向值依临近规则归到0°，45°，90°，135°，180°，225°，270°，315°八个方向上，再对具有相同变化方向的变化幅值进行累加，得到一个8维矢量。则16个子区域得到16个8维矢量：

{m_0°(ψ_n)，m_45°(ψ_n)，…，m_315°(ψ_n)}＝{m_k(ψ_n)|k＝1，…，8}，n＝1，…，16 (7)

④将p_i的主方向定义为：

⑤令其中k′＝1，…，8；n＝1，…，16，则

根据主方向得到调整后的16个8维矢量{m_k′(ψ_n)|k′＝1，2，…，8}。沿约定的子区域顺序，将属于特征点p_i的16个子域ψ_n的16个8维矢量依次拼接形成128维向量：

向量中，主方向所决定的各矢量使SIFT具有了旋转不变性。正因为上述算法是在空间域和尺度域上同时进行特征点的计算与提取的，因此得到的特征点具有尺度不变性，能够正确的提取尺度和视角变化较大的图像序列中存在的特征点，可以有效地防范输入图像中的噪声干扰。图3a、图3b、图4a、图4b、图5a、图5b分别表示了军舰、潜艇和民用船只SIFT特征提取的结果示意图。

可选的，所述步骤(2)包括：

(201)在诸多SIFT特征点坐标位置中，选取其中横纵坐标的最大值和最小值，形成目标外接矩形的四个顶点，并求算目标外接矩形的面积；

设图中有N个SIFT特征点，特征点为P_i，i＝1…N，如图2所示。则最大X坐标(X_max)、最大Y坐标(Y_max)、最小X坐标(X_min)、最小Y坐标(Y_min)分别为如下定义：

则目标外接矩形的长a和宽b分别为：

a＝X_max-X_min (15)

b＝Y_max-Y_min (16)

最后，得到目标外接矩形的面积：

S＝a×b (17)

(202)每一个SIFT特征点具有128维特征，保留每一维特征值，并对所有128个特征值进行累加，将累加的结果作为单点显著性的评价标准；

由于SIFT特征提取时进行了归一化处理，因此在此不再进行归一化。则单点的显著性计算公式如下：

(203)将目标区域内所有特征点的单点显著性累加，并将累加结果与所述目标外接矩形的面积的比值作为目标区域的显著度。

得到单点的显著性特征值以后，需要计算目标区域的显著度。特征点的数量是一个考量标准，但是SIFT的128维向量也应该作为显著性的一个度量因素。特征点足够多，但是其分布状态不一样，引起的关注度也不一样。特征点聚集在一起比分散到周围更能引起人类的注意力。故，定义特征点特征值密度作为目标区域的显著度：

图6～10为针对军舰、潜艇、民用船只、云遥感图像、有海浪的遥感图像所做的显著度计算的示意图。可以看出，针对人造目标，其显著性明显增强，自然现象的显著性较低，这样就可以有效地区分开遥感影像中自然现象对目标的干扰。因此，本方法具有一定的鲁棒性。

此外，为了进一步检验本发明方法的效果，将阈值设为0.03，使用本发明的方法进行舰船检测，并与文献[1](Christina Corbane，Laurent Najman，EmilienPecoul，et al.Acomplete processing chain for ship detection using optical satellite imagery[J].International Journal of Remote Sensing，2010，31(22)：5837-5854.)和文献[2](Yang G，Li B，Ji S，et al.Ship Detection From Optical Satellite Images Based onSea Surface Analysis[J].IEEE Geoscience&Remote Sensing Letters，2014，11(3)：641-645.)的方法进行对比，对比结果如表1所示：

表1实验统计对比

方法	目标总数	检测正确数	虚警数	准确率	虚警率
						方法[1]	200	163	26	81.5％	13％
方法[2]	200	171	19	85.5％	9.5％
						本发明方法	200	190	6	95％	3％

可见，在检测准确率和虚警率方面本发明方法相对于现有技术均取得了显著进步。

总之，本发明的发明人通过观察遥感图像SIFT特征提取的结果发现，人造目标的SIFT特征点十分密集，并呈现出目标应有的轮廓性和规律性，因此提出利用SIFT特征128向量并结合特征点密度量的方式来定义显著度。该方法具有检测性能高的特点，特别适用于可见光遥感影像的海上目标检测、海上目标定位、海上目标识别预处理等应用，是对现有技术的一种重要改进。

所属领域的普通技术人员应当理解：以上任何实施例的讨论仅为示例性的，并非旨在暗示本公开的范围(包括权利要求)被限于这些例子。凡在本发明的精神和原则之内，对以上实施例所做的任何省略、修改、等同替换、改进等，均应包含在本发明的保护范围之内。

Claims (3)

1.一种遥感图像海上目标检测方法，其特征在于，包括以下步骤：

(1)在遥感图像上进行SIFT特征提取；

(2)利用SIFT特征的128维向量，以特征点的特征值密度作为目标区域的显著度；

(3)根据步骤(2)所获得的显著度以及预设阈值，检测海上目标。

2.根据权利要求1所述的遥感图像海上目标检测方法，其特征在于，所述步骤(1)包括：

(101)生成系列不同尺度空间的差分图像；

(102)将各差分图像的局部极值点确定为初始特征点p_i，其中i为特征点顺序号；

(103)去除不稳定的特征点；

(104)通过特征点描述获得SIFT特征点。

3.根据权利要求1所述的遥感图像海上目标检测方法，其特征在于，所述步骤(2)包括：

(201)在诸多SIFT特征点坐标位置中，选取其中横纵坐标的最大值和最小值，形成目标外接矩形的四个顶点，并求算目标外接矩形的面积；

(202)每一个SIFT特征点具有128维特征，保留每一维特征值，并对所有128个特征值进行累加，将累加的结果作为单点显著性的评价标准；

(203)将目标区域内所有特征点的单点显著性累加，并将累加结果与所述目标外接矩形的面积的比值作为目标区域的显著度。