基于自适应尺度图像块加权差测量的红外小目标检测方法和系统技术方案

本发明专利技术公开了一种基于自适应尺度图像块加权差异测量的红外弱小目标检测方法和系统，属于图像目标检测技术领域。本发明专利技术的主要目的是解决传统红外弱小目标检测方法中存在的目标增强能力差和杂波抑制能力差的问题，从而提高目标检测中的检测率并降低虚警率。实现的技术方案是：1)利用Sobel边缘滤波算子获取原红外图像对应的边缘图像；2)根据边缘图像计算以单个像素位置为中心的红外小目标的尺寸信息；3)根据获得的尺寸信息计算自适应尺度图像块的加权差异；4)利用自适应阈值分割提取目标，完成目标检测。该方法充分利用红外弱小目标与背景间存在局部对比度差异的特点，能有效地增强目标并抑制背景杂波，提高红外弱小目标检测的性能。

【技术实现步骤摘要】
基于自适应尺度图像块加权差测量的红外小目标检测方法和系统
本专利技术属于图像目标检测
，具体涉及一种基于自适应尺度图像块加权差测量的红外小目标检测方法和系统。
技术介绍
红外弱小目标检测技术在导弹制导、红外预警和夜间目标检测等多个领域有着广泛的应用。由于红外相机的成像特点以及较远的成像距离，使得红外弱小目标通常没有具体的形状和纹理，且易受到复杂的背景杂波和噪声干扰。因此，在复杂场景下的红外弱小目标检测作为一个富有挑战性的工作受到越来越多的关注。传统的红外弱小目标检测方法是直接对图像进行滤波操作，如最大均值滤波(Max-Mean)，Top-Hat滤波和双边滤波等。近来，基于局部对比度测量(LCM)的方法开始广泛应用于红外弱小目标检测。由于单尺度的LCM方法难以适应不同大小的目标，所以多尺度的LCM被大量研究。例如：申请号(CN201711470286.X)的专利文献提出了一种基于前景加权局部对比度的红外微弱目标检测方法，引入像素的前景概率估计并将其作为局部对比度的加权因子，计算多个尺度图像块的LCM取最大值，然后通过简单的阈值分割检测目标。申请号(CN201910371192.X)的专利文献提出一种基于分级LCM的快速小目标检测方法，先采用单尺度LCM对图像进行第一级滤波，再用阈值分割得到可疑目标区域。然后再用多尺度LCM进行第二级滤波，再用阈值分割得到目标区域增强图像。将可疑目标区域与增强图像进行差值运算，并用简单阈值分割得到最终的检测结果。上述多尺度的LCM检测方法虽然能抑制一些背景杂波，但仍然存在以下不足：(1)多尺度的运算不能有效测量目标与周围背景的实际对比度，限制了目标增强效果；(2)多尺度的LCM难以抑制复杂背景，会造成高的虚警率；(3)多尺度的LCM运算复杂度高，难以实时检测目标。随着深度学习的兴起，有学者提出了用神经网络实现红外弱小目标检测，但红外弱小目标没有固定的形状和纹理特征，且在实际应用中难以获取大量而多样的数据集，所以实际应用价值有限。
技术实现思路
针对现有技术的不足，本专利技术提出一种基于自适应尺度图像块加权差测量的红外小目标检测方法。本方法充分利用红外弱小目标与周围背景构成区域存在信息熵差异的特性，用自适应尺度图像块来准确测量信息熵权值差异。本方法克服了传统的多尺度的LCM的问题，从而有效增强目标和抑制复杂背景杂波，提高目标检测率。为实现上述目的，本专利技术的基本思路是：首先，利用Sobel边缘滤波算子对原红外图像进行滤波处理，得到的对应的边缘图像；然后根据边缘图像计算以单个像素位置为中心的红外小目标的尺寸信息；接着，根据获得的尺寸信息计算自适应尺度图像块的加权差异(ASPWD)大小；最后利用自适应阈值分割提取目标，完成目标检测。本专利技术的技术方案提供一种基于自适应尺度图像块的红外弱小目标检测方法，包括以下步骤：步骤1：输入大小为M×N红外图像I，其中图像I的像素的灰度值表示为I(x,y),x＝1,2,...,M,y＝1,2,...,N；步骤2：通过Sobel算子对所述红外图像I进行边缘检测，获取边缘检测图像F，其中图像F的像素的灰度值表示为F(x,y),x＝1,2,...，M,y＝1，2，...,N；步骤3：从(x,y)＝(1,1)到(x,y)＝(M,N)遍历所述边缘检测图像F所有像素点，在左、上、右、下四个方向上，计算每个像素点(x,y)的λ范围以内的小目标自适应尺寸图L1(x,y)、L2(x,y)、L3(x，y)、L4(x，y)，具体包括以下子步骤：步骤3.1：遍历所述边缘检测图像F，寻找每一个像素(x,y)在左、上、右、下四个方向上像素距离λ以内的局部最大值点的坐标(xk,yk),k＝1,2,3,4。其中λ为常数，(x1,y1),(x2,y2),(x3,y3),(x4,y4)分别表示像素(x,y)在左、上、右、下四个方向上的局部最大值点的坐标，计算公式如下：(xk,yk)＝(i,j)|max(F(i,j)),(i,j)∈Ωk(x,y),k＝1,2,3,4，其中Ω1(x,y)、Ω2(x,y)、Ω3(x,y)、Ω4(x,y)分别表示以像素(x,y)为中心的向左、上、右、下四个方向上像素距离λ以内的搜索区域；步骤3.2：计算像素(x,y)到所述局部最大值点的坐标(xk,yk)的距离，计算公式为：Lk(x,y)＝|x-xk|+|y-yk|,k＝1,2,3,4，其中L1(x,y)、L2(x,y)、L3(x,y)、L4(x,y)表示以(x,y)为中心到对应左、上、右、下四个方向上边缘的距离。步骤4：根据小目标尺寸图L1(x,y)、L2(x,y)、L3(x,y)、L4(x,y)，在红外图像I中，计算每个像素(x,y)的加权差异大小ASPWD(x,y)，具体包括如下子步骤：步骤4.1：在红外图像I中，对每一个像素(x,y)，以所述像素(x,y)为中心构建在左、上、右、下四个方向分别长L1(x,y)、L2(x,y)、L3(x,y)、L4(x,y)的矩形，即目标区域T(x，y)。步骤4.2：以所述像素(x,y)为中心构建1个m×n大小的矩形区域，并排除与目标区域重叠的部分，即得到背景区域B(x,y)。步骤4.3：分别计算所述目标区域T(x,y)的灰度值的均值mT(x,y)和所述背景区域B(x，y)的灰度值的均值mB(x,y)，计算公式为：其中N和NB分别表示区域T(x,y)和区域B(x,y)的像素个数。步骤4.4：分别计算所述目标区域T(x,y)的信息熵IET(x,y)和所述背景区域B(x,y)的灰度值的信息熵IEB(x,y)，计算公式为：其中nk表示区域中存在第k个像素点的灰度值的像素总数，pk表示该灰度值在区域中的概率密度。步骤4.5：根据所述目标区域均值mT(x,y)和信息熵IET(x,y)和与所述背景区域均值mBk(x,y)和信息熵IEB(x,y)，计算该像素点的加权差异ASPWD(x,y)，计算公式为：其中，ASPWD表示自适应尺度图像块加权差测量增强图。步骤5：根据上述求得的各像素点的加权差异ASPWD(x,y)构建一幅ASPWD增强图。对ASPWD增强图进行自适应阈值分割，获取红外弱小目标检测结果，具体包括如下子步骤：步骤5.1：计算ASPWD增强图像的阈值T，计算公式如下：T＝λ×Mx+(1-λ)×Mn，其中Mx和Mn分别表示ASPWD增强图像的最大值灰度值和平均灰度值，λ是一个根据实验调节的常数，推荐设置为0.5。步骤5.2：根据上述阈值T对ASPWD增强图像进行计算，得到最终检测结果图像R，其中图像R的像素的灰度值表示为R(x,y),x＝1,2,...,M,y＝1,2,...,N。计算公式为：检测结果图像R中非零像素的位置即为红外小目标位置。本专利技术还提供一种基于自适应尺度图像块加权差测量的红外小目标检测系统，包括如下模本文档来自技高网...

【技术保护点】
1.基于自适应尺度图像块加权差测量的红外小目标检测方法，其特征在于，包括如下步骤：/n步骤1，输入大小为M×N红外图像I，其中图像I的像素的灰度值表示为I(x,y),x＝1,2,...,M,y＝1,2,...,N；/n步骤2，通过Sobel算子对所述红外图像I进行边缘检测，获取边缘检测图像F，其中图像F的像素的灰度值表示为F(x,y),x＝1,2,...,M,y＝1,2,...,N；/n步骤3，从(x,y)＝(1,1)到(x,y)＝(M,N)遍历所述边缘检测图像F所有像素点，在左、上、右、下四个方向上，计算每个像素点(x,y)的λ范围以内的小目标自适应尺寸图L

【技术特征摘要】
1.基于自适应尺度图像块加权差测量的红外小目标检测方法，其特征在于，包括如下步骤：
步骤1，输入大小为M×N红外图像I，其中图像I的像素的灰度值表示为I(x,y),x＝1,2,...,M,y＝1,2,...,N；
步骤2，通过Sobel算子对所述红外图像I进行边缘检测，获取边缘检测图像F，其中图像F的像素的灰度值表示为F(x,y),x＝1,2,...,M,y＝1,2,...,N；
步骤3，从(x,y)＝(1,1)到(x,y)＝(M,N)遍历所述边缘检测图像F所有像素点，在左、上、右、下四个方向上，计算每个像素点(x,y)的λ范围以内的小目标自适应尺寸图L1(x,y)、L2(x,y)、L3(x,y)、L4(x,y)；
步骤4，根据小目标自适应尺寸图L1(x,y)、L2(x,y)、L3(x,y)、L4(x,y)，在红外图像I中，计算每个像素(x,y)的加权差异ASPWD(x，y)；
步骤5，根据上述求得的各像素点的加权差异ASPWD(x,y)构建一幅ASPWD增强图，对ASPWD增强图进行自适应阈值分割，获取红外弱小目标检测结果。


2.如权利要求1所述的基于自适应尺度图像块加权差测量的红外小目标检测方法，其特征在于：步骤3的具体实现包括以下子步骤；
步骤3.1，遍历所述边缘检测图像F，寻找每一个像素(x,y)在左、上、右、下四个方向上像素距离λ以内的局部最大值点的坐标(xk,yk),k＝1,2,3,4，其中λ为常数，(x1,y1),(x2,y2),(x3,y3),(x4,y4)分别表示像素(x，y)在左、上、右、下四个方向上的局部最大值点的坐标，计算公式如下：
(xk，yk)＝(i，j)|max(F(i，j))，(i，j)∈Ωk(x，y)，k＝1，2，3，4，
其中Ω1(x，y)、Ω2(x，y)、Ω3(x，y)、Ω4(x,y)分别表示以像素(x,y)为中心的向左、上、右、下四个方向上像素距离λ以内的搜索区域；
步骤3.2，计算像素(x,y)到所述局部最大值点的坐标(xk,yk)的距离，计算公式为：
Lk(x,y)＝|x-xk|+|y-yk|,k＝1,2,3,4，
其中L1(x,y)、L2(x,y)、L3(x,y)、L4(x,y)表示以(x,y)为中心到对应左、上、右、下四个方向上边缘的距离，也称为小目标自适应尺寸图。


3.如权利要求1所述的基于自适应尺度图像块加权差测量的红外小目标检测方法，其特征在于：步骤4的具体实现包括如下子步骤；
步骤4.1，在红外图像I中，对每一个像素(x,y)，以所述像素(x,y)为中心构建在左、上、右、下四个方向分别长L1(x，y)、L2(x，y)、L3(x，y)、L4(x，y)的矩形，即目标区域T(x,y)；
步骤4.2，以所述像素(x，y)为中心构建1个m×n大小的矩形区域，并排除与目标区域重叠的部分，即得到背景区域B(x,y)；
步骤4.3，分别计算所述目标区域T(x，y)的灰度值的均值mT(x，y)和所述背景区域B(x，y)的灰度值的均值mB(x，y)，计算公式为：



其中N和NB分别表示区域T(x，y)和区域B(x，y)的像素个数；
步骤4.4，分别计算所述目标区域T(x，y)的信息熵IET(x，y)和所述背景区域B(x，y)的灰度值的信息熵IEB(x，y)，计算公式为：



其中nk表示区域中存在第k个像素点的灰度值的像素总数，pk表示该灰度值在区域中的概率密度；
步骤4.5，根据所述目标区域均值mT(x,y)和信息熵IET(x,y)和与所述背景区域均值mBk(x，y)和信息熵IEB(x，y)，计算该像素点的加权差异ASPWD(x，y)，计算公式为：



其中，ASPWD表示自适应尺度图像块加权差测量增强图。


4.如权利要求1所述的基于自适应尺度图像块加权差测量的红外小目标检测方法，其特征在于：步骤5的具体实现包括如下子步骤；
步骤5.1，计算ASPWD增强图的阈值T，计算公式如下：
T＝λ×Mx+(1-λ)×Mn，
其中Mx和Mn分别表示ASPWD增强图的最大值灰度值和平均灰度值，λ是一个根据实验调节的常数；
步骤5.2，根据上述阈值T对ASPWD增强图进行计算，得到最终检测结果图像R，其中图像R的像素的灰度值表示为R(x，y)，x＝1，2，...，M，y＝1，2，...，N，计算公式为：



检测结果图像R中非零像素的位置即为红外小目标位置。


5.一种基于自适应尺度图像块加...

【专利技术属性】
技术研发人员：马泳，裘兆炳，樊凡，黄珺，张灿，
申请(专利权)人：武汉大学，
类型：发明
国别省市：湖北;42