CN105068054A - 一种船用雷达海杂波的抑制算法 - Google Patents
一种船用雷达海杂波的抑制算法 Download PDFInfo
- Publication number
- CN105068054A CN105068054A CN201510523377.XA CN201510523377A CN105068054A CN 105068054 A CN105068054 A CN 105068054A CN 201510523377 A CN201510523377 A CN 201510523377A CN 105068054 A CN105068054 A CN 105068054A
- Authority
- CN
- China
- Prior art keywords
- curve
- sea
- parameter
- sea clutter
- data
- 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
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/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/28—Details of pulse systems
- G01S7/285—Receivers
- G01S7/292—Extracting wanted echo-signals
- G01S7/2923—Extracting wanted echo-signals based on data belonging to a number of consecutive radar periods
-
- 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/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/35—Details of non-pulse systems
- G01S7/352—Receivers
- G01S7/354—Extracting wanted echo-signals
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明涉及一种船用雷达海杂波的抑制算法,该算法包括有如下步骤:第一步,初始化并建立海情量程数据表;第二步,录入当前海况下的海杂波抑制曲线参数,读取海杂波抑制修正算法参数的设置;第三步,建立海杂波抑制曲线;第四步,依照海杂波抑制算法公式,计算出海杂波抑制曲线,通过海杂波抑制曲线叠加,完成海杂波抑制。本发明的船用雷达海杂波抑制的算法根据实际海情选择海杂波分布模式进行海杂波抑制,在尽可能大的抑制海杂波的同时,把对目标的影响降低到最小。
Description
技术领域
本发明涉及船用雷达探测,特别涉及到一种船用雷达海杂波的抑制算法。
背景技术
船用雷达是船舶的核心导航设备,其关系到国家海防和海事的安全。但是船用雷达的探测和跟踪性能常常会受海浪表面散射的电磁波的影响,海浪表面散射的电磁波也称为海杂波。海杂波的出现会增加虚警概率,或者相应地降低虚警率恒定时的检测概率。对于复杂海况来说,这种情况可能更加严重,海杂波是运动和起伏的,因而很难消除。
现有的船用雷达杂波抑制算法中,一般是采用四种海杂波模型(瑞利分布、对数正态分布、韦布尔分布、K分布)中的一种分布方式建立固定模型进行海杂波抑制,大都没有处理雷达附近的海杂波。这样会导致雷达显示屏的中心位置会有杂波圆环,并且在周围出现大量的干扰信号,影响对真实目标的判断,如图2所示,图中的当前量程3海里的雷达回波图像,可见近程杂波强度很大,真实目标淹没在杂波中。而且在实际应用过程中,大多情况是现有技术四种模型的混合状态,单一固定模型的海杂波抑制也会对近距离目标信号造成很大的衰减,甚至导致漏检的发生。
发明内容
本发明的目的在于克服上述现有技术中存在的不足,提供一种新的船用雷达海杂波的抑制算法。本发明的海杂波抑制算法通过实际海情来选择海杂波分布模式进行海杂波抑制,在尽可能大的抑制海杂波的同时,把对目标的影响降低到最小,呈现出更为清晰真实的雷达探测目标。
为了达到上述发明目的,本发明提供的技术方案如下:
一种船用雷达海杂波的抑制算法,其特征在于,该算法包括有如下步骤:
第一步,基于现有的雷达显控软件进行初始化,建立海情量程数据表;所述的初始化为参数的初始化设置和数据的初始化,建立海情量程数据表分为两个步骤:
(11)读取当前的量程数据,该量程数据包括有start、end、sz和j;
(12)由公式建立当前量程数据表,其中table是建立的量程数据表数组,start为当前量程起始点,end为当前量程终止点,sz为两个距离环之间的差值,j是表格数组步进参数;
第二步,通过当前海况读取海杂波抑制算法的参数设置,PWR:杂波抑制曲线总体放大倍数,P0:手动曲线放大倍数,P1:一次自动曲线放大倍数,P2:二次自动曲线放大倍数,P3:二次自动曲线值域调整参数,:当前海浪抑制参数数据;
第三步,建立海杂波抑制曲线,海杂波抑制曲线建立包括有如下步骤:
(31)建立手动曲线数据表,由公式计算得出,其中sea0是手动曲线数据数组,PWR和P0为第二步中用户依据当前海况和雷达回波图像设置的参数;
(32)建立一次自动曲线,由公式计算得出,其中sea1是一次自动曲线数据数组,PWR、P1和为第二步中用户依据当前海况和雷达回波图像设置的参数,为步骤(31)计算得出的手动曲线数据数组;
(33)建立二次自动曲线,先由公式计算得出二次自动曲线值域调整参数,其中A是二次自动曲线值域参数数据,P3和为第二步中用户设置参数;之后由公式计算得出二次自动曲线数据数组,式中是二次自动曲线数据数组;
第四步,依照上述第三步的海杂波抑制算法公式来计算出海杂波抑制曲线,所述的海杂波抑制曲线由如下公式计算得出:
,
式中为最终杂波抑制曲线数据数组,通过上述海杂波抑制曲线叠加来完成海杂波抑制。
在本发明一种船用雷达海杂波的抑制算法的运行过程中,作为优化的技术方案,还需要断船用雷达的探测距离是否有变化,若船用雷达的探测距离有变化则执行步骤(11),若探测距离无变化则进一步判断海杂波抑制算法参数设置是否有变化,若海杂波抑制算法参数设置有变化则执行步骤(31),若海杂波抑制算法参数设置无变化则继续循环判断船用雷达的探测距离是否有变化。
基于上述技术方案,本发明专利与现有技术相比具有如下技术优点:
本发明船用雷达海杂波抑制算法根据实际海情来选择海杂波分布模式,利用叠加的抑制曲线来进行海杂波抑制,在尽可能大的抑制海杂波的同时,把对目标的影响降低到最小。
附图说明
图1是本发明船用雷达海杂波抑制算法的流程示意图。
图2是采用本发明船用雷达海杂波抑制算法之前实测数据的雷达显示图。
图3是采用本发明船用雷达海杂波抑制算法之后实测数据的雷达显示图。
具体实施方式
下面我们结合附图和具体的实施例来的本发明专利的操作流程和工作方式做进一步的详细阐述,以求更为清楚明了地理解其操作原理和工作流程但不能以此来限制本发明专利的保护范围。
本发明专利通过实际海情调整不同海杂波抑制参数进行海杂波抑制,在尽可能大地抑制海杂波的同时,把对目标的影响降到最小。该船用雷达海杂波抑制算法主要包括有如下步骤:一、初始化,建立海情量程数据表;二、通过当前海况读取海杂波抑制算法参数设置;三、依照海杂波抑制算法公式,计算出杂波抑制曲线;四、杂波抑制曲线叠加,完成海杂波抑制算法。
如图1所示,在具体的操作层面上,本发明的海杂波修正算法包括有如下步骤:
第一步,基于现有的雷达显控软件进行参数和数据初始化,建立海情量程数据表;所述的初始化为参数的初始化设置和数据的初始化,建立海情量程数据表分为两个步骤:
(11)读取当前的量程数据,量程是由用户在现有的雷达显控软件中自行切换,算法根据用户设置的当前量程,在上位机内存中读取;
(12)由公式建立当前量程数据表,其中table是建立的量程数据表数组,start为当前量程起始点,end为当前量程终止点,sz为两个距离环之间的差值,j是表格数组步进参数;
第二步,通过当前海况,来读取海杂波抑制算法参数的设置,PWR:杂波抑制曲线总体放大倍数,P0:手动曲线放大倍数,P1:一次自动曲线放大倍数,P2:二次自动曲线放大倍数,P3:二次自动曲线值域调整参数,:当前海浪抑制参数数据。上述的参数是上位机的雷达显控软件中用户自行设置的参数,在上位机的内存中直接进行读取。也就是说用户根据当前海面情况,依照当前回波图像,自行设置海杂波抑制曲线的4个参数,即PO、P1、P2和P3。
第三步,在船用雷达显控软件中建立海杂波抑制曲线,海杂波抑制曲线建立包括有如下步骤:
(31)建立手动曲线数据表,由公式计算得出,其中sea0是手动曲线数据数组,PWR和P0为第二步中用户依据当前海况和雷达回波图像设置的参数。
(32)建立一次自动曲线,由公式计算得出,其中sea1是一次自动曲线数据数组,PWR、P1和为第二步中用户依据当前海况和雷达回波图像设置的参数,为步骤(31)计算得出的手动曲线数据数组,这里一次曲线的斜率固定,抑制为线性。
(33)建立二次自动曲线,先由公式计算得出二次自动曲线值域调整参数,其中A是二次自动曲线值域参数数据,P3和为第二步中用户设置参数;之后由公式计算得出二次自动曲线数据数组,式中是二次自动曲线数据数组;这里的二次曲线为双曲线,双曲线抑制离雷达中心抑制幅度稿,远离雷达抑制幅度低,针对雷达近程海杂波影像比较大的因素。
第四步,依照上述第三步的三个海杂波抑制算法公式,计算出海杂波抑制曲线,由如下公式计算得出抑制算法曲线:
,
式中为最终杂波抑制曲线数据数组,通过海杂波抑制曲线叠加,完成海杂波抑制。
作为实施例,针对于图2来说,图3是采用本发明船用雷达海杂波抑制算法之后实测数据的雷达显示图。图中设置好杂波抑制曲线4个参数:P0=225,P1=40,P2=91,P3=23,当前量程为39。杂波抑制曲线叠加后的效果,可见近程的雷达杂波环消失,且真实目标清晰可见。从而充分体现了本发明专利海杂波抑制算法的优越性。
在本发明一种船用雷达海杂波的抑制算法的运行过程中,作为优化的技术方案,还需要断船用雷达的探测距离是否有变化,若船用雷达的探测距离有变化则执行步骤(11),若探测距离无变化则进一步判断海杂波抑制算法参数设置是否有变化,若海杂波抑制算法参数设置有变化则执行步骤(31),若海杂波抑制算法参数设置无变化则继续循环判断船用雷达的探测距离是否有变化。
综上所述,海杂波抑制算法采用一次曲线与双曲线数据叠加的方式,并根据当前显示量程,计算得出抑制数据表格,进行海杂波抑制,根据不同海情,从而进行了有效海杂波抑制。
Claims (2)
1.一种船用雷达海杂波的抑制算法,其特征在于,该算法包括有如下步骤:
第一步,基于现有的雷达显控软件进行初始化,建立海情量程数据表;所述的初始化为参数的初始化设置和数据的初始化,建立海情量程数据表分为两个步骤:
(11)读取当前的量程数据,该量程数据包括有start、end、sz和j;
(12)由公式建立当前量程数据表,其中table是建立的量程数据表数组,start为当前量程起始点,end为当前量程终止点,sz为两个距离环之间的差值,j是表格数组步进参数;
第二步,通过当前海况读取海杂波抑制算法的参数设置,PWR:杂波抑制曲线总体放大倍数,P0:手动曲线放大倍数,P1:一次自动曲线放大倍数,P2:二次自动曲线放大倍数,P3:二次自动曲线值域调整参数,:当前海浪抑制参数数据;
第三步,建立海杂波抑制曲线,海杂波抑制曲线建立包括有如下步骤:
(31)建立手动曲线数据表,由公式计算得出,其中sea0是手动曲线数据数组,PWR和P0为第二步中用户依据当前海况和雷达回波图像设置的参数;
(32)建立一次自动曲线,由公式计算得出,其中sea1是一次自动曲线数据数组,PWR、P1和为第二步中用户依据当前海况和雷达回波图像设置的参数,为步骤(31)计算得出的手动曲线数据数组;
(33)建立二次自动曲线,先由公式计算得出二次自动曲线值域调整参数,其中A是二次自动曲线值域参数数据,P3和为第二步中用户设置参数;之后由公式计算得出二次自动曲线数据数组,式中是二次自动曲线数据数组;
第四步,依照上述第三步的海杂波抑制算法公式来计算出海杂波抑制曲线,所述的海杂波抑制曲线由如下公式计算得出:
,
式中为最终杂波抑制曲线数据数组,通过上述海杂波抑制曲线叠加来完成海杂波抑制。
2.根据权利要求1所述的一种船用雷达海杂波的抑制算法,其特征在于,在运行步骤中需判断船用雷达的探测距离是否有变化,若船用雷达的探测距离有变化则执行步骤(11),若探测距离无变化则进一步判断海杂波抑制算法参数设置是否有变化,若海杂波抑制算法参数设置有变化则执行步骤(31),若海杂波抑制算法参数设置无变化则继续循环判断船用雷达的探测距离是否有变化。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510523377.XA CN105068054B (zh) | 2015-08-25 | 2015-08-25 | 一种船用雷达海杂波的抑制算法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510523377.XA CN105068054B (zh) | 2015-08-25 | 2015-08-25 | 一种船用雷达海杂波的抑制算法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105068054A true CN105068054A (zh) | 2015-11-18 |
CN105068054B CN105068054B (zh) | 2017-06-20 |
Family
ID=54497469
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510523377.XA Active CN105068054B (zh) | 2015-08-25 | 2015-08-25 | 一种船用雷达海杂波的抑制算法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105068054B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105717494A (zh) * | 2016-01-27 | 2016-06-29 | 重庆邮电大学 | 一种基于小波变换的船用雷达海杂波抑制曲线设计方法 |
CN107797097A (zh) * | 2017-10-30 | 2018-03-13 | 中船航海科技有限责任公司 | 一种船用雷达中心控制方法 |
CN111323755A (zh) * | 2019-12-23 | 2020-06-23 | 北京海兰信数据科技股份有限公司 | 岸基非相参雷达杂波抑制的方法及装置 |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US3893117A (en) * | 1972-05-10 | 1975-07-01 | Tokyo Keiki Kk | Marine radar transmission and reception system |
CN101008673A (zh) * | 2007-01-24 | 2007-08-01 | 南京大桥机器有限公司 | 二次测风雷达接收控制、数据处理方法 |
CN101806887A (zh) * | 2010-03-19 | 2010-08-18 | 清华大学 | 基于空间跟踪滤波器的海杂波抑制和目标检测方法 |
CN102176000A (zh) * | 2011-01-27 | 2011-09-07 | 电子科技大学 | 一种船用雷达海杂波的抑制方法 |
US20110291877A1 (en) * | 2010-06-01 | 2011-12-01 | Raytheon Company | Methods and apparatus for non-isotropic sea clutter modeling |
CN102540161A (zh) * | 2011-12-01 | 2012-07-04 | 宁波成电泰克电子信息技术发展有限公司 | 一种船用雷达海杂波参数获取方法 |
CN103645467A (zh) * | 2013-11-11 | 2014-03-19 | 北京环境特性研究所 | 海杂波抑制以及海杂波背景中目标检测的方法和系统 |
CN104133195A (zh) * | 2014-07-23 | 2014-11-05 | 成都天奥信息科技有限公司 | 一种海浪干扰抑制方法 |
CN104849707A (zh) * | 2015-04-27 | 2015-08-19 | 中国电子科技集团公司第二十八研究所 | 一种基于多雷达分布式检测的海杂波抑制方法 |
-
2015
- 2015-08-25 CN CN201510523377.XA patent/CN105068054B/zh active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US3893117A (en) * | 1972-05-10 | 1975-07-01 | Tokyo Keiki Kk | Marine radar transmission and reception system |
CN101008673A (zh) * | 2007-01-24 | 2007-08-01 | 南京大桥机器有限公司 | 二次测风雷达接收控制、数据处理方法 |
CN101806887A (zh) * | 2010-03-19 | 2010-08-18 | 清华大学 | 基于空间跟踪滤波器的海杂波抑制和目标检测方法 |
US20110291877A1 (en) * | 2010-06-01 | 2011-12-01 | Raytheon Company | Methods and apparatus for non-isotropic sea clutter modeling |
CN102176000A (zh) * | 2011-01-27 | 2011-09-07 | 电子科技大学 | 一种船用雷达海杂波的抑制方法 |
CN102540161A (zh) * | 2011-12-01 | 2012-07-04 | 宁波成电泰克电子信息技术发展有限公司 | 一种船用雷达海杂波参数获取方法 |
CN103645467A (zh) * | 2013-11-11 | 2014-03-19 | 北京环境特性研究所 | 海杂波抑制以及海杂波背景中目标检测的方法和系统 |
CN104133195A (zh) * | 2014-07-23 | 2014-11-05 | 成都天奥信息科技有限公司 | 一种海浪干扰抑制方法 |
CN104849707A (zh) * | 2015-04-27 | 2015-08-19 | 中国电子科技集团公司第二十八研究所 | 一种基于多雷达分布式检测的海杂波抑制方法 |
Non-Patent Citations (3)
Title |
---|
J LIU等: ""Radar sea clutter suppression and target indication with a spatial tracking filter"", 《TSINGHUA SCIENCE AND TECHNOLOGY》 * |
姚云萍等: ""一种船用雷达的海杂波抑制算法"", 《火控雷达技术》 * |
姚元飞等: ""降低船用连续波雷达虚警概率的海杂波抑制方法"", 《电讯技术》 * |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105717494A (zh) * | 2016-01-27 | 2016-06-29 | 重庆邮电大学 | 一种基于小波变换的船用雷达海杂波抑制曲线设计方法 |
CN105717494B (zh) * | 2016-01-27 | 2018-06-05 | 重庆邮电大学 | 一种基于小波变换的船用雷达海杂波抑制曲线设计方法 |
CN107797097A (zh) * | 2017-10-30 | 2018-03-13 | 中船航海科技有限责任公司 | 一种船用雷达中心控制方法 |
CN107797097B (zh) * | 2017-10-30 | 2020-04-17 | 中船航海科技有限责任公司 | 一种船用雷达中心控制方法 |
CN111323755A (zh) * | 2019-12-23 | 2020-06-23 | 北京海兰信数据科技股份有限公司 | 岸基非相参雷达杂波抑制的方法及装置 |
Also Published As
Publication number | Publication date |
---|---|
CN105068054B (zh) | 2017-06-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP7182869B2 (ja) | 物標検出装置 | |
CN102540176B (zh) | 物标探测方法、物标探测装置及雷达装置 | |
Liu et al. | Wave Height Estimation from Shipborne X‐Band Nautical Radar Images | |
CN104297758B (zh) | 一种基于2d脉冲式激光雷达的辅助靠泊装置及其方法 | |
Park et al. | A study on propeller noise source localization in a cavitation tunnel | |
CN105068054A (zh) | 一种船用雷达海杂波的抑制算法 | |
CN110910448A (zh) | 一种海上溢油漂移轨迹溯源方法 | |
KR20170054175A (ko) | 미식별 수중신호 위치의 추정 방법 및 시스템 | |
Thornhill et al. | Real time local sea state measurement using wave radar and ship motions | |
CN117452391B (zh) | 海上风电桩基的冲刷监测方法、装置、设备、系统及介质 | |
CN112882037B (zh) | 侧扫声呐海底线检测方法及装置 | |
CN107561517A (zh) | 一种雷达海杂波剩余虚假点迹抑制方法 | |
CN105487055A (zh) | 一种航海雷达目标检测及显示系统 | |
CN104568846B (zh) | 基于布里渊散射海水盐跃层的二维扫描探测法 | |
KR101400924B1 (ko) | 원격탐사자료를 이용한 기름유출 감시 신뢰도를 높이는 기름유출 감시방법 | |
Xu et al. | Hydrographic data inspection and disaster monitoring using shipborne radar small range images with electronic navigation chart | |
CN114909610B (zh) | 一种水下油气管道泄漏检测与定位的方法及控制系统 | |
KR101961380B1 (ko) | 비음향 수중 잠수함 및 돌출물 탐지 시스템, 방법, 및 상기 방법을 실행시키기 위한 컴퓨터 판독 가능한 프로그램을 기록한 기록 매체 | |
Anderson | Prospects for submarine detection with hf over-the-horizon radar | |
CN113989691A (zh) | 船只的异常行为检测方法以及相关装置 | |
AU2020309309A1 (en) | Method for determining the spatial coverage of a detection system | |
Fietz et al. | Background Noise Cancellation from Active Sonar Images Utilizing Transmission-free Recordings | |
Ehlers et al. | Advances in signal processing for maritime applications | |
JP2016121970A (ja) | エコー信号処理装置、レーダ装置、エコー信号処理方法、及びプログラム | |
CN116930976B (zh) | 基于小波模极大值的侧扫声呐图像的海底线检测方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |