一种合成孔径雷达图像舰船目标协同增强方法和装置

1.本发明涉及合成孔径雷达(sar)图像处理技术领域，尤其涉及一种合成孔径雷达图像舰船目标协同增强方法和装置。


背景技术：

2.合成孔径雷达(sar)属于一种主动成像传感器，可提供海面舰船等目标的高分辨图像。相比于光学、红外等被动型传感器，sar图像不受光照和天气的影响，有全天候、全时段的监视能力。对sar图像中舰船目标的增强、分割、检测在军事海防、可持续渔业等军事民用方面有着重要的应用。
3.近年来，众多专家学者提出了多种sar图像中目标的增强方法。然而，现有的目标增强方法忽略了由单幅sar图像中分割产生的不同图像切片之间的相关性，导致了有限的目标增强效果。


技术实现要素：

4.本发明旨在至少在一定程度上解决相关技术中的技术问题之一。
5.为此，本发明的目的为克服已有技术的不足之处，提出一种合成孔径雷达图像舰船目标协同增强方法。本发明通过开发不同sar图像切片之间基于密度特征的相关性，协同增强其中的弱目标、抑制强海杂波，有利于后续目标检测性能的提升。
6.本发明的另一个目的在于提出一种合成孔径雷达图像舰船目标协同增强装置。
7.为达上述目的，本发明一方面提出了一种合成孔径雷达图像舰船目标协同增强方法，包括以下步骤：
8.输入原始合成孔径雷达sar图像，并获取所述sar图像中含有舰船目标的所有sar图像切片的切片集合，以获得所述切片集合中的超像素；基于所述切片集合计算出所述切片集合中所有超像素的局部灰度对比度，根据直方图信息计算得到所述局部灰度的超像素直方图密度值；基于所述超像素直方图密度值，计算得到交叉相关特征；基于所述交叉相关特征进行特征融合，获得图像增强结果并输出。
9.本发明实施例的合成孔径雷达图像舰船目标协同增强方法，可用于增大海面sar图像中舰船目标与海杂波的对比度，增强对舰船目标的检测性能。
10.另外，根据本发明上述实施例的合成孔径雷达图像舰船目标协同增强方法还可以具有以下附加的技术特征：
11.进一步地，所述局部灰度对比度由下面公式计算得到：
[0012][0013]
其中，为局部灰度对比度，m表示超像素的索引，为局部灰度对比度，m表示超像素的索引，表示第n个sar图像切片中超像素的数量，δ
n,m
表示第n个sar图像切片中第m个超像素的平均灰度值，表示第n个sar图像切片中第m个超像素的邻接超像素的索引集合，p表示中
的元素。
[0014]
进一步地，所述根据直方图信息计算得到所述局部灰度的超像素直方图密度值，包括：
[0015]
预设表示第n个sar图像切片中的直方图信息，则所述超像素直方图密度值，由下面公式计算得到：
[0016][0017]
其中，表示落在以为中心的直方图各自中的数量，k表示直方图划分的各自的数量,k表示直方图中格子的索引，表示截断距离，ρ
n,m,q
为超像素直方图密度值，indi(.)为当括号内的值小于零时，indi(.)输出1，否则输出0，q表示sar图像切片的索引，q＝1,2,
…
,n。
[0018]
进一步地，所述交叉相关特征由下面公式计算得到：
[0019][0020]
其中，为交叉相关特征，φq为ρ
n,m,q
的权重。
[0021]
进一步地，所述φq由下面公式计算得到：
[0022]
φq＝1/[1+κln(1+εq)]
[0023]
其中，κ＝50为常数，ε为f
ta
[θ
(intra)
]和f
ca
[θ
(intra)
]交叉点曲线下面积。
[0024]
进一步地，预设第n个切片中的目标和杂波的服从两个不同的gamma分布，则所述f
ta
[θ
(intra)
]和f
ca
[θ
(intra)
]由下面公式计算得到：
[0025][0026][0027]
其中，ta表示目标区域，ca表示杂波区域，ξ和μ表示gamma分布的形状参数和尺度参数，γ(.)表示gamma函数。
[0028]
进一步地，预设0《θ
*
《+∞表示所述f
ta
[θ
(intra)
]和f
ca
[θ
(intra)
]的交叉点，则所述ε由下面公式计算得到：
[0029][0030]
其中，γ(.,.)表示下不完备gamma函数。
[0031]
进一步地，所述图像增强结果由下面公式计算得到：
[0032][0033]
其中，n＝1,2,
…
,n，,n，为图像增强结果。
[0034]
进一步地，所述方法还包括，对ρ
n,m,q
进行归一化：
[0035][0036]
为达到上述目的，本发明另一方面提出了一种合成孔径雷达图像舰船目标协同增强装置，包括：
[0037]
像素获取模块，用于输入原始合成孔径雷达sar图像，并获取所述sar图像中含有舰船目标的所有sar图像切片的切片集合，以获得所述切片集合中的超像素；
[0038]
第一计算模块，用于基于所述切片集合计算出所述切片集合中所有超像素的局部灰度对比度，根据直方图信息计算得到所述局部灰度的超像素直方图密度值；
[0039]
第二计算模块，用于基于所述超像素直方图密度值，计算得到交叉相关特征；
[0040]
融合输出模块，用于基于所述交叉相关特征进行特征融合，获得图像增强结果并输出。
[0041]
本发明实施例的合成孔径雷达图像舰船目标协同增强装置，可用于增大海面sar图像中舰船目标与海杂波的对比度，增强对舰船目标的检测性能。
[0042]
本发明的有益效果：
[0043]
现有的sar图像中舰船目标增强方法的不足之处在于，只开发了单个sar图像切片内部的特征，并没有联合考虑不同海面sar图像切片之间的相似性，例如，舰船目标在不同海面sar图像切片中均呈现空域的稀疏性。图像切片相关性的考虑缺乏这导致现有目标增强方法的性能较低。本发明基于不同海面sar图像切片中的舰船目标的低密度特征，提出了密度特征的sar图像舰船目标增强方法，通过开发不同切片之间的目标/杂波在密度上相似性，显著增强了舰船目标，有望提升我国对海面舰船目标的快速响应能力。
[0044]
本发明附加的方面和优点将在下面的描述中部分给出，部分将从下面的描述中变得明显，或通过本发明的实践了解到。
附图说明
[0045]
本发明上述的和/或附加的方面和优点从下面结合附图对实施例的描述中将变得明显和容易理解，其中：
[0046]
图1为根据本发明实施例的合成孔径雷达图像舰船目标协同增强方法的流程图；
[0047]
图2为根据本发明实施例的合成孔径雷达图像舰船目标协同增强方法的效果示意图；
[0048]
图3为根据本发明实施例的合成孔径雷达图像舰船目标协同增强装置的结构示意图。
具体实施方式
[0049]
需要说明的是，在不冲突的情况下，本技术中的实施例及实施例中的特征可以相互组合。下面将参考附图并结合实施例来详细说明本发明。
[0050]
为了使本技术领域的人员更好地理解本发明方案，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是
本发明一部分的实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都应当属于本发明保护的范围。
[0051]
下面参照附图描述根据本发明实施例提出的合成孔径雷达图像舰船目标协同增强方法及装置，首先将参照附图描述根据本发明实施例提出的合成孔径雷达图像舰船目标协同增强方法。
[0052]
图1是本发明一个实施例的合成孔径雷达图像舰船目标协同增强方法的流程图。
[0053]
如图1所示，该合成孔径雷达图像舰船目标协同增强方法包括以下步骤：
[0054]
步骤s1，输入原始合成孔径雷达sar图像，并获取sar图像中含有舰船目标的所有sar图像切片的切片集合，以获得切片集合中的超像素。
[0055]
具体的，先输入原始合成孔径雷达sar图像，下面作为一种实现方式：
[0056]
原始sar图像，尺寸为m
×
n；
[0057]
sar图像切片大小mc
×
nc,满足条件256≤mc《m,256≤nc《n；
[0058]
尺度因子α∈(0,1)，用于计算超像素直方图密度特征。
[0059]
优选的，采用加权稀疏优化方法，获得原始sar图像中含有舰船目标的所有sar图像切片，将获得的所有含有目标的切片集合记为其中，in表示含有目标的第n个切片，n＝1,2,
…
,n,n表示原始sar图像中含有舰船目标切片的sar图像的数量。
[0060]
优选的，采用简单线性迭代聚类算法获得切片集合中的超像素。
[0061]
步骤s2，基于切片集合计算出切片集合中所有超像素的局部灰度对比度，根据直方图信息计算得到局部灰度的超像素直方图密度值。
[0062]
具体地，计算中所有超像素的局部灰度对比度
[0063][0064]
其中，m表示超像素的索引，其中，m表示超像素的索引，表示第n个sar图像切片中超像素的数量，δ
n,m
表示第n个sar图像切片中第m个超像素的平均灰度值，表示第n个sar图像切片中第m个超像素的邻接超像素的索引集合，p表示中的元素。令表示第n个sar图像切片中所有的集合。
[0065]
进一步地，令表示第n个sar图像切片中的直方图信息，其中表示落在以为中心的直方图各自中的数量,k表示直方图划分的各自的数量,k表示直方图中格子的索引。接下来，计算第n个sar图像切片中每一个在第q个(q＝1，2，
……
n)sar图像切片中的超像素直方图密度值ρ
n,m,q
：
[0066]
[0067]
其中，表示截断距离，indi(.)的含义为，当括号内的值小于零时，indi(.)输出1，否则输出0，q表示sar图像切片的索引，q＝1,2,
…
,n。接下来对ρ
n,m,q
进行归一化：
[0068][0069]
步骤s3，基于超像素直方图密度值，计算得到交叉相关特征。
[0070]
具体地，本发明开发不同sar图像切片在密度特征上的相似性，计算交叉相关特征
[0071][0072]
其中，φq表示ρ
n,m,q
的权重。φq的计算方式如下。假设第n个切片中的目标和杂波的服从两个不同的gamma分布：
[0073][0074][0075]
其中，ta表示目标区域，ca表示杂波区域，ξ和μ表示gamma分布的形状参数和尺度参数，γ(.)表示gamma函数，{ξ
ta
,μ
ta
,ξ
ca
,μ
ca
}可以由otsu算法获得。令0《θ
*
《+∞表示f
ta
[θ
(intra)
]和f
ca
[θ
(intra)
]曲线的交叉点，计算f
ta
[θ
(intra)
]和f
ca
[θ
(intra)
]交叉点曲线下面积ε：
[0076][0077]
其中，γ(.,.)表示下不完备gamma函数。则权重的计算方式如下：
[0078]
φq＝1[1+κln(1+εq)]
[0079]
其中，κ＝50为常数。
[0080]
步骤s4，基于交叉相关特征进行特征融合，获得图像增强结果并输出。
[0081]
具体地，根据交叉相关特征，再进行sar图像地特征融合。对每一个sar图像切片中的超像素获取其增强结果
[0082][0083]
则输出的结果为：n＝1,2,
…
,n，
[0084]
进一步地，图2为本发明实施例的技术效果图，如图2所示，图中第一行为sar图像切片；第二行现有的方法所呈现的图像效果；第三行为利用本发明所呈现的图像效果。
[0085]
根据本发明实施例的合成孔径雷达图像舰船目标协同增强方法，通过输入原始合成孔径雷达sar图像，并获取sar图像中含有舰船目标的所有sar图像切片的切片集合，以获得切片集合中的超像素；基于切片集合计算出切片集合中所有超像素的局部灰度对比度，
根据直方图信息计算得到局部灰度的超像素直方图密度值；基于超像素直方图密度值，计算得到交叉相关特征；基于交叉相关特征进行特征融合，获得图像增强结果并输出。本发明通过开发不同sar图像切片之间基于密度特征的相关性，协同增强其中的弱目标、抑制强海杂波，有利于后续目标检测性能的提升。
[0086]
为了实现上述实施例，如图3所示，本实施例中还提供了一种合成孔径雷达图像舰船目标协同增强装置10，该装置10包括：像素获取模块100、第一计算模块200、第二计算模块300和融合输出模块400。
[0087]
像素获取模块100，用于输入原始合成孔径雷达sar图像，并获取sar图像中含有舰船目标的所有sar图像切片的切片集合，以获得切片集合中的超像素；
[0088]
第一计算模块200，用于基于切片集合计算出切片集合中所有超像素的局部灰度对比度，根据直方图信息计算得到局部灰度的超像素直方图密度值；
[0089]
第二计算模块300，用于基于超像素直方图密度值，计算得到交叉相关特征；
[0090]
融合输出模块400，用于基于交叉相关特征进行特征融合，获得图像增强结果并输出。
[0091]
根据本发明实施例的合成孔径雷达图像舰船目标协同增强装置，通过输入原始合成孔径雷达sar图像，并获取sar图像中含有舰船目标的所有sar图像切片的切片集合，以获得切片集合中的超像素；基于切片集合计算出切片集合中所有超像素的局部灰度对比度，根据直方图信息计算得到局部灰度的超像素直方图密度值；基于超像素直方图密度值，计算得到交叉相关特征；基于交叉相关特征进行特征融合，获得图像增强结果并输出。本发明通过开发不同sar图像切片之间基于密度特征的相关性，协同增强其中的弱目标、抑制强海杂波，有利于后续目标检测性能的提升。
[0092]
需要说明的是，前述对合成孔径雷达图像舰船目标协同增强方法的实施例的解释说明也适用于该实施例的合成孔径雷达图像舰船目标协同增强装置，此处不再赘述。
[0093]
此外，术语“第一”、“第二”仅用于描述目的，而不能理解为指示或暗示相对重要性或者隐含指明所指示的技术特征的数量。由此，限定有“第一”、“第二”的特征可以明示或者隐含地包括至少一个该特征。在本发明的描述中，“多个”的含义是至少两个，例如两个，三个等，除非另有明确具体的限定。
[0094]
在本说明书的描述中，参考术语“一个实施例”、“一些实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本发明的至少一个实施例或示例中。在本说明书中，对上述术语的示意性表述不必须针对的是相同的实施例或示例。而且，描述的具体特征、结构、材料或者特点可以在任一个或多个实施例或示例中以合适的方式结合。此外，在不相互矛盾的情况下，本领域的技术人员可以将本说明书中描述的不同实施例或示例以及不同实施例或示例的特征进行结合和组合。
[0095]
尽管上面已经示出和描述了本发明的实施例，可以理解的是，上述实施例是示例性的，不能理解为对本发明的限制，本领域的普通技术人员在本发明的范围内可以对上述实施例进行变化、修改、替换和变型。