CN102338874B - 用于被动式多传感器目标跟踪的全局概率数据关联方法 - Google Patents
用于被动式多传感器目标跟踪的全局概率数据关联方法 Download PDFInfo
- Publication number
- CN102338874B CN102338874B CN 201110172026 CN201110172026A CN102338874B CN 102338874 B CN102338874 B CN 102338874B CN 201110172026 CN201110172026 CN 201110172026 CN 201110172026 A CN201110172026 A CN 201110172026A CN 102338874 B CN102338874 B CN 102338874B
- Authority
- CN
- China
- Prior art keywords
- sensor
- measurement
- target
- probability
- tracking
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 58
- 238000012545 processing Methods 0.000 claims abstract description 13
- 238000012163 sequencing technique Methods 0.000 claims abstract description 4
- 238000005259 measurement Methods 0.000 claims description 86
- 239000011159 matrix material Substances 0.000 claims description 18
- 238000005315 distribution function Methods 0.000 claims description 5
- 238000004422 calculation algorithm Methods 0.000 claims description 3
- 238000001514 detection method Methods 0.000 claims description 3
- 238000009795 derivation Methods 0.000 claims description 2
- 238000013139 quantization Methods 0.000 claims description 2
- 238000010276 construction Methods 0.000 claims 1
- CZZYITDELCSZES-UHFFFAOYSA-N diphenylmethane Chemical compound C=1C=CC=CC=1CC1=CC=CC=C1 CZZYITDELCSZES-UHFFFAOYSA-N 0.000 claims 1
- 230000000694 effects Effects 0.000 abstract description 5
- 230000007547 defect Effects 0.000 description 4
- 238000004364 calculation method Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 238000003672 processing method Methods 0.000 description 2
- 238000007476 Maximum Likelihood Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000010408 sweeping Methods 0.000 description 1
Images
Landscapes
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开用于被动式多传感器目标跟踪的全局概率数据关联方法,该方法能够有效的克服原有方法非线性部分局部线性化造成的误差以及由于量测不完全可观性导致的数据关联中的概率分布不完全,原方法中有偏的概率分布估计也导致了其跟踪门设置的不合理,在本方法中也得到了更正。此外本方法中还提出了一种多传感器顺序处理排序的机制,能够有效的提高目标跟踪效果。本方法尤其适用于量测原始数据可靠性不高的被动式无线传感器网络目标跟踪系统中。
Description
技术领域
本发明涉及用于被动式多传感器目标跟踪的全局概率数据关联方法,尤其涉及一种适用于量测原始数据可靠性不高的被动式无线传感器网络目标跟踪系统中的基于目标位置预测的全局概率数据关联方法。
背景技术
目标跟踪理论和技术在民用和军事领域都有重要的应用价值,该领域一直是学术界、工业界乃至军方重点发展的方向。传统的主动式目标跟踪系统,如雷达,通过发射雷达波并接收反射回波来获得目标方位角和距离的量测,容易被敌方轻易侦测到跟踪系统的位置,进而遭到敌方摧毁。因此,在该应用场景中,被动式传感器体现了很大的优势。但被动式传感器仅能得到目标方位角,而不能有效的获得目标距离的量测,因此单个传感器不能提供足够的量测信息来保证目标位置的可观性,需要多传感器进行协同工作。本发明针对的即是被动式多传感器的目标跟踪系统,尤其适用于被动式无线传感器网络目标跟踪系统。
在目标跟踪系统中,由于跟踪环境的干扰以及传感器自身性能的限制,在测量过程中不可避免的引入量测噪声和杂波干扰,即使目标只有一个,有效量测可能有多个。需要运用数据关联技术来建立目标源和量测的对应关系,并在正确的数据关联前提下进行目标跟踪过程。适用于主动式传感器目标跟踪系统的数据关联方法已趋于成熟,其中概率数据关联方法在杂波环境的单目标跟踪问题中具有很好的理论基础和实用效果。目前并没有专门针对被动式传感器目标跟踪的数据关联方法,现有的方法是针对被动式传感器的量测特点,利用观测方程中量测与目标状态的关系将目标状态协方差矩阵转化为量测协方差矩阵,再将方位角的量测直接作为被关联对象,并沿用原有主动式传感器目标跟踪的概率数据关联方法,实现被动式传感器目标跟踪的数据关联。
现有的用于被动式传感器的概率数据关联的基本方法的主要步骤如下:
4.根据上步得到的关联概率,应用全期望公式,对目标状态和估计协方差矩阵进行更新:
5.以上为单传感器数据关联过程,再利用顺序处理结构序贯处理多传感器量测,即得到多传感器下目标跟踪的结果。
现有的用于被动式传感器的概率数据关联的方法虽然能实现目标跟踪功能,但由于被动式传感器的量测对目标状态不完全可观的缺陷,直接沿用传统概率数据关联的方法存在一些问题:首先,被动式传感器的量测与目标状态的对应关系存在的非线性,局部线性化仅是一种近似的逼近,因此量测残差协方差的计算存在一定的误差;其次,由于被动式传感器的量测仅是1维信号,比2维的目标位置缺失了一部分信息,且量测残差协方差的分布一般并不关于预测量测方向角对称,现有方法中隐含假设了该条件,在跟踪门设置以及关联概率计算时造成了一定的偏差;此外,针对多传感器的应用,顺序处理结构的序贯处理方法需要按照合理的次序对各传感器量测进行处理,若采用随机确定传感器量测的处理顺序,则忽略了顺序处理方法中处理次序对结果的影响,显然效果不理想,若按照传感器量测误差协方差排序,则忽略了被动式传感器量测的不完全性,这种不完全的观测可能会导致对目标状态的更新所作的贡献变得很小,因此该排序方法也不是全局最优的。
发明内容
本发明的目的在于针对现有技术的不足,提供用于被动式多传感器目标跟踪的全局概率数据关联方法。
为了实现上述目的,本发明采取如下技术方案:用于被动式多传感器目标跟踪的全局概率数据关联方法,包括以下步骤:
(1)随机部署若干个被动式传感器,并确定各传感器的位置以及量测误差协方差。
(2)根据各传感器位置和传感器量测误差协方差,并结合设定的椭圆跟踪门区域,量化各传感器量测对目标跟踪更新的贡献程度,并由此对传感器量测的处理顺序进行排序。
(3)根据先验预测的目标位置和对应的误差协方差矩阵,得到跟踪门内目标位置的相对概率分布。
(4)根据当前处理的传感器位置和待关联的量测,得到观测射线方程,并计算观测射线与跟踪门边界的交点,由此判决是否接受量测。
(5)根据得到的目标位置的相对概率分布函数和临界交点,得到跟踪门内量测方向上各点的概率积分,并进一步得到量测源自目标的概率。
(6)根据上步得到的关联概率,应用全期望公式,对目标状态和估计协方差矩阵进行更新。
(7)重复步骤(3)~(6),完成顺序结构的被动式多传感器目标跟踪,得到目标实时位置信息。
本发明的有益效果是,本发明很好地克服了原有直接沿用传统概率数据关联的被动式多传感器目标跟踪方法中的不足;首先,由于直接考察目标位置误差协方差矩阵,避免了原方法中非线性近似的误差;其次,本方法中应用了量测对应的全局先验概率,使得关联概率的计算是无偏的,同时更正了原有方法跟踪门判决规则在该应用场景中存在的偏差;此外,本方法提出的量化被动式传感器对目标跟踪更新贡献程度的方法充分考虑了传感器量测特性以及被动式传感器不完全观测性质,能够有效的对处理顺序进行排序,以得到全局最优的估计结果。
附图说明
图1是被动式多传感器目标跟踪场景示意图;
图2是本发明中量化被动式传感器量测贡献方法示意图;
图3是本发明中跟踪门判决方法示意图;
图4是本发明中量测的先验概率分布和原方法的对比图;
图5是本发明的跟踪效果和原有方法的对比图。
具体实施方式
本发明专利在原有概率数据关联方法的基础上,通过分析被动式传感器的量测特性以及目标跟踪的应用场景,实现了适用于被动式多传感器目标跟踪的全局概率数据关联方法,并能够扩展为多目标跟踪情形。本方法的原理是通过直接分析目标位置的相关信息来实现全局的概率数据关联,克服了原有方法计算关联概率的有偏性,尤其适用于量测原始数据可靠性不高的被动式无线传感器网络目标跟踪系统中。
本发明用于被动式多传感器目标跟踪的全局概率数据关联方法,包括以下步骤:
1、随机部署若干个被动式传感器,并确定各传感器的位置以及量测误差协方差。
2、根据各传感器位置和传感器量测误差协方差,并结合设定的椭圆跟踪门区域,量化各传感器量测对目标跟踪更新的贡献程度,由此对传感器量测的处理顺序进行排序。
如附图2所示,构造量化直线,斜率k为,其中经过为传感器i预测的量测角,可以利用如下公式求得:,并经过点。由此可以求得直线在椭圆内的线段长度,记为。同时传感器i的量测误差的标准差,取2倍标准差范围扫过直线的部分,计算该线段长度:,其中为传感器到目标预测位置的距离,可通过公式计算如下:。
3、根据先验预测的目标位置和对应的误差协方差矩阵,得到跟踪门内目标位置的相对概率分布。
根据先验的目标位置预测和对应的误差协方差矩阵,得到目标位置的相对概率分布:。这里考虑到各传感器的误差协方差,需要更新误差协方差阵:,其中,用 代替概率分布函数中的计算 ,同理得到更新的跟踪门边界椭圆方程。
4、根据当前处理的传感器位置和待关联的量测,得到观测射线方程,并计算观测射线与跟踪门边界的交点,由此判决是否接受量测。
根据当前传感器位置和待关联量测,得到观测射线方程,如附图3所示,并得到观测射线与椭圆方程的交点。若有两个交点,则判决量测在跟踪门内,接受量测;若只有一个交点或是没有交点,则判决量测在跟踪门外,拒绝量测。易得出以下结论:在被动式传感器中,接受的量测并不是关于预测中心量测对称的,这和原方法是不同的,显然本方法的跟踪门设定规则对于位置估计是全局最优的。
5、根据得到的目标位置的相对概率分布函数和临界交点,得到跟踪门内量测方向上各点的概率积分,并进一步得到量测源自目标的概率。
根据上步得到的接受的量测集合,其中为确认量测数目,分别计算量测对应的所有跟踪门内目标位置的先验概率积分,并利用经过预测位置的量测概率积分对结果进行归一化,可表示为,其中和分别是与椭圆方程两交点的横坐标。图4显示了在某一相同参数条件下,原方法和本方法得到的量测的先验概率。
这里考虑到积分运算的算法复杂度较高,给出一种近似的最大似然的估计方法:根据上步得到的与椭圆方程两交点,计算两点的平均值,易证明该点也落在量测射线方程上且是射线方程上先验概率最大值点,并取该值作为量测的先验概率。
再根据概率数据关联算法中的相关推导得到量测源自目标的概率:
;
6、根据上步得到的关联概率,应用全期望公式,对目标状态和估计协方差矩阵进行更新。
利用如下方程计算:
,
7、重复步骤(3)~(6),完成顺序结构的被动式多传感器目标跟踪,到目标实时位置信息。
本发明考虑到被动式传感器目标跟踪的具体应用场景,采用了一种新型的全局概率数据关联方法,无偏地得到了量测对应的先验概率分布,并提出了一种有效的传感器处理顺序排序方法,进一步提高估计精度。本发明的有益之处在于避免了原数据关联方法中线性化误差以及有偏的概率分布估计,同时发现并改进了原方法简单的对称跟踪门规则的不足,从而得到更为可信的数据关联结果,使得目标跟踪效果提高,图5给出了20次相同参数设置下,原有方法和本方法的跟踪误差比较,可以看出本方法具有一定的性能优势,跟踪误差较为平稳且处于较低范围。
Claims (1)
1.一种用于被动式多传感器目标跟踪的全局概率数据关联方法,其特征在于,包括以下步骤:
(1)随机部署若干个被动式传感器,并确定各传感器的位置以及量测误差协方差;
(2)根据各传感器位置和传感器量测误差协方差,并结合设定的椭圆跟踪门区域,量化各传感器量测对目标跟踪更新的贡献程度,由此对传感器量测的处理顺序进行排序;
(3)根据先验预测的目标位置和对应的误差协方差矩阵,得到跟踪门内目标位置的相对概率分布;
(4)根据当前处理的传感器位置和待关联的量测,得到观测射线方程,并计算观测射线与跟踪门边界的交点,由此判决是否接受量测;
(5)根据得到的目标位置的相对概率分布函数和临界交点,得到跟踪门内量测方向上各点的概率积分,并进一步得到量测源自目标的概率;
(6)根据上步得到的关联概率,应用全期望公式,对目标状态和估计协方差矩阵进行更新;
(7)重复步骤(3)~(6),完成顺序结构的被动式多传感器目标跟踪,得到目标实时位置信息;
所述步骤(2)具体为,设第i个传感器的位置为[XNi,YNi]和该传感器的量测误差Ri,同时设先验预测的目标位置为和对应的误差协方差矩阵Ps,则椭圆跟踪门边界方程o1为: 其中γ0为设定的椭圆跟踪门门限,用如下公式计算其中β为回波密度即目标量测与杂波量测的比值,PD为检测概率;构造量化直线li,斜率k为ki=-1/tanθi,其中经过θi为传感器i预测的量测角,利用如下公式求得:
并经过点由此求得直线li在椭圆o1内的线段长度,记为d1,i,同时传感器i的量测误差的标准差取2倍标准差范围扫过直线li的部分,计算该线段长度:d2,i=ditan(2σi),其中di为传感器到目标预测位置的距离,通过公式计算如下:将各传感器的量测贡献程度用αi来量化,公式如下:αi=d1,i/d2,i,其中量化值越大越优先处理;所述步骤(3)中,根据先验的目标位置预测和误差协方差矩阵Ps,得到目标位置[X,Y]的相对概率分布:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201110172026 CN102338874B (zh) | 2011-06-24 | 2011-06-24 | 用于被动式多传感器目标跟踪的全局概率数据关联方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201110172026 CN102338874B (zh) | 2011-06-24 | 2011-06-24 | 用于被动式多传感器目标跟踪的全局概率数据关联方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102338874A CN102338874A (zh) | 2012-02-01 |
CN102338874B true CN102338874B (zh) | 2013-05-08 |
Family
ID=45514701
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN 201110172026 Expired - Fee Related CN102338874B (zh) | 2011-06-24 | 2011-06-24 | 用于被动式多传感器目标跟踪的全局概率数据关联方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102338874B (zh) |
Families Citing this family (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103017771B (zh) * | 2012-12-27 | 2015-06-17 | 杭州电子科技大学 | 一种静止传感器平台的多目标联合分配与跟踪方法 |
CN103096444B (zh) * | 2013-01-29 | 2016-08-10 | 浙江大学 | 一种基于传感器节点策略选择的水下无线传感器网络目标跟踪方法 |
CN103576137B (zh) * | 2013-09-27 | 2015-05-27 | 电子科技大学 | 一种基于成像策略的多传感器多目标定位方法 |
CN104155651A (zh) * | 2014-08-30 | 2014-11-19 | 电子科技大学 | 一种用于极化雷达目标跟踪的概率数据关联方法 |
CN104237845B (zh) * | 2014-09-19 | 2016-08-24 | 中国电子科技集团公司第二十八研究所 | 一种被动传感器目标位置估计误差均值方差分析方法 |
CN105137418B (zh) * | 2015-07-28 | 2017-06-16 | 中国人民解放军海军航空工程学院 | 基于全邻模糊聚类的多目标跟踪与数据互联方法 |
CN108871416A (zh) * | 2018-03-19 | 2018-11-23 | 西安电子科技大学 | 基于虚假点消除的角度冗余数据关联方法、被动定位系统 |
CN109298414B (zh) * | 2018-09-29 | 2020-11-20 | 中国科学院电子学研究所 | 雷达多运动目标实时跟踪方法 |
CN109901152B (zh) * | 2019-03-25 | 2023-04-11 | 西安电子科技大学 | 异步动平台多基地传感器系统目标检测跟踪联合估计方法 |
CN110542885B (zh) * | 2019-08-13 | 2021-09-21 | 北京理工大学 | 一种复杂交通环境下的毫米波雷达目标跟踪方法 |
CN111488812B (zh) * | 2020-04-01 | 2022-02-22 | 腾讯科技(深圳)有限公司 | 障碍物位置识别方法、装置、计算机设备和存储介质 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20050234679A1 (en) * | 2004-02-13 | 2005-10-20 | Evolution Robotics, Inc. | Sequential selective integration of sensor data |
CN101770024B (zh) * | 2010-01-25 | 2011-08-31 | 上海交通大学 | 多目标跟踪方法 |
CN101984362B (zh) * | 2010-08-12 | 2012-10-31 | 中国人民解放军海军航空工程学院 | 基于数据压缩的集中式多源广义相关跟踪方法 |
CN101975575B (zh) * | 2010-10-15 | 2012-03-28 | 西安电子科技大学 | 基于粒子滤波的被动传感器多目标跟踪方法 |
-
2011
- 2011-06-24 CN CN 201110172026 patent/CN102338874B/zh not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
CN102338874A (zh) | 2012-02-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102338874B (zh) | 用于被动式多传感器目标跟踪的全局概率数据关联方法 | |
CN104297748B (zh) | 一种基于轨迹增强的雷达目标检测前跟踪方法 | |
CN106405537B (zh) | 基于位置信息和多普勒信息的雷达航迹起始方法 | |
CN111007495A (zh) | 一种基于双重融合最大熵模糊聚类jpda的目标航迹优化方法 | |
CN110208806B (zh) | 一种航海雷达图像降雨识别方法 | |
CN103217673B (zh) | 一种非均匀Weibull杂波背景下的CFAR检测方法 | |
CN107436434B (zh) | 基于双向多普勒估计的航迹起始方法 | |
CN115685096B (zh) | 一种基于逻辑回归的二次雷达副瓣抑制方法 | |
CN107064865A (zh) | 基于深度聚类的极坐标动态规划无源协同定位方法 | |
CN103729462B (zh) | 一种基于稀疏表示处理遮挡的行人检索方法 | |
CN111830501B (zh) | Hrrp历史特征辅助的信号模糊数据关联方法及系统 | |
CN111259332B (zh) | 一种杂波环境下的模糊数据关联方法及多目标跟踪方法 | |
CN116542912A (zh) | 一种多目标视觉跟踪的柔性体桥梁振动检测模型、应用 | |
CN116859356B (zh) | 基于旋转矩阵优化求解的车载4d毫米波雷达自标定方法 | |
Bie et al. | Online method to impute missing loop detector data for urban freeway traffic control | |
CN112051567B (zh) | 人体目标微多普勒频率估计方法 | |
RU2562616C1 (ru) | Способ получения радиотехнической информации и радиотехнический комплекс для его осуществления | |
CN105652256B (zh) | 一种基于极化信息的高频地波雷达tbd方法 | |
CN105353353B (zh) | 多重搜索粒子概率假设密度滤波的多目标跟踪方法 | |
CN115220002B (zh) | 一种固定单站的多目标数据关联跟踪方法和相关装置 | |
Hamed et al. | An application of K-Nearest Neighbor interpolation on calibrating corrosion measurements collected by two non-destructive techniques | |
CN105046707A (zh) | 基于n阶多项式函数拟合海杂波的SAR图像船只检测方法 | |
CN115144847A (zh) | 针对海杂波中弱目标探测的参数空间多通道目标搜索技术 | |
Ge et al. | Multi-radar hybrid detection algorithm based on information entropy | |
CN108387879A (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 | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20130508 Termination date: 20140624 |
|
EXPY | Termination of patent right or utility model |