CN104931934B - 一种基于pam聚类分析的雷达点迹凝聚方法 - Google Patents

一种基于pam聚类分析的雷达点迹凝聚方法 Download PDF

Info

Publication number
CN104931934B
CN104931934B CN201510360108.6A CN201510360108A CN104931934B CN 104931934 B CN104931934 B CN 104931934B CN 201510360108 A CN201510360108 A CN 201510360108A CN 104931934 B CN104931934 B CN 104931934B
Authority
CN
China
Prior art keywords
cluster
pam
radar
plot
mark
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
Application number
CN201510360108.6A
Other languages
English (en)
Other versions
CN104931934A (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.)
724th Research Institute of CSIC
Original Assignee
724th Research Institute of CSIC
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 724th Research Institute of CSIC filed Critical 724th Research Institute of CSIC
Priority to CN201510360108.6A priority Critical patent/CN104931934B/zh
Publication of CN104931934A publication Critical patent/CN104931934A/zh
Application granted granted Critical
Publication of CN104931934B publication Critical patent/CN104931934B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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

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

本发明涉及一种基于PAM聚类分析的雷达点迹凝聚方法。本发明主要适用于杂波较多、群目标航迹交叉等复杂环境下的雷达点迹提取处理。本发明的处理流程是：首先按方位划分扇区，对每个扇区内的点迹按特征相似度采用PAM聚类算法进行分类，每个划分聚类簇内的点迹认为是同一目标的分裂点迹，然后采用特征贴近度加权计算方法进行目标点迹凝聚处理，并对每个聚类结果建立时间序列模型，进行目标延续判决，判决结果和软件杂波图、跟踪波门、录取波门等信息作为下一扫描周期聚类个数计算的主要依据，迭代执行以上步骤。本发明利用了雷达探测回波的统计特性进行聚类和点迹凝聚处理，有效的提高了雷达点迹提取的准确度和可靠性。

Description

一种基于PAM聚类分析的雷达点迹凝聚方法
技术领域
本发明属于雷达数据处理领域，是一种采用PAM(Partitioning Around Medoid，围绕中心点的划分)聚类分析算法的雷达目标点迹凝聚方法。
背景技术
在雷达探测目标的过程中，受到噪声干扰、目标自身的电磁散射效应等因素导致回波存在一定程度的起伏，一个目标可能在不同位置产生多个点迹，即出现目标分裂现象；同时在复杂杂波环境下，包括常见的气候杂波(云、雨、海浪等)以及干扰虚假杂波与真实目标的量测信息揉合在一起。为完成雷达的数据处理首先要找到同一目标分裂的各个可能点迹并将其凝聚起来。
传统的点迹凝聚方法首先寻找某一区域内幅度最大值点作为中心，然后将凝聚波门内的所有点迹通过质心法或比幅法进行凝聚处理，在背景比较干净的条件下效果尚可，但在杂波较多、航迹交叉或小范围内存在多批目标等复杂环境时凝聚得到的点迹可能存在较大偏差，从而导致后续的目标跟偏、跟丢等跟踪不稳定现象。
因此需要设计一种适用于复杂环境下的雷达点迹凝聚方法，以最大限度地减少点迹凝聚结果的可能偏差，进而提高雷达数据处理系统的跟踪精度和航迹稳定度。
孟凡等提出的“一种基于线性回归和聚类分析理论的高速机动目标跟踪方法及其实现”(申请号201210437109.2)专利申请中，对航迹相关时的分裂航迹使用聚类分析进行聚类，从而对高速机动目标进行跟踪，而本发明在点迹相关层面进行聚类处理，从而得到较高置信度的点迹数据。目前雷达数据处理领域已有的参考文献中：通常采用基于均值聚类分析技术(K-均值等)，其对雷达噪声和异常点迹过于敏感，本发明采用较为健壮的围绕中心点划分方法；区域内可能存在的目标个数的设定对最终的聚类划分有很大影响，一般采用模糊化处理，在实际工程实现中有较大偏差，本发明提出基于时间序列判决的多参数综合方法进行目标个数估计。
发明内容
本发明的目的在于提供一种利用基于PAM聚类和时间序列预测的点迹凝聚方法用于复杂环境下雷达目标点迹的提取处理。实现本发明目的的技术解决方案为：对雷达天线扫描周期内的点迹按方位划分扇区存储，对每个扇区内的点迹采用改进的PAM聚类算法进行聚类处理，PAM聚类算法中聚类个数初始值由当前扇区录取波门数和软件杂波图发现的运动目标数确定，通过PAM聚类后使得划分的每个簇中的点都具有较高的贴近度，作为一个可能目标的分裂点集合进行点迹凝聚处理，并对每个簇的聚点建立时间序列预测模型，后续PAM聚类算法的聚类个数由上一天线扫描周期每个簇聚点的时间序列预测结果、软件杂波图、跟踪波门、录取波门等参数综合判决得到。
本发明结构简单，逻辑清楚，计算复杂度低，易于工程实现。
本发明与现有技术相比，其显著优点为：该方法提高了复杂环境下雷达点迹凝聚的准确度和点迹提取的智能性。
下面结合附图2和附图3对本发明作进一步详细描述。
附图说明
附图1是本发明的功能实现框图。
附图2是本发明的工作流程图。
附图3是PAM聚类分析算法流程图。
具体实施方式
本发明采用一种基于PAM聚类分析的点迹凝聚方法进行雷达点迹的提取，具体实施步骤为(参见附图2)：
步骤1：将经检测处理后过慢门限的全程全方位点迹按方位扇区存储；
步骤2：提取扇区内疑似目标信息，计算聚类个数，方法如下：
(1)对当前扇区原始点迹建立三维软件杂波图，生成运动点迹；
(2)获取跟踪滤波、雷达综合显示、时间序列预测等软件模块输出的跟踪波门、录取波门、时间序列外推预测结果等中间信息；
(3)对(1)和(2)进行目标冗余归并处理；
步骤3：当前扇区内按照雷达的目标点迹的回波特性(包含点迹的方位、距离、仰角、信噪比、幅度、RCS、多普勒速度等)计算点迹相似度；
步骤4：按照点迹相似度进行PAM聚类处理，方法如下：
(1)选择距离步骤2中确定的K个目标位置相似度最贴近的点迹作为初始的簇中心点；
(2)REPEAT；
(3)将N-K个剩余的点迹指派给和它相似度最近的中心点所代表的簇，并计算每个剩余点迹与其簇中心的相似度；
(4)REPEAT；
(5)选择一个未被选择的中心点Oi，以Oi为中心建立相似度波门；
(6)REPEAT；
(7)选择一个中心点波门内未被选择过的非中心点对象Oh；
(8)计算用Oh代替Oi的总代价并记录在S中；
(9)UNTIL所有的中心点波门内的非中心点都被选择过；
(10)UNTIL所有的中心点都被选择过；
(11)IF在S中的所有非中心点代替所有中心点后的计算出的总代价有小于0的存在，THEN找出S中的用非中心点替代中心点后代价最小的一个，并用该非中心点替代对应的中心点，形成一个新的K个中心点的集合；
(12)UNTIL所有的S都大于0，没有再发生簇的重新分配；
步骤5：每个簇内采用贴近度加权法进行点迹凝聚处理，由n维欧几里德公式计算的贴近度作为点迹凝聚的权值，每点贴近度σi计算公式如下：
凝聚处理方法如下：
其中n为参与凝聚点数，为所有参与凝聚点的贴近度之和，αi、βi和γi为各分裂点迹的方位、距离和仰角参数，αs、βs和γs为凝聚处理后点迹的方位、距离和仰角；
步骤6：对每个簇按聚点的目标统计特性进行时间序列建模，并进行目标延续判决；
步骤7：按划分扇区迭代执行步骤2～步骤6。

Claims (2)

1.一种基于PAM聚类分析的雷达点迹凝聚方法，其特征在于：通过PAM聚类算法对扇区内点迹进行聚类，由当前扇区录取波门数和软件杂波图发现的运动目标数设置PAM聚类算法中聚类个数初始值，PAM聚类后使得划分的每个簇中的点都具有满足期望阈值的贴近度，作为一个可能目标的分裂点集合进行点迹凝聚处理，并对每个簇的聚点建立时间序列预测模型进行目标延续判决，由上一天线扫描周期每个簇聚点的时间序列预测结果、软件杂波图、跟踪波门、录取波门参数，综合计算得到后续PAM聚类算法的聚类个数；所述PAM为围绕中心点划分。
2.根据权利要求1所述的基于PAM聚类分析的雷达点迹凝聚方法，其特征在于：利用目标统计特性的特征贴近度进行加权计算，使用多维欧几里德公式计算的贴近度作为点迹凝聚的权值，加权平均后得到方位、距离、仰角的凝聚结果。
CN201510360108.6A 2015-06-25 2015-06-25 一种基于pam聚类分析的雷达点迹凝聚方法 Active CN104931934B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510360108.6A CN104931934B (zh) 2015-06-25 2015-06-25 一种基于pam聚类分析的雷达点迹凝聚方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510360108.6A CN104931934B (zh) 2015-06-25 2015-06-25 一种基于pam聚类分析的雷达点迹凝聚方法

Publications (2)

Publication Number Publication Date
CN104931934A CN104931934A (zh) 2015-09-23
CN104931934B true CN104931934B (zh) 2017-07-14

Family

ID=54119187

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510360108.6A Active CN104931934B (zh) 2015-06-25 2015-06-25 一种基于pam聚类分析的雷达点迹凝聚方法

Country Status (1)

Country Link
CN (1) CN104931934B (zh)

Families Citing this family (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106291503B (zh) * 2016-08-31 2018-12-28 西南交通大学 一种雷达探测系统密度聚类预警方法
CN109839621A (zh) * 2017-11-24 2019-06-04 西安艾索信息技术有限公司 一种改进的tbd算法
CN110426692B (zh) * 2019-07-20 2021-02-26 中国船舶重工集团公司第七二四研究所 参差加抖动中重频pd模式点迹提取方法
CN110542897B (zh) * 2019-08-01 2021-08-13 北京理工大学 一种基于希尔排序的距离差分多扩展目标点迹凝聚方法
CN111289954B (zh) * 2020-03-31 2022-03-15 四川长虹电器股份有限公司 用于毫米波雷达目标跟踪的点云划分与航迹匹配方法
RU2734071C1 (ru) * 2020-04-16 2020-10-12 Акционерное общество «Всероссийский научно-исследовательский институт радиотехники» (АО «ВНИИРТ») Межобзорное устройство картографирования пассивных помех для импульсно-доплеровских радиолокационных станций и способ его осуществления
CN111521983B (zh) * 2020-04-30 2022-03-08 中国民用航空总局第二研究所 一种基于tws雷达数据处理系统的时钟校准方法及系统
CN112363131B (zh) * 2020-09-03 2024-05-14 江西华讯方舟智能技术有限公司 车载毫米波雷达数据的处理方法、装置及计算机存储介质
CN112731296B (zh) * 2020-11-30 2024-03-08 惠州市德赛西威智能交通技术研究院有限公司 一种汽车毫米波雷达点迹凝聚方法及系统
CN114631035A (zh) * 2020-12-29 2022-06-14 深圳市大疆创新科技有限公司 获取点云聚类波门的方法、雷达、可移动平台和存储介质
CN113702964B (zh) * 2021-08-23 2023-09-26 中国北方工业有限公司 一种基于航迹信息的雷达自适应区域凝聚方法
CN113820685B (zh) * 2021-10-19 2024-07-19 中安锐达(北京)电子科技有限公司 一种基于扇块结构的一次点迹凝聚方法
CN114779204B (zh) * 2022-06-23 2022-09-06 北京海兰信数据科技股份有限公司 一种基于雷达目标幅度最小二乘跟踪处理方法及系统
CN116592894B (zh) * 2023-07-11 2023-12-26 天津中远海运散运数字科技有限公司 一种航线规划方法、装置及存储介质

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2011003453A1 (fr) * 2008-07-08 2011-01-13 Thales Traitement de donnees multi-cibles pour radars passifs multi-statiques et multi-canaux
CN104036146A (zh) * 2014-06-26 2014-09-10 中国电子科技集团公司第二十八研究所 一种用于雷达目标点迹凝聚的点迹聚类方法
CN104166135A (zh) * 2014-09-05 2014-11-26 武汉中原电子集团有限公司 一种宽带雷达目标的原始点迹凝聚处理方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2011003453A1 (fr) * 2008-07-08 2011-01-13 Thales Traitement de donnees multi-cibles pour radars passifs multi-statiques et multi-canaux
CN104036146A (zh) * 2014-06-26 2014-09-10 中国电子科技集团公司第二十八研究所 一种用于雷达目标点迹凝聚的点迹聚类方法
CN104166135A (zh) * 2014-09-05 2014-11-26 武汉中原电子集团有限公司 一种宽带雷达目标的原始点迹凝聚处理方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
一种基于三坐标雷达的点迹凝聚方法;周喃;《雷达与对抗》;20131231;第46-50页 *
基于自适应模糊聚类分析的目标冗余信息处理;李为民;《计算机应用》;20050430;第949-951页 *
雷达点迹凝聚处理技术及其数据分析;杨文琳等;《信号处理》;20010430;第130-137页 *
雷达点迹的目标智能特征提取方法研究;吉军等;《信息技术》;20130630;第83-87页 *

Also Published As

Publication number Publication date
CN104931934A (zh) 2015-09-23

Similar Documents

Publication Publication Date Title
CN104931934B (zh) 一种基于pam聚类分析的雷达点迹凝聚方法
CN109002835A (zh) 一种基于最大熵模糊聚类的粒子滤波数据关联方法
CN104714225B (zh) 一种基于广义似然比的动态规划检测前跟踪方法
CN104239712B (zh) 一种雷达抗干扰效能实时评估方法
CN104297748A (zh) 一种基于轨迹增强的雷达目标检测前跟踪方法
CN103353988B (zh) 异源sar景象特征匹配算法性能评估方法
Ni et al. Gait-based person identification and intruder detection using mm-wave sensing in multi-person scenario
CN102879766A (zh) 一种高频地波雷达弱目标检测跟踪方法和装置
CN103366365A (zh) 基于人工免疫多目标聚类的sar图像变化检测方法
CN105301584B (zh) 同时解距离模糊的ipphdf机动多目标跟踪方法
CN106054169A (zh) 基于跟踪信息的多站雷达信号融合检测方法
CN104181524B (zh) 一种自适应粒子数的多目标粒子滤波检测前跟踪方法
CN109242883A (zh) 基于深度sr-kcf滤波的光学遥感视频目标跟踪方法
CN104751469A (zh) 基于核模糊c均值聚类的图像分割方法
CN104881561A (zh) 一种基于Hough变换的多维参数的检测前跟踪方法
CN106934324A (zh) 基于简化多假设算法的雷达数据关联方法
CN103678949A (zh) 基于密度分析和谱聚类的多扩展目标跟踪量测集划分方法
Zhao et al. Mutation grey wolf elite PSO balanced XGBoost for radar emitter individual identification based on measured signals
CN111401168A (zh) 一种无人机的多层雷达特征提取与选择方法
Zhang et al. Weighted likelihood CFAR detection for Weibull background
CN107436434A (zh) 基于双向多普勒估计的航迹起始方法
CN111562569A (zh) 基于加权群稀疏约束的Weibull背景下多目标恒虚警检测方法
Zhang et al. Target detection for HFSWR based on an S 3 D algorithm
CN102621535B (zh) 一种高效的协方差矩阵结构估计方法
CN104463227A (zh) 基于fqpso和目标分解的极化sar图像分类方法

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