CN101325007A - 基于毫米波交通雷达的公路背景识别与公路背景抑制方法 - Google Patents

基于毫米波交通雷达的公路背景识别与公路背景抑制方法 Download PDF

Info

Publication number
CN101325007A
CN101325007A CNA2008100409535A CN200810040953A CN101325007A CN 101325007 A CN101325007 A CN 101325007A CN A2008100409535 A CNA2008100409535 A CN A2008100409535A CN 200810040953 A CN200810040953 A CN 200810040953A CN 101325007 A CN101325007 A CN 101325007A
Authority
CN
China
Prior art keywords
highway
background
power spectrum
highway background
power
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
Application number
CNA2008100409535A
Other languages
English (en)
Other versions
CN101325007B (zh
Inventor
余稳
张�浩
孙晓玮
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shanghai Institute of Microsystem and Information Technology of CAS
Original Assignee
Shanghai Institute of Microsystem and Information Technology of CAS
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Shanghai Institute of Microsystem and Information Technology of CAS filed Critical Shanghai Institute of Microsystem and Information Technology of CAS
Priority to CN2008100409535A priority Critical patent/CN101325007B/zh
Publication of CN101325007A publication Critical patent/CN101325007A/zh
Application granted granted Critical
Publication of CN101325007B publication Critical patent/CN101325007B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Radar Systems Or Details Thereof (AREA)
  • Traffic Control Systems (AREA)

Abstract

本发明涉及一种基于毫米波交通雷达的公路背景识别与公路背景抑制方法,利用同时含有车辆信息和背景信息的交通雷达回波信号实现了公路背景功率谱自动识别和实时更新,并完成了公路背景抑制。该方法以侧向毫米波交通雷达的时域数字回波信号作为信号输入,应用周期图法得到随时间变化的雷达回波功率谱时间序列;在每个频点处将功率谱时间序列按照功率谱幅度升序排列;选择部分幅度有序功率谱序列作为公路背景功率谱估计的有效数据;对有效幅度有序功率谱序列进行相干平均;利用自学习方法实现公路背景功率谱的实时更新;进行公路背景抑制。本发明能准确实时地识别公路背景功率谱,并能够很好的抑制公路背景,大幅度提高了雷达回波信号的信噪比。

Description

基于毫米波交通雷达的公路背景识别与公路背景抑制方法
技术领域
本发明设计智能交通中的交通信息采集领域,具体涉及了一种基于毫米波交通雷达的公路背景识别与公路背景抑制方法,适用于各种交通雷达的背景识别和背景抑制。
背景技术
交通信息采集系统是智能交通系统不可缺少的基本组成部分。交通信息采集系统依赖于交通信息传感器,传感器实时地准确地获取交通数据的能力对交通控制、公共安全和交通规划都是至关重要的。随着经济和社会的发展,现代城市却反而变得更加拥挤,交通问题引起了广泛的关注。生活节奏的加快使人们更加期望拥有更为良好的交通状况。
许多种交通信息采集设备已经得到了广泛的运用,如掩埋式线圈检测器已经有了五十多年的应用历史,并应用至今。但是,此类掩埋式交通信息检测器通常都有较高的失效率,并有较高的安装和维护成本。更为糟糕的是,此类掩埋式交通信息检测器的安装与维护会导致略面破坏和车道关闭。近几十年来,许多种非掩埋式的交通信息检测设备得到了广泛的应用,这其中包括声学和光学交通信息传感器。但是,这些交通信息传感器有着明显的缺点。首先,这些传感器只能检测单条行车道,而对检测整条公路的交通状况无能为力;第二,这类传感器通常对外部环境变化非常敏感,比如天气变化,昼夜变化等等。
毫米波交通雷达的出现在一定程度上解决了这些问题,通过分析不同车道对应的不同频率段的功率变化,毫米波交通雷达可以实现多车道交通信息检测的功能。毫米波交通雷达有许多其他交通信息检测器难以比拟的优点,如全天候工作,无检测盲区,低安装维护成本,不对公路路面造成损害等等。毫米波交通雷达的这些优良特性使得它成为了掩埋式线圈和声光交通信息检测器的良好的替代品。
但在实际应用中,公路背景会导致功率混叠,严重影响了毫米波交通信息检测雷达的检测性能。这种功率混叠在公路上存在大型静态背景时尤为严重,如公路上存在金属防护栏或绿化带等等。在这种情况下,背景功率谱的准确识别和背景抑制就成为了提高毫米波交通雷达检测精度和性能的首要问题。在应用毫米波交通信息检测器时,准确的识别公路背景并进行有效的背景抑制是非常必要的。
申请号为02136541.5,公开号为CN 1400475,名称为“智能交通中车流量监控的雷达射频装置及测定方法”的发明专利申请公开说明书,公开了一种毫米波车流量检测雷达的硬件系统以及交通信息检测方法。该方法中并没有涉及背景识别与背景抑制。
发明内容
本发明的目的在于提供一种毫米波交通雷达的公路背景识别与公路背景抑制方法,其利用同时含有车辆信息和背景信息的数据实时的估计公路背景,并实现公路背景抑制,提高毫米波交通信息检测雷达的准确性。
本发明具体是这样实现的:
基于毫米波交通雷达的公路背景识别与公路背景抑制方法,包括以下步骤:
步骤1、利用毫米波交通雷达的时域数字回波信号求得雷达回波功率谱时间序列;
步骤2、在每个频点处将雷达回波功率值时间序列按照功率谱幅度升序排列;
步骤3、滤除干扰序列,保留幅度低于有效阀值的有序功率谱序列作为公路背景功率谱的有效数据;
步骤4、对有效幅度有序功率谱序列进行相干平均,得到当前公路背景功率谱;
步骤5、公路背景功率谱通过自学习方法进行实时更新;
步骤6、利用谱相减法对公路背景进行抑制。
所述步骤1进一步包括:
采用N=2n点矩形窗w(n)截取雷达数字回波信号SIF,求得快速傅立叶变换为:
Ssamp(k)=SIF(k)*W(k)。
由于雷达回波信号功率谱的正频部分与负频部分严格对称,可舍去负频部分,得到的雷达回波功率幅度为:
P R ( k ) = 1 N | S samp ( k ) | 2 .
交通雷达回波信号功率通常包含车辆回波信号功率Pv,背景回波信号功率Pb和白噪声功率Pn。频率点k处的雷达回波信号功率可以表示为:
PR(k)=Pv(k)+Pb(k)+Pn(k)。
白噪声功率符合N(0,σ2)的高斯分布,则白噪声功率Pn(k)=σ2为很小的恒定值,可以忽略。这样,雷达回波信号功率仅由运动着的车辆产生的动态功率和相对静止的公路背景产生的静态功率组成。
所述步骤3进一步包括:
回波信号幅度有序功率序列的数据空间归一化;
归一化有效数据的选择。
所述步骤4具体如下:
在有效幅度有序功率序列的归一化数据空间内进行相干平均,如下式表示:
1 M Σ i = 0 M P Ros ( k , i ) ≈ 1 M Σ i = 0 M P vos ( k , i ) + P b 0 ( k ) ,
式中PRos为幅度有序功率,Pvos为幅度有序功率的车辆回波部分,Pb0为当前公路背景功率谱,M为归一化有效数据容量。在归一化有效数据空间内,车辆产生的动态功率 Σ i = 0 M P v ( k , i ) → 0 , 当前公路背景功率谱如下式表示:
P b 0 ( k ) ≈ 1 M Σ i = 0 M P Ros ( k , i ) .
所述步骤5具体如下:
在每个频点处,实时更新公路背景功率谱的自学习过程如下式表示:
P b t ( k ) = ( 1 - r ) P 0 t - 1 ( k ) + r P 0 t ( k ) .
式中Pb为公路背景功率谱,r为学习速率。
所述步骤6具体如下:
在每个频点处,用于完成公路背景抑制的谱相减过程可如下式表示:
Figure A20081004095300085
式中χ为谱相减系数,γ1、γ2为谱相减指数。
所述的谱相减系数χ满足:χ∈[1,2];
所述的谱相减指数满足γ1、γ2:γ1=1/γ2=2。
本发明能够准确实时地识别公路背景功率谱,并能够很好的抑制公路背景,大幅度提高了雷达回波信号的信噪比。
附图说明
图1为本发明的一个实施例的背景识别和背景抑制的流程图。
图2为本发明的一个实施例中的交通雷达安装和公路环境示意图。
图3为本发明的一个实施例中的公路背景识别结果。
图4为本发明的一个实施例中出现的公路背景功率与车辆功率发生功率混叠时的公路回波功率谱。
图5为本发明的一个实施例中的公路背景抑制结果。
具体实施方式
以下结合附图及实施例对本发明作进一步描述。
一种基于毫米波交通雷达的公路背景识别与公路背景抑制方法,包括以下步骤:
步骤1、利用毫米波交通雷达的时域数字回波信号求得雷达回波功率谱时间序列;
步骤2、在每个频点处将雷达回波功率值时间序列按照功率谱幅度升序排列;
步骤3、滤除车辆所产生的干扰序列(即滤除幅度大于有效阀值的有序功率谱序列),保留幅度低于有效阀值的有序功率谱序列作为公路背景功率谱的有效数据;
步骤4、对有效幅度有序功率谱序列进行相干平均,得到当前公路背景功率谱;
步骤5、公路背景功率谱通过自学习方法进行实时更新;
步骤6、利用谱相减法对公路背景进行抑制。
所述步骤1进一步包括:
采用N=2n点矩形窗w(n)截取雷达数字回波信号SIF,求得快速傅立叶变换为:
Ssamp(k)=SIF(k)*W(k)。
由于雷达回波信号功率谱的正频部分与负频部分严格对称,可舍去负频部分,得到的雷达回波功率幅度为:
P R ( k ) = 1 N | S samp ( k ) | 2 .
交通雷达回波信号功率通常包含车辆回波信号功率Pv,背景回波信号功率Pb和白噪声功率Pn。频率点k处的雷达回波信号功率可以表示为:
PR(k)=Pv(k)+Pb(k)+Pn(k)。
白噪声功率符合N(0,σ2)的高斯分布,则白噪声功率Pn(k)=σ2为很小的恒定值,可以忽略。这样,雷达回波信号功率仅由运动着的车辆产生的动态功率和相对静止的公路背景产生的静态功率组成。
所述步骤3进一步包括:
回波信号幅度有序功率序列的数据空间归一化;
归一化有效数据的选择。
所述步骤4具体如下:
在有效幅度有序功率序列的归一化数据空间内进行相干平均,如下式表示:
1 M Σ i = 0 M P Ros ( k , i ) ≈ 1 M Σ i = 0 M P vos ( k , i ) + P b 0 ( k ) ,
式中PRos为幅度有序功率,Pvos为幅度有序功率的车辆回波部分,Pb0为当前公路背景功率谱,M为归一化有效数据容量。在归一化有效数据空间内,车辆产生的动态功率 Σ i = 0 M P v ( k , i ) → 0 , 当前公路背景功率谱如下式表示:
P b 0 ( k ) ≈ 1 M Σ i = 0 M P Ros ( k , i ) .
所述步骤5具体如下:
在每个频点处,实时更新公路背景功率谱的自学习过程如下式表示:
P b t ( k ) = ( 1 - r ) P b t - 1 ( k ) + r P b 0 t ( k ) .
式中Pb为公路背景功率谱,r为学习速率。
所述步骤6具体如下:
在每个频点处,用于完成公路背景抑制的谱相减过程可如下式表示:
Figure A20081004095300115
式中χ为谱相减系数,γ1、γ2为谱相减指数。
所述的谱相减系数χ满足:χ∈[1,2];
所述的谱相减指数满足γ1、γ2:γ1=1/γ2=2。
图1描述了根据本发明的一个具体实施例的背景识别和背景抑制的流程图。首先利用同时包含车辆和背景信息的回波信号,由周期图法求得雷达回波信号功率谱,并对回波信号功率谱时间序列按照功率谱幅度排序。然后应用数据滤波器选择出有效数据,并对有效数据应用相干平均,求得当前公路背景功率谱。然后,应用自学习过程求得最终的公路背景功率谱。最后进行公路背景抑制。
图2描述了根据本发明的一个具体实施例中的交通雷达安装和公路环境示意图。在已有系统样机的基础上,我们在上海市长宁路上进行了实地测试。毫米波交通雷达201安装在公路旁的支架202上,并保证雷达的波速方向垂直于公路方向。毫米波交通雷达201可以同时检测正向行车道204和反向行车道206的交通信息。公路为双向八车道,在第一车道与非机动车道之间存在高为1.5m的金属防护栏203。并在正向行车道207和反向行车道208之间存在防护绿化带205,该绿化带由灌木和草坪组成。在一次测试中,仅正向行车道204的第四车道存在车辆207。
图3描述了根据本发明的一个具体实施例中的公路背景识别结果。功率谱峰303是由第一车道与非机动车道之间的金属防护栏203造成的。由灌木和草坪组成的防护绿化带205导致了功率谱包络305。
图4描述了根据本发明的一个具体实施例中出现的公路背景功率与车辆功率发生功率混叠时的公路回波功率谱。第一车道与非机动车道之间的金属防护栏203造成的功率谱峰303可比拟与正向行车道204第四车道的车辆207造成的功率谱峰407,严重影响了交通雷达的检测性能。由防护绿化带205产生的功率谱包络305与车辆207造成的功率谱峰407处于同一频段内,造成了频率混叠,这同样对交通雷达的检测精度造成了不良影响。
图5描述了根据本发明的一个具体实施例中的公路背景抑制结果。测试中,谱相减系数χ=1.1,谱相减指数γ1=1/γ2=2。

Claims (8)

1、基于毫米波交通雷达的公路背景识别与公路背景抑制方法,其特性在于,包括如下步骤:
步骤1、利用毫米波交通雷达的时域数字回波信号求得雷达回波功率谱时间序列;
步骤2、在每个频点处将雷达回波功率值时间序列按照功率谱幅度升序排列;
步骤3、滤除干扰序列,保留幅度低于有效阀值的有序功率谱序列作为公路背景功率谱的有效数据;
步骤4、对有效幅度有序功率谱序列进行相干平均,得到当前公路背景功率谱;
步骤5、公路背景功率谱通过自学习方法进行实时更新;
步骤6、利用谱相减法对公路背景进行抑制。
2、如权利要求1所述的基于毫米波交通雷达的公路背景识别与公路背景抑制方法,其特性在于,所述步骤1具体如下:
采用N=2n点矩形窗w(n)截取雷达数字回波信号SIF,求得快速傅立叶变换为:
Ssamp(k)=SIF(k)*W(k);
舍去雷达回波信号功率谱中与正频部分对称的负频部分,得到的雷达回波功率幅度为:
P R ( k ) = 1 N | S samp ( k ) | 2 .
3、如权利要求2所述的基于毫米波交通雷达的公路背景识别与公路背景抑制方法,其特征在于:
交通雷达回波信号功率包含车辆回波信号功率Pv,背景回波信号功率Pb和白噪声功率Pn;频率点k处的雷达回波信号功率为:
PR(k)=Pv(k)+Pb(k)+Pn(k);
忽略其中符合高斯分布的微小恒定值Pn(k)。
4、如权利要求1所述的基于毫米波交通雷达的公路背景识别与公路背景抑制方法,其特征在于,所述步骤3进一步包括:
回波信号幅度有序功率序列的数据空间归一化;
归一化有效数据的选择。
5、如权利要求1所述的基于毫米波交通雷达的公路背景识别与公路背景抑制方法,其特征在于,所述步骤4具体如下:
在有效幅度有序功率序列的归一化数据空间内进行相干平均,如下式表示:
1 M Σ i = 0 M P Ros ( k , i ) ≈ 1 M Σ i = 0 M P vos ( k , i ) + P b 0 ( k ) ,
式中PRos为幅度有序功率,Pvos为幅度有序功率的车辆回波部分,Pb0为当前公路背景功率谱,M为归一化有效数据容量;
在归一化有效数据空间内,车辆产生的动态功率 Σ i = 0 M P v ( k , i ) → 0 , 当前公路背景功率谱如下式表示:
P b 0 ( k ) ≈ 1 M Σ i = 0 M P Ros ( k , i ) .
6、如权利要求1所述的基于毫米波交通雷达的公路背景识别与公路背景抑制方法,其特征在于,所述步骤5具体如下:
在每个频点处,实时更新公路背景功率谱的自学习过程如下式表示:
P b t ( k ) = ( 1 - r ) P b t - 1 ( k ) + r P b 0 t ( k ) .
式中Pb为公路背景功率谱,r为学习速率。
7、如权利要求1所述的基于毫米波交通雷达的公路背景识别与公路背景抑制方法,其特征在于,所述步骤6具体如下:
在每个频点处,用于完成公路背景抑制的谱相减过程如下式表示:
Figure A2008100409530004C1
式中χ为谱相减系数,γ1、γ2为谱相减指数。
8、如权利要求7所述的基于毫米波交通雷达的公路背景识别与公路背景抑制方法,其特征在于:
所述的谱相减系数χ满足: χ ⋐ [ 1,2 ]
所述的谱相减指数满足γ1、γ2:γ1=1/γ2=2。
CN2008100409535A 2008-07-24 2008-07-24 基于毫米波交通雷达的公路背景识别与公路背景抑制方法 Expired - Fee Related CN101325007B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN2008100409535A CN101325007B (zh) 2008-07-24 2008-07-24 基于毫米波交通雷达的公路背景识别与公路背景抑制方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN2008100409535A CN101325007B (zh) 2008-07-24 2008-07-24 基于毫米波交通雷达的公路背景识别与公路背景抑制方法

Publications (2)

Publication Number Publication Date
CN101325007A true CN101325007A (zh) 2008-12-17
CN101325007B CN101325007B (zh) 2011-04-06

Family

ID=40188529

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2008100409535A Expired - Fee Related CN101325007B (zh) 2008-07-24 2008-07-24 基于毫米波交通雷达的公路背景识别与公路背景抑制方法

Country Status (1)

Country Link
CN (1) CN101325007B (zh)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105527624A (zh) * 2014-09-29 2016-04-27 北京敏视达雷达有限公司 一种雷达回波动态估算噪声的方法和气象雷达系统
CN108020838A (zh) * 2016-11-02 2018-05-11 惠州市德赛西威汽车电子股份有限公司 一种毫米波雷达信号在自适应巡航中的处理方法
CN110494762A (zh) * 2017-04-10 2019-11-22 株式会社电装 周边监视雷达装置
CN113050056A (zh) * 2021-03-17 2021-06-29 森思泰克河北科技有限公司 一种雷达工作环境识别方法、装置及终端设备

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5798983A (en) * 1997-05-22 1998-08-25 Kuhn; John Patrick Acoustic sensor system for vehicle detection and multi-lane highway monitoring
CN1291238C (zh) * 2002-08-16 2006-12-20 中国科学院上海微系统与信息技术研究所 智能交通中车流量监控的雷达射频装置及测定方法
US7091901B2 (en) * 2004-05-14 2006-08-15 Kustom Signals, Inc. Traffic radar system with improved patrol speed capture
CN1912949A (zh) * 2005-08-12 2007-02-14 上海雷弗电子有限公司 一种车流量检测雷达车道识别方法

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105527624A (zh) * 2014-09-29 2016-04-27 北京敏视达雷达有限公司 一种雷达回波动态估算噪声的方法和气象雷达系统
CN108020838A (zh) * 2016-11-02 2018-05-11 惠州市德赛西威汽车电子股份有限公司 一种毫米波雷达信号在自适应巡航中的处理方法
CN108020838B (zh) * 2016-11-02 2021-08-20 惠州市德赛西威汽车电子股份有限公司 一种毫米波雷达信号在自适应巡航中的处理方法
CN110494762A (zh) * 2017-04-10 2019-11-22 株式会社电装 周边监视雷达装置
CN110494762B (zh) * 2017-04-10 2022-11-11 株式会社电装 周边监视雷达装置
CN113050056A (zh) * 2021-03-17 2021-06-29 森思泰克河北科技有限公司 一种雷达工作环境识别方法、装置及终端设备

Also Published As

Publication number Publication date
CN101325007B (zh) 2011-04-06

Similar Documents

Publication Publication Date Title
CN102915637B (zh) 一种信号控制交叉口的交通状态评价方法及系统
CN101923782B (zh) 基于地磁传感技术的交通车辆速度获取方法
CN103761889B (zh) 一种防追尾碰撞的前方车辆有效目标确定方法
CN103852756B (zh) 一种连续波雷达目标检测跟踪方法
CN110531336A (zh) 一种物体检测识别方法及系统
CN105374212B (zh) 基于智能终端传感的高速公路车辆车道识别方法及系统
CN107591002A (zh) 一种基于分布式光纤的高速公路交通参数实时估计方法
CN102707285A (zh) 车载毫米波防撞雷达系统的频域恒虚警检测方法
Liu et al. The calculation method of road travel time based on license plate recognition technology
CN109324317A (zh) 毫米波雷达系统及其定位测速方法
CN101325007B (zh) 基于毫米波交通雷达的公路背景识别与公路背景抑制方法
CN103473932B (zh) 一种结合振荡标线的音频信号车型识别系统
CN103593649A (zh) 一种用于车道偏离预警的车道线检测方法
Chen et al. Road vehicle detection and classification using magnetic field measurement
CN102608606A (zh) 汽车盲区监控雷达准确检测多目标的有效方法
CN104299417B (zh) 基于波形检测的车辆识别方法
US20210350699A1 (en) Method for Vehicle Classification Using Multiple Geomagnetic Sensors
CN110133632B (zh) 一种基于cwd时频分析的复合调制信号识别方法
CN103680138A (zh) 一种车辆路况监测方法
CN106054169A (zh) 基于跟踪信息的多站雷达信号融合检测方法
CN105092467A (zh) 一种高速铁路车轮踏面擦伤快速监测装置及方法
CN105467428A (zh) 一种基于短时能量检测和频谱特征分析的地震波预警方法
CN108344982A (zh) 基于长时间相参积累的小型无人机目标雷达检测方法
CN104933871B (zh) 一种交通车辆检测方法及系统
CN104503432A (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
C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20110406

Termination date: 20150724

EXPY Termination of patent right or utility model