CN106324577B - 一种基于标准差椭圆的高分辨雷达检测点聚集方法 - Google Patents
一种基于标准差椭圆的高分辨雷达检测点聚集方法 Download PDFInfo
- Publication number
- CN106324577B CN106324577B CN201610685807.2A CN201610685807A CN106324577B CN 106324577 B CN106324577 B CN 106324577B CN 201610685807 A CN201610685807 A CN 201610685807A CN 106324577 B CN106324577 B CN 106324577B
- Authority
- CN
- China
- Prior art keywords
- standard deviation
- ellipse
- target
- detection points
- detection
- 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
- 238000000034 method Methods 0.000 title claims abstract description 26
- 238000001514 detection method Methods 0.000 claims abstract description 86
- 230000002776 aggregation Effects 0.000 claims abstract description 17
- 238000004220 aggregation Methods 0.000 claims abstract description 17
- 238000005070 sampling Methods 0.000 claims description 6
- 230000009466 transformation Effects 0.000 claims description 3
- 230000004931 aggregating effect Effects 0.000 claims description 2
- 238000010586 diagram Methods 0.000 description 8
- 238000004364 calculation method Methods 0.000 description 2
- 238000002592 echocardiography Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 238000009825 accumulation Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000001427 coherent effect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000012732 spatial analysis Methods 0.000 description 1
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
-
- 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/66—Radar-tracking systems; Analogous systems
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-均值的检测点聚集方法需要场景中目标个数的先验信息,并且初始聚类中心的选取对聚集结果存在较大影响,这使复杂场景下目标的检测点聚集较为困难。
发明内容
针对上述现有技术,本发明提供了一种基于标准差椭圆的高分辨雷达检测点聚集方法,能够在复杂场景下能够更准确地聚集检测点,具体步骤如下:
步骤一、建立标准差椭圆:在距离多普勒平面图上,针对包含所有高分辨雷达检测点的集合,建立标准差椭圆;
步骤二、判断集合内的目标个数:计算标准差椭圆内过圆心且平行于距离采样单元轴的线段mn的长度lmn,并将其与目标在雷达径向上的最大分布长度Dr进行比较;如果Dr<lmn,则判断检测点集合里包含多个目标,执行步骤三;否则,判断检测点集合里只包含一个目标,执行步骤四;
步骤三、子集合划分:将横坐标小于或等于圆心横坐标的检测点划分为一个子集合,将横坐标大于圆心横坐标的检测点划分为另一个子集合;然后针对每个子集合,采用步骤一~步骤三的方法,分别建立标准差椭圆模型、目标个数判断和子集合划分,直到每个子集合都只包含一个目标,执行步骤四;
步骤四、聚集检测点:针对每个只包含1个目标的集合或子集合,对该集合或子集合进行聚集,选取集合或子集合内幅值最大的点作为聚集结果或将集合或子集合内的检测点按幅值进行加权平均得到聚集结果。
其中,建立标准差椭圆模型的步骤如下:
步骤1.1:假设在距离多普勒平面目标经过检测后形成数量为n的检测点,这些检测点集合为G:
G={p1(r1,d1),p2(r2,d2),...,pi(ri,di),..,pn(rn,dn)},1≤i≤n (1)
其中,ri和di分别为检测点pi所在的距离采样单元和多普勒通道号;
计算n个检测点的平均中心
步骤1.2:将所有的检测点进行坐标变换为:
步骤1.3:计算标准差椭圆圆心(C1,C2):
步骤1.4:计算标准差椭圆转角θ:
其中,
其中,转角θ为椭圆长轴与距离多普勒平面横轴的夹角;
步骤1.5:计算椭圆的长轴σ1和短轴σ2:
步骤1.6:建立目标检测点椭圆方程:
λ1(r-C1)2+λ2(r-C1)(d-C2)+λ3(d-C2)2+λ4=0 (8)
其中,
优选的,在步骤四中,针对集合中椭圆内的检测点进行聚集。
有益效果:
该方法不需要场景目标的先验信息,也不需要给定初始的聚类中心,在任何复杂场景下,均可根据目标的检测点的分布特征进行点迹聚类,因此对检测点的聚集更加准确;
在强杂波环境下,此算法能进一步减少“干扰杂波点”对目标点迹聚集的影响。
附图说明
图1为标准差椭圆模型示意图。
图2为子算法-建立检测点的标准差椭圆的流程图。
图3为主算法-基于标准差椭圆的检测点聚集方法的流程图。
图4为基于标准差椭圆的检测点聚集方法的工作流程示意图。
图4(a)为检测点集合G示意图。
图4(b)为标准差椭圆示意图。
图4(c)为将集合G划分为两个子集合G1和G2示意图。
图4(d)为将检测点分为左右两个部分G1和G2的结果示意图。
图4(e)为两个子集合的标准差椭圆示意图。
图4(f)和图4(g)分别为剔除前和剔除后的效果图。
图4(h)为两个目标检测点区域G1'和G2'的示意图。
具体实施方式
本发明提供了一种基于标准差椭圆的高分辨雷达检测点聚集方法。目标检测点的分布特征能够展现目标的空间结构信息,而目标的空间结构不随时间变化,并且不同目标的结构特征差异较大,因此,本发明利用目标的检测点分布特征进行检测点聚集。在空间分析领域,标准差椭圆是点数据分布特征分析的一个重要工具,利用标准差椭圆可识别一组数据的方向及分布趋势,并得到该组数据的一些特性。基于此,提出一种基于标准差椭圆的检测点聚集方法,根据目标检测点的分布特征建立椭圆模型,并利用椭圆的特征参数聚集不同目标的检测点,判断目标个数并排除干扰点。
下面结合附图对本发明的实施方式作详细说明。具体流程如图3所示,具体包括如下步骤:
步骤1、建立检测点的标准差椭圆
高分辨雷达目标回波经相参积累和检测后形成检测点通常具有一定的方向偏离性并且检测点具有高斯分布的特点,所以采用标准差椭圆模型来描述检测点的分布。
假设目标经过检测后在距离多普勒平面上(以距离采样单元为横坐标,多普勒通道信号为纵坐标的构成的平面)形成数量为n的检测点,这些检测点集合为
G={p1(r1,d1),p2(r2,d2),...,pi(ri,di),..,pn(rn,dn)},1≤i≤n (10)
其中,ri和di分别为检测点pi所在的距离采样单元和多普勒通道号。采用如下步骤建立建立检测点标准差椭圆模型,如图1所示:
步骤1.1:计算n个检测点平均中心
步骤1.2:将所有的检测点进行坐标变换为(ri′,di′):
步骤1.3:计算标准差椭圆圆心(C1,C2):
步骤1.4:计算标准差椭圆转角θ:
其中,
其中,转角θ为椭圆长轴与r轴的夹角,r轴为距离多普勒平面的横轴,如图1所示。
步骤1.5:计算椭圆的长轴σ1和短轴σ2:
步骤1.6:建立目标检测点椭圆方程:
λ1(r-C1)2+λ2(r-C1)(d-C2)+λ3(d-C2)2+λ4=0 (17)
其中,
λ1=σ1 2cos2θ+σ2 2sin2θ
λ2=2(σ1 2-σ2 2)sinθcosθ
λ3=σ1 2sin2θ+σ2 2cos2θ
λ4=-σ1 2σ2 2
(18)
步骤1.7:计算椭圆内过圆心且平行于距离采样单元轴的线段mn的长度lmn:
步骤1.8:得标准差椭圆的特征参数集合T:
T={C1,C2,λ1,λ2,λ3,λ4,lmn} (20)
步骤2、判断检测点集合内含目标个数
针对图4(a)的检测点集合G,经计算后形成如图4(b)所示的标准差椭圆,因单个目标所占距离单元数有限,所以可以通过线段mn的长度来判断当前集合内的检测点是否属于同一个目标,实现检测点的聚集,具体流程:
判断线段mn长度是否满足如下条件:
lmn≤Dr (21)
其中,Dr为目标在雷达径向上的最大分布长度;如果不满足,则判定集合G内包含多个目标的检测点,执行步骤3;否则判定集合G内只包含一个目标的检测点,执行步骤4。
步骤3、划分检测点集合为两个子集合
如果集合G内包含多个目标的检测点,则通过以下判断条件
ri≤C1 (22)
将集合G划分为两个子集合,即以图4(c)中的虚线将检测点分为左右两个部分G1和G2,结果如图4(d);然后针对每个子集合重复步骤1~2,计算两个子集合的标准差椭圆,结果如图4(e)所示,直到每个子集合的检测点建立的标准差椭圆参数lmn都满足式(12);
步骤4、排除干扰检测点
由于在复杂情况下,干扰回波经检测后形成干扰检测点,干扰检测点一般分布在目标检测点的周围,针对只包含一个目标的检测点集合或子集,通过以下判断条件排除所述标准差椭圆外的干扰点,即认为不满足式(14)的检测点为干扰点
λ1(ri-C1)2+λ2(ri-C1)(di-C2)+λ3(di-C2)2+λ4<0 (23)
图4(f)和图4(g)分别代表剔除前和剔除后的效果图。最后形成两个目标检测点区域G1'和G2',如图4(h)所示。通过此方法能够减少干扰检测点对后续量测计算的影响。
步骤5、聚集检测点
椭圆内的多个检测点均为同一目标的检测结果,选取标准差椭圆内幅值最大的点作为聚集结果或将椭圆圆内的检测点按幅值进行加权平均得到聚集结果,以进行后续的跟踪处理。
基于标准差椭圆的检测点聚集方法利用了检测点分布特征,建立检测点的标准差椭圆模型,通过判断椭圆的特征参数将不同目标的检测点聚集起来,同时进一步减少了干扰噪声对目标量测的影响。
综上所述,以上仅为本发明的较佳实施例而已,并非用于限定本发明的保护范围。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (3)
1.一种基于标准差椭圆的高分辨雷达检测点聚集方法,其特征在于,包括以下步骤:
步骤一、建立标准差椭圆:在距离多普勒平面图上,针对包含所有高分辨雷达检测点的集合,建立标准差椭圆;
步骤二、判断集合内的目标个数:计算标准差椭圆内过圆心且平行于距离采样单元轴的线段mn的长度lmn,并将其与目标在雷达径向上的最大分布长度Dr进行比较;如果Dr<lmn,则判断检测点集合里包含多个目标,执行步骤三;否则,判断检测点集合里只包含一个目标,执行步骤四;
步骤三、子集合划分:将横坐标小于或等于圆心横坐标的检测点划分为一个子集合,将横坐标大于圆心横坐标的检测点划分为另一个子集合;然后针对每个子集合,采用步骤一~步骤三的方法,分别建立标准差椭圆模型、目标个数判断和子集合划分,直到每个子集合都只包含一个目标,执行步骤四;
步骤四、聚集检测点:针对每个只包含1个目标的集合或子集合,对该集合或子集合进行聚集,选取集合或子集合内幅值最大的点作为聚集结果或将集合或子集合内的检测点按幅值进行加权平均得到聚集结果。
2.如权利要求1所述的基于标准差椭圆的高分辨雷达检测点聚集方法,其特征在于,采用如下步骤建立标准差椭圆模型:
步骤1.1:假设在距离多普勒平面目标经过检测后形成数量为n的检测点,这些检测点集合为G:
G={p1(r1,d1),p2(r2,d2),...,pi(ri,di),..,pn(rn,dn)},1≤i≤n (1)
其中,ri和di分别为检测点pi所在的距离采样单元和多普勒通道号;
计算n个检测点的平均中心
步骤1.2:将所有的检测点进行坐标变换为:
步骤1.3:计算标准差椭圆圆心(C1,C2):
步骤1.4:计算标准差椭圆转角θ:
其中,
其中,转角θ为椭圆长轴与距离多普勒平面横轴的夹角;
步骤1.5:计算椭圆的长轴σ1和短轴σ2:
步骤1.6:建立目标检测点椭圆方程:
λ1(r-C1)2+λ2(r-C1)(d-C2)+λ3(d-C2)2+λ4=0 (8)
其中,
3.如权利要求1所述的基于标准差椭圆的高分辨雷达检测点聚集方法,其特征在于,在步骤四中,针对集合中椭圆内的检测点进行聚集。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610685807.2A CN106324577B (zh) | 2016-08-18 | 2016-08-18 | 一种基于标准差椭圆的高分辨雷达检测点聚集方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610685807.2A CN106324577B (zh) | 2016-08-18 | 2016-08-18 | 一种基于标准差椭圆的高分辨雷达检测点聚集方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106324577A CN106324577A (zh) | 2017-01-11 |
CN106324577B true CN106324577B (zh) | 2019-03-12 |
Family
ID=57743403
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610685807.2A Active CN106324577B (zh) | 2016-08-18 | 2016-08-18 | 一种基于标准差椭圆的高分辨雷达检测点聚集方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106324577B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP3518001B1 (en) * | 2018-01-25 | 2020-09-16 | Aptiv Technologies Limited | Method for increasing the reliability of determining the position of a vehicle on the basis of a plurality of detection points |
EP3546983B1 (en) * | 2018-03-28 | 2021-12-08 | Aptiv Technologies Limited | Method for identifying objects in a traffic space |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101894269A (zh) * | 2010-07-16 | 2010-11-24 | 西安电子科技大学 | 基于多分类器系统的合成孔径雷达自动目标识别方法 |
CN101963662A (zh) * | 2010-09-20 | 2011-02-02 | 北京理工大学 | 基于短时分数阶傅里叶域滤波的自聚焦预处理方法 |
CN105372656A (zh) * | 2015-08-31 | 2016-03-02 | 首都师范大学 | 一种表述地面沉降空间演化趋势的方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20070018890A1 (en) * | 2005-07-22 | 2007-01-25 | Kulyukin Vladimir A | Multi-sensor wayfinding device |
-
2016
- 2016-08-18 CN CN201610685807.2A patent/CN106324577B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101894269A (zh) * | 2010-07-16 | 2010-11-24 | 西安电子科技大学 | 基于多分类器系统的合成孔径雷达自动目标识别方法 |
CN101963662A (zh) * | 2010-09-20 | 2011-02-02 | 北京理工大学 | 基于短时分数阶傅里叶域滤波的自聚焦预处理方法 |
CN105372656A (zh) * | 2015-08-31 | 2016-03-02 | 首都师范大学 | 一种表述地面沉降空间演化趋势的方法 |
Non-Patent Citations (2)
Title |
---|
前向散射雷达地表杂波物理建模及频谱扩展分析;胡程 等;《中国科学: 信息科学》;20101231;第40卷(第12期);第1646-1659页 |
姜元 等.基于CLEAN的宽带多散射点目标微动提取.《计算机工程与应用》.2015,第51卷(第S1期), |
Also Published As
Publication number | Publication date |
---|---|
CN106324577A (zh) | 2017-01-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110187318B (zh) | 一种雷达数据处理方法 | |
CN108520023B (zh) | 一种基于混合聚类算法的雷暴核识别及追踪方法 | |
CN104635203B (zh) | 一种基于粒子滤波算法的无线电干扰源测向定位方法 | |
CN106680806A (zh) | 一种多雷达点迹融合方法 | |
CN108168564A (zh) | 一种基于lhd灰色关联度的航迹关联方法 | |
CN108181620B (zh) | 一种三坐标雷达点迹质量评估方法 | |
CN112051568B (zh) | 一种两坐标雷达的俯仰测角方法 | |
JP6773606B2 (ja) | レーダ装置及びそのレーダ信号処理方法 | |
CN102338874A (zh) | 用于被动式多传感器目标跟踪的全局概率数据关联方法 | |
CN106526549A (zh) | 联合两坐标雷达和三坐标雷达统计量的假目标鉴别方法 | |
CN106680783B (zh) | 基于站址误差融合算法对抗欺骗式假目标的方法 | |
CN106324577B (zh) | 一种基于标准差椭圆的高分辨雷达检测点聚集方法 | |
CN103728615B (zh) | 相控阵二次雷达多目标检测方法及系统 | |
CN111505649B (zh) | 一种拖曳式被动阵列声呐低信噪比舰船运动目标检测方法 | |
CN103675808A (zh) | 一种单脉冲雷达导引头的不可分辨多目标检测方法 | |
CN103809161B (zh) | 雷达网抗距离欺骗+soj复合干扰方法 | |
CN110988856B (zh) | 一种基于密度聚类的目标检测点迹凝聚算法 | |
Pengcheng et al. | Radar group target recognition based on HRRPs and weighted mean shift clustering | |
CN110261828B (zh) | 基于距离-角误差二维聚类的横向干扰判定方法 | |
CN111580105B (zh) | 一种用于太赫兹雷达高分辨成像的自适应处理方法 | |
CN105652256B (zh) | 一种基于极化信息的高频地波雷达tbd方法 | |
CN106093891A (zh) | 基于多普勒速度检验的雷达网抗密集距离假目标干扰方法 | |
CN116381611A (zh) | 一种基于区间扩展合并的雷达点迹凝聚方法 | |
CN107728121B (zh) | 一种基于可变窗的局部拟合优度检验方法 | |
CN106443618B (zh) | 一种基于标准差椭圆参数辅助的高分辨雷达航迹关联方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CB03 | Change of inventor or designer information | ||
CB03 | Change of inventor or designer information |
Inventor after: Liu Haibo Inventor after: Wang E Inventor after: Zhou Chao Inventor after: Long Teng Inventor after: Wang Yanhua Inventor after: Wang Junfu Inventor before: Liu Haibo Inventor before: Wang E Inventor before: Zhou Chao Inventor before: Long Teng |