一种主动声呐回波图中小目标运动状态估计方法技术

本发明专利技术公开了一种主动声呐回波图中小目标运动状态估计方法，首先，通过高斯拉普拉斯变换(Laplacian of Gaussian，LoG)来描述声呐回波图中小目标的边缘恒常性，以此解决小目标的非刚性特性造成的目标回波不稳定的问题；然后，针对舰船运动产生的尾流和空化噪声问题，引入基于上下文信息的自适应正则化参数对目标函数中的非局部项进行正则化，最后将经过升采样的声流场图像与边缘检测图像进行warping处理，得到连续的声呐回波图像中的小目标运动状态估计结果。本发明专利技术有效地减轻了水下运动小目标非刚性特性给运动估计精度带来的影响，同时提高了当声呐回波中存在船舶尾流噪声和空化噪声时的声流估计精度。化噪声时的声流估计精度。化噪声时的声流估计精度。

【技术实现步骤摘要】
一种主动声呐回波图中小目标运动状态估计方法


[0001]本专利技术属于图像处理领域，具体涉及一种目标运动状态估计方法。

技术介绍

[0002]长期以来，利用光流等运动流场的运动估计方法被广泛地应用于智能监控、遥感、计算机视觉等多个领域。在最近30年，光流及类似的运动流场估计技术被越来越多地应用于水下目标的运动估计中。其中，可查到的文献包括：《IEEE Journal of OceanicEngineering》中1998年发表的《Robust tracking of multiple objects in sector scan sonarimage sequences using optical flow motion estimation》，该技术使用从前视图像声呐回波图中估计的2D运动场作为目标跟踪和识别的先验信息；以及，《IEEE Winter Conferenceon Applications of Computer Vision》中2014年发表的《Determining underwater vehiclemovement from sonar data in relatively featureless seafloor tracking missions》，该技术使用了基于特征点的稀疏光流来减少稠密光流算法的巨大计算量。
[0003]上述技术所采用的方法均为基于变分法的经典光流估计算法(提出者为Horn和 Schunck，通常简称为HS方法)，HS方法的核心框架于1981年在《Articial Intelligence》中发表的《Determining optical flow》一文中被提出。但是，近年来越来越多的研究表明，HS方法中使用的二次惩罚函数影响了算法在实际使用中的鲁棒性。其中，在《International Journal of Computer Vision》中2014年发表的《A quantitative analysis ofcurrent practices in optical flow estimation and the principles behind them》文中定量分析了目前已有的光流估计算法性能，并在光流目标函数中增加了加权的非局部项来应对噪声干扰。但是，从目前公开可查的资料来看，现有的基于运动流场的目标运动估计方法主要应用于侧扫声呐、高分辨图像声呐等声呐系统中。这是由于在水下小目标探测常用的高频主动声呐的回波图像中，目标回波强度的时变性让传统的光流方程约束条件难以实现。因此，迄今为止，尚未见到将该技术应用于目标状态估计中的报告。

技术实现思路

[0004]为了克服现有技术的不足，本专利技术提供了一种主动声呐回波图中小目标运动状态估计方法，首先，通过高斯拉普拉斯变换(Laplacian of Gaussian，LoG)来描述声呐回波图中小目标的边缘恒常性，以此解决小目标的非刚性特性造成的目标回波不稳定的问题；然后，针对舰船运动产生的尾流和空化噪声问题，引入基于上下文信息的自适应正则化参数对目标函数中的非局部项进行正则化，最后将经过升采样的声流场图像与边缘检测图像进行warping处理，得到连续的声呐回波图像中的小目标运动状态估计结果。本专利技术有效地减轻了水下运动小目标非刚性特性给运动估计精度带来的影响，同时提高了当声呐回波中存在船舶尾流噪声和空化噪声时的声流估计精度。
[0005]本专利技术解决其技术问题所采用的技术方案包括以下步骤：
[0006]步骤1：利用LoG算子对声呐回波图像序列中的所有声呐回波图像里的小目标进行边缘检测，得到第二声呐回波图像序列；
[0007]步骤2：由下式确定正则化参数λ
NL
：
[0008][0009]其中，FG
i,j
为声呐回波图像中非零像素点，H和L分别为声呐回波图像中的高度和长度，th为阈值，Num(FG
i,j
)表示计算声呐回波图像中非零像素点个数；
[0010]定义t＝1，定义图像分辨率参数c*d；
[0011]步骤3：从第二声呐回波图像序列取出第t帧和第t+1帧声呐回波图像，进行M 次降采样，直到第t帧和第t+1帧声呐回波图像的分辨率都变为c*d；此时，第t帧声呐回波图像和第t帧声呐回波图像在M次降采样过程中生成的M帧不同分辨率的声呐回波图像构成第t帧声呐回波图金字塔，第t帧声呐回波图金字塔共有M+1层，其中第t帧声呐回波图金字塔的第1层为第t帧声呐回波图像降采样最终得到的分辨率为c*d的声呐回波图像，按分辨率从小到大，第t帧声呐回波图金字塔的第M+1层为第二声呐回波图像序列的第t帧声呐回波图像；
[0012]第t+1帧声呐回波图像和第t+1帧声呐回波图像在M次降采样过程中生成的M帧不同分辨率的声呐回波图像构成第t+1帧声呐回波图金字塔，第t+1帧声呐回波图金字塔共有M+1层，其中第t+1帧声呐回波图金字塔的第1层为第t+1帧声呐回波图像降采样最终得到的分辨率为c*d的声呐回波图像，按分辨率从小到大，第t+1帧声呐回波图金字塔的第M+1层为第二声呐回波图像序列的第t+1帧声呐回波图像；
[0013]定义k＝1；
[0014]步骤4：使用第t帧声呐回波图金字塔的第k层声呐回波图像和第t+1帧声呐回波图金字塔的第k层声呐回波图像定义第k层声流目标函数E
Ak
(u
k
,v
k
)：
[0015][0016]其中，LoG
t
(
·
)为第t帧声呐回波图像边缘检测模块输出，LoG
t+1
(
·
)为第t+1帧声呐回波图像边缘检测模块输出；(u
k
,v
k
)为声流场图像，u
k
和v
k
分别是第k层声呐回波图像水平方向声流u和垂直方向声流v的集合，u
i,j
表示像素点(i,j)水平方向声流，u
i+1,j
表示像素点(i+1,j)水平方向声流，u
i,j+1
表示像素点(i,j+1)水平方向声流，v
i,j
表示像素点(i,j)垂直方向声流，v
i+1,j
表示像素点(i+1,j)垂直方向声流，v
i,j+1
表示像素点 (i,j+1)垂直方向声流，u
i',j'
表示像素点(i',j')水平方向声流，v
i',j'
表示像素点(i',j')垂直方向声流；λ为正则化参数；N
i,j
为点(i,j)的邻域点集合；为权重，表征点(i',j') 属于运动目标边界点的概率；ρ
D
(
·
)和ρ
S
(
·
)分别为数据项惩罚函数和空间项惩罚函数；
[0017]采用交替最小化法对式(3)和式(4)交替最小化实现对式(2)的最小化，本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种主动声呐回波图中小目标运动状态估计方法，其特征在于，包括以下步骤：步骤1：利用LoG算子对声呐回波图像序列中的所有声呐回波图像里的小目标进行边缘检测，得到第二声呐回波图像序列；步骤2：由下式确定正则化参数λ
NL
：其中，FG
i,j
为声呐回波图像中非零像素点，H和L分别为声呐回波图像中的高度和长度，th为阈值，Num(FG
i,j
)表示计算声呐回波图像中非零像素点个数；定义t＝1，定义图像分辨率参数c*d；步骤3：从第二声呐回波图像序列取出第t帧和第t+1帧声呐回波图像，进行M次降采样，直到第t帧和第t+1帧声呐回波图像的分辨率都变为c*d；此时，第t帧声呐回波图像和第t帧声呐回波图像在M次降采样过程中生成的M帧不同分辨率的声呐回波图像构成第t帧声呐回波图金字塔，第t帧声呐回波图金字塔共有M+1层，其中第t帧声呐回波图金字塔的第1层为第t帧声呐回波图像降采样最终得到的分辨率为c*d的声呐回波图像，按分辨率从小到大，第t帧声呐回波图金字塔的第M+1层为第二声呐回波图像序列的第t帧声呐回波图像；第t+1帧声呐回波图像和第t+1帧声呐回波图像在M次降采样过程中生成的M帧不同分辨率的声呐回波图像构成第t+1帧声呐回波图金字塔，第t+1帧声呐回波图金字塔共有M+1层，其中第t+1帧声呐回波图金字塔的第1层为第t+1帧声呐回波图像降采样最终得到的分辨率为c*d的声呐回波图像，按分辨率从小到大，第t+1帧声呐回波图金字塔的第M+1层为第二声呐回波图像序列的第t+1帧声呐回波图像；定义k＝1；步骤4：使用第t帧声呐回波图金字塔的第k层声呐回波图像和第t+1帧声呐回波图金字塔的第k层声呐回波图像定义第k层声流目标函数E
Ak
(u
k
,v
k
)：其中，LoG
t
(
·
)为第t帧声呐回波图像边缘检测模块输出，LoG
t+1
(
·
)为第t+1帧声呐回波图像边缘检测模块输出；(u
k
,v
k
)为声流场图像，u
k
和v
k
分别是第k层声呐回波图像水平方向声流u和垂直方向声流v的集合，u
i,j
表示像素点(i,j)水平方向声流，u
i+1,j...

【专利技术属性】
技术研发人员：韩一娜，魏卓群，赵爽，刘清宇，宋俊，
申请(专利权)人：西北工业大学，
类型：发明
国别省市：