CN106569208B - 基于杂波图的机场跑道异物检测方法 - Google Patents
基于杂波图的机场跑道异物检测方法 Download PDFInfo
- Publication number
- CN106569208B CN106569208B CN201611004864.6A CN201611004864A CN106569208B CN 106569208 B CN106569208 B CN 106569208B CN 201611004864 A CN201611004864 A CN 201611004864A CN 106569208 B CN106569208 B CN 106569208B
- Authority
- CN
- China
- Prior art keywords
- clutter map
- unit
- detection
- clutter
- detection unit
- 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.)
- Expired - Fee Related
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
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/88—Radar or analogous systems specially adapted for specific applications
- G01S13/885—Radar or analogous systems specially adapted for specific applications for ground probing
-
- 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
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V3/00—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
- G01V3/12—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation operating with electromagnetic waves
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- General Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Computer Networks & Wireless Communication (AREA)
- Life Sciences & Earth Sciences (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geophysics (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明涉及一种基于杂波图的机场跑道异物检测方法,首先给定线性调频连续波雷达的参数,并给定无FOD跑道回波;然后,根据性能要求对回波进行下变频、混频、低通滤波、放大快速傅里叶变换后得到频率‑幅度谱,建立相应杂波图;检测时先根据杂波图特性生成相应阀值;再以波束为单位进行信号处理将处理,得到检测单元回波频谱特性强度;根据GO‑CFAR检测判断是否存在异物,如果存在异物输出该检测单元的坐标,如果没有异物则更新杂波图。本发明不会因杂物而影响杂波图的刷新,保证了杂波图的完整性,提高了雨天雾天等气象条件的适应性。
Description
技术领域
本发明属于雷达技术领域,特别涉及一种基于杂波图的机场跑道异物检测识别方法。
技术背景
机场跑道水泥地面和边灯等固有物体的回波反射会形成杂波,和目标回波混叠在一起进入接收通道,影响信号检测,造成虚警。
杂波图处理技术通过对检测单元进行多次扫描测量值进行迭代和平均得到检测单元处的杂波背景的强度。杂波图存储每个杂波单元的背景强度估计值,并进行迭代更新。一般情况下,杂波虽在距离和方位上变化十分剧烈,但同一距离单元上的杂波强度随时间的变化是缓慢的,故而可以使用杂波图作为检测基准与待检测信号进行比对,降低虚警发生。
杂波图最早由Nitzberg等人提出,是一种时域上的自适应处理方法。这种方法把雷达空间分成杂波图单元进行工作。杂波图技术首先将整个检测区域划分成若干个杂波图单元,对于每个单元,利用其自身有限次的回波输入进行反复迭代,根据恒虚警检测技术决定阈值,最后对于某空域的目标回波可用其所在杂波图单元的存储的检测阀值进行检测。
现有静态更新杂波图检测方式应用于机场跑道时,由于静态更新杂波图检测方式的遗忘因子是固定的,无法保证不同检测单元的不同时刻遗忘因子均是较好的,易造成杂波图边缘不稳定。
发明内容
本发明的目的在于提出一种基于杂波图的机场跑道异物检测方法。
实现本发明目的的技术方案为:一种基于杂波图的机场跑道异物检测方法,包括以下步骤:
步骤1,给定跑道扇形检测区域范围,采用线性调频连续波雷达对扇形检测区域进行异物检测,通过混频、低通滤波、频率灵敏控制和快速傅里叶变换,提取回波信号中的频率-幅度谱,即距离-幅度谱;
步骤2,根据步骤1给定扇形检测区域范围内无异物的回波频率幅度谱,得到各个位置的雷达回波幅度值,建立雷达杂波图;设定波束单元有N个,杂波图距离单元有M个,则位于(m,n)的杂波图单元的估计值为1≤m≤M,1≤n≤N,lm为杂波图单元距离发射天线的距离,为发射天线的转动角度;
步骤3,根据杂波图的距离-幅度谱选取杂波图的分布模型,对于给定的虚警率,由公式得到GO-CFAR门限因子;
步骤4,检测时连续波雷达回波采用步骤1方式提取频率-幅度谱,读取位于(m,n)的杂波图观测值为采用步骤3得到的GO-CFAR门限因子,结合GO-CFAR判决法判定当前检测的杂波图单元中是否存在目标,如果有目标则其位置为转到步骤3继续检测下一个波束宽度,若没有目标则转到步骤5;
步骤5,根据第x次检测的杂波图观测值和第x-1次杂波图的估计值对相应检测单元的遗忘因子进行动态更新;
步骤6,根据相应检测单元的遗忘因子、第x次检测的杂波图观测值和第x-1次杂波图的估计值对相应检测单元的杂波图进行更新。
与现有技术相比,本发明的显著优点为:(1)本发明在一定范围内遗忘因子动态更新,减缓杂波图刷新速度;(2)本发明不会因杂物而影响杂波图的刷新,保证了杂波图的完整性;(3)本发明提高了杂波图边缘的稳定性,实现对不同信噪比的相对较好的杂波因子,实现各个检测单元均有较好的遗忘因子。
附图说明
图1为本发明基于杂波图的机场跑道异物检测方法的流程图。
图2为杂波图单元示意图。
图3为GO-CFAR判决示意图。
图4为遗忘因子更新流程图。
图5为杂波图更新流程图。
具体实施方式
结合图1,一种基于杂波图的机场跑道异物检测方法,包括以下步骤:
步骤1,给定跑道扇形检测区域范围,采用线性调频连续波雷达对扇形检测区域进行异物检测,通过混频、低通滤波、频率灵敏控制和快速傅里叶变换,提取回波信号中的频率-幅度谱,即距离-幅度谱;
步骤2,根据步骤1给定扇形检测区域范围内无异物的回波频率幅度谱,得到各个位置的雷达回波幅度值,建立雷达杂波图;设定波束单元有N个,杂波图距离单元有M个,则位于(m,n)的杂波图单元的估计值为1≤m≤M,1≤n≤N,lm为杂波图单元距离发射天线的距离,为发射天线的转动角度;具体为:
步骤2.1,假设雷达波束宽度夹角为雷达信号处理时采用2M点FFT,检测距离为S,则杂波图的检测单元为夹角为边长为S/M的拱形;
步骤2.2,所述检测单元杂波图的估计值是角度为距离为lm的雷达回波幅度,首次检测的估计值所述检测单元杂波图的观测值是连续波雷达回波信号处理后得到的相应位置的幅值。
步骤3,根据杂波图的距离-幅度谱选取杂波图的分布模型,对于给定的虚警率,由公式得到GO-CFAR门限因子;具体为:
M个样本数据组成的矢量的联合概率密度函数为:
xi为第i个相邻单元的观测值,β2为噪声功率;
设观测数据x的似然函数,记为ν,则有:
关于β2的导数为零,得
得到最大似然估计为:
得到GO-CFAR门限因子:
步骤4,检测时连续波雷达回波采用步骤1方式提取频率-幅度谱,读取位于(m,n)的杂波图观测值为采用步骤3得到的GO-CFAR门限因子,结合GO-CFAR判决法判定当前检测的杂波图单元中是否存在目标,如果有目标则其位置为转到步骤3继续检测下一个波束宽度,若没有目标则转到步骤5;
判定当前检测的杂波图单元中是否存在目标的具体过程为:
若检测信号的相应单元的回波强度X1,X2分别为检测单元左侧的门限估计和检测单元右侧的门限估计,则认为相应单元对应位置含有目标,检测单元以极坐标表示,转化为直角坐标后FOD坐标为转到步骤3继续检测下一个波束宽度;否则,认为相应检测单元中无目标,执行步骤5。
步骤5,如图4所示,根据第x次检测的杂波图观测值和第x-1次杂波图的估计值对相应检测单元的遗忘因子进行动态更新;遗忘因子为:
式中Λ为遗忘因子的上限,α为GO-CFAR门限因子,
步骤6,根据相应检测单元的遗忘因子、第x次检测的杂波图观测值和第x-1次杂波图的估计值对相应检测单元的杂波图进行更新;如图5所示,具体为
当雷达扫过一个波束宽度,通过对回波的检测及信号处理得到相应检测单元的杂波图观测值杂波图的估计值为得到当前周期内杂波图的估计值:
其中,为该次更新的遗忘因子。
进一步的,上述步骤1具体为:
信号处理时采用2M点FFT,测量距离为S,发射信号载频F,调频斜率k;发射信号为:
ff=Acos(F+kt)
A为雷达发射波强度,t为当前时间;
第j个距离单元回波信号为:
fc=ajcos(F+k(t-τj))
aj为第j个距离单元的回波强度,τj为第j个距离单元的时延;
则接收信号发射信号混频得:
ff*fc=A*aj/2(cos(2F+2kt-kτj)+cos(kτj))
经过低通滤波器可得,第j个距离单元的下变频Gj:
Gj=A*aj/2coskτj
故Δf=kτ,τ为时延,电磁波的传播速度为c,故频差与距离成正比;经过2M点快速傅里叶变换可得,单位脉冲的杂波图中的检测单元个数Mφ为:
(θ1,θ2)为给定探测范围,如图2所示为一个杂波图检测单元。
进一步的,上述步骤4具体为:
对雷达回波采用步骤1中方法进行信号处理到频率幅度谱,得到杂波图观测值由于噪声特性在时间和空间上剧烈变化时,不太可能出现相邻目标但是干扰严重不均匀的情况,故更关注杂波边缘的虚警,采用GO-CFAR判决准则:
X=max(X1,X2)
GO-CFAR的原理如图3所示,在GO-CFAR中,门限与发现概率Pd的关系为:
SNR为信噪比,T为门限值,m为单个参考单元的长度,即m=M/2;
当没有信号,也就是SNR=0时,门限与虚警概率Pfa的关系为
得到GO-CFAR进行判决则判决公式为:
当判决为H1时,则由
得出FOD出现的位置,然后检测下一组波束的回波;
当判决为H0,则转入步骤5。
Claims (3)
1.一种基于杂波图的机场跑道异物检测方法,其特征在于,包括以下步骤:
步骤1,给定跑道扇形检测区域范围,采用线性调频连续波雷达对扇形检测区域进行异物检测,通过混频、低通滤波、频率灵敏控制和快速傅里叶变换,提取回波信号中的频率-幅度谱,即距离-幅度谱;
步骤2,根据步骤1给定扇形检测区域范围内无异物的回波频率幅度谱,得到各个位置的雷达回波幅度值,建立雷达杂波图;设定波束单元有N个,杂波图距离单元有M个,则位于(m,n)的杂波图单元的估计值为1≤m≤M,1≤n≤N,lm为杂波图单元距离发射天线的距离,为发射天线的转动角度;
步骤3,根据杂波图的距离-幅度谱选取杂波图的分布模型,对于给定的虚警率,由公式得到GO-CFAR门限因子;具体过程为:
M个样本数据组成的矢量的联合概率密度函数为:
xi为第i个相邻单元的观测值,β2为噪声功率;
设观测数据x的似然函数,记为ν,则有:
关于β2的导数为零,得
得到最大似然估计为:
得到GO-CFAR门限因子:
步骤4,检测时连续波雷达回波采用步骤1方式提取频率-幅度谱,读取位于(m,n)的杂波图观测值为采用步骤3得到的GO-CFAR门限因子,结合GO-CFAR判决法判定当前检测的杂波图单元中是否存在目标,如果有目标则其位置为转到步骤3继续检测下一个波束宽度,若没有目标则转到步骤5;判定当前检测的杂波图单元中是否存在目标的具体过程为:
若检测信号的相应单元的回波强度X1,X2分别为检测单元左侧的门限估计和检测单元右侧的门限估计,则认为相应单元对应位置含有目标,检测单元以极坐标表示,转化为直角坐标后FOD坐标为转到步骤3继续检测下一个波束宽度;否则,认为相应检测单元中无目标,执行步骤5;
步骤5,根据第x次检测的杂波图观测值和第x-1次杂波图的估计值对相应检测单元的遗忘因子进行动态更新;所述遗忘因子为:
式中Λ为遗忘因子的上限,α为GO-CFAR门限因子,
步骤6,根据相应检测单元的遗忘因子、第x次检测的杂波图观测值和第x-1次杂波图的估计值对相应检测单元的杂波图进行更新。
2.根据权利要求书1所述的基于杂波图的机场跑道异物检测方法,其特征在于,步骤2的具体过程为:
步骤2.1,假设雷达波束宽度夹角为雷达信号处理时采用2M点FFT,检测距离为S,则杂波图的检测单元为夹角为边长为S/M的拱形;
步骤2.2,所述检测单元杂波图的估计值是角度为距离为lm的雷达回波幅度,首次检测的估计值
3.根据权利要求书1所述的基于杂波图的机场跑道异物检测方法,其特征在于,所述步骤6中杂波图更新的具体过程为:
当雷达扫过一个波束宽度,通过对回波的检测及信号处理得到相应检测单元的杂波图观测值杂波图的估计值为得到当前周期内杂波图的估计值:
其中,为该次更新的遗忘因子。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611004864.6A CN106569208B (zh) | 2016-11-15 | 2016-11-15 | 基于杂波图的机场跑道异物检测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611004864.6A CN106569208B (zh) | 2016-11-15 | 2016-11-15 | 基于杂波图的机场跑道异物检测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106569208A CN106569208A (zh) | 2017-04-19 |
CN106569208B true CN106569208B (zh) | 2019-04-12 |
Family
ID=58542295
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201611004864.6A Expired - Fee Related CN106569208B (zh) | 2016-11-15 | 2016-11-15 | 基于杂波图的机场跑道异物检测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106569208B (zh) |
Families Citing this family (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107247262B (zh) * | 2017-05-31 | 2019-11-19 | 中国电子科技集团公司第二十九研究所 | 一种机场跑道异物分层检测方法 |
CN108196238B (zh) * | 2018-01-12 | 2021-07-27 | 西安电子科技大学 | 高斯背景下基于自适应匹配滤波的杂波图检测方法 |
CN108445477B (zh) * | 2018-01-19 | 2022-02-01 | 成都赛英科技有限公司 | 机场场面异物检测雷达的高精度测距方法 |
CN108226890B (zh) * | 2018-01-19 | 2021-08-13 | 成都赛英科技有限公司 | 基于时间向统计的机场异物雷达检测方法 |
CN109143181A (zh) * | 2018-07-25 | 2019-01-04 | 南京理工大学 | 一种便于编程的雷达点迹航迹互联方法 |
CN108761400B (zh) * | 2018-08-24 | 2020-05-15 | 中科宇达(北京)科技有限公司 | 机场跑道fod检测雷达系统及处理方法 |
CN109884638A (zh) * | 2019-03-13 | 2019-06-14 | 苏州理工雷科传感技术有限公司 | 一种信号处理机、雷达系统及信号处理方法 |
CN110161473B (zh) * | 2019-05-08 | 2021-04-30 | 南京莱斯电子设备有限公司 | 一种基于基元杂波图与图层的杂波图检测与更新方法 |
CN110320514A (zh) * | 2019-06-12 | 2019-10-11 | 南京理工大学 | 基于车载侧视探测雷达的fod检测方法 |
CN110441766B (zh) * | 2019-07-02 | 2023-02-17 | 中国航空工业集团公司雷华电子技术研究所 | 一种机场道面fod探测雷达变门限检测方法 |
CN110286373B (zh) * | 2019-07-19 | 2022-05-24 | 湖南华诺星空电子技术有限公司 | 一种复杂气候条件下的fod雷达雨雪杂波抑制方法 |
CN111060880B (zh) * | 2019-12-30 | 2021-09-07 | 无锡市雷华科技有限公司 | 一种基于恒虚警检测原理的气象杂波抑制方法 |
CN112213695B (zh) * | 2020-10-15 | 2023-04-11 | 四川大学 | 一种基于单元平均杂波图的机场跑道异物检测方法 |
CN112731386B (zh) * | 2020-12-28 | 2024-02-20 | 西南技术物理研究所 | 基于杂波图匹配的机场跑道异物检测方法 |
CN115575946B (zh) * | 2022-10-20 | 2024-06-04 | 中国人民解放军海军航空大学 | 雷达目标斜对称子空间智能融合检测方法 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4468638A (en) * | 1982-02-11 | 1984-08-28 | The Bendix Corporation | Linear sweep frequency modulator for FM/CW radio altimeter |
US6757668B1 (en) * | 1999-11-05 | 2004-06-29 | General Electric Company | Information fusion of classifiers in systems with partial redundant information |
CN103149563B (zh) * | 2013-01-28 | 2014-10-08 | 北京航空航天大学 | 机场跑道异物二维探测系统及方法 |
CN103675923B (zh) * | 2013-12-17 | 2016-07-13 | 无锡市雷华科技有限公司 | 基于静态杂波图的异物检测方法及异物检测设备 |
CN105389921B (zh) * | 2015-12-09 | 2019-03-26 | 中国民用航空总局第二研究所 | 一种机场跑道异物的监测系统和方法 |
-
2016
- 2016-11-15 CN CN201611004864.6A patent/CN106569208B/zh not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
CN106569208A (zh) | 2017-04-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106569208B (zh) | 基于杂波图的机场跑道异物检测方法 | |
US11650286B2 (en) | Method for separating targets and clutter from noise, in radar signals | |
CN105445734B (zh) | 具有基于相位的多目标检测的雷达系统 | |
CN104898103B (zh) | 基于多通道杂波图的低速目标检测方法 | |
CN109490874A (zh) | 确定雷达目标作为位置地标的适合性的方法 | |
CN104502907A (zh) | 一种机载雷达地面动静目标稳定跟踪方法 | |
CN107576959B (zh) | 一种基于区域映射解模糊的高重频雷达目标检测前跟踪方法 | |
JP2018205174A (ja) | レーダ装置及びそのレーダ信号処理方法 | |
CN108693523A (zh) | 基于锯齿波线性调频雷达多目标测距测速的方法及系统 | |
CN108535730A (zh) | 一种多普勒气象雷达解速度模糊方法和系统 | |
CN107942305A (zh) | 双偏振雷达系统初始差分相位的在线标定方法 | |
Onunka et al. | Autonomous marine craft navigation: On the study of radar obstacle detection | |
Jose et al. | Millimetre wave radar spectra simulation and interpretation for outdoor slam | |
KR20150059864A (ko) | 엑스밴드 이중편파 레이더 원시자료를 이용한 강수체 구분 방법 | |
Liang | Waveform design and diversity in radar sensor networks: theoretical analysis and application to automatic target recognition | |
KR101885922B1 (ko) | 레이더 비디오 신호를 이용하는 클러터 모델 생성 장치, 시스템 및 방법 | |
CN113608193A (zh) | 一种基于UNet的雷达多目标距离和速度估计方法 | |
CN113419238A (zh) | 基于毫米波雷达的山体滑坡监测方法、电子设备、存储介质 | |
Laws et al. | Ship tracking by HF radar in coastal waters | |
Yoshikawa et al. | Dual-directional radar observation for preliminary assessment of the Ku-band broadband radar network | |
CN114280612B (zh) | 一种绝缘体目标的毫米波雷达恒虚警检测方法 | |
Abeynayake et al. | Ground penetrating radar applications in buried improvised explosive device detection | |
McCann et al. | Visualising the aspect-dependent radar cross section of seabirds over a tidal energy test site using a commercial marine radar system | |
CN113093174B (zh) | 基于phd滤波雷达起伏微弱多目标的检测前跟踪方法 | |
CN104297735A (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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20190412 Termination date: 20201115 |
|
CF01 | Termination of patent right due to non-payment of annual fee |