CN111413702B - 用于宽带探鱼仪的高效目标分割方法 - Google Patents
用于宽带探鱼仪的高效目标分割方法 Download PDFInfo
- Publication number
- CN111413702B CN111413702B CN202010401601.9A CN202010401601A CN111413702B CN 111413702 B CN111413702 B CN 111413702B CN 202010401601 A CN202010401601 A CN 202010401601A CN 111413702 B CN111413702 B CN 111413702B
- Authority
- CN
- China
- Prior art keywords
- target
- wave value
- fish finder
- value data
- data set
- 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.)
- Active
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
- G01S15/00—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
- G01S15/88—Sonar systems specially adapted for specific applications
- G01S15/96—Sonar systems specially adapted for specific applications for locating fish
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Acoustics & Sound (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Abstract
本发明提供一种用于宽带探鱼仪的高效目标分割方法,包括步骤:S1:根据一背景噪声值确定一检测阈值;S2:设定一滑动窗的宽度,滑动窗的宽度小于等于一探鱼仪的距离分辨率;S3:利用一匹配滤波器对探鱼仪的输出信号进行匹配滤波,获得一滤波输出信号的离散滤波数据集,所述离散滤波数据集包括多个波值数据;S4:利用滑动窗和检测阈值对离散滤波数据集分割目标,并输出目标的位置。本发明的一种用于宽带探鱼仪的高效目标分割方法,运算量小,当两个目标回波重叠时,只要目标间隔大于探鱼仪的最小距离分辨率,就可将两个目标区分开来,并快速获得目标的数量和位置。
Description
技术领域
本发明涉及探鱼技术领域,尤其涉及一种用于宽带探鱼仪的高效目标分割方法。
背景技术
渔业资源声学评估中会使用到回波计数和回波积分方法。在单体间距大于探鱼仪距离分辨力的情况下可以直接使用回波计数的方法,当密度较高时使用回波积分方法,由于回波积分需要获得单体的平均目标强度,因此这两种方法都需要对目标进行距离上的分割。
目前距离分割使用的方法是设置一个检测阈值PLDL,从回波信号强度高于检测阈值时开始到强度低于阈值时结束,统计该脉冲的长度,如果该长度在设定的范围则认为找到了1个目标,目标的距离为峰值出现的位置。但现有方法当两个目标回波重叠时即使目标间隔大于探鱼仪分辨率,探鱼仪也无法区分。
发明内容
针对上述现有技术中的不足,本发明提供一种用于宽带探鱼仪的高效目标分割方法,运算量小,当两个目标回波重叠时,只要目标间隔大于探鱼仪的最小距离分辨率,就可将两个目标区分开来,并快速获得目标的数量和位置。
为了实现上述目的,本发明提供一种用于宽带探鱼仪的高效目标分割方法,包括步骤:
S1:根据一背景噪声值确定一检测阈值,所述检测阈值大于所述背景噪声值;
S2:设定一滑动窗的宽度,所述滑动窗的宽度小于等于一探鱼仪的距离分辨率;
S3:利用一匹配滤波器对所述探鱼仪的输出信号进行匹配滤波,获得一滤波输出信号的离散滤波数据集,所述离散滤波数据集包括多个波值数据;
S4:利用所述滑动窗和所述检测阈值对所述离散滤波数据集分割目标,并输出所述目标的位置。
优选地,所述背景噪声值通过停止所述探鱼仪发射并记录所述探鱼仪的回波强度估算获得,当无法停止所述探鱼仪发射时利用超出海底范围或远距离收到的回波强度进行估算。
优选地,所述滑动窗的宽度M满足公式(1):
其中,ΔR表示所述探鱼仪的距离分辨率;Fs表示所述探鱼仪的采样率;C表示水中声速;
且所述滑动窗的宽度M为满足公式(1)的最大整奇数,M=2N+1,N为自然数。
优选地,所述探鱼仪的输出信号s(t)的解析表达式为:
其中,t为时间,τ为信号脉宽,rect表示矩形函数;
所述匹配滤波器的时域冲击响应函数h(t)的表达式为:
h(t)=w(t)s*(-t) (3);
其中,w(t表示加权函数;
所述滤波输出信号r(t)的表达式为:
优选地,所述S5步骤进一步包括步骤:
S51:设置一变量n的初值值为1;
S52:判断所述变量n是否大于等于所述离散滤波数据集的波值数据的总数L;如是,结束步骤;否则,继续后续步骤;
S53:获取所述离散滤波数据集的第n个所述波值数据至第n+M-1个所述波值数据中最大的一个所述波值数据作为目标波值数据,所述目标波值数据为所述离散滤波数据集中第p个所述波值数据,设定所述目标波值数据的位置为p;设定一目标波值e的数值为所述目标波值数据的数值;
S54:判断所述目标波值e是否大于所述检测阈值;如是,继续后续步骤;否则,将变量n赋值为n+2N+1,并返回步骤S52;
S55:判断所述目标波值数据的位置是否等于n+N;如是,继续后续步骤;否则,将变量n赋值为n+1,并返回步骤S52;
S56:将所述目标波值数据作为所述目标并输出所述目标的位置;
S57:将变量n赋值为p+1,并返回步骤S52。
优选地,所述S5步骤进一步包括步骤:
S51:设置一变量n的初值值为1;
S52:判断所述变量n是否大于等于所述离散滤波数据集的波值数据的总数L;如是,结束步骤;否则,继续后续步骤;
S53:获取所述离散滤波数据集的第n个所述波值数据至第n+M-1个所述波值数据中最大的一个所述波值数据作为目标波值数据,所述目标波值数据为所述离散滤波数据集中第p个所述波值数据,设定所述目标波值数据的位置为p;设定一目标波值e的数值为所述目标波值数据的数值;
S54:判断所述目标波值e是否大于所述检测阈值;如是,继续后续步骤;否则,将变量n赋值为n+2N+1,并返回步骤S52;
S55:判断所述目标波值数据的位置是否小于n+N;如是,将变量n赋值为n+N,并返回步骤S52;否则,继续后续步骤;
S56:判断所述目标波值数据的位置是否等于n+N;如是,将所述目标波值数据作为所述目标并输出所述目标的位置,将变量n赋值为p+1,并返回步骤S52;如否,继续后续步骤;
S57:将变量n赋值为n+1;
S58:判断N是否大于L;
如是,结束步骤;
否则,判断所述离散滤波数据集的第n+M-1个波值数据的数值是否大于所述目标波值e;如是,将所述目标波值数据的位置赋值为n+M-1,将所述目标波值e赋值为所述离散滤波数据集的第n+M-1个波值数据的数值,返回步骤S57;否则,返回步骤S56。
本发明由于采用了以上技术方案,使其具有以下有益效果:
本发明的一种用于宽带探鱼仪的高效目标分割方法,当两个目标回波重叠时,只要目标间隔大于探鱼仪的最小距离分辨率,就可以将两个目标区分开来;且运算量小,可快速获得目标的数量和位置。
附图说明
图1为本发明实施例的用于宽带探鱼仪的高效目标分割方法的流程图;
图2为本发明实施例一的分割目标波值数据步骤的流程图;
图3为本发明实施例一的探鱼仪的回波强度图;
图4为本发明实施例一的5个目标的探鱼仪的输出信号图;
图5为本发明实施例一的目标分割结果图;
图6为本发明实施例二的分割目标波值数据步骤的流程图。
具体实施方式
下面根据附图1~图6,给出本发明的较佳实施例,并予以详细描述,使能更好地理解本发明的功能、特点。
请参阅图1,本发明实施例一的一种用于宽带探鱼仪的高效目标分割方法,包括步骤:
S1:根据一背景噪声值确定一检测阈值,检测阈值大于背景噪声值。
其中,背景噪声值通过停止探鱼仪发射并记录探鱼仪的回波强度估算获得;当无法停止探鱼仪发射时利用超出海底范围或远距离收到的回波强度进行估算。
S2:设定一滑动窗的宽度,滑动窗的宽度小于等于一探鱼仪的距离分辨率。
其中,滑动窗的宽度M满足公式(1):
其中,ΔR表示探鱼仪的距离分辨率;Fs表示探鱼仪的采样率;C表示水中声速;
且滑动窗的宽度M为满足公式(1)的最大整奇数,M=2N+1,N为自然数。
S3:利用一匹配滤波器对探鱼仪的输出信号进行匹配滤波,获得一滤波输出信号的离散滤波数据集r(n),离散滤波数据集r(n)包括L个波值数据。
其中,探鱼仪的输出信号s(t)的解析表达式为:
其中,t为时间,τ为信号脉宽,rect表示矩形函数;
匹配滤波器的时域冲击响应函数h(t)的表达式为:
h(t)=w(t)s*(-t) (3);
其中,w(t表示加权函数;
滤波输出信号r(t)的表达式为:
S4:利用滑动窗和检测阈值对离散滤波数据集r(n)分割目标,并输出目标的位置。
从离散滤波数据集的第1个波值数据开始,查找滑动窗的宽度M内波值数据的最大值,当波值数据的最大值大于检测阈值且该波值数据的位置处在滑动窗中心时,认为找到1个目标波值数据,移动滑动窗,重复步骤,直至所有波值数据查找完成。
请参阅图2,S5步骤进一步包括步骤:
S51:设置一变量n的初值值为1;
S52:判断变量n是否大于等于离散滤波数据集的波值数据的总数L;如是,结束步骤;否则,继续后续步骤;
S53:获取离散滤波数据集的第n个波值数据至第n+M-1个波值数据中最大的一个波值数据作为目标波值数据,目标波值数据为离散滤波数据集中第p个波值数据,设定目标波值数据的位置为p;设定一目标波值e的数值为目标波值数据的数值;
S54:判断目标波值e是否大于检测阈值;如是,继续后续步骤;否则,将变量n赋值为n+2N+1,并返回步骤S52;
S55:判断目标波值数据的位置是否等于n+N;如是,继续后续步骤;否则,将变量n赋值为n+1,并返回步骤S52;
S56:将目标波值数据作为目标并输出目标的位置;
S57:将变量n赋值为p+1,并返回步骤S52。
例如:使用信号为带宽B=60kHz,脉宽τ=1ms的线性调频(LFM)信号,采样率为2MHz,使用的加权函数w(t)为汉明窗,匹配滤波后输出探鱼仪的回波强度图,请参阅图3。
可见探鱼仪的最小的距离分辨能力约为40mm,则M=53。
模拟在100m距离处有5个目标,距离分别为:100.0m 100.08m 100.12m100.15m100.17m,5个目标的回波状态请参阅图4。
使用本发明实施例一进行目标分割,获得目标分割结果图,请参阅图5。
前三个目标间隔距离大于距离分辨率,被识别出来,最后两个目标的间距小于最小分辨率被识别为1个目标。
请参阅图6,本发明实施例二的一种用于宽带探鱼仪的高效目标分割方法,其步骤与实施例一基本相同,其区别在于,S5步骤进一步包括步骤:
S51:设置一变量n的初值值为1;
S52:判断变量n是否大于等于离散滤波数据集的波值数据的总数L;如是,结束步骤;否则,继续后续步骤;
S53:获取离散滤波数据集的第n个波值数据至第n+M-1个波值数据中最大的一个波值数据作为目标波值数据,目标波值数据为离散滤波数据集中第p个波值数据,设定目标波值数据的位置为p;设定一目标波值e的数值为目标波值数据的数值;
S54:判断目标波值e是否大于检测阈值;如是,继续后续步骤;否则,将变量n赋值为n+2N+1,并返回步骤S52;
S55:判断目标波值数据的位置是否小于n+N;如是,将变量n赋值为n+N,并返回步骤S52;否则,继续后续步骤;
S56:判断目标波值数据的位置是否等于n+N;如是,将目标波值数据作为目标并输出目标的位置,将变量n赋值为p+1,并返回步骤S52;如否,继续后续步骤;
S57:将变量n赋值为n+1;
S58:判断N是否大于L;
如是,结束步骤;
否则,判断离散滤波数据集的第n+M-1个波值数据的数值是否大于目标波值e;如是,将目标波值数据的位置赋值为n+M-1,将目标波值e赋值为离散滤波数据集的第n+M-1个波值数据的数值,返回步骤S57;否则,返回步骤S56。
以上结合附图实施例对本发明进行了详细说明,本领域中普通技术人员可根据上述说明对本发明做出种种变化例。因而,实施例中的某些细节不应构成对本发明的限定,本发明将以所附权利要求书界定的范围作为本发明的保护范围。
Claims (3)
1.一种用于宽带探鱼仪的高效目标分割方法,包括步骤:
S1:根据一背景噪声值确定一检测阈值,所述检测阈值大于所述背景噪声值;
S2:设定一滑动窗的宽度,所述滑动窗的宽度小于等于一探鱼仪的距离分辨率;
S3:利用一匹配滤波器对所述探鱼仪的输出信号进行匹配滤波,获得一滤波输出信号的离散滤波数据集,所述离散滤波数据集包括多个波值数据;
S4:利用所述滑动窗和所述检测阈值对所述离散滤波数据集分割目标,并输出所述目标的位置;
所述S4步骤进一步包括步骤:
S51:设置一变量n的初值值为1;
S52:判断所述变量n是否大于等于所述离散滤波数据集的波值数据的总数L;如是,结束步骤;否则,继续后续步骤;
S53:获取所述离散滤波数据集的第n个所述波值数据至第n+M-1个所述波值数据中最大的一个所述波值数据作为目标波值数据,所述目标波值数据为所述离散滤波数据集中第p个所述波值数据,设定所述目标波值数据的位置为p;设定一目标波值e的数值为所述目标波值数据的数值;
S54:判断所述目标波值e是否大于所述检测阈值;如是,继续后续步骤;否则,将变量n赋值为n+2N+1,并返回步骤S52;
S55:判断所述目标波值数据的位置是否等于n+N;如是,继续后续步骤;否则,将变量n赋值为n+1,并返回步骤S52;
S56:将所述目标波值数据作为所述目标并输出所述目标的位置;
S57:将变量n赋值为p+1,并返回步骤S52;
所述滑动窗的宽度M满足公式(1):
其中,ΔR表示所述探鱼仪的距离分辨率;Fs表示所述探鱼仪的采样率;C表示水中声速;
且所述滑动窗的宽度M为满足公式(1)的最大整奇数,M=2N+1,N为自然数。
2.根据权利要求1所述的用于宽带探鱼仪的高效目标分割方法,其特征在于,所述背景噪声值通过停止所述探鱼仪发射并记录所述探鱼仪的回波强度估算获得,当无法停止所述探鱼仪发射时利用超出海底范围或远距离收到的回波强度进行估算。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010401601.9A CN111413702B (zh) | 2020-05-13 | 2020-05-13 | 用于宽带探鱼仪的高效目标分割方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010401601.9A CN111413702B (zh) | 2020-05-13 | 2020-05-13 | 用于宽带探鱼仪的高效目标分割方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111413702A CN111413702A (zh) | 2020-07-14 |
CN111413702B true CN111413702B (zh) | 2022-04-05 |
Family
ID=71493746
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010401601.9A Active CN111413702B (zh) | 2020-05-13 | 2020-05-13 | 用于宽带探鱼仪的高效目标分割方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111413702B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112084886B (zh) * | 2020-08-18 | 2022-03-15 | 眸芯科技(上海)有限公司 | 提升检测神经网络目标检测性能的方法及装置 |
CN112884803B (zh) * | 2020-08-18 | 2023-11-24 | 眸芯科技(上海)有限公司 | 基于dsp的实时智能监控目标检测方法及装置 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4089002A (en) * | 1977-04-05 | 1978-05-09 | The Bendix Corporation | CFAR threshold adaptive target coordinate extractor |
JPH0436680A (ja) * | 1990-05-31 | 1992-02-06 | Furuno Electric Co Ltd | シラス漁用魚群探知装置 |
US5311189A (en) * | 1991-09-19 | 1994-05-10 | Deutsche Aerospace Ag | Method for distinguishing between at least two targets |
CN103969639A (zh) * | 2014-05-09 | 2014-08-06 | 哈尔滨工程大学 | 五波束探鱼仪的信号处理系统及其信号处理方法 |
CN106970360A (zh) * | 2017-05-25 | 2017-07-21 | 哈尔滨工业大学 | 一种导航雷达多次反射假回波抑制方法 |
CN108009544A (zh) * | 2017-12-13 | 2018-05-08 | 北京小米移动软件有限公司 | 目标检测方法及装置 |
CN108614257A (zh) * | 2011-04-14 | 2018-10-02 | 古野电气株式会社 | 水中探测装置、鱼群探测机以及水中探测方法 |
CN109061625A (zh) * | 2018-06-22 | 2018-12-21 | 西安空间无线电技术研究所 | 一种距离速度匹配的运动目标检测方法 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7760587B2 (en) * | 2007-01-04 | 2010-07-20 | Ocean Acoustical Services and Instrumentation Systems (OASIS), Inc. | Methods of and systems for monitoring the acoustic transmission conditions in underwater areas using unmanned, mobile underwater vehicles |
US7916933B2 (en) * | 2007-01-18 | 2011-03-29 | Northrop Grumman Systems Corporation | Automatic target recognition system for detection and classification of objects in water |
CN101231342A (zh) * | 2007-07-24 | 2008-07-30 | 哈尔滨工程大学 | 拖曳式水平渔探仪 |
-
2020
- 2020-05-13 CN CN202010401601.9A patent/CN111413702B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4089002A (en) * | 1977-04-05 | 1978-05-09 | The Bendix Corporation | CFAR threshold adaptive target coordinate extractor |
JPH0436680A (ja) * | 1990-05-31 | 1992-02-06 | Furuno Electric Co Ltd | シラス漁用魚群探知装置 |
US5311189A (en) * | 1991-09-19 | 1994-05-10 | Deutsche Aerospace Ag | Method for distinguishing between at least two targets |
CN108614257A (zh) * | 2011-04-14 | 2018-10-02 | 古野电气株式会社 | 水中探测装置、鱼群探测机以及水中探测方法 |
CN103969639A (zh) * | 2014-05-09 | 2014-08-06 | 哈尔滨工程大学 | 五波束探鱼仪的信号处理系统及其信号处理方法 |
CN106970360A (zh) * | 2017-05-25 | 2017-07-21 | 哈尔滨工业大学 | 一种导航雷达多次反射假回波抑制方法 |
CN108009544A (zh) * | 2017-12-13 | 2018-05-08 | 北京小米移动软件有限公司 | 目标检测方法及装置 |
CN109061625A (zh) * | 2018-06-22 | 2018-12-21 | 西安空间无线电技术研究所 | 一种距离速度匹配的运动目标检测方法 |
Non-Patent Citations (2)
Title |
---|
Real-Time Wideband Acoustic Signal Processing for Fish Counting;Gabriel Vasile et al.;《OCEANS 2018 MTS》;20181025;全文 * |
多波束渔用声呐作用距离预报建模及性能分析;李国栋 等;《渔业现代化》;20200229;第47卷(第1期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN111413702A (zh) | 2020-07-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Roch et al. | Automated extraction of odontocete whistle contours | |
CN111413702B (zh) | 用于宽带探鱼仪的高效目标分割方法 | |
Filho et al. | Preprocessing passive sonar signals for neural classification | |
CN106249218B (zh) | 一种基于多频点回波信息融合的目标跟踪方法 | |
CN109977724B (zh) | 一种水下目标分类方法 | |
CN113050071B (zh) | 激光雷达数据处理方法、装置、设备及存储介质 | |
US11835649B2 (en) | Method and apparatus for radar signal processing using convolutional neural network | |
KR102011959B1 (ko) | 펄스 압축 과정에서 간섭신호를 탐지하는 레이더 수신신호 처리 방법 및 그를 위한 장치 | |
CN109507675B (zh) | 基于频分系统实现水下多目标时延估计处理的方法 | |
CN110632563B (zh) | 一种基于短时傅里叶变换的脉内频率编码信号参数测量方法 | |
CN110808791B (zh) | 基于分数阶傅里叶变换的线性调频信号的检测方法 | |
CN110515044A (zh) | 低慢小雷达空频域四维杂波抑制方法 | |
EP1609001A1 (fr) | Procede de traitement de signaux, et sonar actif le mettant en oeuvre | |
CN110299926B (zh) | 一种面向低信噪比环境的水声信号检测方法 | |
CN116794634A (zh) | 一种激光制导雷达的抗高重频干扰方法及装置 | |
CN110673118A (zh) | 一种主动声纳单频脉冲串波形设计及检测算法 | |
Gubnitsky et al. | Inter-pulse estimation for sperm whale click detection | |
RU2646857C1 (ru) | Способ селекции движущихся целей | |
CN108459316B (zh) | 一种多途信道下的回波信号参数估计方法 | |
CN112487887B (zh) | 一种提高主动目标检测信噪比的方法 | |
CN111123274A (zh) | 一种水下声呐成像系统目标检测方法 | |
CN112929053B (zh) | 一种跳频信号特征提取与参数估计方法 | |
CN114325671B (zh) | 改进的基于水下声信标信号多脉冲累积的抗多普勒方法 | |
CN113138372B (zh) | 一种基于改进型m/n检测器的多径环境下雷达目标检测方法 | |
CN115685214A (zh) | 一种主动声纳梳状谱波形检测算法 |
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 |