一种极地声信号增强和识别方法

1.本发明属于水声信号处理、水声学和水声目标检测领域，涉及一种极地声信号增强和识别方法。


背景技术：

2.极地地理位置独特气候寒冷，海洋表面的浮冰不仅影响卫星通信，同时也隔绝了大气和海水的热交换，因此形成了特殊的声场环境，极地声信号中不仅包含自然的冰裂信号，还包括人类科考船航行、通信的信号以及动物活动的信号。但是在一次科考活动由于不同研究人员的研究目的不同，导致采集的信号往往是多种活动交织在一起，从中快速检测并提取到不同类别高质量的有效信号目前还有很大的难度，主要在于，包含不同研究人员所需的信号并不相同，且对于其需要的有效信号的拾取目前主要通过人工手段对海量的信号进行筛选，对于人工成本和时间成本都是极大的开销。因此，如果能使用自动手段对极地采集的声信号进行分析，从中消除自然冰裂信号造成的脉冲噪声进而识别出人类活动信号以及动物活动信号的片段，对于研究效率将具有极大的提升。


技术实现要素：

3.针对上述现有技术，本发明要解决的技术问题是提供一种极地声信号增强和识别方法，主要包括对采集到的极地声信号进行脉冲信号时频联合提取并降噪来提高信号质量，对采集到的极地声信号中检测出对研究有利有效的信号片段，从中进行时频联合提取并降噪并自动提取出由于不同活动而导致的信号片段，提升科研工作效率，克服目前仅通过人工手段对极地声信号进行筛选的问题。
4.本发明的目的是这样实现的：包括以下步骤：
5.步骤1：对所有极地信号进行小波变换得到二维时频图，然后所有信号片段进行时频联合提取，以得到极地信号中的脉冲噪声出现位置。
6.步骤2：首先对步骤1生成的二维时频图中出现的脉冲噪声位置使用空间平滑滤波。针对动物信号和通信信号持续时间有一定范围，将信号分割为若干个时间片段，然后对全部的时间片段执行如下操作：
7.步骤a：计算互相关积分列向量。对各时间片段的二维时频图，设置相关样本长度corrnum和时间窗长度segnum，满足segnum＞corrnum，corrnum表示每一个采样点计算与之相关的后续点个数，segnum表示每个时间片段被切分为若干个时间窗的长度，将时间片段按时间窗长度segnum划分为若干个时间窗，对每个时间窗执行以下操作：对时间窗内的所有采样点对应的时频图进行互相关计算，得到互相关矩阵seccoord
sup
，对矩阵seccoord
sup
沿行方向求和，得到时间窗内的积分列向量；
8.步骤b：求取时间片段的互相关积分曲线；
9.步骤c：根据设定的阈值判断时间片段的互相关积分曲线中是否存在动物信号或者通信信号，若互相关积分值在一段时间内大于设定的阈值，则判定为存在动物信号，并记
录信号持续时间；
10.步骤3：在步骤1基础上得到的脉冲噪声出现的位置附近一个窗长为w的时域片段进行中值滤波以消除脉冲噪声，然后将信号分为若干的时间片段，对所有片段进行带通滤波滤除低频段海洋环境噪声得到较为干净的中频端信号，然后对滤波后结果取绝对值计算以获得信号的包络，接下来对包络信号进行快速傅里叶变换，得到信号中的线谱成分，然后将线谱成分的经过梳状峰值滤波器后进行净化，然后计算净化后的线谱的方差值，若线谱方差值大于设定的阈值则说明该信号片段包含船舶信号，记录下该时频片段的位置，否则继续处理后续的信号片段。
11.本发明还包括这样一些结构特征：
12.1.步骤1中对极地脉冲声信号时频联合提取，其包括时频图片段中，选定fl-fh的频段fp，设定一个门限值pausethreshold，在fp频段中大于pausethreshold的值置为1，小于pausethreshold的功率值置为0，然后统计时频图片段在fp频段内每一列中1所占的比例为pexist，在所有信号片段中用所有时刻的pexist与时域幅度绝对值相乘得到的结果confidence作为判别脉冲信号是否存在的依据，若confidence值大于设定的阈值0.5则说明该位置处存在脉冲噪声，然后使用相同的方法得到其余所有信号片段的脉冲噪声位置。
13.2.步骤2中使用由步骤1得到的极地脉冲噪声位置对时频图对应位置进行使用水平方向空间平滑滤波器来滤除极地声信号时频图中的脉冲噪声，滤波器核长度为k，该滤波器核为
[0014][0015]
对于原时频图中的点(i,j),滤波器输出为
[0016][0017]
对于边缘部分，使用0进行填充。
[0018]
3.步骤a所述对经过脉冲噪声去除的信号时间窗内的所有采样点对应的时频图进行互相关计算，得到互相关矩阵seccoord
sup
包括：
[0019]
对时间窗内的所有采样点对应的功率谱进行互相关计算：
[0020][0021]
其中，ρ
xy
为第x列和第y列的互相关值，cov(x,y)代表时频图第x列和第y列的协方差，σ
x
代表时频图第x列的标准差，σy代表时频图第y列的标准差；ρ
xy
作为互相关矩阵seccoord第x行、第y列的元素，得到互相关矩阵seccoord，所述互相关矩阵seccoord为
segnum维方阵；取seccoord的上三角矩阵seccoord
utri
，每行保留从seccoord
utri
主对角线元素开始的前corrnum个元素取值，其他位置全部置0，令q＝segnum-corrnum+1，从第q+i开始，每行补充i个元素ρ
(q+a)(segnum+a)
得到矩阵seccoord
sup
，其中i＝1,2,3,
…
,corrnum-1，a＝1,2，
…
，i。
[0022]
4.步骤b所述求取时间片段的互相关积分曲线包括：
[0023]
将时间片段内全部时间窗的积分列向量依次连接组成一个列向量矩阵，所述列向量矩阵中第m行元素表示时间片段内第m个采样点对应的互相关值，然后将列向量矩阵转置，建立以时间为横坐标轴、以互相关值为纵坐标轴的二维坐标系，以采样点的采样时间为横坐标，以采样点的互相关值为纵坐标，得到时间片段各个采样点的互相关积分曲线。
[0024]
5.步骤3所述再对船舶噪声线谱进行提取之前，以步骤1中得到的脉冲噪声出现的位置为中心一个窗长为w的片段进行中值滤波以消除脉冲噪声，具体为对长度为w的信号片段进行排序，令出现脉冲位置的值为排序后的中位数。
[0025]
6.步骤3所述通过将线谱成分的经过梳状峰值滤波器后进行净化，然后计算净化后的线谱的方差值与设定阈值的比较，若该方差值大于设定的阈值则判定该信号片段为船舶信号，否则判断该片段不包含船舶信号。
[0026]
与现有技术相比，本发明的有益效果是：(1)本发明的方法能有效对极地脉冲声信号进行时频联合降噪，提高信号质量；(2)本发明的方法利用互相关积分算法对通信和动物信号进行提取，算法执行效率高；(3)本发明的方法利用船舶信号线谱的谐波特征来进行提取和净化，计算复杂度低。
附图说明
[0027]
图1是本发明的整体流程图；
[0028]
图2是降噪后对动物/通信信号的识别提取结果；
[0029]
图3是降噪后提取的舰船噪声线谱结果。
具体实施方式
[0030]
下面结合附图与具体实施方式对本发明作进一步详细描述。
[0031]
结合图1，本发明包含以下步骤：
[0032]
步骤1：对所有极地信号进行小波变换得到二维时频图，然后所有时间片段进行时频联合提取，以准确捕获极地信号中的脉冲噪声，具体步骤如下：
[0033]
步骤a：在步骤1得到的时频图片段中，选定fl-fh的频段fp，设定一个门限值pausethreshold，在fp频段中大于pausethreshold的值置为1，小于pausethreshold的功率值置为0，然后统计时频图片段在fp频段内每一列中1所占的比例为pexist；
[0034]
步骤b：在所有信号片段中用所有时刻的pexist与时域幅度绝对值相乘得到的结果confidence作为判别脉冲信号存在与否的依据。
[0035]
步骤c：根据confidence与设定的阈值做比较，来确定脉冲信号在所有时间点中出现的位置，若confidence大于设定的阈值则判断该位置出现脉冲噪声，否则判断该位置不存在脉冲噪声。
[0036]
步骤2：对步骤1生成的二维时频图中出现的脉冲噪声进行空间平滑滤波，以剔除
脉冲噪声的干扰,具体为使用水平方向大小为k的运动滤波器。
[0037]
针对动物信号和通信信号持续时间有一定范围，将信号分割为若干个时间片段，然后对全部的时间片段执行如下操作：
[0038]
步骤a：计算互相关积分列向量。对各时间片段的二维时频图，设置相关样本长度corrnum和时间窗长度segnum，满足segnum＞corrnum，corrnum表示每一个采样点计算与之相关的后续点个数，表示每个时间片段被切分为若干个时间窗的长度，将时间片段按时间窗长度划分为若干个时间窗，对每个时间窗执行以下操作：根据动物信号和通信信号在一个周期内的各帧功率谱是连续的且彼此具有强相关性的特点，对时间窗内的所有采样点对应的时频图进行互相关计算，得到互相关矩阵seccoord
sup
，对矩阵seccoord
sup
沿行方向求和，得到时间窗内的互相关积分列向量；
[0039]
步骤b：求取时间片段的互相关积分曲线；
[0040]
步骤c：根据设定的阈值与互相关积分曲线的值大小判断是否存在动物信号或者通信信号，若该阈值大于互相关积分曲线的值则判断该位置存在动物或通信信号，并记录该满足该阈值的持续时间，即为信号的持续时间；
[0041]
步骤3：对步骤1基础上得到的脉冲噪声出现的位置附近一个窗长为w的片段进行中值滤波以消除脉冲噪声，具体为对长度为w的信号片段进行排序，令出现脉冲位置的值为排序后的中位数，然后以一定窗长window进行处理，对该窗长内信号片段进行带通滤波滤除低频段海洋环境噪声得到较为干净的中频端信号，然后对滤波后结果取绝对值计算以获得信号的包络，接下来对包络信号进行快速傅里叶变换，得到信号中的线谱成分，然后再将线谱成分经过梳状滤波器后进行净化得到其基频以及各次谐波，然后计算净化后的线谱的方差值，并通过与设定的阈值比较来判断船舶信号的存在与否，若计算的方差值大于设定阈值则判定该信号片段为船舶信号，否则滑动到下一个window窗口继续执行步骤3，直到将所有数据处理完毕。
[0042]
下面结合具体参数数值给出本发明的实施例如下：
[0043]
步骤1：极地信号脉冲噪声提取。
[0044]
对信号进行小波变换得到二维时频图，频率分辨率为25hz。对于采集的极地信号脉冲噪声所在频段大多分布在5khz到15khz之间，而因此设置fl为5khz，fh为15khz,在该频段内对二维时频图上的使用脉冲功率门限pausethreshold为-90来对信号进行划分，小于门限的数据置为0，大于门限的数据置为1，此时再对该频段内所有的列进行统计，求出每一列中1所占的比例为pexist，然后再使用时频联合的方式进一步进行筛选，使用每个时域点上的幅度绝对值与时频图上对应的pexist相乘，得到confidence作为判断脉冲信号存在与否的依据，再用confidence与设定的阈值0.5进行比较，大于该阈值则说明该时间点存在脉冲噪声信号，小于该阈值说明该时间点不存在脉冲噪声信号。
[0045]
步骤2：对步骤1生成的二维时频图中出现的脉冲噪声进行空间平滑滤波，具体为使用长度为7的水平运动滤波器，以剔除脉冲噪声的干扰。滤波器核为
[0046][0047]
对于原时频图中的每点滤波器输出为
[0048][0049]
对于边缘部分，使用0进行填充。
[0050]
然后将信号分割为若干个时间片段，信号长度为n s，采样率为fs，然后对全部的时间片段执行如下操作：
[0051]
步骤a：
[0052]
计算互相关积分列向量。对每次取出的时间片段内的时频图,设置相关样本长度corrnum和时间窗长度segnum(segnum＞corrnum0，corrnum用以决定对于每一个采样点需要计算与之相关的后续点的个数，segnum用以判断时间窗的长度。依次取segnum个采样点集合secmeta(j0,j＝1,...,m,m＝2fs/segnum；之后对secmeta(j)中的所有采样点对应的功率谱进行互相关计算：
[0053][0054]
其中，ρ
xy
代表时频图的第x列和第y列的互相关值即pearson相关系数，cov(x,y)代表时频图第x列和第y列的协方差，σ
x
代表时频图第x列的标准差。之后得到该secmeta集合的互相关矩阵seccoord，其中seccoord为segnum维方阵，seccoord矩阵为：
[0055][0056]
互相关矩阵简化处理。seccoord为对称矩阵，为了简化计算，取其上三角矩阵，记为seccoord
utri
：
[0057][0058]
由于每一次计算只关心与每一个采样点临近的corrnum个采样点的相关性，因此，矩阵seccoord
utri
沿对角线方向每一行只保留前corrnum个值，其后的全部位置的点均被置
为0，记为seccoord
trap
：
[0059][0060]
但当corrnum+q-1＝segnum时，即q＝segnum-corrnum+1时，第q行以后(不含q行)的每行不足corrnum个值(共corrnum-1行不足corrnum个值)，故第q+1行补一个值为ρ
(q+1)(segnum+1)
，第q+2行补两个值为ρ
(q+2)(segnum+1)
和ρ
(q+2)(segnum+2)
，依此类推，第segnum行(最后一行)补corrnum-1个值。补值后的矩阵记为seccoord
sup
:
[0061][0062]
求取时间片段的互相关积分曲线。对于矩阵seccoord
sup
，沿行方向对其求和，得到积分向量
[0063]
之后对每个时间窗重复上述步骤，即对每一个secmeta(j),j＝1,2,3,...,2fs/segnum均计算积分向量ccavg(j)，并将所有ccavg(j)按顺序连接起来，得到第i个时间片段的总相关积分向量allccavg(i)，allccavg(i)的维度是1列，2fs行。该向量是当前时间片段内所有采样点对应的互相关值，即第一行代表第一个采样点的互相关值，然后将该向量转置，按照每个采样点的顺序依次将互相关值画在二维坐标系上，得。
[0064]
到当前时间片段各个时间点的互相关积分曲线
[0065]
步骤c：对于由总相关积分向量allccavg(i)按时间点形成的互相关积分曲线，通过设置合理的阈值shd对其进行判决，根据设置好的阈值，自动从互相关积分曲线中判断是否存在动物信号和通信信号，并计算其起止时间，一个时间点的互相关积分结果大于阈值，则判断为当前时间点有动物信号和通信信号，在大于阈值之前的那个重合点即为动物信号和通信信号的起始时间，记为t
start
；一个时间点的互相关积分结果小于阈值，则判断为当前时间点动物信号和通信信号消失，在小于阈值之前的那个重合点即为动物信号和通信信号的终止时间，记为t
end
，由于每次切片时间为2秒，当前时间片段动物信号和通信信号会出现
有始无终的情况，此时，仅记录动物信号和通信信号起点，当处理下个时间片段时记录动物信号和通信信号终点。之后对每个时间片段进行以上的步骤，得到全部时间信号的起止时间，并记录保存为信号提取结果如图2。
[0066]
步骤3：对步骤1基础上得到的脉冲噪声出现的位置左右长度为10点的片段进行中值滤波以消除脉冲噪声，然后以长度为5s，重叠度为50％的窗口进行滑动，对所有信号片段进行5-10khz带通滤波滤除低频段海洋环境噪声得到较为干净的中频端信号，然后对滤波后结果取绝对值计算以获得信号的包络，接下来对包络信号进行快速傅里叶变换，得到信号中的线谱成分，然后再将线谱成分经过10阶梳状峰值滤波器后进行净化得到其基频以及各次谐波，梳状滤波器的中心频率ω0为线谱的基频，带宽bw设置为归一化频率0.013，品质因数为35，最后计算净化后的线谱的方差值，若该值大于设定的阈值0.7则判定该段信号为船舶信号，线谱的提取的结果如图3所示。
[0067]
综上，本发明提供一种极地声信号增强和识别方法，对采集到的极地声信号的数据进行时频联合降噪并提取信号片段。本发明属于水声信号处理领域。本发明是这样实现的：步骤1：对极地声信号中的脉冲干扰噪声进行时频联合提取，得到所有脉冲噪声出现的位置；步骤2：对步骤1中获得的脉冲噪声位置在信号二维时频图上进行空间平滑滤波以去除脉冲噪声，然后用在时频图上做互相关得到互相关积分曲线，若互相关曲线积分值大于设定的阈值则判断动物/通信信号存在；步骤3：以对步骤1获得的脉冲噪声位置为中心进行时域窗口的中值滤波，然后对信号提取线谱成分，并通过梳状滤波器进行净化，最后通过比较净化后线谱的方差值，若大于设定的方差阈值则判断船舶信号存在。本发明实现针对在极地中进行长期声呐信号采集获得的海量数据的自动增强和提取，降低对人工干预的依赖，提升效率。本发明的优点在于：(1)针对极地脉冲信号进行时频联合降噪，增强信号质量；(2)利用互相关算法对通信和动物信号进行提取，算法执行效率高；(3)利用船舶信号线谱的谐波特征来进行提取和净化，计算复杂度低。