CN106646416A - 一种机场跑道异物检测方法 - Google Patents
一种机场跑道异物检测方法 Download PDFInfo
- Publication number
- CN106646416A CN106646416A CN201611122439.7A CN201611122439A CN106646416A CN 106646416 A CN106646416 A CN 106646416A CN 201611122439 A CN201611122439 A CN 201611122439A CN 106646416 A CN106646416 A CN 106646416A
- Authority
- CN
- China
- Prior art keywords
- frequency
- target
- spectrum
- component
- signal
- 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/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
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、对中频信号进行有限长度的FFT变换,得到中频信号的频谱;S2、根据频谱估计第一频谱分量的频率,判断第一频谱分量的频率是否为FFT变换的频率分辨率的整数倍:若是则将中频信号的第一频谱分量的频率作为其对应的目标的实际频率;若否则对第一频谱分量的频率进行频率插值,将频率插值结果作为其对应的目标的实际频率;S3、根据目标的实际频率得到目标的幅度和相位,并计算目标的距离;S4、将中频信号减去目标的实际频率对应的余弦信号得到的剩余信号作为新的中频信号,转入步骤S1,直到得到所有目标的距离。本发明提高提高目标检测精度。
Description
技术领域
本发明涉及机场安全防护技术领域。更具体地,涉及一种机场跑道异物检测方法。
背景技术
机场跑道异物(Foreign Object Debris,FOD),指可能损伤航空器的某种外来物质、碎屑或物体,如金属零件、碎石块、纸屑、树叶等。FOD严重影响飞机的安全,据保守估计,每年全球因FOD造成的损失约为30-40亿美元。近年来,随着民航事业的不断发展,FOD的自动监测也引起了世界范围内的关注。
目前存在的FOD监测***,如英国的Tarsier,以色列的FODetect,美国的FODFinder等***主要采用毫米波雷达或毫米波雷达/光学复合模式。毫米波雷达具有高分辨率和小尺寸的优点,可全天候、全天时自动监测机场跑道上的异物,比较适合安装于机场跑道附近,能够准确地探测到机场跑道上的外来异物。上述FOD监测***均采用线性调频连续波雷达体制,这种体制的雷达发射功率低,接收机灵敏度高,探测距离远,无距离盲区,能够检测到小目标。
现有技术中的FOD监测***检测静止目标是将目标检测转化为与目标距离相关的频率参数估计。频率参数估计的方法可分为时域和频域两类,时域方法估计涉及多参数计算,而且计算量较大,无法实现实时性;相比而言,频域方法主要基于对时域信号进行离散傅里叶变换,这可通过快速傅里叶变换(FFT)来实现加速,同时不涉及多参数估计问题,因而检测速度更快,更能满足机场跑道异物快速检测要求。
尽管基于FFT的频域估计方法具有快速计算的优点,但FFT方法自身存在两点不足:1)有限采样时间导致的频谱泄露;2)频率离散化带来的栅栏效应。这都会影响频率估计的准确度,进而影响目标检测率。此外,该方法无法检测出强目标干扰下的弱目标。
因此,需要提供一种机场跑道异物检测方法。
发明内容
本发明的目的在于提供一种机场跑道异物检测方法,以提高强目标干扰下的弱目标检测概率。
为达到上述目的,本发明采用下述技术方案:
一种机场跑道异物检测方法,该方法包括如下步骤:
S1、对目标反射信号经混频后得到的中频信号进行有限长度的FFT变换,得到中频信号的频谱;
S2、根据中频信号的频谱估计中频信号的第一频谱分量的频率,判断第一频谱分量的频率是否为FFT变换的频率分辨率的整数倍:若是则将中频信号的第一频谱分量的频率作为第一频率分量对应的目标的实际频率;若否则对第一频谱分量的频率进行频率插值,将频率插值后得到的第一频谱分量的频率作为第一频率分量对应的目标的实际频率,从而得到更准确的第一频率分量对应的目标的实际频率,提高目标检测精度;
S3、根据第一频率分量对应的目标的实际频率得到第一频率分量对应的目标的幅度和相位,并计算第一频率分量对应的目标的距离;
S4、将中频信号减去第一频率分量对应的目标的实际频率对应的余弦信号得到的剩余信号作为新的中频信号,转入步骤S1,直到得到所有目标的距离。
优选地,步骤S2中所述根据中频信号的频谱估计中频信号的第一频谱分量的频率的公式为:
其中,k1为第一频谱分量最大值的频率,fs为采样频率,N为FFT变换的长度。
优选地,FFT变换的频率分辨率为:Δf=fs/N。
优选地,步骤S2中对第一频谱分量的频率进行频率插值的具体过程为:
目标的实际频率的取值公式为:
其中,δ为目标的实际频率偏离第一频谱分量最大值的频率k1的偏离值;
第一频谱分量的频率与相邻频谱分量的频率比例参数如下:
α1=R{SIF(k1-1)/SIF(k1)}
α2=R{SIF(k1+1)/SIF(k1)}
其中,R{·}为代表取复数的实部;
计算偏离值δ的两个估计值如下:
δ1=α1/(1-α1)
δ2=-α2/(1-α2)
通过偏离值δ的两个估计值计算得到偏离值δ的取值为:
将偏离值δ代入目标的实际频率的取值公式计算得到目标的实际频率
优选地,步骤S3中计算第一频率分量对应的目标的距离的计算公式为:
其中,α为调频斜率,c为光速。
本发明的有益效果如下:
本发明所述技术方案该采用频域插值与时域窗函数的方法来降低FFT技术自身的缺陷,通过插值利用离散点频谱信息获取目标的实际频率(实际频谱位置),获得目标的距离,进而可获得目标的位置信息;同时采用频率分量分离的方法提高强目标干扰下的弱目标检测概率,通过多次迭代得到所有目标的距离和位置信息。
附图说明
下面结合附图对本发明的具体实施方式作进一步详细的说明;
图1示出机场跑道异物检测方法的流程图。
图2示出频率插值示意图。
具体实施方式
为了更清楚地说明本发明,下面结合优选实施例和附图对本发明做进一步的说明。附图中相似的部件以相同的附图标记进行表示。本领域技术人员应当理解,下面所具体描述的内容是说明性的而非限制性的,不应以此限制本发明的保护范围。
如图1所示,本发明公开的一种机场跑道异物检测方法,包括如下步骤:
S1、对目标反射信号经混频后得到的中频信号进行有限长度的FFT变换,得到中频信号的频谱;
S2、根据中频信号的频谱估计中频信号的第一频谱分量的频率,判断第一频谱分量的频率是否为FFT变换的频率分辨率的整数倍:若是则将中频信号的第一频谱分量的频率作为第一频率分量对应的目标的实际频率;若否则对第一频谱分量的频率进行频率插值,将频率插值后得到的第一频谱分量的频率作为第一频率分量对应的目标的实际频率,从而得到更准确的第一频率分量对应的目标的实际频率,提高目标检测精度;
S3、根据第一频率分量对应的目标的实际频率得到第一频率分量对应的目标的幅度和相位,并计算第一频率分量对应的目标的距离;
S4、将中频信号减去第一频率分量对应的目标的实际频率对应的余弦信号得到的剩余信号作为新的中频信号,转入步骤S1,直到得到所有目标的距离。
本方案中,步骤S1的具体过程为:
FOD监测***雷达的发射信号(线性调频连续波)的信号模型为:
其中,f0为发射信号的起始频率,α=B/T为发射信号的调频斜率,B为扫频带宽,T为扫频周期,为发射信号的初相。
FOD监测***雷达的发射信号照射到目标后反射回来的目标反射信号的信号模型为:
其中,p为目标个数,i=1,2,…p为第i个目标的回波时延,ri为第i个目标到FOD监测***雷达的距离,c为光速。
目标反射信号经过混频后得到的中频信号的信号模型为:
从式(3)可以看出,中频信号是由p个余弦信号叠加而成的,每个余弦信号对应的频率为:
这样,目标的距离与频率存在如下关系
对式(3)得到的中频信号进行长度为N的FFT变换,得到中频信号的频谱SIF(K)。
本方案中,步骤S2的具体过程为:
根据中频信号的频谱估计中频信号的第一频谱分量的频率
其中,k1为第一频谱分量最大值的频率(频谱位置),fs为采样频率;
判断第一频谱分量的频率是否为FFT变换的频率分辨率Δf=fs/N的整数倍:
若是则将中频信号的第一频谱分量的频率作为第一频率分量对应的目标的实际频率;
若否,则目标的实际频率会因FFT的栅栏效应而位于离散频谱相邻两谱线的中间,设目标的实际频率为:
其中,δ为目标的实际频率(实际频谱位置)偏离第一频谱分量最大值的频率k1的偏离值,如图2所示。
下面给出偏离值δ的估计方法,主要思路是利用第一频谱分量两边的频谱信息通过插值来获取目标的实际频率的优化值:
第一频谱分量的频率与相邻频谱分量的频率比例参数如下:
α1=R{SIF(k1-1)/SIF(k1)} (8)
α2=R{SIF(k1+1)/SIF(k1)} (9)
其中,R{·}为代表取复数的实部;
由式(8)、(9)计算偏离值δ的两个估计值如下:
δ1=α1/(1-α1) (10)
δ2=-α2/(1-α2) (11)
通过偏离值δ的两个估计值计算得到偏离值δ的取值为
将偏离值δ代入式(7)计算得到目标的实际频率
本方案中,步骤S4的具体过程为:
将第一频率分量对应的目标的实际频率代入中频信号的频谱中得到第一频率分量对应的目标的幅度和相位;
并计算第一频率分量对应的目标的距离,由式(5)和(6)可得,第一频率分量对应的目标的距离的计算公式为:
因此,对中频信号进行FFT变化就能得到目标的频谱分布,如(6)所示,再结合式(13)便可得到目标的距离,进而可通过距离得到目标的位置信息。
本方案中,步骤S5的具体过程为:
从式(3)中发现,中频信号是由多个余弦信号叠加而成的,对其进行FFT后会得到若干频率分量,这些频率分量的幅度值各不相同,对于弱目标,其对应的频谱分量就会受到来自强目标频谱分量的干扰,增加检测难度,从而产生漏检。为提高机场跑道弱目标检测概率,提升FOD监测***检测能力,对强弱目标同时存在的场景,采用频率分量分离方法来完成检测,具体为:
将中频信号减去2)减去步骤S2得到的第一频率分量对应的目标的实际频率对应的余弦信号得到剩余信号,转入步骤S1,即,将剩余信号作为新的中频信号,对该剩余信号再进行上述计算频谱、估计中频信号的第一频谱分量的频率、频率插值、得到第一频率分量对应的目标的幅度和相位并计算距离等步骤,直到得到所有目标的距离。
显然,本发明的上述实施例仅仅是为清楚地说明本发明所作的举例,而并非是对本发明的实施方式的限定,对于所属领域的普通技术人员来说,在上述说明的基础上还可以做出其它不同形式的变化或变动,这里无法对所有的实施方式予以穷举,凡是属于本发明的技术方案所引伸出的显而易见的变化或变动仍处于本发明的保护范围之列。
Claims (5)
1.一种机场跑道异物检测方法,其特征在于,该方法包括如下步骤:
S1、对目标反射信号经混频后得到的中频信号进行有限长度的FFT变换,得到中频信号的频谱;
S2、根据中频信号的频谱估计中频信号的第一频谱分量的频率,判断第一频谱分量的频率是否为FFT变换的频率分辨率的整数倍:若是则将中频信号的第一频谱分量的频率作为第一频率分量对应的目标的实际频率;若否则对第一频谱分量的频率进行频率插值,将频率插值后得到的第一频谱分量的频率作为第一频率分量对应的目标的实际频率;
S3、根据第一频率分量对应的目标的实际频率得到第一频率分量对应的目标的幅度和相位,并计算第一频率分量对应的目标的距离;
S4、将中频信号减去第一频率分量对应的目标的实际频率对应的余弦信号得到的剩余信号作为新的中频信号,转入步骤S1,直到得到所有目标的距离。
2.根据权利要求1所述的机场跑道异物检测方法,其特征在于,步骤S2中所述根据中频信号的频谱估计中频信号的第一频谱分量的频率的公式为:
其中,k1为第一频谱分量最大值的频率,fs为采样频率,N为FFT变换的长度。
3.根据权利要求2所述的机场跑道异物检测方法,其特征在于,FFT变换的频率分辨率为:Δf=fs/N。
4.根据权利要求3所述的机场跑道异物检测方法,其特征在于,步骤S2中对第一频谱分量的频率进行频率插值的具体过程为:
目标的实际频率的取值公式为:
其中,δ为目标的实际频率偏离第一频谱分量最大值的频率k1的偏离值;
第一频谱分量的频率与相邻频谱分量的频率比例参数如下:
α1=R{SIF(k1-1)/SIF(k1)}
α2=R{SIF(k1+1)/SIF(k1)}
其中,R{·}为代表取复数的实部;
计算偏离值δ的两个估计值如下:
δ1=α1/(1-α1)
δ2=-α2/(1-α2)
通过偏离值δ的两个估计值计算得到偏离值δ的取值为:
将偏离值δ代入目标的实际频率的取值公式计算得到目标的实际频率
5.根据权利要求4所述的机场跑道异物检测方法,其特征在于,步骤S3中计算第一频率分量对应的目标的距离的计算公式为:
其中,α为调频斜率,c为光速。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611122439.7A CN106646416B (zh) | 2016-12-08 | 2016-12-08 | 一种机场跑道异物检测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611122439.7A CN106646416B (zh) | 2016-12-08 | 2016-12-08 | 一种机场跑道异物检测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106646416A true CN106646416A (zh) | 2017-05-10 |
CN106646416B CN106646416B (zh) | 2019-04-05 |
Family
ID=58819289
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201611122439.7A Active CN106646416B (zh) | 2016-12-08 | 2016-12-08 | 一种机场跑道异物检测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106646416B (zh) |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0886149B1 (en) * | 1997-06-18 | 2004-09-08 | Honda Giken Kogyo Kabushiki Kaisha | Frequency-modulated continuous-wave radar system |
CN103777199A (zh) * | 2014-02-24 | 2014-05-07 | 中国科学院电子学研究所 | 一种调频连续波雷达***的测距方法 |
CN105301591A (zh) * | 2015-10-22 | 2016-02-03 | 上海无线电设备研究所 | 一种道路交通监控雷达及其实现方法 |
CN105549110A (zh) * | 2015-12-08 | 2016-05-04 | 北京无线电计量测试研究所 | 一种机场跑道异物检测设备和方法 |
-
2016
- 2016-12-08 CN CN201611122439.7A patent/CN106646416B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0886149B1 (en) * | 1997-06-18 | 2004-09-08 | Honda Giken Kogyo Kabushiki Kaisha | Frequency-modulated continuous-wave radar system |
CN103777199A (zh) * | 2014-02-24 | 2014-05-07 | 中国科学院电子学研究所 | 一种调频连续波雷达***的测距方法 |
CN105301591A (zh) * | 2015-10-22 | 2016-02-03 | 上海无线电设备研究所 | 一种道路交通监控雷达及其实现方法 |
CN105549110A (zh) * | 2015-12-08 | 2016-05-04 | 北京无线电计量测试研究所 | 一种机场跑道异物检测设备和方法 |
Non-Patent Citations (2)
Title |
---|
P. FEIL ET AL.: "Foreign Objects Debris Detection (FOD) on Airport Runways Using a Broadband 78 GHz Sensor", 《2008 38TH EUROPEAN MICROWAVE CONFERENCE》 * |
吴静 等: "多目标环境下跑道异物监测雷达的CFAR检测技术", 《光电工程》 * |
Also Published As
Publication number | Publication date |
---|---|
CN106646416B (zh) | 2019-04-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101718870B (zh) | 图像域的高速、微弱目标航迹检测方法 | |
CN101963662B (zh) | 基于短时分数阶傅里叶域滤波的自聚焦预处理方法 | |
CN104898103B (zh) | 基于多通道杂波图的低速目标检测方法 | |
CN104730518B (zh) | 一种基于高斯拟合的雷达多普勒谱估计海面流场的方法 | |
CN104166135A (zh) | 一种宽带雷达目标的原始点迹凝聚处理方法 | |
CN103278820B (zh) | 临近空间慢速平台sar动目标检测方法及成像方法 | |
CN106338727A (zh) | 一种车载辅助驾驶雷达的目标检测方法 | |
CN107861117A (zh) | 一种适用于连续波周界监视雷达的多目标参数测量方法 | |
CN105487060B (zh) | 一种双通道四斜率调制的多目标提取方法 | |
CN103293521B (zh) | 一种利用x波段雷达探测近海海域水深的方法 | |
JP2010197185A (ja) | 降水分布の推定システムおよび降水分布の推定方法 | |
CN103344949A (zh) | 基于Radon-线性正则模糊函数的雷达微弱动目标检测方法 | |
CN107358250A (zh) | 基于双波段雷达微多普勒融合的人体步态识别方法及*** | |
CN104316936B (zh) | 一种综合的dme脉冲干扰抑制方法 | |
CN102798855B (zh) | 一种基于数字电视信号的直升机目标识别方法 | |
CN106814359B (zh) | 一种机场跑道异物检测雷达信号预处理方法 | |
CN105093215A (zh) | 基于多普勒信息的雷达对低空慢速小目标的跟踪方法 | |
CN104483668A (zh) | 一种高精度雷达信号检测和跟踪***及其方法 | |
CN101498788A (zh) | 一种逆合成孔径雷达的目标转角估计和横向定标方法 | |
CN105044691A (zh) | 一种海杂波背景下的快速雷达性能评估方法 | |
Nai et al. | On the mitigation of wind turbine clutter for weather radars using range‐Doppler spectral processing | |
CN104714225A (zh) | 一种基于广义似然比的动态规划检测前跟踪方法 | |
CN106125088A (zh) | 基于激光雷达测深***确定海水深度的方法 | |
CN106872969A (zh) | 基于mtd脉冲积累及滑动处理的雷达目标角度估计方法 | |
CN105259537A (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 |