CN112924956A - 一种降低低频噪声级的背景均衡方法 - Google Patents
一种降低低频噪声级的背景均衡方法 Download PDFInfo
- Publication number
- CN112924956A CN112924956A CN202110161326.2A CN202110161326A CN112924956A CN 112924956 A CN112924956 A CN 112924956A CN 202110161326 A CN202110161326 A CN 202110161326A CN 112924956 A CN112924956 A CN 112924956A
- Authority
- CN
- China
- Prior art keywords
- data
- frequency
- time
- valley
- peak
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/52—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
- G01S7/537—Counter-measures or counter-counter-measures, e.g. jamming, anti-jamming
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Abstract
本发明属于声纳信号显示技术领域,涉及一种降低低频噪声级的背景均衡方法,该方法设计了一种新的背景均衡模式,通过分段搜索时频图图像中的峰值和谷值并与预设门限进行比较,然后通过逐点最小值处理输出背景均衡后的时频图。本发明的有益结果是:1.本发明所提方法不需要人工设置经验参数;2.本发明最终输出的背景均衡后的时频图相比S3PM法、OTA法和波束扫描特征法在各频点具有更低的平均噪声级,尤其在低频段对噪声的抑制效果更好,使得时频图中的线谱能够被更有效地进行检测分析。
Description
技术领域
本发明属于声纳信号显示技术领域,涉及一种降低低频噪声级的背景均衡方法。
背景技术
声纳装备检测水下目标的一种重要手段是将水下潜艇、鱼雷等目标的辐射噪声作为接收信号,然后通过检测分析接收信号的时间—频率图像(以下简称时频图)中各个频率分量上的“线谱”(即在某段时间内出现的该频率分量上的亮线,表示该频率分量在这段时间的幅度较高),来对目标的有无和类型进行判断。通常来说,水下目标辐射噪声的线谱分量主要集中在较低频段(若无特别说明,后面所指的低频段一般为1000Hz以下频段),分布于低频段不同频率分量的线谱是不同类型目标区别于环境噪声的主要特征之一。然而由于水下环境噪声等因素的影响,时频图中各个频率分量中均会存在幅度强弱不一的噪声干扰,形成非均匀的干扰背景,严重影响低频段的线谱检测分析。另一方面,声纳接收设备中通常存在较强的低频自噪声,这使得时频图中低频段的噪声干扰相对其他频段更加明显。
背景均衡作为一种声纳信号显示技术,能降低背景噪声强度的动态变化程度,实现所谓的恒虚警检测。其主要目的是为了减弱非均匀干扰背景的影响,使数据的动态范围适应图形显示要求。具体是通过对时频图图像进行处理,使非均匀干扰背景变成均匀干扰背景,即各个频率分量中的噪声干扰幅度变得强弱差不多。这样时频图中就不会出现成片的强噪声干扰,得到较为干净的显示画面,从而有效检测分析各频率分量中的线谱,提高对目标的检测能力。
目前常用的背景均衡方法包括双边分离窗法(two-pass split-windownormalizer,S3PM)(S.Stergiopoulo.Noise normalization technique for beamformedtowed array data[J].The Journal of the Acoustical Society of America,1995,97(4):2334-2345.)、排序截断平均法(order truncate average,OTA)(郑佼.被动声纳显示优化技术研究[D].哈尔滨工程大学硕士学位论文,2013.)和波束扫描特征法(郑佼.被动声纳显示优化技术研究[D].哈尔滨工程大学硕士学位论文,2013.)等。
上述方法中,S3PM法需要用到三个经验参数:双边分离窗的绝对值上下限G与E,以及门限常数p;OTA法需要用到两个经验参数:窗口长度K和门限常数α;波束扫描特征法不需要用到经验参数。当经验参数设置合适时,这三种方法对同一幅时频图图像进行处理后的平均噪声级排序为:OTA法最好(平均噪声级最低),波束扫描特征法其次,S3PM法再次。但是,尽管这三种方法可以有效减少时频图中的非均匀干扰背景,但它们均不能较好地抑制由声纳接收设备产生的低频噪声;此外,OTA法和S3PM法的性能较多地依赖于人工设置的经验参数,如果经验参数设置不合适,这两种方法的效果会大打折扣。
发明内容
本发明的目的是:提出一种降低低频噪声级的背景均衡方法,该方法设计了一种新的背景均衡模式,通过分段搜索时频图图像中的峰值和谷值并将获得的两组新序列进行逐点比较,然后通过逐点最小值处理输出背景均衡后的时频图。
本发明的技术方案包括以下步骤:
第一步,取出时频图中某一时刻t对应的数据行向量,t=1,2,...,T,T表示时频图在时间维上的像素点数;若对应数据为列向量则将其转置为行向量。
第二步,搜索得到数据行向量中的所有峰值和所有谷值。
第三步,对于某一峰值,搜索找到它左右相邻距离其第二近的谷值,然后截取这两个谷值之间的数据,记为x(1),x(2),...,x(M),共M点,并记这段数据中最大的谷值和最小的谷值分别为xmin1和xmin2,最大的峰值为xmax。
第六步,重复第三步到第五步直至遍历所有峰值,对所有新数据序列拼接得到完整的新数据序列z1(1),z1(2),...,z1(N)(N为时频图在频率维上的数据点数)。
第七步,对于某一谷值,搜索找到它左右相邻距离其第二近的峰值,然后截取这两个峰值之间的数据,记为y(1),y(2),...,y(Q),共Q点,并记这段数据中最大的峰值和最小的峰值分别为ymax1和ymax2,最小的谷值为ymin。
第十步,重复第七步到第九步直至遍历所有谷值,对所有新数据序列拼接得到完整的新数据序列z2(1),z2(2),...,z2(N),N为时频图在频率维上的数据点数。
第十一步,将第六步和第十步得到的两组完整的新数据序列进行逐点最小值处理,得到处理后的完整新数据序列:
z(n)=min[z1(n),z2(n)],n=1,2,...,N;
第十二步,重复第一步到第十一步直至遍历所有时刻,对所有的处理后的完整新数据序列沿时刻维进行堆叠得到背景均衡后的时频图。
本发明的有益结果是:
1.本发明所提方法不需要人工设置经验参数,减少了所提方法的人工干预程度,使方法性能更加稳定。
2.本发明最终输出的背景均衡后的时频图相比S3PM法、OTA法和波束扫描特征法在各频点具有更低的平均噪声级,尤其在低频段对噪声的抑制效果更好,使得时频图中的线谱能够被更有效地进行检测分析。
附图说明
图1是本发明的实现流程图;
图2是利用本发明的具体实施方式进行实测数据试验的结果:(a)为背景均衡处理前的时频图;(b)为利用S3PM法对图2(a)的处理结果;(c)为利用OTA法对图2(a)的处理结果;(d)为利用波束扫描特征法对图2(a)的处理结果;(e)为利用本发明所提方法对图2(a)的处理结果;(f)为背景均衡处理前和采用4种方法(S3PM法,OTA法,波束扫描特征法,本发明方法)进行背景均衡处理后各频点上平均噪声级的对比结果;(g)为图2(f)在0到500Hz的放大图。
具体实施方式
图1是本发明的实现流程图。具体实施方式包括以下步骤:
第一步,取出时频图中某一时刻t(t=1,2,...,T,T表示时频图在时间维上的像素点数)对应的数据行向量(若对应数据为列向量则将其转置为行向量)。
第二步,搜索得到数据行向量中的所有峰值和所有谷值。
第三步,对于某一峰值,搜索找到它左右相邻距离其第二近的谷值,然后截取这两个谷值之间的数据,记为x(1),x(2),...,x(M),共M点,并记这段数据中最大的谷值和最小的谷值分别为xmin1和xmin2,最大的峰值为xmax。
如果该峰值的左边仅搜索到1个或0个谷值,则x(1)取第一步中数据行向量的第1个值;如果该峰值的右边仅搜索到1个或0个谷值,则x(M)取第一步中数据行向量的最后1个值。
第六步,重复第三步到第五步直至遍历所有峰值,对所有新数据序列拼接得到完整的新数据序列z1(1),z1(2),...,z1(N)(N为时频图在频率维上的数据点数)。
第七步,对于某一谷值,搜索找到它左右相邻距离其第二近的峰值,然后截取这两个峰值之间的数据,记为y(1),y(2),...,y(Q),共Q点,并记这段数据中最大的峰值和最小的峰值分别为ymax1和ymax2,最小的谷值为ymin。
如果该谷值的左边仅搜索到1个或0个峰值,则y(1)取第一步中数据行向量的第1个值;如果该谷值的右边仅搜索到1个或0个峰值,则y(Q)取第一步中数据行向量的最后1个值。
第十步,重复第七步到第九步直至遍历所有谷值,对所有新数据序列拼接得到完整的新数据序列z2(1),z2(2),...,z2(N)(N为时频图在频率维上的数据点数)。
第十一步,将第六步和第十步得到的两组完整的新数据序列进行逐点最小值处理,得到处理后的完整新数据序列:
z(n)=min[z1(n),z2(n)],n=1,2,...,N
第十二步,重复第一步到第十一步直至遍历所有时刻,对所有的处理后的完整新数据序列进行堆叠得到背景均衡后的时频图。
图2是利用本发明的具体实施方式进行实测数据试验的结果,其中图2(a)到图2(e)的横坐标表示频率,单位为“赫兹(Hz)”,纵坐标表示时间,单位为“秒(s)”,色条颜色表示图中像素点的噪声级,单位为“dB”;图2(f)和图2(g)的横坐标表示频率,单位为“赫兹(Hz)”,纵坐标表示平均噪声级,单位为“dB”。
图2(a)为背景均衡处理前的时频图,可以发现图中500Hz以下频段存在若干条线谱,但有不少被淹没在了低频段(尤其是在50Hz以下频段)的噪声干扰当中;此外,图中在第50到200s之间还存在3条由贯穿全频段的噪声干扰形成的亮线。我们的目的是通过背景均衡方法,使得时频图中的线谱可以被更有效地检测分析,同时降低低频段的噪声干扰。图2(b)到图2(e)分别是利用S3PM法、OTA法、波束扫描特征法和本发明所提方法对图2(a)的处理结果。经验参数按照上述文献中设置的值进行取值,即S3PM法中,G=16,E=3,p=1.14;OTA法中,K=25,α=1.05。可以发现4种方法的处理结果均比处理前的时频图更容易观察到被低频噪声干扰淹没的线谱(如图中第210到240s处频率为120Hz的线谱,以及若干条1000Hz及更高频率的线谱),同时4种方法均基本消除了图中在第50到200s之间存在的3条由贯穿全频段的噪声干扰形成的亮线。图2(f)对比了背景均衡处理前和采用4种方法处理后各频点上的平均噪声级,图2(g)是其在0到500Hz的放大图,结果表明本发明所提方法具有比前述3种方法更低的平均噪声级,尤其是在50Hz以下频段有明显降低(以30Hz处的平均噪声级为例,本发明方法相比于前述3种方法中最优的OTA法进一步降低了10dB左右)。
Claims (3)
1.一种降低低频噪声级的背景均衡方法,其特征在于,该方法包括以下步骤:
第一步,取出时频图中某一时刻t对应的数据行向量,t=1,2,...,T,T表示时频图在时间维上的像素点数;若对应数据为列向量则将其转置为行向量;
第二步,搜索得到数据行向量中的所有峰值和所有谷值;
第三步,对于某一峰值,搜索找到它左右相邻距离其第二近的谷值,然后截取这两个谷值之间的数据,记为x(1),x(2),...,x(M),共M点,并记这段数据中最大的谷值和最小的谷值分别为xmin1和xmin2,最大的峰值为xmax;
第六步,重复第三步到第五步直至遍历所有峰值,对所有新数据序列拼接得到完整的新数据序列z1(1),z1(2),...,z1(N),N为时频图在频率维上的数据点数;
第七步,对于某一谷值,搜索找到它左右相邻距离其第二近的峰值,然后截取这两个峰值之间的数据,记为y(1),y(2),...,y(Q),共Q点,并记这段数据中最大的峰值和最小的峰值分别为ymax1和ymax2,最小的谷值为ymin;
第十步,重复第七步到第九步直至遍历所有谷值,对所有新数据序列拼接得到完整的新数据序列z2(1),z2(2),...,z2(N),N为时频图在频率维上的数据点数;
第十一步,将第六步和第十步得到的两组完整的新数据序列进行逐点最小值处理,得到处理后的完整新数据序列:
z(n)=min[z1(n),z2(n)],n=1,2,...,N;
第十二步,重复第一步到第十一步直至遍历所有时刻,对所有的处理后的完整新数据序列进行堆叠得到背景均衡后的时频图。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110161326.2A CN112924956B (zh) | 2021-02-05 | 2021-02-05 | 一种降低低频噪声级的背景均衡方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110161326.2A CN112924956B (zh) | 2021-02-05 | 2021-02-05 | 一种降低低频噪声级的背景均衡方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112924956A true CN112924956A (zh) | 2021-06-08 |
CN112924956B CN112924956B (zh) | 2022-05-13 |
Family
ID=76170775
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110161326.2A Active CN112924956B (zh) | 2021-02-05 | 2021-02-05 | 一种降低低频噪声级的背景均衡方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112924956B (zh) |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101430759A (zh) * | 2008-12-04 | 2009-05-13 | 上海大学 | 优化的人脸识别预处理方法 |
WO2013140910A1 (ja) * | 2012-03-22 | 2013-09-26 | 古野電気株式会社 | 信号サーチ方法、信号サーチプログラム、信号サーチ装置、gnss信号受信装置、および情報機器端末 |
CN105353373A (zh) * | 2015-12-16 | 2016-02-24 | 武汉大学 | 一种基于Hough变换探地雷达目标提取方法和装置 |
US20180174274A1 (en) * | 2016-05-19 | 2018-06-21 | Shenzhen University | Method and System for Eliminating Stripe Noise in Infrared Images |
CN109102479A (zh) * | 2018-06-29 | 2018-12-28 | 中国船舶重工集团公司第七〇五研究所 | 一种新图像算子的声纳目标增强方法 |
CN109117698A (zh) * | 2017-12-27 | 2019-01-01 | 南京世海声学科技有限公司 | 一种基于最小均方误差准则的噪声背景估计方法 |
CN111693954A (zh) * | 2020-05-08 | 2020-09-22 | 中国人民解放军空军工程大学 | 空间旋转群目标微多普勒分离方法 |
-
2021
- 2021-02-05 CN CN202110161326.2A patent/CN112924956B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101430759A (zh) * | 2008-12-04 | 2009-05-13 | 上海大学 | 优化的人脸识别预处理方法 |
WO2013140910A1 (ja) * | 2012-03-22 | 2013-09-26 | 古野電気株式会社 | 信号サーチ方法、信号サーチプログラム、信号サーチ装置、gnss信号受信装置、および情報機器端末 |
CN105353373A (zh) * | 2015-12-16 | 2016-02-24 | 武汉大学 | 一种基于Hough变换探地雷达目标提取方法和装置 |
US20180174274A1 (en) * | 2016-05-19 | 2018-06-21 | Shenzhen University | Method and System for Eliminating Stripe Noise in Infrared Images |
CN109117698A (zh) * | 2017-12-27 | 2019-01-01 | 南京世海声学科技有限公司 | 一种基于最小均方误差准则的噪声背景估计方法 |
CN109102479A (zh) * | 2018-06-29 | 2018-12-28 | 中国船舶重工集团公司第七〇五研究所 | 一种新图像算子的声纳目标增强方法 |
CN111693954A (zh) * | 2020-05-08 | 2020-09-22 | 中国人民解放军空军工程大学 | 空间旋转群目标微多普勒分离方法 |
Non-Patent Citations (5)
Title |
---|
JIAHUA ZHU等: "Range sidelobe suppression for using Golay complementary waveforms in multiple moving target detection", 《SIGNAL PROCESSING》 * |
李伟 等: "基于组合滤波设计的背景噪声归一化方法", 《火力与指挥控制》 * |
王晓宇 等: "一种新的宽带声呐波束域背景均衡方法", 《航海工程》 * |
王疆一 等: "基于多频段滤波的舰船噪声DEMON谱分析", 《声学技术》 * |
邱家兴 等: "一种自适应门限时间方位历程显示背景均衡算法", 《舰船科学技术》 * |
Also Published As
Publication number | Publication date |
---|---|
CN112924956B (zh) | 2022-05-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107609488B (zh) | 一种基于深度卷积网络的舰船噪声识别分类方法 | |
CN114937055B (zh) | 基于人工智能的图像自适应分割方法与系统 | |
CN111476809A (zh) | 一种侧扫声呐图像目标识别方法 | |
CN116012346B (zh) | 一种输电线路监测方法及系统 | |
CN109102479A (zh) | 一种新图像算子的声纳目标增强方法 | |
CN109785260B (zh) | 一种针对侧扫声呐图像的非线性增强方法 | |
CN110675410A (zh) | 基于选择性搜索算法的侧扫声呐沉船目标非监督探测方法 | |
CN112528775A (zh) | 一种水下目标的分类方法 | |
CN115661649B (zh) | 一种基于bp神经网络的船载微波雷达图像溢油探测方法及系统 | |
CN113567969B (zh) | 一种基于水下声信号的非法采砂船自动监控方法及系统 | |
CN109782261A (zh) | 对水下入侵目标引起接收信号变化的无监督学习检测方法 | |
CN112380961A (zh) | 一种基于人工智能的检测气泡流型评估气密性方法及系统 | |
CN112924956B (zh) | 一种降低低频噪声级的背景均衡方法 | |
CN107783096A (zh) | 一种用于方位历程图显示的二维背景均衡方法 | |
Han et al. | A novel sonar image denoising algorithm based on block matching | |
CN115187855A (zh) | 一种海底底质声呐图像分类方法 | |
CN114841218A (zh) | 一种语音啸叫检测算法 | |
RU2736188C1 (ru) | Способ отображения гидроакустической информации | |
CN114298094A (zh) | 一种基于主成分分析的线谱自动提取方法 | |
CN113516187A (zh) | 一种采用局部特征对比度的红外弱小目标检测算法 | |
Jia et al. | A piecewise nonlinear enhancement method of side scan sonar images | |
Ye et al. | A Gray Scale Correction Method for Side-Scan Sonar Images Based on GAN | |
CN109886899A (zh) | 一种受动态干扰的海天线图像检测方法 | |
Marta et al. | The Comparison of Image Enhancement Methods for A Realtime Under Water Vision System of The Break Water Implementation | |
Manek et al. | A soft weighted median filter for removing general purpose image noise |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |