基于曲率和的主动声呐干涉条纹特征提取方法

专利检索2022-05-10  11



1.本发明属于水声学和水声信号处理领域,涉及主动声呐干涉条纹的特征提取方法,特别涉及基于曲率和的主动声呐干涉条纹特征提取方法。


背景技术:

2.声波在浅水域中传播是一种复杂的现象,受水中噪声及海底底部属性的影响,其传播模式之间相互干扰,声波与信道边界之间相互作用导致信号衰减和多径传播,这种传播产生的时频或时空结构可生成恒定强度的干涉条纹。在主动声呐研究中,由于目标散射回波的多路径效应,不同简正波相互叠加,给干涉条纹的提取造成困难,原因在于条纹之间的宽度与频率、声源距离等有关,此外,由于干涉条纹谱图是由接收器处的声场强度得到的,因此还要求声源频率带为均匀的,实验条件下条纹生成与特征提取更为困难。
3.主动声呐干涉条纹的研究可作为一种水声探测的重要研究方向,是水声领域的一个关键技术。目前关于主动声呐干涉条纹特征提取理论基于jorge e.quijano的实验分析,其根据均方误差的方法统计出条纹特征,不足之处在于只能得到一个范围内的干涉条纹特征值,无法精确感知特征变化。


技术实现要素:

4.针对现有技术的不足,本发明提出了基于曲率和的主动声呐干涉条纹特征提取方法,通过hough曲线检测法得到条纹谱图中曲线,然后基于最小二乘拟合曲线法得到曲线模型,最后通过累加曲率和的方式,精确计算每束条纹的特征值,解决现有技术中主动声呐的条纹特征表征的问题。
5.基于曲率和的主动声呐干涉条纹特征提取方法,具体包括以下步骤:
6.步骤一、采集主动声呐干涉条纹谱图
7.在深度为z
s
的位置放置主动声呐的声源与接收器,声源发射出的宽带脉冲信号经过浅海波导传播后,经位于深度z
t
的目标散射后,被位于深度z
r
的接收器接收;建立接收器位置的声压场模型:
[0008][0009]
其中,c为常数,r1、r2分别为声源到目标、目标到接收器的距离;ω是声波频率,表示散射矩阵,a
m
、a
n
分别表示入射波的入射角度和反射波的反射角度,分别表示入射声波的方位角和反射声波的方位角;m为入射信号的信道个数,n为目标反射后的信道个数,m>n;k
m
是第m个信道传播的水平波束,k
n
是第n个信道传播的水平波束;ψ
m
是目标相关的深度相关模式函数,其中ψ
m
(z
s
)是目标深度相关模式函数,ψ
m
(z
t
)是发射声源深度相关模式函数;ψ
n
是反射模态相关的模式函数,其中ψ
n
(z
t
)是反射模态中目标深度相关模式函数,ψ
n
(z
r
)是反射模态中接收器深度相关模式函数。
[0010]
步骤二、hough变换曲线检测
[0011]
对步骤一得到的主动声呐干涉条纹谱图进行hough变换,通过边缘检测算法,将主动声呐干涉条纹谱图中的曲线转换到参数空间中的hough矩阵,在参数空间中设置门限,提取大于门限的点,得到空间中的拟合曲线。
[0012]
步骤三、最小二乘法曲线拟合
[0013]
对步骤二通过hough变换得到的检测曲线,运用最小二乘法建立一个新的函数模型,来逼近检测曲线上的已知样本点得到拟合曲线。{(x
i
,y
i
):i=1,2,3,....,n}为检测曲线上的一组数据点,拟合曲线的模型为两条曲线在第i个数据点处的距离误差为则第i个实际给定数据点y
i
与拟合曲线模型之差的平方和s为:
[0014][0015]
其中为拟合曲线上的第i个点。求解令平方和s取值最小的参数,得到拟合曲线
[0016]
步骤四、基于曲率和的干涉条纹特征提取
[0017]
对于步骤三得到的拟合曲线求解对应的曲率和,作为干涉条纹的特征值。
[0018]
曲率和的特征提取方法为:定义法曲率k
n
,得到主曲率的方程为:
[0019]
(eg

f2)k
n2

(en

2fm gl)k
n
(ln

m2)=0
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(3)
[0020]
其中e、f、g为第一基本形式i的系数,l、m、n为第二基本形式ii的系数。
[0021]
设du:dv是曲面上一点p(u,v)处的任意一个切方向,该切方向与u线的夹角记为θ,根据曲面上两条曲线夹角公式得到:
[0022][0023]
又因为p点处各方向的法曲率公式为:
[0024][0025]
将式子(4)带入式子(5)可得:
[0026][0027]
式子(6)为欧拉公式,反映了法曲率随方向变化的规律,因为高斯曲率k=k1k2,平均曲率为根据式子(3)和式子(6),由韦达定理可得:
[0028][0029]
[0030]
根据曲率和公式求解条纹特征值:
[0031][0032]
m为拟合曲线上点的个数,m为采样最大点数,g(x,y,σ)为高斯标准差的高斯函数:
[0033][0034]
其中σ为高斯标准差。
[0035]
本发明具有以下有益效果:
[0036]
通过计算拟合条纹曲线的曲率得到条纹特征,计算过程中用到一阶和二阶高斯导数滤波器,具有良好的稳定性和鲁棒性,充分的利用了微分几何中曲率不变的特性,能够同时提现条纹谱图的宏观和微观特征信息,最后利用曲率和的概念表示条纹特征信息,可作为主动声呐波导不变量的一种衡量标准,为主动声呐干涉条纹特征提取打下基础。
附图说明
[0037]
图1是实施例使用的模型布放装置;
[0038]
图2是采集到的主动声呐干涉条纹谱图;
[0039]
图3是经过hough变换的曲线检测结果图;
[0040]
图4最小二乘拟合曲线模型函数;
[0041]
图5是通过曲率和求解得到的条纹特征值。
具体实施方式
[0042]
以下结合附图对本发明作进一步的解释说明;
[0043]
本实例在理想波导条件下进行仿真实验,目标及水听器布放位置如图1所示,其中,水深h=15m,水中传播声速与密度分别为c=1500m/s,ρ=1.0g/cm3,声源及水听器深度为9m,目标距水面深度为12m,在距声源水平距离3m的平面上移动,目标移动过程中,记录接收器处声压场变化情况。根据公式得出波导条件下条纹谱图,如图2所示,由弯曲的曲线构成,要得到条纹特征需要对条纹曲线进行处理。
[0044]
hough变换可检测曲线条纹,使主动声呐干涉条纹谱图更为清晰,其中图3为hough变换检测结果,原理在于利用图像空间的点和参数空间的线之间存在的对偶关系,将图像空间的问题转化到参数空间来解决,直线上每个坐标点(x,y)对应到参数坐标系中的一条曲线。而两个不同的点对应两条相交曲线,曲线的交点可以唯一确定一条直线的参数,因此参数坐标系中最多曲线经过的交点坐标就是直线所具有的参数。获取该点的方法是通过对每个不同的(θ,ρ)采用一个累加器进行统计,将最终得到统计量最大的(θ,ρ)作为直线参数直角坐标系中的一点,所以取直角坐标系上曲线的切点进行计算,在极坐标系可形成类似曲线的映射,通过设置门限来对极坐标系中的值进行提取,可得到hough变换后的计算值。
[0045]
图4是最小二乘拟合曲线的结果,选取条纹谱曲线散点,通过最小二乘法进行曲线拟合,借助于matlab的cftool工具箱,可以自适应拟合曲线函数,在曲线拟合强度极强相关,其中拟合函数为:
[0046]
y=814.4*x
(

0.168)
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(10)
[0047]
拟合结果的相关系数为0.983,偏离程度较低,选取多束曲线进行拟合,得到一组拟合曲线,在通过曲率求和的方法得到条纹特征值。
[0048]
由拟合曲线函数模型可进行曲率求解,如图5所示为曲率和求解结果,可作为条纹特征,由公式(10)可知,通过求解一阶导数和二阶导数,得到曲线函数一点处的曲率,其中高斯标准差σ=1将曲线分成很多细小的部分,每个部分可看做一段圆弧,圆弧的弯曲程度越大,其曲率半径越小,而曲率就是曲率半径的倒数,本方法采用的曲率求解方法具有良好的稳定性和鲁棒性,可以将不同参数下的图像信息进行融合,以形成多分辨率信息提取,根据推导公式过程可知,高斯二阶导数滤波器滤波后的结果与高斯标准差σ相关联,这里利用多个不同的值和σ达到提高分辨率的效果,采用的高斯滤波器在x轴和y轴方向是可分离的,因此二阶导数能直接突出曲率变化快的区域。
[0049]
通过求取曲线曲率和作为条纹曲线特征,曲率是用来反映几何体弯曲程度的量,可作为条纹特征的衡量指标,提取出的条纹特征值可作为主动声呐波导不变量的重要参考,本发明通过计算拟合条纹曲线的曲率得到条纹特征,计算过程中利用hough变换检测干涉条纹曲线,最小二乘拟合得到曲线模型,最后求得条纹曲线的曲率值,用到一阶和二阶高斯导数滤波器,具有良好的稳定性和鲁棒性,充分的利用了微分几何中曲率不变的特性,能够同时提现条纹谱图的宏观和微观特征信息,最后利用曲率和的概念表示条纹特征信息,可作为主动声呐波导不变量的一种衡量标准,为主动声呐干涉条纹特征提取打下基础。
转载请注明原文地址:https://win.8miu.com/read-550165.html

最新回复(0)