CN105893621B - 基于多维航迹聚类的目标行为规律挖掘方法 - Google Patents
基于多维航迹聚类的目标行为规律挖掘方法 Download PDFInfo
- Publication number
- CN105893621B CN105893621B CN201610280669.XA CN201610280669A CN105893621B CN 105893621 B CN105893621 B CN 105893621B CN 201610280669 A CN201610280669 A CN 201610280669A CN 105893621 B CN105893621 B CN 105893621B
- Authority
- CN
- China
- Prior art keywords
- track
- cluster
- neighbour
- clusterid
- goal behavior
- 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
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F16/00—Information retrieval; Database structures therefor; File system structures therefor
- G06F16/20—Information retrieval; Database structures therefor; File system structures therefor of structured data, e.g. relational data
- G06F16/28—Databases characterised by their database models, e.g. relational or object models
- G06F16/284—Relational databases
- G06F16/285—Clustering or classification
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F16/00—Information retrieval; Database structures therefor; File system structures therefor
- G06F16/20—Information retrieval; Database structures therefor; File system structures therefor of structured data, e.g. relational data
- G06F16/24—Querying
- G06F16/245—Query processing
- G06F16/2458—Special types of queries, e.g. statistical queries, fuzzy queries or distributed queries
- G06F16/2462—Approximate or statistical queries
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F16/00—Information retrieval; Database structures therefor; File system structures therefor
- G06F16/20—Information retrieval; Database structures therefor; File system structures therefor of structured data, e.g. relational data
- G06F16/28—Databases characterised by their database models, e.g. relational or object models
- G06F16/284—Relational databases
- G06F16/285—Clustering or classification
- G06F16/287—Visualization; Browsing
Abstract
本发明公开了一种基于多维航迹聚类的目标行为规律挖掘方法,所述方法包括:步骤1,设置目标的属性和类型标签;步骤2,计算目标航迹间的相似性度量;步骤3,选取每条航迹的近邻航迹;步骤4,把当前核心航迹分配到簇中;步骤5,把当前核心航迹的近邻中不属于其他簇的航迹分配的当前簇中;步骤6,把当前核心航迹的近邻中分配到当前簇的航迹中的核心航迹添加到当前核心航迹的近邻中;步骤7,重复步骤4,5,6;步骤8,设置目标行为规律编标签;步骤9,将目标行为规律可视化。该方法充分利用目标的属性、类型、位置、速度和航向信息,参数设置简单,准确率高,工程易实现,在模式识别和智能情报处理领域有广阔的应用前景。
Description
技术领域
本发明涉及数据挖掘中的聚类分析技术和信息融合中的高层融合技术,属于模式识别和智能情报处理领域。
背景技术
随着目标探测手段的不断丰富和多传感器信息融合技术的不断完善,各种各样的目标被检测、跟踪和识别,乃至形成连续稳定的目标航迹。海量的目标航迹数据在各种目标情报处理系统、预警监视系统、航行管制系统和视频监控系统中存储和积累,这些航迹数据中蕴含着大量的信息和知识。数据挖掘技术依靠跨学科的特点已经在各个领域得到了广泛应用,轨迹数据挖掘技术也成为数据挖掘领域的一个热点。在信息融合领域,基于电子信息系统中积累的大量历史航迹数据,利用数据挖掘技术中的无监督聚类技术,可以通过挖掘目标的多维特征来发现目标的行为规律,进而对其行为和意图进行分析和预测,这对于目标任务分类、异常情报检测、态势评估和指挥决策都具有非常重要的意义。
通过对目标航迹进行聚类分析,可以挖掘目标的行为规律。现有目标航迹聚类方法没有充分利用目标的属性、类型、位置、速度和航向等多维特征,只能反映目标空间位置变化的规律,在挖掘目标行为规律时具有局限性。
发明内容
本发明针对现有目标航迹聚类方法的局限性,提供了一种基于多维航迹聚类的目标行为规律挖掘方法。利用目标的属性、类型、位置、速度和航向信息,通过对多维航迹数据的聚类分析,实现对目标行为规律的挖掘。具体包括以下步骤:
步骤1,设置目标的属性和类型标签。本方法通过设置目标属性、类型和簇号标签来组成目标行为规律标签,从而表示目标的规律行为。在预警监视领域,将目标的属性分为我方,敌方和友方,分别用标签1,2和3表示,将目标类型分为军用飞机、民用飞机、军舰和民船等,分别用标签1,2,3和4等来表示。
步骤2,利用多维航迹数据中的位置、速度和航向信息,计算目标航迹间的多因素Hausdorff距离,作为目标航迹间的相似性度量,并存放在距离矩阵M中。
步骤3,输入近邻阈值ε,计算每条航迹TRi的ε近邻Nε(TRi)。
步骤4,设置簇号clusterId=1,输入最小近邻航迹数量MinTRs,从目标航迹数据库TD中随意选取一条航迹TRi,如果航迹TRi的ε近邻数量|Nε(TRi)|≥MinTRs,则该航迹为核心航迹,把它分配到簇CclusterId中。
步骤5,判断当前核心航迹的ε近邻Nε(TRi)中的航迹属不属于其他簇,将不属于其他簇的航迹分配到簇CclusterId中。
步骤6,判断当前核心航迹的ε近邻Nε(TRi)中分配到簇CclusterId中的航迹是不是核心航迹,如果是核心航迹,则把它添加到Nε(TRi)中。
步骤7,重复步骤5和步骤6,直到CclusterId不能再扩展。然后,令CclusterId=CclusterId+1,选择另外一条未分配到簇中的目标航迹进行步骤4,步骤5和步骤6。以此类推,直到C不能再扩展。
步骤8,设置目标行为规律标签。在预警监视领域,对敌方军用飞机的多维航迹数据聚类后得到的簇集为C={C1,C2,…,Cm},则目标行为规律的标签可以依次设置为211,213,...,21m。
步骤9,将目标行为规律可视化。为了直观的表示出目标的行为规律,我们可以选取每个簇中ε近邻数量最多的航迹作为行为规律的特征航迹。在可视化显示时只画出特征航迹,并且标出特征航迹对应的目标行为规律标签,用箭头表示出目标航向,用相邻航迹点间隔的长短表示目标的速度大小。
本发明的有益效果在于:本发明提供的一种基于多维航迹聚类的目标行为规律挖掘方法,充分利用了目标的属性、类型、位置、速度和航向信息,可以对不同属性、类型的目标行为规律进行区分,不仅可以挖掘出运动位置相近的目标行为规律,还可以区分速度和航向差异较大的目标行为。
附图说明
图1为一个人为构造的数据集中的目标航迹示意图。
图2为目标行为规律示意图。
具体实施方式
下面以一个人为构建的多维航迹数据集为例,根据附图详细阐述本发明,使本发明的技术路线和操作步骤更加清晰。
数据集:这个人为构建的多维航迹数据集包括代表8种目标行为规律的1200条航迹数据和300条不规律的航迹数据,TD={TR1,TR2,…,TR1500},目标航迹的顺序是随机的。多维航迹数据包含目标的属性、类型、位置、速度和航向信息。假设目标为敌方军用飞机。图1为数据集中的目标航迹示意图。
步骤1,设置目标的属性和类型标签。因为假设目标为敌方军用飞机,所以将属性标签设置为2,类型标签设置为1。
步骤2,利用多维航迹数据中的位置、速度和航向信息,计算目标航迹间的多因素Hausdorff距离,作为目标航迹间的相似性度量,并存放在距离矩阵M中,M为对称矩阵,并且其中的元素Mi,j代表第i条航迹TRi与第j条航迹TRj的相似性度量。多因素Hausdorff距离的具体计算方法如下:
1)考虑两个目标的位置信息、速度信息和运动方向信息,两个目标之间的多因素距离为:
其中表示点Pa与点Pa的速度信息,表示点Pa与点Pa的航向信息,wd表示位置信息的权重因子,wv表示速度信息的权重因子,wθ表示运动方向信息的权重因子,权重因子的取值取决于多因素距离的应用场景,满足wd≥0,wv≥0,wθ≥0和wd+wv+wθ=1。
2)基于多因素距离mfdist(Pa,Pb),航迹TRi到TRj的定向多因素Hausdorff距离为:
3)基于定向多因素Hausdorff距离,航迹TRi与TRj之间的多因素Hausdorff距离为:
步骤3,输入近邻阈值ε=0.03,计算每条航迹TRi的ε近邻Nε(TRi)。
Nε(TRi)={TRj∈TD|δM(TRi,TRj)≤ε} (4)
步骤4,设置簇号clusterId=1,输入最小近邻航迹数量MinTRs=4,从目标航迹数据库TD中随意选取一条航迹TRi,如果航迹TRi的ε近邻数量|Nε(TRi)|≥4,则该航迹为核心航迹,把它分配到簇CclusterId中。
步骤5,判断当前核心航迹的ε近邻Nε(TRi)中的航迹属不属于其他簇,将不属于其他簇的航迹分配到簇CclusterId中。
步骤6,判断当前核心航迹的ε近邻Nε(TRi)中的分配到簇CclusterId中的航迹是不是核心航迹,如果是核心航迹,则把它添加到Nε(TRi)中。
步骤7,重复步骤5和步骤6,直到CclusterId不能再扩展。然后,令CclusterId=CclusterId+1,选择另外一条未分配到簇中的目标航迹进行步骤4,步骤5和步骤6。以此类推,直到C不能再扩展。
步骤8,设置目标行为规律标签。对多维航迹数据聚类后得到的簇集为C={C1,C2,…,C8},则目标行为规律标签可依次设置为211,212,213,214,215,216,217,218。
步骤9,将目标行为规律可视化。为了直观的表示出目标的行为规律,选取每个簇中ε近邻数量最多的航迹作为行为规律的特征航迹。在可视化显示时只画出特征航迹,并且标出特征航迹对应的目标行为规律标签,用箭头表示出目标航向,用相邻航迹点间隔的长短表示目标的速度大小。图2为目标行为规律示意图。
Claims (3)
1.一种基于多维航迹聚类的目标行为规律挖掘方法,其特征在于,包括以下步骤:
步骤1,设置目标的属性和类型标签;
步骤2,利用多维航迹数据中的位置、速度和航向信息,计算目标航迹间的相似性度量,并存放在距离矩阵M中;
步骤3,输入近邻阈值ε,计算每条航迹TRi的ε近邻Nε(TRi);
步骤4,设置簇号clusterId=1,输入最小近邻航迹数量MinTRs,从目标航迹数据库TD中随意选取一条航迹TRi,如果航迹TRi的ε近邻数量|Nε(TRi)|≥MinTRs,则该航迹为核心航迹,把它分配到簇CclusterId中;
步骤5,判断当前核心航迹的ε近邻Nε(TRi)中的航迹属不属于其他簇,将不属于其他簇的航迹分配到簇CclusterId中;
步骤6,判断当前核心航迹的ε近邻Nε(TRi)中分配到簇CclusterId中的航迹是不是核心航迹,如果是核心航迹,则把该分配到簇CclusterId中的航迹的ε近邻添加到Nε(TRi)中;
步骤7,重复步骤5和步骤6,直到CclusterId不能再扩展,然后,令clusterId=clusterId+1,选择另外一条未分配到簇中的目标航迹进行步骤4,步骤5和步骤6,以此类推,直到簇集C不能再扩展;
步骤8,设置目标行为规律标签;
步骤9,将目标行为规律可视化;
所述步骤2中目标航迹间相似性度量的具体计算方法如下:
1)考虑两个目标的位置信息、速度信息和运动方向信息,两个目标之间的多因素距离为:
其中表示点Pa与点Pb的速度信息,表示点Pa与点Pb的航向信息,wd表示位置信息的权重因子,wv表示速度信息的权重因子,wθ表示运动方向信息的权重因子,权重因子的取值取决于多因素距离的应用场景,满足wd≥0,wv≥0,wθ≥0和wd+wv+wθ=1;
2)基于多因素距离mfdist(Pa,Pb),航迹TRi到TRj的定向多因素Hausdorff距离为:
3)基于定向多因素Hausdorff距离,航迹TRi与TRj之间的多因素Hausdorff距离为:
将δM(TRi,TRj)作为航迹TRi与TRj之间的相似性度量。
2.根据权利要求1所述的基于多维航迹聚类的目标行为规律挖掘方法,其特征在于,步骤3中航迹TRi的ε近邻Nε(TRi)的计算方法为:
Nε(TRi)={TRj∈TD|δM(TRi,TRj)≤ε}
TD={TR1,TR2,…,TRn}为目标多维航迹数据集。
3.根据权利要求1所述的基于多维航迹聚类的目标行为规律挖掘方法,其特征在于,步骤9中将目标规律可视化的具体方法为:
选取每个簇中ε近邻数量最多的航迹作为行为规律的特征航迹,在可视化显示时只画出特征航迹,并且标出特征航迹对应的目标行为规律标签,用箭头表示出目标航向,用相邻航迹点间隔的长短表示目标的速度大小。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610280669.XA CN105893621B (zh) | 2016-04-29 | 2016-04-29 | 基于多维航迹聚类的目标行为规律挖掘方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610280669.XA CN105893621B (zh) | 2016-04-29 | 2016-04-29 | 基于多维航迹聚类的目标行为规律挖掘方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105893621A CN105893621A (zh) | 2016-08-24 |
CN105893621B true CN105893621B (zh) | 2019-11-05 |
Family
ID=56703084
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610280669.XA Active CN105893621B (zh) | 2016-04-29 | 2016-04-29 | 基于多维航迹聚类的目标行为规律挖掘方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105893621B (zh) |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107480722B (zh) * | 2017-08-22 | 2020-03-17 | 中国人民解放军海军航空大学 | 基于归纳式一致性多类分类的目标行为模式在线分类方法 |
CN110135451B (zh) * | 2019-03-27 | 2020-06-26 | 中电莱斯信息系统有限公司 | 一种基于点到线段集距离的航迹聚类方法 |
CN109977546B (zh) * | 2019-03-27 | 2020-08-04 | 北京航空航天大学 | 一种基于无监督学习的四维航迹在线异常检测方法 |
CN110210537B (zh) * | 2019-05-22 | 2021-06-25 | 中国电子科技集团公司第二十八研究所 | 一种基于无监督聚类的船舶历史航迹规律提取方法 |
CN111881125B (zh) * | 2020-08-04 | 2023-06-02 | 中国人民解放军海军航空大学 | 一种海上非作战目标实时清洗方法及系统 |
CN113157800B (zh) * | 2021-03-31 | 2023-03-28 | 西南电子技术研究所(中国电子科技集团公司第十研究所) | 实时发现空中动态目标识别方法 |
CN114047508B (zh) * | 2022-01-13 | 2022-04-26 | 中国人民解放军海军工程大学 | 一种对海雷达探测范围统计评估方法及系统 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7263472B2 (en) * | 2004-06-28 | 2007-08-28 | Mitsubishi Electric Research Laboratories, Inc. | Hidden markov model based object tracking and similarity metrics |
CN104408401A (zh) * | 2014-10-28 | 2015-03-11 | 中国科学院自动化研究所 | 一种时敏目标的在轨检测方法 |
CN104715154A (zh) * | 2015-03-23 | 2015-06-17 | 西安交通大学 | 基于kmdl准则判据的核k-均值航迹关联方法 |
-
2016
- 2016-04-29 CN CN201610280669.XA patent/CN105893621B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7263472B2 (en) * | 2004-06-28 | 2007-08-28 | Mitsubishi Electric Research Laboratories, Inc. | Hidden markov model based object tracking and similarity metrics |
CN104408401A (zh) * | 2014-10-28 | 2015-03-11 | 中国科学院自动化研究所 | 一种时敏目标的在轨检测方法 |
CN104715154A (zh) * | 2015-03-23 | 2015-06-17 | 西安交通大学 | 基于kmdl准则判据的核k-均值航迹关联方法 |
Non-Patent Citations (1)
Title |
---|
时空约束下的热点路径空间分布检测算法;段炼 等;《计算机工程与设计》;20140316;第35卷(第3期);第861-866页 * |
Also Published As
Publication number | Publication date |
---|---|
CN105893621A (zh) | 2016-08-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105893621B (zh) | 基于多维航迹聚类的目标行为规律挖掘方法 | |
Etienne et al. | Spatio-temporal trajectory analysis of mobile objects following the same itinerary | |
CN111179638B (zh) | 一种基于时间序列的船舶ais目标航行监控方法 | |
CN110188093A (zh) | 一种基于大数据平台针对ais信息源的数据挖掘系统 | |
CN109000645A (zh) | 复杂环境目标经典航迹提取方法 | |
CN105654139A (zh) | 一种采用时间动态表观模型的实时在线多目标跟踪方法 | |
CN110175633A (zh) | 一种基于北斗卫星传输数据的渔船轨迹异常检测方法 | |
CN106022372B (zh) | 基于多维特征的目标异常行为检测方法 | |
CN111475544B (zh) | 船舶轨迹数据中离群点的检测方法和装置 | |
CN108710637A (zh) | 基于时空关系的出租车异常轨迹实时检测方法 | |
CN111582380A (zh) | 一种基于时空特征的船舶轨迹密度聚类方法及装置 | |
CN104715154B (zh) | 基于kmdl准则判据的核k‑均值航迹关联方法 | |
CN108804539A (zh) | 一种时间与空间双视角下的轨迹异常检测方法 | |
CN104091078B (zh) | 基于d‑s证据理论的产品多信息融合标识失效补救方法 | |
Toloue et al. | Anomalous behavior detection of marine vessels based on Hidden Markov Model | |
CN111291098A (zh) | 一种航迹模式挖掘方法与装置 | |
CN115147790A (zh) | 一种基于图神经网络的车辆未来轨迹预测方法 | |
Dästner et al. | Machine learning techniques for enhancing maritime surveillance based on GMTI radar and AIS | |
Ma et al. | Risk prediction for ship encounter situation awareness using long short-term memory based deep learning on intership behaviors | |
CN104008403A (zh) | 一种svm(矢量机)模式的多目标识别判定方法 | |
CN106203526B (zh) | 基于多维特征的目标行为模式在线分类方法 | |
CN115457300A (zh) | 一种基于距离度量和隔离机制的船舶异常行为检测方法 | |
CN114139642A (zh) | 联合船舶行程语义对象与图论理论的船舶航线提取方法 | |
CN105894014A (zh) | 基于多因素不一致度量的异常行为序贯检测方法 | |
Wang et al. | Classification of ship trajectories by using naive Bayesian algorithm |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
TA01 | Transfer of patent application right | ||
TA01 | Transfer of patent application right |
Effective date of registration: 20191011 Address after: 264001 Research and Academic Department, 188 Erma Road, Zhifu District, Yantai City, Shandong Province Applicant after: Naval Aviation University of PLA Address before: 264001 Yantai City, Zhifu Province, No. two road, No. 188, Department of research, Applicant before: Naval Aeronautical Engineering Institute PLA |
|
GR01 | Patent grant | ||
GR01 | Patent grant |