CN104036146B - 一种用于雷达目标点迹凝聚的点迹聚类方法 - Google Patents
一种用于雷达目标点迹凝聚的点迹聚类方法 Download PDFInfo
- Publication number
- CN104036146B CN104036146B CN201410295455.0A CN201410295455A CN104036146B CN 104036146 B CN104036146 B CN 104036146B CN 201410295455 A CN201410295455 A CN 201410295455A CN 104036146 B CN104036146 B CN 104036146B
- Authority
- CN
- China
- Prior art keywords
- point
- mark
- video signal
- test point
- radar
- 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
Landscapes
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开了一种用于雷达目标点迹凝聚的点迹聚类方法,可应用于低分辨率脉冲搜索雷达的目标检测录取,属于雷达信号处理领域。本发明对雷达检测视频进行操作,充分考虑了雷达目标回波形态特征,以点迹为基本单位,遵循“若点迹属于已知目标,则下一相邻方位相同距离的点迹与之属于同一目标”这一基本原则,根据测试点周边点迹的分布来处理新目标、同一目标点迹断裂、特殊回波形状等一系列情况,判别测试点是否属于已知目标。本发明降低了传统点迹聚类时的算法复杂度,同时可以将粘连的两个目标回波辨别为两个点迹,提高了处理速度和精度。
Description
技术领域
本发明涉及雷达信号处理领域,特别是一种用于低分辨率脉冲搜索雷达点迹凝聚的点迹聚类方法。
背景技术
雷达原始视频信号分方位和距离两个维度。低分辨率脉冲搜索雷达扫掠过目标时,天线将接收到一串回波。典型的目标回波如图1所示,在方位上有一定展宽,在距离上有一定延伸。量化采样后的雷达原始视频信号进行检测,大于检测阈值保留的信号单元称为点迹,小于检测阈值清零的信号单元称为空点。检测后的雷达视频信号称为检测视频。理论上,检测视频去除了噪声和干扰,只保留了目标的点迹。但实际上可能依然存在无用的孤立点。进一步的,同一目标点迹经凝聚后发送给雷达终端用于跟踪、录取。
点迹凝聚分为聚类和凝聚两个阶段。聚类是指逐方位逐距离的读入雷达检测视频,判断测试点是点迹,还是空点或是孤立点;若测试点是点迹,则继续判断当前点迹属于哪个目标,标记聚类状态。所有雷达视频信号处理结束时,可根据聚类状态找到目标所包含的所有点迹,然后根据下面两个公式完成凝聚:
公式(1)和(2)中,A是回波幅值,θ与r分别是点迹的方位与距离,与是点迹凝聚后的方位与距离,i与j分别是点迹在方位与距离上的序号。
传统上,点迹聚类分为两步:第一步,统计同一方位上的连续点迹,归纳为线段;第二步,计算方位间相邻线段的相关性,根据设定的相关性阈值,判断两条线段是否属于同一目标,从而完成点迹聚类。
传统点迹聚类存在这样的问题:对于如图2所示的目标回波粘连,会将同一方位上多个目标的点迹统计为一条线段,如图3所示。从而在方位间相邻线段相关性计算时,将两个目标的点迹聚类为同一个目标的点迹。凝聚效果如图4所示,回波信号中央的加号即为凝聚结果,显然,两个相邻的目标回波只有一个加号,被凝聚为同一个目标。
传统点迹聚类的基本单位是线段。如果将基本单位改为点迹,则可以将传统点迹聚类方法归纳为两条特征:
特征1:距离上相邻的点迹属于同一目标;
特征2:方位上相邻的点迹属于同一目标;
发明内容
发明目的:本发明目的在于提供一种更有效的用于雷达目标回波点迹凝聚的点迹聚类方法,克服传统点迹聚类过程中无法区分目标粘连的缺陷。
为了解决上述技术问题,本发明公开了一种用于雷达目标回波点迹凝聚的点迹聚类方法基于雷达目标回波形态特征,以点迹为基本单位,遵循“若点迹属于已知目标,则下一相邻方位相同距离的点迹与之属于同一目标”这一基本原则,根据测试点周边点迹的分布来处理新目标、同一目标点迹断裂、特殊回波形状等一系列情况,完成同一目标的点迹聚类。具体步骤为:
步骤1:逐个方位、逐个距离的顺序读入R×W个雷达检测视频信号处理,R是雷达检测视频信号的方位总数,W是距离总数。定位当前正在处理的视频信号A(i,j)为测试点,i与j分别是视频信号在方位与距离上的序号;
步骤2:若测试点A(i,j)不是点迹,沿当前方位i逐点顺序查询视频信号A(i,j-1)、A(i,j-2)…A(i,j-L),L是小于10的自然数;
若L点视频信号均不是点迹,则测试点是空点,转入步骤4;
若顺序查询L点视频信号时,有视频信号是点迹且已聚类给某目标,且A(i+1,j)也是点迹,则测试点与查询到的点迹应聚类给同一目标,不执行瀑布规程,转入步骤4;
若顺序查询L点视频信号时,有视频信号是点迹且已聚类给某目标,但A(i+1,j)不是点迹,则跳过测试点,转到步骤4;
步骤3:若测试点是点迹,则判断其是否完成聚类:
若测试点已完成聚类,转到步骤4;
若测试点未完成聚类,转到步骤3-1;
步骤3-1:查询视频信号A(i,j-1)的处理状态;
若A(i,j-1)已完成聚类,则测试点与A(i,j-1)应聚类给同一目标;进一步的,若
A(i,j-1)也是点迹,则执行瀑布规程,否则不执行瀑布规程;转到步骤4;
若A(i,j-1)未完成聚类,则转到步骤3-1-1:
步骤3-1-1:顺序查询视频信号A(i,j+1)、A(i,j+2)…A(i,j+M)点的处理状态,M是小于10的自然数:
若M点范围内有点迹已完成聚类,则测试点与该点迹应聚类给同一目标,执行瀑布规程,转到步骤4;
若M点范围内没有点迹完成聚类,且测试点符合鲁棒条件,则应聚类给
新目标并执行瀑布规程;否则,测试点属于孤立点,不做处理转入步骤4;
步骤4:判断是否已处理完R×W个雷达检测视频信号的处理:
若是,则算法结束;
若否,则继续读入新的视频信号作为测试点,返回步骤2。
前述步骤中,鲁棒条件是指当前测试点A(i,j)沿同一方位A(i,j+1)、A(i,j+1)…A(i,j+N)都是点迹。瀑布规程是指当前测试点A(i,j)与同一距离下一方位的点迹A(i+1,j)聚类为同一目标。
本方法在传统聚类方法基础上,以点迹为基本单位对检测视频进行操作,根据测试点周边空点和点迹的分布条件,判别其是否完成聚类。点迹聚类的基本原则是:若点迹已完成聚类,则其下方的点迹应聚类给同一目标。具体操作过程中,需对新点、同一目标点迹断裂、特殊回波形状等一系列情况进行更加细化的处理。其显著优点是:(1)充分利用了雷达目标回波形状特点(2)处理时基本单位更小,处理更为精细(3)可依据雷达分辨率自由设置L、M、N,从而使得算法鲁棒性更强。(4)相较于传统算法,可区分头部分裂,腹部粘连的雷达目标;(5)算法便于在硬件实现加速。(6)算法中的约束条件便于升级优化。
附图说明
下面结合附图和具体实施方式对本发明做更进一步的具体说明,本发明的上述和/或其他方面的优点将会变得更加清楚。
图1典型目标回波形态。
图2典型目标回波粘连。
图3传统点迹聚类缺陷。
图4是传统方法实施效果图。
图5是实施例用于雷达目标回波点迹凝聚的点迹聚类方法的流程图。
图6是实施例处理目标头部部分断裂的示意图。
图7是实施例处理目标粘连的示意图。
图8是实施例实施中“左侧临近”原则示意图。
图9是实施例处理右侧首先聚类时的目标形状示意图。
图10是实施例实施效果图。
具体实施方式
本发明公开了一种用于雷达目标点迹凝聚的点迹聚类方法,可应用于低分辨率脉冲搜索雷达的目标检测录取,属于雷达信号处理领域。本发明对雷达检测视频进行操作,充分考虑了雷达目标回波形态特征,以点迹为基本单位,遵循“若点迹属于已知目标,则下一相邻方位相同距离的点迹与之属于同一目标”这一基本原则,根据测试点周边点迹的分布来处理新目标、同一目标点迹断裂、特殊回波形状等一系列情况,判别测试点是否属于已知目标。本发明降低了传统点迹聚类时的算法复杂度,同时可以将粘连的两个目标回波辨别为两个点迹,提高了处理速度和精度。
本发明可基于目标回波的特性,将点迹视为基本单位,对聚类特征进行扩充,将聚类方法细化,实现目标点迹聚类,从而部分避免传统点迹聚类过程中存在的缺陷,并提高处理效率。
实施例
如图5所示,本实施例以R×W尺寸的检测视频为例进行处理。除了目标点迹和空点,还可能存在无用的孤立点。处理时,逐方位、逐距离的读入检测视频A(i,j),定位其为测试点,具体实施步骤是:
对于读入的测试点A(i,j),首先判断其是点迹还是空点。
若当前测试点A(i,j)是空点,则可能是因检测造成的断裂,如图6所示,空白框的是空点,反斜线框是点迹。如不处理,则同一目标可能处理为2个目标。对此情形处理:沿当前方位i逐点顺序查找A(i,j-1)、A(i,j-2)…A(i,j-L)是否存在点迹。若不存在点迹,则A(i,j)不是因检测造成的目标断裂,而是属于连续的空点,跳过当前测试点处理下一点。
若L点范围内查到点迹且完成聚类,A(i+1,j)也是点迹,则认为当前空点是图6所示情况,空点与所查到点迹聚类给同一目标,不执行瀑布规程;若L点范围内存在点迹完成聚类,但A(i+1,j)不是点迹,此时认为对应图7所示的目标粘连情况,跳过当前测试点,开始处理下一点。
若当前读入的测试点A(i,j)是点迹且完成聚类,那么该点必然由A(i-1,j)执行瀑布规程完成了目标属性传递,直接跳过处理下一点。
若当前读入的测试点A(i,j)是点迹但未完成聚类,此时该目标处于头部或内部,如图8所示,交叉线框是新目标的第一个点迹,后续点迹只能通过第一个点迹传递。对于目标断裂下方的点迹也存在类似问题。此时,遵循的原则是“左侧临近”,即判断A(i,j-1)是否为点迹:
(1)若A(i,j-1)是点迹且完成聚类,则测试点A(i,j)与A(i,j-1)聚类给同一目标,执行瀑布规程;
(2)若A(i,j-1)不是点迹但完成聚类,则测试点A(i,j)与A(i,j-1)聚类给同一目标,不执行瀑布规程。
(3)若A(i,j-1)未完成聚类,此时既可能是新目标,也可能属于图9所示情况:目标自右侧开始聚类,但左侧依然没有开始聚类。对此进行判断,沿当前方位i逐点顺序查找A(i,j+1)、A(i,j+2)…A(i,j+M)是否存在点迹:(3.a)若M点范围内有点迹完成聚类,则当前测试点A(i,j)与该点迹应聚类给同一目标,执行瀑布规程,转到步骤4。(3.b)若M点内没有点迹完成聚类,当前测试点符合鲁棒条件,则该点属于新目标,执行瀑布规程;当前测试点不符合鲁棒条件,则该点属于孤立点,不做任何操作,读入后1点处理。
前述步骤中,鲁棒条件是指当前测试点A(i,j)沿同一方位A(i,j+1)、A(i,j+1)…A(i,j+N)都是点迹。瀑布规程是指当前测试点A(i,j)与同一距离下一方位的点迹A(i+1,j)应聚类给同一目标。至此,对于检测视频中一个点的聚类处理结束,可以读入下一点处理。直至所有的R×W视频信号处理结束。
上述过程中,虽然L、M、N都是变量,但对于低分辨率雷达,其在距离上的信号点较少(一般小于10个),而且目标回波呈“眉毛”状,两头较少,中间最大。因此,L、M、N值设置较小。进行点迹聚类时,虽然特征条件多,但可以采用滑窗的形式一次性读入所有有关点迹进行逻辑判断,一次性完成单个点迹的处理,便于硬件实现。进一步的,算法复杂度为ο(XY),其中X为方位单元数,Y为距离单元数。原有的算法则是首先要进行线段统计,再进行线段相关计算,其复杂度至少为本实施例的2倍。
本实施例通过以点迹为基本单位,通过相邻方位同一距离上两个点迹的聚类,使得两个粘连目标能分别聚类,提高了原有算法的性能。效果图如图10所示。
本发明提供了一种用于雷达目标点迹凝聚的点迹聚类方法,具体实现该技术方案的方法和途径很多,以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也应视为本发明的保护范围。本实施例中未明确的各组成部分均可用现有技术加以实现。
Claims (2)
1.一种用于雷达目标点迹凝聚的点迹聚类方法,其特征在于,基于雷达目标回波的形态,判断雷达检测视频信号是否为点迹,以点迹为基本单位进行聚类,包括以下步骤:
步骤1:逐个方位、逐个距离的顺序读入R×W个雷达检测视频信号处理,R是雷达检测视频信号的方位总数,W是距离总数;定位当前正在处理的视频信号A(i,j)为测试点,i与j分别是视频信号在方位与距离上的序号;
步骤2:若测试点A(i,j)不是点迹,沿当前方位i逐点顺序查询L点视频信号A(i,j-1)、A(i,j-2)…A(i,j-L),L是小于10的自然数;
若L点视频信号均不是点迹,则判定测试点是空点,转入步骤4;
若顺序查询L点视频信号时,有视频信号是点迹且已聚类给一个目标,且视频信号A(i+1,j)也是点迹,则判定测试点与查询到的点迹应聚类给同一目标,不执行瀑布规程,转入步骤4;
若顺序查询L点视频信号时,有视频信号是点迹且已聚类给某目标,但视频信号A(i+1,j)不是点迹,则跳过测试点,转到步骤4;
步骤3:若测试点A(i,j)是点迹,则判断其是否完成聚类:
若测试点A(i,j)已完成聚类,转到步骤4;
若测试点A(i,j)未完成聚类,转到步骤3-1;
步骤3-1:查询视频信号A(i,j-1)的处理状态;
若视频信号A(i,j-1)已完成聚类,则判定测试点A(i,j)与A(i,j-1)应聚类给同一目标;若A(i,j-1)也是点迹,则执行瀑布规程,否则不执行瀑布规程;转到步骤4;
若视频信号A(i,j-1)未完成聚类,则转到步骤3-1-1:
步骤3-1-1:顺序查询M点视频信号A(i,j+1)、A(i,j+2)…A(i,j+M)点的处理状态,M是小于10的自然数:
若M点范围内有点迹已完成聚类,则判定测试点A(i,j)与该点迹应聚类给同一目标,执行瀑布规程,转到步骤4;
若M点范围内没有点迹完成聚类,且测试点A(i,j)符合鲁棒条件,则判定测试点A(i,j)应聚类给新目标并执行瀑布规程;否则,如果测试点A(i,j)不符合鲁棒条件,判定测试点A(i,j)属于孤立点,转入步骤4;
步骤4:判断是否已处理完R×W个雷达检测视频信号的处理:
若是,则结束;
若否,则继续读入下一个的视频信号作为测试点,返回步骤2;
所述瀑布规程是指:将当前测试点A(i,j)与同一距离下一方位的视频信号A(i+1,j)应聚类给同一目标。
2.根据权利要求1所述的用于雷达目标点迹凝聚的点迹聚类方法,其特征在于,步骤3-3-1中的鲁棒条件是指:当前测试点A(i,j)沿同一方位的N点视频信号A(i,j+1)、A(i,j+1)…A(i,j+N)都是点迹,N是小于10的自然数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410295455.0A CN104036146B (zh) | 2014-06-26 | 2014-06-26 | 一种用于雷达目标点迹凝聚的点迹聚类方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410295455.0A CN104036146B (zh) | 2014-06-26 | 2014-06-26 | 一种用于雷达目标点迹凝聚的点迹聚类方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104036146A CN104036146A (zh) | 2014-09-10 |
CN104036146B true CN104036146B (zh) | 2017-04-26 |
Family
ID=51466916
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410295455.0A Active CN104036146B (zh) | 2014-06-26 | 2014-06-26 | 一种用于雷达目标点迹凝聚的点迹聚类方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104036146B (zh) |
Families Citing this family (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104459664B (zh) * | 2014-11-28 | 2017-02-22 | 大连海大船舶导航国家工程研究中心有限责任公司 | 一种基于点迹合并的雷达标绘目标识别方法 |
CN104569950B (zh) * | 2015-01-27 | 2017-04-19 | 西安电子科技大学 | 基于相对功率和相位线性度的海面目标凝聚方法 |
CN104931934B (zh) * | 2015-06-25 | 2017-07-14 | 中国船舶重工集团公司第七二四研究所 | 一种基于pam聚类分析的雷达点迹凝聚方法 |
CN105116394A (zh) * | 2015-06-29 | 2015-12-02 | 中国船舶重工集团公司第七二四研究所 | 一种基于多维细微特征分析的目标点迹提取方法 |
CN106918807B (zh) * | 2017-02-28 | 2019-07-09 | 西安电子科技大学 | 一种雷达回波数据的目标点迹凝聚方法 |
CN107703503B (zh) * | 2017-09-22 | 2020-10-23 | 电子科技大学 | 一种基于gpu加速的点迹凝聚方法 |
CN109100696B (zh) * | 2018-07-06 | 2022-07-01 | 西安电子科技大学 | 基于点状杂波图的慢速运动目标消除方法 |
CN110542897B (zh) * | 2019-08-01 | 2021-08-13 | 北京理工大学 | 一种基于希尔排序的距离差分多扩展目标点迹凝聚方法 |
CN110988856B (zh) * | 2019-12-19 | 2021-08-03 | 电子科技大学 | 一种基于密度聚类的目标检测点迹凝聚算法 |
CN111366924B (zh) * | 2020-05-09 | 2022-11-25 | 苏州雷森电子科技有限公司 | 一种基于多维连通域处理的雷达目标点迹凝聚方法 |
CN112731296B (zh) * | 2020-11-30 | 2024-03-08 | 惠州市德赛西威智能交通技术研究院有限公司 | 一种汽车毫米波雷达点迹凝聚方法及系统 |
CN113253228B (zh) * | 2021-05-11 | 2024-07-09 | 成都西科微波通讯有限公司 | 一种基于雷达凝聚算法的多尺寸目标检测方法 |
CN113406591B (zh) * | 2021-06-15 | 2022-08-12 | 成都西科微波通讯有限公司 | 一种应用于场面监视雷达的目标凝聚方法 |
CN115166677B (zh) * | 2022-07-27 | 2024-09-24 | 中国船舶集团有限公司第七二四研究所 | 一种基于区域优先的雷达目标点迹智能提取方法 |
CN118444279B (zh) * | 2024-07-08 | 2024-09-13 | 南京奥联智驾科技有限公司 | 一种雷达回波信号快速点迹凝聚方法、设备及介质 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103852756A (zh) * | 2012-11-30 | 2014-06-11 | 中国科学院沈阳自动化研究所 | 一种连续波雷达目标检测跟踪方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR100973039B1 (ko) * | 2009-12-07 | 2010-07-29 | 엘아이지넥스원 주식회사 | 레이더 펄스의 클러스터링 방법 |
-
2014
- 2014-06-26 CN CN201410295455.0A patent/CN104036146B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103852756A (zh) * | 2012-11-30 | 2014-06-11 | 中国科学院沈阳自动化研究所 | 一种连续波雷达目标检测跟踪方法 |
Non-Patent Citations (1)
Title |
---|
雷达点迹凝聚处理技术及其数据分析;杨文琳;《信号处理》;20010430;第17卷(第2期);第130-138页 * |
Also Published As
Publication number | Publication date |
---|---|
CN104036146A (zh) | 2014-09-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104036146B (zh) | 一种用于雷达目标点迹凝聚的点迹聚类方法 | |
CN112505648B (zh) | 基于毫米波雷达回波的目标特征提取方法 | |
WO2020078249A1 (zh) | 一种基于雷达信号的多目标生命探测方法及探测雷达 | |
CN103064074B (zh) | 一种脉冲多普勒雷达在强杂波背景下的弱目标检测方法 | |
CN108490410B (zh) | 一种两坐标雷达对海目标联合检测跟踪方法 | |
CN106918807A (zh) | 一种雷达回波数据的目标点迹凝聚方法 | |
CN106896407B (zh) | 一种基于近似负熵的微地震信号初至拾取方法 | |
CN104166135A (zh) | 一种宽带雷达目标的原始点迹凝聚处理方法 | |
CN104714225B (zh) | 一种基于广义似然比的动态规划检测前跟踪方法 | |
CN109188430B (zh) | 一种基于地面监视雷达系统的目标提取方法 | |
CN109919862B (zh) | 雷达图像去噪系统、方法及计算机设备 | |
CN106597365A (zh) | 一种基于时域聚类的复杂电子信号时差定位方法 | |
CN105487061A (zh) | 目标数据关联的多特征信息融合方法 | |
CN106019254A (zh) | 一种uwb冲击生物雷达多人体目标距离向分离辨识方法 | |
CN104407346A (zh) | 一种基于信息融合的移动式跑道异物监测方法 | |
CN107146238A (zh) | 基于特征块优选的运动目标跟踪方法 | |
CN103149552A (zh) | 一种多普勒天气雷达径向速度场的退模糊方法 | |
CN109839619A (zh) | 基于自适应分桶的雷达信号粗分选方法、系统及存储介质 | |
CN109100697B (zh) | 一种基于地面监视雷达系统的目标凝聚方法 | |
CN102914302A (zh) | 一种无人机着陆视觉导航合作目标鲁棒检测方法 | |
CN102567726B (zh) | 极地冰盖边缘区域浮冰自动提取技术 | |
CN105445761A (zh) | 基于多普勒频移的抗欺骗式干扰方法 | |
CN110031854A (zh) | 一种实时的高精度激光多回波距离提取方法 | |
CN109389612A (zh) | 一种城轨列车受电弓滑板边缘检测方法 | |
CN104199009B (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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |