一种基于形态学滤波和SVD的弱小目标检测方法技术

本发明专利技术公开了一种基于形态学滤波和SVD的弱小目标检测方法，步骤1：通过形态学滤波目标增强算法进行背景抑制、噪声去除，得到预处理后图像序列；步骤2：读入Nmax幅所得图像序列，进行帧数估计，得到需要处理的帧数N；步骤3：将N+1幅图像组成的图像合并成二维数据，求其自相关矩阵并对其自相关矩阵进行SVD；步骤4：选择合适的特征向量重构图像序列，得到新的特征图像序列；步骤5：对重构图像序列进行阈值分割，从背景中分离得到原图像中弱小目标的位置；步骤6：对序列中的每幅图像分别进行修正；步骤7：将N代替Nmax后，重复步骤2～7。本发明专利技术将形态学滤波与奇异值分解的方法有效的结合对视频中的弱小目标进行检测，计算时间短，检测效率高，准确性和鲁棒性都比较好。

【技术实现步骤摘要】

本方法属于视频分析领域，具体涉及一种基于形态学滤波和SVD的弱小目标检测方法。
技术介绍
弱小目标的检测在现代战争中的地位是不言而喻的，目前已成为卫星遥感、高能物理、低空预警以及精确制导等领域信息处理的核心技术。由于弱小目标的像元个数很少，缺乏目标的结构信息，可供分割与检测算法利用的信息很少。而传感器接受的目标强度较弱，噪声和背景杂波干扰较强，使图像的信噪比降低，故而我们应当利用好序列图像中目标的连续性和规则性来检测目标。一直以来，如何更好的利用弱小目标帧间信息，提高检测的可靠性和效率，是弱小目标检测的重点。对于背景抑制与弱小目标检测问题，早期的一些工作主要是集中于动态规划和状态估计技术来增加目标的可检测性，但是，在低信噪比情况下可能呈现比较差的性能。目前，已经有人将形态学滤波、遗传算法、神经网络算法、小波变换等方法用于弱小目标检测。但是，在复杂背景下，目标点极易被噪声淹没，实现目标的可靠性检测和识别难度较大。此外，在数据吞吐量大、实时性要求高的条件下难以满足很好的检测性能。
技术实现思路
专利技术目的：针对现有技术存在的问题，本专利技术提供了一种计算时间短，检测精准的基于形态学滤波和奇异值分解(下文简称SVD)的弱小目标检测方法。
技术实现思路
：本专利技术提供了一种基于形态学滤波和SVD的弱小目标检测方法，包括以下步骤：步骤1：输入待检测的视频序列，通过形态学滤波目标增强算法进行背景抑制、噪声去除，得到预处理后图像序列；步骤2：从预处理后的图像序列中，读入Nmax幅图像组成的图像序列，进行帧数估计，得到需要处理的帧数N；步骤3：读入N+1幅图像，包括需要处理的N帧图像和N帧图像的后一幅图像，并将N+1幅图像组成的图像合并成二维数据，求其自相关矩阵并对其自相关矩阵进行SVD；步骤4：选择中间特征值对应特征向量重构图像序列，得到重构图像序列；步骤5：对重构图像序列进行阈值分割，从背景中分离得到原图像中弱小目标的位置；步骤6：对步骤4中得到的重构图像序列中的每幅图像分别进行帧间位置修正与帧内位置修正；步骤7：将N代替Nmax后，重复步骤2～7，直到Nmax幅图像组成的图像序列中最后一幅图像处理完成后输出结果。进一步，所述步骤1中形态学滤波目标增强算法采用圆形作为结构元素。这对于恢复噪声污染图像会产生较好的滤波效果。这样能够得到更好的滤波效果。进一步，所述步骤1中利用形态学滤波目标增强算法先对待检测视频中的每一帧图像进行闭运算，然后对有能填入砂眼噪声之间的图像内部或不会形成退化矩形的区域进行开运算。进一步，所述步骤2中帧数估计的方法为：读入Nmax幅图像组成的图像序列，对读入的所有图像进行两两差分操作得到Nmax-1帧图像序列，然后对Nmax-1帧图像序列中的每一幅图像总像素求和，再用求得的Nmax-1个像素和分别除以对应图像的宽和高得到Nmax-1个评估值，将Nmax-1个评估值求平均到Nmax帧图像变化剧烈程度的值X，最后将图像变化剧烈程度X代入公式N=Nmin,X<0.5Nmax+Nmin2+Nmax-Nmin2cos(πX-π2),X∈[0.5,1.5]Nmax,X>1.5]]>中，得到本次需要处理的图像帧数N，其中，Nmax表示最大处理帧数，Nmin表示最小处理帧数。进一步，所述最大处理帧数Nmax设为25，最小处理帧数Nmin设为5。进一步，所述步骤5中采用最大熵方法对图像进行阈值分割。进一步，所述步骤6中修正方法为：在当前搜索窗口内寻找灰度最大值，然后将灰度最大值的坐标作为下次搜索窗口的中心位置迭代搜索，直到最终搜索窗口不变为止，这时搜索窗口中心的位置即为目标位置。有益效果：与现有技术相比，本专利技术将形态学滤波与奇异值分解的方法有效的结合对视频中的弱小目标进行检测，不仅计算时间短，检测效率高，而且准确性和鲁棒性都比较好。附图说明图1是本专利技术的工作流程图。图2(a)和图2(b)分别为检测视频中第20帧和第50帧图像；图3(a)和图3(b)分别为基于形态学滤波方法对检测视频中第20帧和第50帧图像中弱小目标检测的结果；图4(a)和图4(b)分别为基于SVD对检测视频中第20帧和第50帧图像中弱小目标检测的结果；图5(a)和图5(b)分别为本专利技术对检测视频中第20帧和第50帧图像中弱小目标检测的结果。具体实施方式下面结合附图，对本专利技术进行详细说明。如图1所示，本专利技术提供了一种基于形态学滤波和SVD的弱小目标检测方法，包括以下步骤：步骤1：输入待检测的视频序列，通过形态学滤波目标增强算法进行背景抑制、噪声去除，得到预处理后图像序列；考虑到弱小目标的检测，或多或少会遇到噪声图像与非噪声图像发生重叠形成结团或者某些噪声粒子的半径超过了非噪声粒子的半径，那么这种情况可选择合适半径的圆形作为结构元素，这对于恢复噪声污染图像会产生较好的滤波效果，是因为：(1)圆形的圆化作用可以起到低通滤波的效果；(2)采用圆形滤波可以不必考虑旋转的影响。在确定圆形结构元素的半径时,可采用优化方法,将图象和噪声视为随机过程,通过统计分析,对被噪声污染的颗粒图象进行数量分析,求取统计分布参数,获得出现概率最大的噪声颗粒和未被噪声污染颗粒的半径,选取未被噪声污染颗粒的半径为结构元素半径,得到优化结果。此外，弱小目标的序列中可能既存在胡椒状噪声(差噪声)，也存在砂眼噪声(并噪声)，这样单纯的使用开运算或者闭运算效果都不会好，为了兼容平滑噪声和保留图像边缘及其他有意义的特征，本专利技术先对待检测视频中的每一帧图像进行闭运算，然后对可以进行开运算的区域开运算，能够得到更好的滤波效果，其中可以进行开运算的区域为有可填入砂眼噪声之间的图像内部,不会形成严重退化矩形的区域。根据公式和定义，我们不难得知，该形态学滤波器具有良好的平移不变性、递增性、幂等性和对偶性。这些性质使得此种开-闭或者闭-开滤波器实用性和可行性更强，取得的效果也更佳。步骤2：从预处理后的图像序列中，读入Nmax幅图像组成的图像序列，进行帧数估计，得到需要处理的帧数N；在弱小目标序列中，不同目标需要一次性处理的帧数并不相同，所以需要通过自适应确定来得到需要处理的帧数N。在选择帧数的时候，首先要考虑的是图像变化快慢，即变化快时选择较少帧数来处理，反之选择较多帧数。故而，帧数估计的时候需要确定最大处理帧数Nmax本文档来自技高网...

【技术保护点】
一种基于形态学滤波和SVD的弱小目标检测方法，其特征在于：包括以下步骤：步骤1：输入待检测的视频序列，通过形态学滤波目标增强算法进行背景抑制、噪声去除，得到预处理后图像序列；步骤2：从预处理后的图像序列中，读入Nmax幅图像组成的图像序列，进行帧数估计，得到需要处理的帧数N；步骤3：读入N+1幅图像，包括需要处理的N帧图像和N帧图像的后一幅图像，并将N+1幅图像组成的图像合并成二维数据，求其自相关矩阵并对其自相关矩阵进行SVD；步骤4：选择中间特征值对应特征向量重构图像序列，得到重构图像序列；步骤5：对重构图像序列进行阈值分割，从背景中分离得到原图像中弱小目标的位置；步骤6：对步骤4中得到的重构图像序列中的每幅图像分别进行帧间位置修正与帧内位置修正；步骤7：将N代替Nmax后，重复步骤2～7，直到Nmax幅图像组成的图像序列中最后一幅图像处理完成后输出结果。

【技术特征摘要】
1.一种基于形态学滤波和SVD的弱小目标检测方法，其特征在于：包括以下步骤：
步骤1：输入待检测的视频序列，通过形态学滤波目标增强算法进行背景抑制、噪
声去除，得到预处理后图像序列；
步骤2：从预处理后的图像序列中，读入Nmax幅图像组成的图像序列，进行帧数估
计，得到需要处理的帧数N；
步骤3：读入N+1幅图像，包括需要处理的N帧图像和N帧图像的后一幅图像，并
将N+1幅图像组成的图像合并成二维数据，求其自相关矩阵并对其自相关矩阵进行
SVD；
步骤4：选择中间特征值对应特征向量重构图像序列，得到重构图像序列；
步骤5：对重构图像序列进行阈值分割，从背景中分离得到原图像中弱小目标的位
置；
步骤6：对步骤4中得到的重构图像序列中的每幅图像分别进行帧间位置修正与帧
内位置修正；
步骤7：将N代替Nmax后，重复步骤2～7，直到Nmax幅图像组成的图像序列中最后
一幅图像处理完成后输出结果。
2.根据权利要求1所述的基于形态学滤波和SVD的弱小目标检测方法，其特征在
于：所述步骤1中形态学滤波目标增强算法采用圆形作为结构元素。
3.根据权利要求1所述的基于形态学滤波和SVD的弱小目标检测方法，其特征在
于：所述步骤1中利用形态学滤波目标增强算法先对待检测视频中的每一帧图像进行闭
运算，然后对有能填入砂眼噪声之间的图像内部或不会形成退化矩形的区域进行开运
算。
4.根据权利要求1所述的基于形态学滤波...

【专利技术属性】
技术研发人员：王敏，
申请(专利权)人：河海大学，
类型：发明
国别省市：江苏;32