CN113376612B - 一种基于地形矩阵化及探测的雷达杂波生成方法 - Google Patents
一种基于地形矩阵化及探测的雷达杂波生成方法 Download PDFInfo
- Publication number
- CN113376612B CN113376612B CN202110924385.0A CN202110924385A CN113376612B CN 113376612 B CN113376612 B CN 113376612B CN 202110924385 A CN202110924385 A CN 202110924385A CN 113376612 B CN113376612 B CN 113376612B
- Authority
- CN
- China
- Prior art keywords
- radar
- ray
- height
- calculation
- terrain
- 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
- 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/41—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
- G01S7/414—Discriminating targets with respect to background clutter
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开了一种基于地形矩阵化及探测的雷达杂波生成方法,所述方法包括如下步骤:S1:在待检测区域内划定一个预设区域范围,并对预设区域范围进行坐标系处理,获取处理后的区域内所有高程数据;S2:计算出雷达射线空间直线方程,并随后通过雷达射线空间直线方程遍历区域矩阵的每个计算节点;S3:遍历每个计算节点后,获取每个节点的地形高度是否高于当前投影的了雷达射线高度;若高于雷达射线高度,则判断为雷达射线被遮挡,将遮挡点作为雷达回波点进行反射计算;若低于雷达射线高度,则判断为雷达射线没有被遮挡,用原计算节点作为雷达回波点进行反射计算;S4:设立雷达接收灵敏度,将计算结果与雷达接收灵敏度进行阈值判断。
Description
技术领域
本发明涉及一种探测雷达杂波生成的方法,具体涉及一种基于地形矩阵化及探测的雷达杂波生成方法。
背景技术
雷达通过发射和接收无线电波,可以实现对目标物体的探测,例如, 探测目标物体的距离、方位、高度、速度、形状和目标类型等。在实施探 测时,雷达通过发射雷达波束并接收雷达波束被目标物体的回波信号进行 探测。回波信号中,除了目标物体反射雷达波束的回波信号(无线电波) 之外,由于雷达波束在传播过程中,一部分雷达波束还会触及地物,例如, 地表、山丘及建筑物。因而,可能被地表、山丘及建筑物等阻挡并形成反 射,使得雷达射线空间直线方程部分雷达波束不能有效传播到目标物体,从而不可避免地接收 到地物的反射回波,即地杂波,对从目标物体反射回的回波信号形成干扰, 即地形起伏成了对雷达波束的遮挡,造成雷达的探测误差,降低了雷达 对目标物体的探测精度。
雷达系统的信号处理,通常都是先进行距离维的脉冲压缩,将回波按距离单元进行划分,然后进行杂波抑制。因此,在对杂波抑制算法进行研究与验证时,需要能够产生出相应的距离环的杂波回波,实际中,大部分地区都不是海拔为零,而是具有一定的高度,杂波对应的地表在具有一定高度后将影响其到达雷达的斜距。因此,在进行等距离环杂波仿真时,必须考虑实际的地表高程。
目前所采用的技术中是有将区域进行栅格化,再对每个栅格区域的电磁场进行计算,但这样处理对于大型区域栅格点会过多,计算量会非常大。且对于其中的不同小区域,用户的关注度可能不一样,对计算精度的要求也有所差异。因而现有的栅格化做法难以满足实际应用的需求,极大的限制了雷达射线空间直线方程技术的应用范围。
发明内容
本发明所要解决的技术问题是在一定区域中,地形差异过大后,在进行雷达杂波探测时,往往无法获取山脉或者沟壑遮蔽后的杂波,并且没有能够调整雷达接收灵敏度,以满足不同条件上回波反射探测,本申请文件目的在于提供一种基于地形矩阵化及探测的雷达杂波生成方法,解决上述的问题。
本发明通过下述技术方案实现:
一种基于地形矩阵化及探测的雷达杂波生成方法,所述方法包括如下步骤:S1:在待检测区域内划定一个预设区域范围,并对预设区域范围进行坐标系处理,获取处理后的区域内所有高程数据;S2:计算出雷达射线空间直线方程,并随后通过雷达射线空间直线方程遍历区域矩阵的每个计算节点;S3:遍历每个计算节点后,获取每个节点的地形高度是否高于当前投影的了雷达射线高度;若高于雷达射线高度,则判断为雷达射线被遮挡,将遮挡点作为雷达回波点进行反射计算;若低于雷达射线高度,则判断为雷达射线没有被遮挡,用原计算节点作为雷达回波点进行反射计算;S4:设立雷达接收灵敏度,将计算结果与雷达接收灵敏度进行阈值判断,若计算结果超过雷达接收灵敏度,则记录当前计算节点的位置信息;若计算结果低于雷达接收灵敏度,则不进行记录,略过当前节点位置信息;S5:在遍历整个区域内的计算节点后,将结果输出。
进一步地,所述步骤S1中进行空间坐标系处理方式是采用栅格化处理方式将所有坐标点进行预处理,获取区域内所有坐标点的坐标数据。
进一步地,所述步骤S2中计算雷达射线空间直线方程时,将雷达射线分成若干段,投影每段段点至地形高程上,在通过计算每段段点的经纬高度来计算出雷达射线当前高度。
进一步地,所述步骤S3中的反射计算公式为:
本发明与现有技术相比,具有如下的优点和有益效果:
1、本发明一种基于地形矩阵化及探测的雷达杂波生成方法,通过将区域内进行空间坐标系的栅格化处理,并获取每个点的高程数据,通过流程计算使其能在任意区域中生成满足雷达接收灵敏度以及被山脉或者沟壑遮蔽后的杂波。
附图说明
此处所说明的附图用来提供对本发明实施例的进一步理解,构成本申请的一部分,并不构成对本发明实施例的限定。在附图中:
图1为本发明流程图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚明白,下面结合实施例和附图,对本发明作进一步的详细说明,本发明的示意性实施方式及其说明仅用于解释本发明,并不作为对本发明的限定。
实施例
如图1所示,本发明一种基于地形矩阵化及探测的雷达杂波生成方法,首先对区域进行一定精度的矩阵栅格化后得到矩阵内每个精度点坐标,并且取得区域所有高程数据;再计算处雷达射线的空间直线方程,将雷达射线分成若干段,投影每段段点至地形高程上,在通过计算每段段点的经纬高度来计算出雷达射线当前高度;两个高度相互比较即可得出雷达射线是否被山脉或者沟壑所遮挡的结论;如果雷达射线被遮挡了,那么雷达回波点就是此遮挡点;此后通过自由空间下的雷达方程计算出雷达回波能量,将其与雷达接收灵敏度作比较,大于灵敏度则记录当前杂波点的坐标,反之不做处理。其中使用的雷达方程:
仿真中建立的杂波模型不应过分复杂,应以数据精度要求为限。为了杂波模型能够以可接收的精度来仿真实现,因此需要一个相对复杂的杂波模型,雷达射线空间直线方程雷达方程公式为:
其中Pt为发射天线功率,Gr为发射天线增益,Rmax为最大传播距离,Ae为接收天线的
有效面积,为目标雷达散射截面积,为最小可检测信号;部分辐射的能量被一定距
离的目标截获,目标截获的能量重新辐射到许多方向。其中目标RCS为;一部分重新辐
射的能量返回雷达天线,被雷达天线接收,其中接收天线的有效面积为Ae;在被接收机放大
和信号处理后,在接收机的输出端做出目标回波信号是否存在的判决(最小可检测信噪比
由判决条件确定)。其中最小可检测信号为,大多数雷达发射与接收为同一天线,有天
线理论得,其中λ为雷达信号的波长。其中,天线增益的计算方式通常取决
于的天线的类型,如果天线为全向天线即取天线最大增益为值;如果天线为单向天线,则天
线的计算方式取决于天线方向图函数。
上述的天线函数中,现举例常见的几种如下所示:
高斯函数:
单向余弦函数 :
双向余弦函数:
单向辛克函数:
双向辛克函数 :
相关参数解释:θ为目标点相对于天线的方位角(单位度),θ1天线部署的方位角
(单位度),为天线主波束的水平方向的3db波束宽度(单位度),为目标点相对于天
线的俯仰角(单位度),天线部署的俯仰角(单位度),为天线主波束的垂直方向
3db波束宽度(单位度),az=θ-θ1方位角之差(单位度),el = - 为俯仰角之差(单位
度)。
以上所述的具体实施方式,对本发明的目的、技术方案和有益效果进行了进一步详细说明,所应理解的是,以上所述仅为本发明的具体实施方式而已,并不用于限定本发明的保护范围,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (4)
1.一种基于地形矩阵化及探测的雷达杂波生成方法,其特征在于,所述方法包括如下步骤:
S1:在待检测区域内划定一个预设区域范围,并对预设区域范围进行栅格化坐标系处理,获取处理后的区域内所有高程数据;
S2:计算出雷达射线空间直线方程,并随后通过雷达射线空间直线方程遍历区域矩阵的每个计算节点;
S3:遍历每个计算节点后,获取每个节点的地形高度是否高于当前投影的雷达射线高度;若高于雷达射线高度,则判断为雷达射线被遮挡,将遮挡点作为雷达回波点进行反射计算;若低于雷达射线高度,则判断为雷达射线没有被遮挡,用原计算节点作为雷达回波点进行反射计算;
S4:设立雷达接收灵敏度,将计算结果与雷达接收灵敏度进行阈值判断,若计算结果超过雷达接收灵敏度,则记录当前计算节点的位置信息;若计算结果低于雷达接收灵敏度,则不进行记录,略过当前节点位置信息;
S5:在遍历整个区域内的计算节点后,将结果输出。
2.根据权利要求1所述的一种基于地形矩阵化及探测的雷达杂波生成方法,其特征在于,所述步骤S1中进行空间坐标系处理方式是采用栅格化处理方式将所有坐标点进行预处理,获取区域内所有坐标点的坐标数据。
3.根据权利要求1所述的一种基于地形矩阵化及探测的雷达杂波生成方法,其特征在于,所述步骤S2中计算雷达射线空间直线方程时,将雷达射线分成若干段,投影每段段点至地形高程上,在通过计算每段段点的经纬高度来计算出雷达射线当前高度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110924385.0A CN113376612B (zh) | 2021-08-12 | 2021-08-12 | 一种基于地形矩阵化及探测的雷达杂波生成方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110924385.0A CN113376612B (zh) | 2021-08-12 | 2021-08-12 | 一种基于地形矩阵化及探测的雷达杂波生成方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113376612A CN113376612A (zh) | 2021-09-10 |
CN113376612B true CN113376612B (zh) | 2021-11-23 |
Family
ID=77576976
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110924385.0A Active CN113376612B (zh) | 2021-08-12 | 2021-08-12 | 一种基于地形矩阵化及探测的雷达杂波生成方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113376612B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113866741A (zh) * | 2021-09-30 | 2021-12-31 | 成都众享天地网络科技有限公司 | 一种模拟地杂波点迹的算法 |
CN115032605A (zh) * | 2022-08-10 | 2022-09-09 | 中国航天科工集团八五一一研究所 | 一种基于目标能力画像的雷达目标体系分析方法 |
CN116224260A (zh) * | 2023-05-06 | 2023-06-06 | 成都众享天地网络科技有限公司 | 一种基于干扰的雷达三维威力范围计算方法 |
Citations (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5563601A (en) * | 1985-08-16 | 1996-10-08 | Northrop Grumman Corporation. | Two-port synthetic aperature radar system for radar detection of targets |
CN104345305A (zh) * | 2013-07-30 | 2015-02-11 | 北京敏视达雷达有限公司 | 一种地物杂波残留滤除方法 |
CN104408292A (zh) * | 2014-11-02 | 2015-03-11 | 北京环境特性研究所 | 一种动态目标的快速射线追踪方法 |
CN104656078A (zh) * | 2013-11-25 | 2015-05-27 | 北京环境特性研究所 | 一种基于射线追踪的海面舰船的方法 |
WO2018122282A1 (de) * | 2016-12-29 | 2018-07-05 | Avl List Gmbh | Verfahren und system zur simulationsgestützten bestimmung von echopunkten sowie verfahren zur emulation und emulationsvorrichtung |
CN108267722A (zh) * | 2018-01-23 | 2018-07-10 | 航天建筑设计研究院有限公司 | 地质雷达回波信号物性解构与探测目标数字重构智能化识取方法 |
JP2018146548A (ja) * | 2017-03-09 | 2018-09-20 | 三菱重工業株式会社 | 反射特性評価システムおよび反射特性評価方法 |
CN109683160A (zh) * | 2019-03-19 | 2019-04-26 | 长沙深之瞳信息科技有限公司 | 一种适用于穿墙雷达的目标检测方法 |
CN110082738A (zh) * | 2019-06-18 | 2019-08-02 | 西安电子科技大学 | 基于高斯混合和张量循环神经网络的雷达目标识别方法 |
CN110145969A (zh) * | 2018-08-06 | 2019-08-20 | 蓝箭航天空间科技股份有限公司 | 导弹拦截方法及服务器 |
CN112098958A (zh) * | 2020-10-22 | 2020-12-18 | 中国电波传播研究所(中国电子科技集团公司第二十二研究所) | 一种基于数字地图和气象水文信息的雷达杂波预测方法 |
CN113189562A (zh) * | 2021-07-02 | 2021-07-30 | 成都众享天地网络科技有限公司 | 一种基于高程的地形探测算法 |
CN113191027A (zh) * | 2021-06-28 | 2021-07-30 | 成都众享天地网络科技有限公司 | 基于地形二维、三维矩阵化辐射源信号强度分布计算方法 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
ITTO20070620A1 (it) * | 2007-08-31 | 2009-03-01 | Giancarlo Capaccio | Sistema e metodo per presentare dati visuali rilevati a distanza in immagini multi-spettrali, di fusione, ed a tre dimensioni spaziali. |
DE102008001467A1 (de) * | 2008-04-30 | 2009-11-05 | Robert Bosch Gmbh | Mehrstrahlradarsensor |
CN103728604B (zh) * | 2013-11-19 | 2016-01-13 | 中国国土资源航空物探遥感中心 | 一种宽带合成孔径雷达子带干涉数据处理方法 |
CN106842165B (zh) * | 2017-03-16 | 2020-02-18 | 电子科技大学 | 一种基于不同距离角度分辨率雷达集中式异步融合方法 |
-
2021
- 2021-08-12 CN CN202110924385.0A patent/CN113376612B/zh active Active
Patent Citations (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5563601A (en) * | 1985-08-16 | 1996-10-08 | Northrop Grumman Corporation. | Two-port synthetic aperature radar system for radar detection of targets |
CN104345305A (zh) * | 2013-07-30 | 2015-02-11 | 北京敏视达雷达有限公司 | 一种地物杂波残留滤除方法 |
CN104656078A (zh) * | 2013-11-25 | 2015-05-27 | 北京环境特性研究所 | 一种基于射线追踪的海面舰船的方法 |
CN104408292A (zh) * | 2014-11-02 | 2015-03-11 | 北京环境特性研究所 | 一种动态目标的快速射线追踪方法 |
WO2018122282A1 (de) * | 2016-12-29 | 2018-07-05 | Avl List Gmbh | Verfahren und system zur simulationsgestützten bestimmung von echopunkten sowie verfahren zur emulation und emulationsvorrichtung |
JP2018146548A (ja) * | 2017-03-09 | 2018-09-20 | 三菱重工業株式会社 | 反射特性評価システムおよび反射特性評価方法 |
CN108267722A (zh) * | 2018-01-23 | 2018-07-10 | 航天建筑设计研究院有限公司 | 地质雷达回波信号物性解构与探测目标数字重构智能化识取方法 |
CN110145969A (zh) * | 2018-08-06 | 2019-08-20 | 蓝箭航天空间科技股份有限公司 | 导弹拦截方法及服务器 |
CN109683160A (zh) * | 2019-03-19 | 2019-04-26 | 长沙深之瞳信息科技有限公司 | 一种适用于穿墙雷达的目标检测方法 |
CN110082738A (zh) * | 2019-06-18 | 2019-08-02 | 西安电子科技大学 | 基于高斯混合和张量循环神经网络的雷达目标识别方法 |
CN112098958A (zh) * | 2020-10-22 | 2020-12-18 | 中国电波传播研究所(中国电子科技集团公司第二十二研究所) | 一种基于数字地图和气象水文信息的雷达杂波预测方法 |
CN113191027A (zh) * | 2021-06-28 | 2021-07-30 | 成都众享天地网络科技有限公司 | 基于地形二维、三维矩阵化辐射源信号强度分布计算方法 |
CN113189562A (zh) * | 2021-07-02 | 2021-07-30 | 成都众享天地网络科技有限公司 | 一种基于高程的地形探测算法 |
Non-Patent Citations (1)
Title |
---|
"天发船收高频雷达坐标配准与多径数据处理研究";严万泉;《中国优秀博硕士学位论文全文数据库(硕士)信息科技辑》;20200215;第I136-2005页 * |
Also Published As
Publication number | Publication date |
---|---|
CN113376612A (zh) | 2021-09-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN113376612B (zh) | 一种基于地形矩阵化及探测的雷达杂波生成方法 | |
US8358239B2 (en) | Iterative clutter calibration with phased array antennas | |
Sevgi et al. | An integrated maritime surveillance system based on high-frequency surface-wave radars. 1. Theoretical background and numerical simulations | |
US11391836B2 (en) | Liveliness detection using radar | |
US7136012B2 (en) | Approach radar with array antenna having rows and columns skewed relative to the horizontal | |
KR20190035895A (ko) | 3d 프린트된 루네부르크 렌즈를 사용하는 새로운 차량 레이더 | |
CN104535996B (zh) | 图像/激光测距/低空调频连续波雷达一体化系统 | |
CN114781190B (zh) | 一种雷达探测能力仿真方法及装置 | |
CN107918115B (zh) | 基于多径利用的雷达目标定位方法 | |
GB2496030A (en) | Monopulse radar | |
US20230057370A1 (en) | Method and device for determining an operational geographical zone observed by a sensor | |
CN113191027A (zh) | 基于地形二维、三维矩阵化辐射源信号强度分布计算方法 | |
CN113189562B (zh) | 一种基于高程的地形探测算法 | |
CN108445461A (zh) | 一种多径条件下雷达目标检测方法 | |
Foessel-Bunting | Radar sensor model for three-dimensional map building | |
Ji et al. | Ship monitoring with bistatic compact HFSWR of small aperture | |
Gustafsson | Positioning of objects behind corners using X-band radar | |
Honda et al. | Estimation of target detection rate in aircraft surveillance system | |
Oreshkina et al. | Digital Earth surface maps for radar ground clutter simulation | |
JP6381856B2 (ja) | レーダ信号処理装置及びレーダ信号処理方法 | |
Mecocci et al. | Radar image processing for ship-traffic control | |
CN113504536A (zh) | 一种复杂探测环境下的外辐射源雷达基站直达波获取方法 | |
KR101800281B1 (ko) | 수치 지형 표고 데이터를 이용한 초고주파 영상의 표적 탐지 장치 및 방법 | |
CN113093112B (zh) | 一种面向地面装备在有障碍遮挡下的侦察包络计算方法、装置及存储介质 | |
Wang et al. | Passive distributed arrays cross-location algorithm with windowed deconvolution |
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 | ||
PE01 | Entry into force of the registration of the contract for pledge of patent right | ||
PE01 | Entry into force of the registration of the contract for pledge of patent right |
Denomination of invention: A Method of Radar Clutter Generation Based on Terrain Matrix and Detection Effective date of registration: 20230217 Granted publication date: 20211123 Pledgee: Chengdu SME financing Company Limited by Guarantee Pledgor: Chengdu zhongxiangtiandi Network Technology Co.,Ltd. Registration number: Y2023980032953 |