CN106022372A - 基于多维特征的目标异常行为检测方法 - Google Patents

基于多维特征的目标异常行为检测方法 Download PDF

Info

Publication number
CN106022372A
CN106022372A CN201610330313.2A CN201610330313A CN106022372A CN 106022372 A CN106022372 A CN 106022372A CN 201610330313 A CN201610330313 A CN 201610330313A CN 106022372 A CN106022372 A CN 106022372A
Authority
CN
China
Prior art keywords
flight path
neighbour
target
multifactor
multidimensional
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
CN201610330313.2A
Other languages
English (en)
Other versions
CN106022372B (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.)
Naval Aeronautical University
Original Assignee
Naval Aeronautical Engineering Institute of PLA
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 Naval Aeronautical Engineering Institute of PLA filed Critical Naval Aeronautical Engineering Institute of PLA
Priority to CN201610330313.2A priority Critical patent/CN106022372B/zh
Publication of CN106022372A publication Critical patent/CN106022372A/zh
Application granted granted Critical
Publication of CN106022372B publication Critical patent/CN106022372B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00Pattern recognition
    • G06F18/20Analysing
    • G06F18/25Fusion techniques
    • G06F18/254Fusion techniques of classification results, e.g. of results related to same input data
    • G06F18/256Fusion techniques of classification results, e.g. of results related to same input data of results relating to different input data, e.g. multimodal recognition
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00Pattern recognition
    • G06F18/20Analysing
    • G06F18/25Fusion techniques
    • G06F18/251Fusion techniques of input or preprocessed data

Landscapes

  • Engineering & Computer Science (AREA)
  • Data Mining & Analysis (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Artificial Intelligence (AREA)
  • Evolutionary Biology (AREA)
  • Evolutionary Computation (AREA)
  • Physics & Mathematics (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

本发明公开了一种基于多维特征的目标异常行为检测方法。该方法充分利用目标的属性、类型、位置、速度和航向特征,通过对多维航迹数据的异常检测,实现对目标异常行为的挖掘,具体包括:步骤1,输入目标的多维航迹数据集,设置目标的属性和类型标签;步骤2,计算目标航迹间的多因素定向Hausdorff距离;步骤3,确定每条航迹的近邻航迹;步骤4,计算每条航迹的近邻密度;步骤5,计算每条航迹的多维度局部异常因子;步骤6,对每条航迹进行异常检测判决;步骤7,设置目标异常行为标签。所述方法参数设置简单,准确率高,工程易实现,在模式识别和智能情报处理领域有广阔的应用前景。

Description

基于多维特征的目标异常行为检测方法
技术领域
本发明涉及数据挖掘中的异常检测技术和信息融合中的高层融合技术,属于模式识别和智能情报处理领域。
背景技术
目标航迹数据是由多维数据点组成的多维序列。根据应用场景,可将航迹数据分为预警监视航迹数据,航行管制航迹数据和视频监控航迹数据等;根据目标的类型,可将航迹数据分为飞机航迹数据、船舶航迹数据、车辆航迹数据、行人航迹数据、动物航迹数据和龙卷风航迹数据等。不同的应用场景中,目标航迹数据的多维特征也不一样,例如,在广播式自动相关监视系统中,航迹数据通常包括飞机编号、时间、经度、纬度、高度、速度、航向等多维特征;在预警监视情报处理系统中,航迹数据通常包括目标的批号、属性、类别、数量、型号、机/舷号、时间、经度、维度、高度、速度、航向等多维特征。在信息融合领域,利用数据挖掘中的异常检测技术,可以基于目标的多维航迹特征来挖掘目标的异常行为,实现对异常情报的智能处理,这对于态势评估、威胁估计和指挥决策都具有非常重要的意义。
通过对目标航迹的异常检测,可以挖掘目标的异常行为。一些学者在航迹异常检测方面进行了研究,但是现有方法主要检测目标航迹的位置异常,没有充分利用目标的属性、类型、位置、速度和航向等多维特征,在挖掘目标的异常行为时具有局限性。
发明内容
本发明针对现有航迹异常检测方法的局限性,提供了一种基于多维特征的目标异常行为检测方法。利用目标的属性、类型、位置、速度和航向特征,通过对多维航迹数据的异常检测,实现对目标异常行为的挖掘。具体包括以下步骤:
步骤1,输入目标的多维航迹数据集TD={TR1,TR2,…,TRn},设置目标的属性和类型标签。本方法通过设置目标属性、类型和异常行为序号标签来组成目标异常行为标签,从而表示目标的异常行为。例如,在预警监视领域,将目标的属性分为我方,敌方和友方,分别用标签1,2和3表示,将目标类型分为军用飞机、民用飞机、军舰和民船等,分别用标签1,2,3和4等来表示。
步骤2,利用多维航迹数据中的位置、速度和航向特征,计算目标航迹间的多因素定向Hausdorff距离。
1)考虑目标的位置、速度和航向特征,两个目标数据点之间的多因素距离为:
m f d i s t ( P a , P b ) = w d · d i s t ( P a , P b ) + w v · d i s t ( v P a , v P b ) + w θ · d i s t ( θ P a , θ P b ) - - - ( 1 )
其中dist(Pa,Pb)表示点Pa与点Pb之间位置特征的欧式距离,表示点Pa与点Pb的速度,表示点Pa与点Pb的航向,表示点Pa与点Pb之间速度特征的欧式距离,表示点Pa与点Pb之间航向特征的欧式距离,wd表示位置特征的权重因子,wv表示速度特征的权重因子,wθ表示航向特征的权重因子,权重因子的取值取决于多因素距离的应用场景,满足wd≥0,wv≥0,wθ≥0和wd+wv+wθ=1。
2)基于多因素距离mfdist(Pa,Pb),航迹TRi到TRj的多因素定向Hausdorff距离为:
δ M → ( TR i , TR j ) = m a x P a ∈ TR i { m i n P b ∈ TR j { m f d i s t ( P a , P b ) } } - - - ( 2 )
多因素定向Hausdorff距离用来度量航迹TRi与TRj的相似程度。
步骤3,输入近邻数k,利用目标航迹间的多因素定向Hausdorff距离,确定每条航迹TRi的近邻航迹
1)利用近邻数k和目标航迹间的多因素定向Hausdorff距离,确定多维航迹数据集TD中与TRi的多因素定向Hausdorff距离第k近的航迹NN(TRi,k);
2)航迹TRi的多因素近邻边界距离为:
δ M → b o r d e r ( TR i , k ) = δ M → ( TR i , N N ( TR i , k ) ) - - - ( 3 )
3)每条航迹TRi的近邻航迹为:
当近邻数为k时,航迹TRi的近邻航迹为
步骤4,利用步骤3确定的近邻航迹计算每条航迹TRi的近邻密度ρ(TRi,k)。
1)航迹TRi到TRj的多因素可达距离为:
δ M → r e a c h ( TR i , TR j , k ) = m a x { δ M → b o r d e r ( TR j , k ) , δ M → ( TR j , TR i ) } - - - ( 5 )
2)航迹TRi的近邻密度为:
当近邻数为k时,航迹TRi的近邻密度为ρ(TRi,k)。
步骤5,利用每条航迹TRi和TRi近邻航迹的近邻密度,计算TRi的多维度局部异常因子MDLOF(TRi,k):
航迹TRi的异常程度越大,多维度局部异常因子的取值越大。
步骤6,输入异常阈值ε,对每条航迹TRi进行异常检测判决,并将结果存放到异常指示向量中。
如果MDLOF(TRi,k)>ε,否则,
步骤7,设置目标异常行为标签。例如,在预警监视领域,对民用航班的多维航迹数据进行异常检测后得到的异常指示向量为:则目标的异常行为标签可以依次设置为3201、3202、...。
本发明的有益效果在于:本发明提供的一种基于多维特征的目标异常行为检测方法,充分利用了目标的属性、类型、位置、速度和航向信息,可以对不同属性、类型的目标异常行为进行区分,而且可以挖掘目标位置、速度和运动方向与周围目标运动规律差异较大的目标异常行为。
附图说明
图1为民用航班在三维空间中的航迹图。
图2为民用航班在二维平面中的航迹图。
图3为民用航班的异常行为在二维平面中的示意图。
具体实施方式
下面以一个广播式自动相关监视系统(ADS-B)接收到的部分民用航班的飞行航迹数据集为例,根据附图详细阐述本发明,使本发明的技术路线和操作步骤更加清晰。
该ADS-B数据集包括2015年5月接收到的237条民航飞行航迹,每条航迹包括若干个多维数据点。我们可以直接读取目标的维度、经度、高度组成的位置特征,速度特征和航向特征,目标属性为友方,类型为民用飞机。为了计算目标航迹之间的多因素定向Hausdorff距离,我们将目标航迹数据点的位置特征由地理坐标中的维度、经度和高度坐标转换到局部直角坐标系中。在转换时,选择北京首都国际机场的地理坐标作为局部直角坐标系的原点,然后将这237条目标航迹的所有数据点都转换到这个局部直角坐标系中。图1为这些民用航班在三维空间中的航迹图,图2为这些民用航班在二维平面中的航迹图。
步骤1,输入多维航迹数据集TD={TR1,TR2,…,TR237},设置目标的属性和类型标签。因为目标为民用航班,所以将属性标签设置为3,类型标签设置为2。
步骤2,利用ADS-B多维航迹数据中的位置、速度和航向特征,由式(1)计算目标数据点之间的多因素距离mfdist(Pa,Pb),由式(2)计算目标航迹间的多因素定向Hausdorff距离
步骤3,输入近邻数k=4,利用目标航迹间的多因素定向Hausdorff距离,由式(3)和式(4)确定每条航迹TRi的近邻航迹
步骤4,利用步骤3确定的近邻航迹由式(5)和式(6)计算每条航迹TRi的近邻密度ρ(TRi,k)。
步骤5,利用每条航迹TRi和TRi近邻航迹的近邻密度,由式(7)计算TRi的多维度局部异常因子MDLOF(TRi,k)。
步骤6,输入异常阈值ε=1,对每条航迹TRi进行异常检测判决,并将结果存放到异常指示向量中,共得到6条异常航迹。
步骤7,设置目标异常行为标签。因为目标的属性标签为3,类型标签为2,异常行为序号为1、2、3、4、5、6,所以目标的异常行为标签可以依次设置为3201、3202、3203、3204、3205、3206。
图3为民用航班的异常行为在二维平面中的示意图,箭头表示飞机的运动方向。

Claims (5)

1.一种基于多维特征的目标异常行为检测方法,其特征在于,包括以下步骤:
步骤1,输入目标的多维航迹数据集TD={TR1,TR2,…,TRn},设置目标的属性和类型标签;
步骤2,利用多维航迹数据中的位置、速度和航向特征,计算目标航迹间的多因素定向Hausdorff距离;
步骤3,输入近邻数k,利用目标航迹间的多因素定向Hausdorff距离,确定每条航迹TRi的近邻航迹
步骤4,利用步骤3确定的近邻航迹计算每条航迹TRi的近邻密度ρ(TRi,k);
步骤5,利用每条航迹TRi和TRi近邻航迹的近邻密度,计算TRi的多维度局部异常因子MDLOF(TRi,k);
步骤6,输入异常阈值ε,对每条航迹TRi进行异常检测判决,并将结果存放到异常指示向量中;
步骤7,设置目标异常行为标签。
2.根据权利要求1所述的基于多维特征的目标异常行为检测方法,其特征在于,步骤2具体为:
1)考虑目标的位置、速度和航向特征,两个目标数据点之间的多因素距离为:
m f d i s t ( P a , P b ) = w d · d i s t ( P a , P b ) + w v · d i s t ( v P a , v P b ) + w θ · d i s t ( θ P a , θ P b )
其中dist(Pa,Pb)表示点Pa与点Pb之间位置特征的欧式距离,表示点Pa与点Pb的速度,表示点Pa与点Pb的航向,表示点Pa与点Pb之间速度特征的欧式距离,表示点Pa与点Pb之间航向特征的欧式距离,wd表示位置特征的权重因子,wv表示速度特征的权重因子,wθ表示航向特征的权重因子,权重因子的取值取决于多因素距离的应用场景,满足wd≥0,wv≥0,wθ≥0和wd+wv+wθ=1;
2)基于多因素距离mfdist(Pa,Pb),航迹TRi到TRj的多因素定向Hausdorff距离为:
δ M → ( TR i , TR j ) = m a x P a ∈ TR i { m i n P b ∈ TR j { m f d i s t ( P a , P b ) } }
多因素定向Hausdorff距离用来度量航迹TRi与TRj的相似程度。
3.根据权利要求1所述的基于多维特征的目标异常行为检测方法,其特征在于,步骤3具体为:
1)输入近邻数k,利用目标航迹间的多因素定向Hausdorff距离,确定多维航迹数据集TD中与TRi的多因素定向Hausdorff距离第k近的航迹NN(TRi,k);
2)航迹TRi的多因素近邻边界距离为:
δ M → b o r d e r ( TR i , k ) = δ M → ( TR i , N N ( TR i , k ) ) ;
3)每条航迹TRi的近邻航迹为:
当近邻数为k时,航迹TRi的近邻航迹为
4.根据权利要求1所述的基于多维特征的目标异常行为检测方法,其特征在于,步骤4具体为:
1)航迹TRi到TRj的多因素可达距离为:
δ M → r e a c h ( TR i , TR j , k ) = max { δ M → b o r d e r ( TR j , k ) , δ M → ( TR j , TR i ) }
2)航迹TRi的近邻密度为:
当近邻数为k时,航迹TRi的近邻密度为ρ(TRi,k)。
5.根据权利要求1所述的基于多维特征的目标异常行为检测方法,其特征在于,步骤5具体为:
通过比较TRi与TRi近邻航迹的近邻密度来计算TRi的多维度局部异常因子:
航迹TRi的异常程度越大,多维度局部异常因子的取值越大。
CN201610330313.2A 2016-05-17 2016-05-17 基于多维特征的目标异常行为检测方法 Active CN106022372B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610330313.2A CN106022372B (zh) 2016-05-17 2016-05-17 基于多维特征的目标异常行为检测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610330313.2A CN106022372B (zh) 2016-05-17 2016-05-17 基于多维特征的目标异常行为检测方法

Publications (2)

Publication Number Publication Date
CN106022372A true CN106022372A (zh) 2016-10-12
CN106022372B CN106022372B (zh) 2019-05-10

Family

ID=57098082

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610330313.2A Active CN106022372B (zh) 2016-05-17 2016-05-17 基于多维特征的目标异常行为检测方法

Country Status (1)

Country Link
CN (1) CN106022372B (zh)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107480722A (zh) * 2017-08-22 2017-12-15 中国人民解放军海军航空工程学院 基于归纳式一致性多类分类的目标行为模式在线分类方法
CN110175633A (zh) * 2019-04-29 2019-08-27 杭州电子科技大学 一种基于北斗卫星传输数据的渔船轨迹异常检测方法
CN110288855A (zh) * 2019-06-05 2019-09-27 上海鹰觉科技有限公司 针对海上船只行为的分类方法及其分类系统
CN113538974A (zh) * 2021-07-14 2021-10-22 电子科技大学 一种基于多源数据融合的飞行目标异常检测方法
CN116071667A (zh) * 2023-04-07 2023-05-05 北京理工大学 基于历史数据的指定区域内异常飞机目标检测方法及系统

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20090005985A1 (en) * 2007-06-29 2009-01-01 Gm Global Technology Operations, Inc. GPS-based in-vehicle sensor calibration algorithm
CN104899263A (zh) * 2015-05-22 2015-09-09 华中师范大学 一种基于特定区域的船舶轨迹挖掘分析与监控方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20090005985A1 (en) * 2007-06-29 2009-01-01 Gm Global Technology Operations, Inc. GPS-based in-vehicle sensor calibration algorithm
CN104899263A (zh) * 2015-05-22 2015-09-09 华中师范大学 一种基于特定区域的船舶轨迹挖掘分析与监控方法

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107480722A (zh) * 2017-08-22 2017-12-15 中国人民解放军海军航空工程学院 基于归纳式一致性多类分类的目标行为模式在线分类方法
CN110175633A (zh) * 2019-04-29 2019-08-27 杭州电子科技大学 一种基于北斗卫星传输数据的渔船轨迹异常检测方法
CN110175633B (zh) * 2019-04-29 2020-12-15 杭州电子科技大学 一种基于北斗卫星传输数据的渔船轨迹异常检测方法
CN110288855A (zh) * 2019-06-05 2019-09-27 上海鹰觉科技有限公司 针对海上船只行为的分类方法及其分类系统
CN113538974A (zh) * 2021-07-14 2021-10-22 电子科技大学 一种基于多源数据融合的飞行目标异常检测方法
CN116071667A (zh) * 2023-04-07 2023-05-05 北京理工大学 基于历史数据的指定区域内异常飞机目标检测方法及系统
CN116071667B (zh) * 2023-04-07 2023-07-18 北京理工大学 基于历史数据的指定区域内异常飞机目标检测方法及系统

Also Published As

Publication number Publication date
CN106022372B (zh) 2019-05-10

Similar Documents

Publication Publication Date Title
CN106022372B (zh) 基于多维特征的目标异常行为检测方法
CN110309383B (zh) 基于改进的dbscan算法的船舶轨迹聚类分析方法
Liu et al. Self-adaptive dynamic obstacle avoidance and path planning for USV under complex maritime environment
CN104062973B (zh) 一种基于图像标志物识别的移动机器人slam方法
CN105893621B (zh) 基于多维航迹聚类的目标行为规律挖掘方法
CN102081801B (zh) 多特征自适应融合船舶跟踪和航迹检测方法
CN106355194A (zh) 一种基于激光成像雷达的无人艇水面目标处理方法
CN103729861B (zh) 一种多目标跟踪方法
CN104715154B (zh) 基于kmdl准则判据的核k‑均值航迹关联方法
CN103810382B (zh) 机载分布式多传感器两级数据融合策略选择方法
CN109144060A (zh) 一种船舶航线的危险识别方法及系统
CN102997911A (zh) 一种被动传感器组网探测多目标方法
CN107738852A (zh) 定位方法、定位地图构建方法及机器人
CN108106617A (zh) 一种无人机自动避障方法
Shan et al. LiDAR-based stable navigable region detection for unmanned surface vehicles
CN104050368A (zh) 系统误差下基于误差补偿的群航迹精细关联算法
CN111458702B (zh) 一种机动目标跟踪方法、系统、电子设备及存储介质
CN114462519A (zh) 一种基于航行轨迹特征的船舶异常行为识别方法
CN114139373B (zh) 一种面向无人机反制车的多传感器自动协同管理方法
Patterson et al. Utilizing geographic information system data for unmanned aerial vehicle position estimation
CN105353352B (zh) 改进搜索策略的mm‑pphdf机动多目标跟踪方法
Teacy et al. Observation modelling for vision-based target search by unmanned aerial vehicles
Hofstetter et al. On ambiguities in feature-based vehicle localization and their a priori detection in maps
Cheng et al. Graph-based proprioceptive localization using a discrete heading-length feature sequence matching approach
Veneruso et al. Analysis of ground infrastructure and sensing strategies for all-weather approach and landing in Urban Air Mobility

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: 20190416

Address after: 264001 scientific research office, two Road 188, Zhifu District, Yantai, Shandong

Applicant after: Naval Aviation University of PLA

Address before: 264001 scientific research department, 188 road two, Zhifu District, Yantai, Shandong

Applicant before: Naval Aeronautical Engineering Institute PLA

GR01 Patent grant
GR01 Patent grant