[go: up one dir, main page]

CN116049759A - Multi-sensor based monitoring target fusion method, device, and electronic equipment - Google Patents

Multi-sensor based monitoring target fusion method, device, and electronic equipment Download PDF

Info

Publication number
CN116049759A
CN116049759A CN202211649534.8A CN202211649534A CN116049759A CN 116049759 A CN116049759 A CN 116049759A CN 202211649534 A CN202211649534 A CN 202211649534A CN 116049759 A CN116049759 A CN 116049759A
Authority
CN
China
Prior art keywords
target
value
current
contour
sensor
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.)
Pending
Application number
CN202211649534.8A
Other languages
Chinese (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.)
Chitic Control Engineering Co ltd
Original Assignee
Chitic Control Engineering Co ltd
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 Chitic Control Engineering Co ltd filed Critical Chitic Control Engineering Co ltd
Priority to CN202211649534.8A priority Critical patent/CN116049759A/en
Publication of CN116049759A publication Critical patent/CN116049759A/en
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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
    • G01S13/00Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
    • G01S13/86Combinations of radar systems with non-radar systems, e.g. sonar, direction finder
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/41Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/18Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Data Mining & Analysis (AREA)
  • Mathematical Physics (AREA)
  • Mathematical Optimization (AREA)
  • Theoretical Computer Science (AREA)
  • Computational Mathematics (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Probability & Statistics with Applications (AREA)
  • Operations Research (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Evolutionary Biology (AREA)
  • Algebra (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

The invention discloses a monitoring target fusion method based on multiple sensors, a device and electronic equipment thereof, wherein the fusion method comprises the following steps: and (3) sampling the obtained state estimation values of the target states of the monitoring targets and the obtained contour estimation values of the target contours at the previous moment to obtain weighted particles of the predicted measurement values of the monitoring targets at the current moment, processing the measurement data of each sensor to obtain the current state estimation values and the current contour estimation values of all the monitoring targets corresponding to the sensors at the current moment, screening a target sensor set corresponding to each monitoring target, and fusing the current state estimation values and the current contour estimation values corresponding to each target sensor in the target sensor set to obtain the target state values and the target contour values of the monitoring targets at the current moment. The invention solves the technical problems of poor timeliness and low accuracy of calculating the state value and the contour value of the monitoring target in the related technology.

Description

基于多传感器的监测目标融合方法及其装置、电子设备Monitoring target fusion method based on multi-sensor and its device and electronic equipment

技术领域Technical Field

本发明涉及多传感器信息融合技术领域,具体而言,涉及一种基于多传感器的监测目标融合方法及其装置、电子设备。The present invention relates to the technical field of multi-sensor information fusion, and in particular to a multi-sensor based monitoring target fusion method and device, and electronic equipment.

背景技术Background Art

随着“智慧城市”的建设,越来越多种类的传感器被应用于各类场合。为了避免视觉传感器容易侵犯隐私的缺点,可以使用诸如毫米波雷达等的传感器,在避免侵犯隐私的同时能够进行快速准确的定位与跟踪。而又由于传感器的分辨率高,导致一个监测目标(例如,汽车、行人等)会对应产生多个量测,又因辨别监测目标的形状有利于进一步判断,因此考虑将其作为扩展目标进行处理。With the construction of "smart cities", more and more types of sensors are being used in various occasions. In order to avoid the disadvantage of visual sensors that easily infringe privacy, sensors such as millimeter wave radars can be used to quickly and accurately locate and track while avoiding privacy violations. However, due to the high resolution of the sensor, one monitoring target (for example, cars, pedestrians, etc.) will generate multiple measurements. Since identifying the shape of the monitoring target is conducive to further judgment, it is considered to be treated as an extended target.

相关技术中,对于扩展目标的定位跟踪,往往需要进行诸如聚类等的预处理操作,再进行跟踪,该方法对于计算力有限的嵌入式平台来说计算负担过大,难以满足实时性要求。In the related art, the positioning and tracking of extended targets often requires preprocessing operations such as clustering before tracking. This method has too much computational burden for embedded platforms with limited computing power and is difficult to meet real-time requirements.

针对上述的问题,目前尚未提出有效的解决方案。To address the above-mentioned problems, no effective solution has been proposed yet.

发明内容Summary of the invention

本发明实施例提供了一种基于多传感器的监测目标融合方法及其装置、电子设备,以至少解决相关技术中计算监测目标的状态值以及轮廓值的时效性较差以及准确率较低的技术问题。The embodiment of the present invention provides a monitoring target fusion method based on multiple sensors, a device thereof, and an electronic device, so as to at least solve the technical problems of poor timeliness and low accuracy in calculating the state value and profile value of the monitoring target in the related art.

根据本发明实施例的一个方面,提供了一种基于多传感器的监测目标融合方法,包括:获取在前一时刻下各监测目标的目标状态的状态估计值以及目标轮廓的轮廓估计值;对所述监测目标的所述状态估计值以及所述轮廓估计值进行粒子采样,得到在当前时刻下所述监测目标的预测量测值的带权粒子;基于在所述当前时刻下所述监测目标的预测量测值的带权粒子,对每个传感器的量测数据进行处理,得到在所述当前时刻下与所述传感器对应的所有所述监测目标的当前状态估计值以及当前轮廓估计值;基于所述当前状态估计值以及所述当前轮廓估计值,筛选与每个所述监测目标对应的目标传感器集合;融合所述目标传感器集合中与每个目标传感器对应的所述当前状态估计值以及所述当前轮廓估计值,得到在所述当前时刻下所述监测目标的目标状态值以及目标轮廓值。According to one aspect of an embodiment of the present invention, a monitoring target fusion method based on multiple sensors is provided, comprising: obtaining a state estimation value of a target state and a contour estimation value of a target contour of each monitoring target at a previous moment; performing particle sampling on the state estimation value and the contour estimation value of the monitoring target to obtain weighted particles of a predicted measurement value of the monitoring target at a current moment; processing measurement data of each sensor based on the weighted particles of the predicted measurement value of the monitoring target at the current moment to obtain a current state estimation value and a current contour estimation value of all the monitoring targets corresponding to the sensor at the current moment; screening a set of target sensors corresponding to each monitoring target based on the current state estimation value and the current contour estimation value; fusing the current state estimation value and the current contour estimation value corresponding to each target sensor in the target sensor set to obtain a target state value and a target contour value of the monitoring target at the current moment.

可选地,在获取在前一时刻下各监测目标的目标状态的状态估计值以及目标轮廓的轮廓估计值之后,还包括:基于在所述前一时刻下所述监测目标的所述状态估计值,确定所述监测目标在所述当前时刻下的状态预测值;基于在所述前一时刻下所述监测目标的所述轮廓估计值,确定所述监测目标在所述当前时刻下的轮廓预测值;基于所述状态预测值以及所述轮廓预测值,得到在所述当前时刻下所述监测目标的预测量测值,其中,所述预测量测值是所述监测目标在所述当前时刻下的预测状态所对应的量测值。Optionally, after obtaining the state estimation value of the target state of each monitoring target at the previous moment and the contour estimation value of the target contour, it also includes: determining the state prediction value of the monitoring target at the current moment based on the state estimation value of the monitoring target at the previous moment; determining the contour prediction value of the monitoring target at the current moment based on the contour estimation value of the monitoring target at the previous moment; and obtaining the predicted measurement value of the monitoring target at the current moment based on the state prediction value and the contour prediction value, wherein the predicted measurement value is the measurement value corresponding to the predicted state of the monitoring target at the current moment.

可选地,对所述监测目标的所述状态估计值以及所述轮廓估计值进行粒子采样,得到在当前时刻下所述监测目标的预测量测值的带权粒子的步骤,包括:对所述状态估计值以及所述轮廓估计值进行粒子采样,生成多个粒子;基于预设参数,计算每个所述粒子的初始粒子权重;基于所有所述初始粒子权重,对各所述初始粒子权重进行归一化处理,得到每个粒子的目标粒子权重,其中,将携带有所述目标粒子权重的所述粒子表征为带权粒子;更新所述带权粒子,得到所述监测目标在所述当前时刻下关联所述预测量测值的带权粒子。Optionally, the step of performing particle sampling on the state estimation value and the profile estimation value of the monitoring target to obtain weighted particles of the predicted measurement value of the monitoring target at the current moment includes: performing particle sampling on the state estimation value and the profile estimation value to generate multiple particles; calculating the initial particle weight of each particle based on preset parameters; normalizing each initial particle weight based on all the initial particle weights to obtain a target particle weight for each particle, wherein the particles carrying the target particle weights are characterized as weighted particles; updating the weighted particles to obtain weighted particles associated with the predicted measurement value of the monitoring target at the current moment.

可选地,基于在所述当前时刻下所述监测目标的预测量测值的带权粒子,对每个传感器的量测数据进行处理之前,还包括:获取每个所述传感器的所述量测数据,其中,所述量测数据至少包括:量测向量集合、量测总数;基于所述量测数据,构建面向监测目标的目标关联变量,其中,所述目标关联变量用于表示预设监测目标在所述当前时刻下产生的第一预设位数的量测指示的量测索引;基于所述量测数据,构建面向量测的量测关联变量,其中,所述量测关联变量用于表示在所述当前时刻下由所述量测索引所映射的监测目标产生的第二预设位数的量测。Optionally, based on the weighted particles of the predicted measurement value of the monitoring target at the current moment, before processing the measurement data of each sensor, it also includes: obtaining the measurement data of each of the sensors, wherein the measurement data at least includes: a measurement vector set and a total number of measurements; based on the measurement data, constructing a target-associated variable for the monitoring target, wherein the target-associated variable is used to represent a measurement index indicating a first preset number of bits of measurement generated by a preset monitoring target at the current moment; based on the measurement data, constructing a measurement-associated variable for vector measurement, wherein the measurement-associated variable is used to represent a second preset number of bits of measurement generated by the monitoring target mapped by the measurement index at the current moment.

可选地,基于在所述当前时刻下所述监测目标的预测量测值的带权粒子,对每个传感器的量测数据进行处理,得到在所述当前时刻下与所述传感器对应的所有所述监测目标的当前状态估计值以及当前轮廓估计值的步骤,包括:基于在所述当前时刻下所述监测目标的预测量测值的带权粒子以及所述量测向量集合,对所述目标关联变量进行量测评估,得到似然概率;在所述似然概率属于预设概率阈值范围的情况轮廓估计值下,基于所述目标关联变量以及所述量测关联变量,构建关联因子图;基于所述关联因子图,进行迭代计算,得到第一参数值以及第二参数值;基于所述第一参数值以及所述第二参数值,更新所述传感器的所述监测目标的量测数据,得到更新结果;基于所述更新结果,得到与所述传感器对应的在所述当前时刻下所述监测目标的当前轮廓估计值,并基于所述更新结果,对所述带权粒子进行加权计算,得到与所述传感器对应的所述监测目标的当前状态估计值。Optionally, based on the weighted particles of the predicted measurement values of the monitored target at the current moment, the measurement data of each sensor is processed to obtain the current state estimation values and current profile estimation values of all the monitored targets corresponding to the sensor at the current moment, including: based on the weighted particles of the predicted measurement values of the monitored target at the current moment and the measurement vector set, the target-associated variables are measured and evaluated to obtain a likelihood probability; under the profile estimation value when the likelihood probability belongs to a preset probability threshold range, an association factor graph is constructed based on the target-associated variables and the measurement-associated variables; based on the association factor graph, an iterative calculation is performed to obtain a first parameter value and a second parameter value; based on the first parameter value and the second parameter value, the measurement data of the monitored target of the sensor is updated to obtain an updated result; based on the updated result, the current profile estimation value of the monitored target corresponding to the sensor at the current moment is obtained, and based on the updated result, the weighted particles are weighted calculated to obtain the current state estimation value of the monitored target corresponding to the sensor.

可选地,基于所述当前状态估计值以及所述当前轮廓估计值,筛选与每个所述监测目标对应的目标传感器集合的步骤,包括:将所述当前状态估计值表征为预设高斯分布的均值,并将所述当前轮廓估计值表征为所述预设高斯分布的协方差矩阵;基于所述均值以及所述协方差矩阵,计算每个所述传感器对所述监测目标的估计概率;对所有所述估计概率进行归一化处理,得到目标估计概率;在所述目标估计概率大于预设概率阈值的情况下,将与所述目标估计概率对应的所述传感器加入至所述监测目标对应的所述目标传感器集合。Optionally, based on the current state estimation value and the current profile estimation value, the step of screening the target sensor set corresponding to each of the monitored targets includes: characterizing the current state estimation value as the mean of a preset Gaussian distribution, and characterizing the current profile estimation value as the covariance matrix of the preset Gaussian distribution; based on the mean and the covariance matrix, calculating the estimated probability of each sensor for the monitored target; normalizing all the estimated probabilities to obtain a target estimated probability; and when the target estimated probability is greater than a preset probability threshold, adding the sensor corresponding to the target estimated probability to the target sensor set corresponding to the monitored target.

可选地,融合所述目标传感器集合中与每个目标传感器对应的所述当前状态估计值以及所述当前轮廓估计值,得到在所述当前时刻下所述监测目标的目标状态值以及目标轮廓值的步骤,包括:计算所有与每个所述目标传感器对应的所述当前状态估计值的平均值,得到所述目标状态值;基于所述当前轮廓估计值,确定所述监测目标的当前目标轮廓的第一半轴长值、第二半轴长值以及偏转角度值;基于所有所述目标传感器的所述协方差矩阵,确定目标协方差矩阵;基于所述目标协方差矩阵,确定目标偏转角度值;基于所述目标偏转角度值,计算每个所述目标传感器的目标第一半轴长值以及目标第二半轴长值;基于所有所述目标第一半轴长值以及所述目标第二半轴长值,得到目标半轴长值;基于所述目标半轴长值以及所述目标偏转角度值,得到所述目标轮廓值。Optionally, the step of fusing the current state estimation value and the current contour estimation value corresponding to each target sensor in the target sensor set to obtain the target state value and the target contour value of the monitored target at the current moment includes: calculating the average value of all the current state estimation values corresponding to each of the target sensors to obtain the target state value; determining the first semi-axis length value, the second semi-axis length value and the deflection angle value of the current target contour of the monitored target based on the current contour estimation value; determining the target covariance matrix based on the covariance matrix of all the target sensors; determining the target deflection angle value based on the target covariance matrix; calculating the target first semi-axis length value and the target second semi-axis length value of each target sensor based on the target deflection angle value; obtaining the target semi-axis length value based on all the target first semi-axis length values and the target second semi-axis length values; obtaining the target contour value based on the target semi-axis length value and the target deflection angle value.

根据本发明实施例的另一方面,还提供了一种基于多传感器的监测目标融合装置,包括:获取单元,用于获取在前一时刻下各监测目标的目标状态的状态估计值以及目标轮廓的轮廓估计值;采样单元,用于对所述监测目标的所述状态估计值以及所述轮廓估计值进行粒子采样,得到在当前时刻下所述监测目标的预测量测值的带权粒子;处理单元,用于基于在所述当前时刻下所述监测目标的预测量测值的带权粒子,对每个传感器的量测数据进行处理,得到在所述当前时刻下与所述传感器对应的所有所述监测目标的当前状态估计值以及当前轮廓估计值;筛选单元,用于基于所述当前状态估计值以及所述当前轮廓估计值,筛选与每个所述监测目标对应的目标传感器集合;融合单元,用于融合所述目标传感器集合中与每个目标传感器对应的所述当前状态估计值以及所述当前轮廓估计值,得到在所述当前时刻下所述监测目标的目标状态值以及目标轮廓值。According to another aspect of an embodiment of the present invention, a monitoring target fusion device based on multiple sensors is also provided, including: an acquisition unit, used to acquire a state estimation value of a target state of each monitoring target and a contour estimation value of a target contour at a previous moment; a sampling unit, used to perform particle sampling on the state estimation value and the contour estimation value of the monitoring target to obtain weighted particles of a predicted measurement value of the monitoring target at a current moment; a processing unit, used to process the measurement data of each sensor based on the weighted particles of the predicted measurement value of the monitoring target at the current moment to obtain a current state estimation value and a current contour estimation value of all the monitoring targets corresponding to the sensor at the current moment; a screening unit, used to screen a set of target sensors corresponding to each monitoring target based on the current state estimation value and the current contour estimation value; and a fusion unit, used to fuse the current state estimation value and the current contour estimation value corresponding to each target sensor in the target sensor set to obtain a target state value and a target contour value of the monitoring target at the current moment.

可选地,所述融合装置还包括:第一确定模块,用于在获取在前一时刻下各监测目标的目标状态的状态估计值以及目标轮廓的轮廓估计值之后,基于在所述前一时刻下所述监测目标的所述状态估计值,确定所述监测目标在所述当前时刻下的状态预测值;第二确定模块,用于基于在所述前一时刻下所述监测目标的所述轮廓估计值,确定所述监测目标在所述当前时刻下的轮廓预测值;第一输出模块,用于基于所述状态预测值以及所述轮廓预测值,得到在所述当前时刻下所述监测目标的预测量测值,其中,所述预测量测值是所述监测目标在所述当前时刻下的预测状态所对应的量测值。Optionally, the fusion device also includes: a first determination module, which is used to determine the state prediction value of the monitoring target at the current moment based on the state estimation value of the monitoring target at the previous moment after obtaining the state estimation value of the target state of each monitoring target and the contour estimation value of the target contour at the previous moment; a second determination module, which is used to determine the contour prediction value of the monitoring target at the current moment based on the contour estimation value of the monitoring target at the previous moment; and a first output module, which is used to obtain the predicted measurement value of the monitoring target at the current moment based on the state prediction value and the contour prediction value, wherein the predicted measurement value is the measurement value corresponding to the predicted state of the monitoring target at the current moment.

可选地,所述采样单元包括:第一生成模块,用于对所述状态估计值以及所述轮廓估计值进行粒子采样,生成多个粒子;第一计算模块,用于基于预设参数,计算每个所述粒子的初始粒子权重;第一处理模块,用于基于所有所述初始粒子权重,对各所述初始粒子权重进行归一化处理,得到每个粒子的目标粒子权重,其中,将携带有所述目标粒子权重的所述粒子表征为带权粒子;第一更新模块,用于更新所述带权粒子,得到所述监测目标在所述当前时刻下关联所述预测量测值的带权粒子。Optionally, the sampling unit includes: a first generation module, used to perform particle sampling on the state estimation value and the profile estimation value to generate multiple particles; a first calculation module, used to calculate the initial particle weight of each particle based on preset parameters; a first processing module, used to normalize each initial particle weight based on all the initial particle weights to obtain the target particle weight of each particle, wherein the particle carrying the target particle weight is characterized as a weighted particle; a first update module, used to update the weighted particles to obtain the weighted particles of the monitored target associated with the predicted measurement value at the current moment.

可选地,所述融合装置还包括:第一获取模块,用于基于在所述当前时刻下所述监测目标的预测量测值的带权粒子,对每个传感器的量测数据进行处理之前,获取每个所述传感器的所述量测数据,其中,所述量测数据至少包括:量测向量集合、量测总数;第一构建模块,用于基于所述量测数据,构建面向监测目标的目标关联变量,其中,所述目标关联变量用于表示预设监测目标在所述当前时刻下产生的第一预设位数的量测指示的量测索引;第二构建模块,用于基于所述量测数据,构建面向量测的量测关联变量,其中,所述量测关联变量用于表示在所述当前时刻下由所述量测索引所映射的监测目标产生的第二预设位数的量测。Optionally, the fusion device also includes: a first acquisition module, used to acquire the measurement data of each sensor before processing the measurement data of each sensor based on the weighted particles of the predicted measurement value of the monitoring target at the current moment, wherein the measurement data at least includes: a measurement vector set and a total number of measurements; a first construction module, used to construct a target-associated variable for the monitoring target based on the measurement data, wherein the target-associated variable is used to represent a measurement index indicating a first preset number of bits of measurement generated by a preset monitoring target at the current moment; a second construction module, used to construct a measurement-associated variable for vector measurement based on the measurement data, wherein the measurement-associated variable is used to represent a second preset number of bits of measurement generated by the monitoring target mapped by the measurement index at the current moment.

可选地,所述处理单元包括:第一评估模块,用于基于在所述当前时刻下所述监测目标的预测量测值的带权粒子以及所述量测向量集合,对所述目标关联变量进行量测评估,得到似然概率;第三构建模块,用于在所述似然概率属于预设概率阈值范围的情况轮廓估计值下,基于所述目标关联变量以及所述量测关联变量,构建关联因子图;第二计算模块,用于基于所述关联因子图,进行迭代计算,得到第一参数值以及第二参数值;第二更新模块,用于基于所述第一参数值以及所述第二参数值,更新所述传感器的所述监测目标的量测数据,得到更新结果;第三计算模块,用于基于所述更新结果,得到与所述传感器对应的在所述当前时刻下所述监测目标的当前轮廓估计值,并基于所述更新结果,对所述带权粒子进行加权计算,得到与所述传感器对应的所述监测目标的当前状态估计值。Optionally, the processing unit includes: a first evaluation module, which is used to perform measurement evaluation on the target-associated variables based on the weighted particles of the predicted measurement values of the monitored target at the current moment and the measurement vector set to obtain a likelihood probability; a third construction module, which is used to construct a correlation factor graph based on the target-associated variables and the measurement-associated variables when the likelihood probability belongs to a profile estimation value within a preset probability threshold range; a second calculation module, which is used to perform iterative calculation based on the correlation factor graph to obtain a first parameter value and a second parameter value; a second update module, which is used to update the measurement data of the monitored target of the sensor based on the first parameter value and the second parameter value to obtain an updated result; and a third calculation module, which is used to obtain a current profile estimation value of the monitored target corresponding to the sensor at the current moment based on the updated result, and to perform weighted calculation on the weighted particles based on the updated result to obtain a current state estimation value of the monitored target corresponding to the sensor.

可选地,所述筛选单元包括:第一表征模块,用于将所述当前状态估计值表征为预设高斯分布的均值,并将所述当前轮廓估计值表征为所述预设高斯分布的协方差矩阵;第四计算模块,用于基于所述均值以及所述协方差矩阵,计算每个所述传感器对所述监测目标的估计概率;第二处理模块,用于对所有所述估计概率进行归一化处理,得到目标估计概率;第一加入模块,用于在所述目标估计概率大于预设概率阈值的情况下,将与所述目标估计概率对应的所述传感器加入至所述监测目标对应的所述目标传感器集合。Optionally, the screening unit includes: a first characterization module, used to characterize the current state estimation value as the mean of a preset Gaussian distribution, and to characterize the current profile estimation value as the covariance matrix of the preset Gaussian distribution; a fourth calculation module, used to calculate the estimated probability of each sensor for the monitoring target based on the mean and the covariance matrix; a second processing module, used to normalize all the estimated probabilities to obtain the target estimated probability; and a first adding module, used to add the sensor corresponding to the target estimated probability to the target sensor set corresponding to the monitoring target when the target estimated probability is greater than a preset probability threshold.

可选地,所述融合单元包括:第五计算模块,用于计算所有与每个所述目标传感器对应的所述当前状态估计值的平均值,得到所述目标状态值;第三确定模块,用于基于所述当前轮廓估计值,确定所述监测目标的当前目标轮廓的第一半轴长值、第二半轴长值以及偏转角度值;第四确定模块,用于基于所有所述目标传感器的所述协方差矩阵,确定目标协方差矩阵;第五确定模块,用于基于所述目标协方差矩阵,确定目标偏转角度值;第六计算模块,用于基于所述目标偏转角度值,计算每个所述目标传感器的目标第一半轴长值以及目标第二半轴长值;第二输出模块,用于基于所有所述目标第一半轴长值以及所述目标第二半轴长值,得到目标半轴长值;第三输出模块,用于基于所述目标半轴长值以及所述目标偏转角度值,得到所述目标轮廓值。Optionally, the fusion unit includes: a fifth calculation module, used to calculate the average value of all the current state estimation values corresponding to each of the target sensors to obtain the target state value; a third determination module, used to determine the first semi-axis length value, the second semi-axis length value and the deflection angle value of the current target contour of the monitored target based on the current contour estimation value; a fourth determination module, used to determine the target covariance matrix based on the covariance matrix of all the target sensors; a fifth determination module, used to determine the target deflection angle value based on the target covariance matrix; a sixth calculation module, used to calculate the target first semi-axis length value and the target second semi-axis length value of each of the target sensors based on the target deflection angle value; a second output module, used to obtain the target semi-axis length value based on all the target first semi-axis length values and the target second semi-axis length values; a third output module, used to obtain the target contour value based on the target semi-axis length value and the target deflection angle value.

根据本发明实施例的另一方面,还提供了一种计算机可读存储介质,所述计算机可读存储介质包括存储的计算机程序,其中,在所述计算机程序运行时控制所述计算机可读存储介质所在设备执行上述所述的基于多传感器的监测目标融合方法。According to another aspect of an embodiment of the present invention, a computer-readable storage medium is also provided, wherein the computer-readable storage medium includes a stored computer program, wherein when the computer program is running, the device where the computer-readable storage medium is located is controlled to execute the above-mentioned multi-sensor based monitoring target fusion method.

根据本发明实施例的另一方面,还提供了一种电子设备,包括一个或多个处理器和存储器,所述存储器用于存储一个或多个程序,其中,当所述一个或多个程序被所述一个或多个处理器执行时,使得所述一个或多个处理器实现上述所述的基于多传感器的监测目标融合方法。According to another aspect of an embodiment of the present invention, there is also provided an electronic device, comprising one or more processors and a memory, wherein the memory is used to store one or more programs, wherein when the one or more programs are executed by the one or more processors, the one or more processors implement the above-mentioned multi-sensor based monitoring target fusion method.

在本公开中,获取在前一时刻下各监测目标的目标状态的状态估计值以及目标轮廓的轮廓估计值,对监测目标的状态估计值以及轮廓估计值进行粒子采样,得到在当前时刻下监测目标的预测量测值的带权粒子,基于在当前时刻下监测目标的预测量测值的带权粒子,对每个传感器的量测数据进行处理,得到在当前时刻下与传感器对应的所有监测目标的当前状态估计值以及当前轮廓估计值,基于当前状态估计值以及当前轮廓估计值,筛选与每个监测目标对应的目标传感器集合,融合目标传感器集合中与每个目标传感器对应的当前状态估计值以及当前轮廓估计值,得到在当前时刻下监测目标的目标状态值以及目标轮廓值。在本公开中,可以先对前一时刻下各监测目标的目标状态的状态估计值以及目标轮廓的轮廓估计值进行粒子采样,以得到当前时刻下监测目标的预测量测值的带权粒子,然后根据带权粒子对每个传感器的量测数据进行处理,以得到当前时刻下每个传感器对应的所有监测目标的当前状态估计值以及当前轮廓估计值,再筛选每个监测目标的目标传感器集合,之后融合每个目标传感器对应的当前状态估计值以及当前轮廓估计值,以得到当前时刻下监测目标的目标状态值以及目标轮廓值,能够实时计算监测目标的状态值以及轮廓值,并且利用筛选结果能够高效得到各监测目标的状态值以及轮廓值的准确估计,进而解决了相关技术中计算监测目标的状态值以及轮廓值的时效性较差以及准确率较低的技术问题。In the present disclosure, a state estimation value of a target state and a contour estimation value of a target contour of each monitored target at a previous moment are obtained, particle sampling is performed on the state estimation value and the contour estimation value of the monitored target to obtain weighted particles of predicted measurement values of the monitored target at the current moment, and based on the weighted particles of predicted measurement values of the monitored target at the current moment, measurement data of each sensor is processed to obtain current state estimation values and current contour estimation values of all monitored targets corresponding to the sensor at the current moment, and based on the current state estimation values and the current contour estimation values, a set of target sensors corresponding to each monitored target is screened, and the current state estimation values and the current contour estimation values corresponding to each target sensor in the target sensor set are fused to obtain the target state value and the target contour value of the monitored target at the current moment. In the present disclosure, particle sampling can be first performed on the state estimation value of the target state and the contour estimation value of the target contour of each monitored target at the previous moment to obtain weighted particles of the predicted measurement value of the monitored target at the current moment, and then the measurement data of each sensor is processed according to the weighted particles to obtain the current state estimation value and current contour estimation value of all monitored targets corresponding to each sensor at the current moment, and then the target sensor set of each monitored target is screened, and then the current state estimation value and current contour estimation value corresponding to each target sensor are fused to obtain the target state value and target contour value of the monitored target at the current moment, which can calculate the state value and contour value of the monitored target in real time, and the screening results can be used to efficiently obtain accurate estimates of the state value and contour value of each monitored target, thereby solving the technical problems of poor timeliness and low accuracy in calculating the state value and contour value of the monitored target in the related art.

附图说明BRIEF DESCRIPTION OF THE DRAWINGS

此处所说明的附图用来提供对本发明的进一步理解,构成本申请的一部分,本发明的示意性实施例及其说明用于解释本发明,并不构成对本发明的不当限定。在附图中:The drawings described herein are used to provide a further understanding of the present invention and constitute a part of this application. The exemplary embodiments of the present invention and their descriptions are used to explain the present invention and do not constitute an improper limitation of the present invention. In the drawings:

图1是根据本发明实施例的一种可选的基于多传感器的监测目标融合方法的流程图;FIG1 is a flow chart of an optional multi-sensor based monitoring target fusion method according to an embodiment of the present invention;

图2是根据本发明实施例的一种可选的面向生命线安全的多传感器多扩展目标融合算法原理的流程图;2 is a flow chart of an optional multi-sensor multi-extended target fusion algorithm principle for lifeline safety according to an embodiment of the present invention;

图3是根据本发明实施例的一种可选的基于多传感器的监测目标融合装置的示意图;FIG3 is a schematic diagram of an optional multi-sensor based monitoring target fusion device according to an embodiment of the present invention;

图4是根据本发明实施例的一种用于基于多传感器的监测目标融合方法的电子设备(或移动设备)的硬件结构框图。FIG4 is a hardware structure block diagram of an electronic device (or mobile device) for a multi-sensor based monitoring target fusion method according to an embodiment of the present invention.

具体实施方式DETAILED DESCRIPTION

为了使本技术领域的人员更好地理解本发明方案,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分的实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都应当属于本发明保护的范围。In order to enable those skilled in the art to better understand the scheme of the present invention, the technical scheme in the embodiments of the present invention will be clearly and completely described below in conjunction with the drawings in the embodiments of the present invention. Obviously, the described embodiments are only part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without creative work should fall within the scope of protection of the present invention.

需要说明的是,本发明的说明书和权利要求书及上述附图中的术语“第一”、“第二”等是用于区别类似的对象,而不必用于描述特定的顺序或先后次序。应该理解这样使用的数据在适当情况下可以互换,以便这里描述的本发明的实施例能够以除了在这里图示或描述的那些以外的顺序实施。此外,术语“包括”和“具有”以及他们的任何变形,意图在于覆盖不排他的包含,例如,包含了一系列步骤或单元的过程、方法、系统、产品或设备不必限于清楚地列出的那些步骤或单元,而是可包括没有清楚地列出的或对于这些过程、方法、产品或设备固有的其它步骤或单元。It should be noted that the terms "first", "second", etc. in the specification and claims of the present invention and the above-mentioned drawings are used to distinguish similar objects, and are not necessarily used to describe a specific order or sequence. It should be understood that the data used in this way can be interchanged where appropriate, so that the embodiments of the present invention described herein can be implemented in an order other than those illustrated or described herein. In addition, the terms "including" and "having" and any variations thereof are intended to cover non-exclusive inclusions, for example, a process, method, system, product or device that includes a series of steps or units is not necessarily limited to those steps or units clearly listed, but may include other steps or units that are not clearly listed or inherent to these processes, methods, products or devices.

需要说明的是,本公开所涉及的相关信息(包括但不限于用户设备信息、用户个人信息等)和数据(包括但不限于用于展示的数据、分析的数据等),均为经用户授权或者经过各方充分授权的信息和数据。例如,本系统和相关用户或机构间设置有接口,在获取相关信息之前,需要通过接口向前述的用户或机构发送获取请求,并在接收到前述的用户或机构反馈的同意信息后,获取相关信息。It should be noted that the relevant information (including but not limited to user device information, user personal information, etc.) and data (including but not limited to data for display, data for analysis, etc.) involved in this disclosure are all information and data authorized by the user or fully authorized by all parties. For example, an interface is set between this system and the relevant user or organization. Before obtaining relevant information, it is necessary to send an acquisition request to the aforementioned user or organization through the interface, and obtain relevant information after receiving the consent information fed back by the aforementioned user or organization.

本发明提出了一种面向生命线安全(如城市生命线安全)的多传感器多扩展目标融合算法(A mult i-sensor and mult i-extended target fus ion algorithm forurban l ifel ine safety),该算法基于因子图与和积算法,并使用投票机制优选传感器估计结果,实现多传感器多扩展目标的快速融合。可以先对上一时刻的各监测目标的状态以及轮廓估计进行粒子采样,然后将对这些带权粒子进行更新操作,以得到当前时刻状态的预测值。再由各传感器并行进行迭代数据关联,使得各自得到各目标状态以及轮廓的估计。之后通过投票机制优选估计结果,最后进行融合,以得到各目标状态以及轮廓的估计。该算法具有低计算复杂度、高计算速度的优点,同时,利用优选估计结果可以实现快速融合,能够高效得到各目标状态以及轮廓的准确估计。The present invention proposes a multi-sensor and multi-extended target fusion algorithm for urban lifeline safety, which is based on factor graph and sum-product algorithm, and uses voting mechanism to optimize sensor estimation results, so as to realize fast fusion of multi-sensor multi-extended targets. First, particle sampling can be performed on the state and contour estimation of each monitoring target at the previous moment, and then these weighted particles can be updated to obtain the predicted value of the state at the current moment. Then, each sensor performs iterative data association in parallel, so that each obtains an estimate of the state and contour of each target. After that, the estimation results are optimized through voting mechanism, and finally fused to obtain an estimate of the state and contour of each target. The algorithm has the advantages of low computational complexity and high computational speed. At the same time, fast fusion can be realized by using the optimized estimation results, and accurate estimates of the state and contour of each target can be obtained efficiently.

下面结合各个实施例来详细说明本发明。The present invention is described in detail below in conjunction with various embodiments.

实施例一Embodiment 1

根据本发明实施例,提供了一种基于多传感器的监测目标融合方法的实施例,需要说明的是,在附图的流程图示出的步骤可以在诸如一组计算机可执行指令的计算机系统中执行,并且,虽然在流程图中示出了逻辑顺序,但是在某些情况下,可以以不同于此处的顺序执行所示出或描述的步骤。According to an embodiment of the present invention, an embodiment of a monitoring target fusion method based on multiple sensors is provided. It should be noted that the steps shown in the flowchart of the accompanying drawings can be executed in a computer system such as a set of computer executable instructions, and although a logical order is shown in the flowchart, in some cases, the steps shown or described can be executed in an order different from that shown here.

图1是根据本发明实施例的一种可选的基于多传感器的监测目标融合方法的流程图,如图1所示,该方法包括如下步骤:FIG. 1 is a flow chart of an optional multi-sensor based monitoring target fusion method according to an embodiment of the present invention. As shown in FIG. 1 , the method includes the following steps:

步骤S101,获取在前一时刻下各监测目标的目标状态的状态估计值以及目标轮廓的轮廓估计值。Step S101, obtaining a state estimation value of a target state of each monitored target and a contour estimation value of a target contour at a previous moment.

步骤S102,对监测目标的状态估计值以及轮廓估计值进行粒子采样,得到在当前时刻下监测目标的预测量测值的带权粒子。Step S102, performing particle sampling on the state estimation value and the profile estimation value of the monitoring target to obtain weighted particles of the predicted measurement value of the monitoring target at the current moment.

步骤S103,基于在当前时刻下监测目标的预测量测值的带权粒子,对每个传感器的量测数据进行处理,得到在当前时刻下与传感器对应的所有监测目标的当前状态估计值以及当前轮廓估计值。Step S103, based on the weighted particles of the predicted measurement values of the monitoring targets at the current moment, the measurement data of each sensor is processed to obtain the current state estimation values and the current contour estimation values of all monitoring targets corresponding to the sensor at the current moment.

步骤S104,基于当前状态估计值以及当前轮廓估计值,筛选与每个监测目标对应的目标传感器集合。Step S104: based on the current state estimation value and the current profile estimation value, filter the target sensor set corresponding to each monitoring target.

步骤S105,融合目标传感器集合中与每个目标传感器对应的当前状态估计值以及当前轮廓估计值,得到在当前时刻下监测目标的目标状态值以及目标轮廓值。Step S105 , fusing the current state estimation value and the current contour estimation value corresponding to each target sensor in the target sensor set to obtain the target state value and the target contour value of the monitored target at the current moment.

通过上述步骤,可以获取在前一时刻下各监测目标的目标状态的状态估计值以及目标轮廓的轮廓估计值,对监测目标的状态估计值以及轮廓估计值进行粒子采样,得到在当前时刻下监测目标的预测量测值的带权粒子,基于在当前时刻下监测目标的预测量测值的带权粒子,对每个传感器的量测数据进行处理,得到在当前时刻下与传感器对应的所有监测目标的当前状态估计值以及当前轮廓估计值,基于当前状态估计值以及当前轮廓估计值,筛选与每个监测目标对应的目标传感器集合,融合目标传感器集合中与每个目标传感器对应的当前状态估计值以及当前轮廓估计值,得到在当前时刻下监测目标的目标状态值以及目标轮廓值。在本发明实施例中,可以先对前一时刻下各监测目标的目标状态的状态估计值以及目标轮廓的轮廓估计值进行粒子采样,以得到当前时刻下监测目标的预测量测值的带权粒子,然后根据带权粒子对每个传感器的量测数据进行处理,以得到当前时刻下每个传感器对应的所有监测目标的当前状态估计值以及当前轮廓估计值,再筛选每个监测目标的目标传感器集合,之后融合每个目标传感器对应的当前状态估计值以及当前轮廓估计值,以得到当前时刻下监测目标的目标状态值以及目标轮廓值,能够实时计算监测目标的状态值以及轮廓值,并且利用筛选结果能够高效得到各监测目标的状态值以及轮廓值的准确估计,进而解决了相关技术中计算监测目标的状态值以及轮廓值的时效性较差以及准确率较低的技术问题。Through the above steps, the state estimation value of the target state and the contour estimation value of the target contour of each monitored target at the previous moment can be obtained, and the state estimation value and the contour estimation value of the monitored target can be particle sampled to obtain the weighted particles of the predicted measurement value of the monitored target at the current moment. Based on the weighted particles of the predicted measurement value of the monitored target at the current moment, the measurement data of each sensor is processed to obtain the current state estimation value and the current contour estimation value of all monitored targets corresponding to the sensor at the current moment. Based on the current state estimation value and the current contour estimation value, the target sensor set corresponding to each monitored target is screened, and the current state estimation value and the current contour estimation value corresponding to each target sensor in the target sensor set are fused to obtain the target state value and the target contour value of the monitored target at the current moment. In an embodiment of the present invention, particle sampling can be first performed on the state estimation value of the target state and the contour estimation value of the target contour of each monitored target at the previous moment to obtain weighted particles of the predicted measurement value of the monitored target at the current moment, and then the measurement data of each sensor is processed according to the weighted particles to obtain the current state estimation value and current contour estimation value of all monitored targets corresponding to each sensor at the current moment, and then the target sensor set of each monitored target is screened, and then the current state estimation value and current contour estimation value corresponding to each target sensor are fused to obtain the target state value and target contour value of the monitored target at the current moment, so that the state value and contour value of the monitored target can be calculated in real time, and the screening result can be used to efficiently obtain accurate estimates of the state value and contour value of each monitored target, thereby solving the technical problems of poor timeliness and low accuracy in calculating the state value and contour value of the monitored target in the related art.

下面结合上述各步骤对本发明实施例进行详细说明。The embodiment of the present invention is described in detail below in combination with the above steps.

在本实施例中,假设使用S个传感器对K个监测目标进行监测,其中,S,K是固定且已知的。使用xk(t)表示监测目标k∈{1,2,…,K}在t时刻的状态。In this embodiment, it is assumed that S sensors are used to monitor K monitoring targets, where S and K are fixed and known. x k (t) is used to represent the state of the monitoring target k∈{1,2,…,K} at time t.

首先,设置监测目标的模型为:First, set the model of the monitoring target as:

xk(t)=Fxk(t-1)+vk(t-1)(1);x k (t)=Fx k (t-1)+v k (t-1)(1);

其中,vk(t-1)是服从零均值的高斯白噪声,即vk(t-1)~N(0,Qk(t-1)),Qk(t-1)表示协方差,F为状态转移矩阵,xk(t-1)表示监测目标在t-1时刻的状态。Among them, vk (t-1) is Gaussian white noise with zero mean, that is, vk (t-1)~N(0, Qk (t-1)), Qk (t-1) represents the covariance, F is the state transfer matrix, and xk (t-1) represents the state of the monitored target at time t-1.

传感器s∈{1,2,…,S}在t时刻产生

Figure BDA0004011310020000081
个量测,将第
Figure BDA0004011310020000082
个测量记为
Figure BDA0004011310020000091
定义量测向量
Figure BDA0004011310020000092
Sensor s∈{1,2,…,S} generates
Figure BDA0004011310020000081
A measurement,
Figure BDA0004011310020000082
The measurement is recorded as
Figure BDA0004011310020000091
Defining the measurement vector
Figure BDA0004011310020000092

对目标轮廓进行描述,目标轮廓可以采用长短半轴长以及偏转角度(l1,l2,θ)进行描述,本实施例中假设目标轮廓的偏转角度与目标运动方向一致。The target contour can be described by using the major and minor semi-axis lengths and the deflection angle (l 1 , l 2 , θ). In this embodiment, it is assumed that the deflection angle of the target contour is consistent with the target motion direction.

设置目标在传感器s处的观测模型为:The observation model of the target at sensor s is set as:

Figure BDA0004011310020000093
Figure BDA0004011310020000093

其中,

Figure BDA0004011310020000094
表示t时刻目标k在传感器s上产生的某个量测,μk(t)是服从零均值的高斯白噪声,即μk(t)~N(0,Rk(t)),Rk(t)表示协方差,Xk(t)表示监测目标在t时刻的轮廓,h(·)是一个非线性函数:in,
Figure BDA0004011310020000094
represents a measurement of target k on sensor s at time t, μ k (t) is Gaussian white noise with zero mean, that is, μ k (t) ~ N (0, R k (t)), R k (t) represents the covariance, X k (t) represents the profile of the monitored target at time t, and h (·) is a nonlinear function:

Figure BDA0004011310020000095
Figure BDA0004011310020000095

其中,h1,h2是均值为0,在区间[-1,1]之间服从均匀分布的随机变量。Among them, h 1 and h 2 are random variables with mean 0 and uniform distribution in the interval [-1,1].

本实施例中,所有监测目标的初始目标状态xk(0)和初始目标轮廓Xk(0)是已知的。In this embodiment, the initial target states x k (0) and initial target profiles X k (0) of all monitored targets are known.

步骤S101,获取在前一时刻下各监测目标的目标状态的状态估计值以及目标轮廓的轮廓估计值。Step S101, obtaining a state estimation value of a target state of each monitored target and a contour estimation value of a target contour at a previous moment.

在本发明实施例中,可以先获得t-1时刻(即前一时刻)各监测目标的状态以及轮廓参数的估计(即获取在前一时刻下各监测目标的目标状态的状态估计值以及目标轮廓的轮廓估计值)。In an embodiment of the present invention, the state of each monitored target at time t-1 (i.e., the previous moment) and the estimation of the profile parameters can be obtained first (i.e., the state estimation value of the target state of each monitored target and the profile estimation value of the target profile at the previous moment are obtained).

可选地,在获取在前一时刻下各监测目标的目标状态的状态估计值以及目标轮廓的轮廓估计值之后,还包括:基于在前一时刻下监测目标的状态估计值,确定监测目标在当前时刻下的状态预测值;基于在前一时刻下监测目标的轮廓估计值,确定监测目标在当前时刻下的轮廓预测值;基于状态预测值以及轮廓预测值,得到在当前时刻下监测目标的预测量测值,其中,预测量测值是监测目标在当前时刻下的预测状态所对应的量测值。Optionally, after obtaining the state estimation value of the target state of each monitoring target at the previous moment and the contour estimation value of the target contour, it also includes: determining the state prediction value of the monitoring target at the current moment based on the state estimation value of the monitoring target at the previous moment; determining the contour prediction value of the monitoring target at the current moment based on the contour estimation value of the monitoring target at the previous moment; and obtaining the predicted measurement value of the monitoring target at the current moment based on the state prediction value and the contour prediction value, wherein the predicted measurement value is the measurement value corresponding to the predicted state of the monitoring target at the current moment.

在本发明实施例中,对每个监测目标k,有:In the embodiment of the present invention, for each monitoring target k, there are:

Figure BDA0004011310020000101
Figure BDA0004011310020000101

Figure BDA0004011310020000102
Figure BDA0004011310020000102

其中,

Figure BDA0004011310020000103
分别表示t-1时刻下监测目标k的状态以及轮廓的估计。in,
Figure BDA0004011310020000103
They respectively represent the state and contour estimation of the monitored target k at time t-1.

在本实施例中,可以根据在前一时刻下监测目标的状态估计值

Figure BDA0004011310020000104
以及轮廓估计值
Figure BDA0004011310020000105
确定监测目标在当前时刻下的状态预测值以及轮廓预测值(即可以根据前一时刻下监测目标的状态以及轮廓的估计,通过公式(4)、公式(5)以及公式(1)得到监测目标在当前时刻下的状态预测值以及轮廓预测值)。In this embodiment, the state estimation value of the monitoring target at the previous moment can be used.
Figure BDA0004011310020000104
and the silhouette estimate
Figure BDA0004011310020000105
Determine the state prediction value and profile prediction value of the monitoring target at the current moment (that is, the state prediction value and profile prediction value of the monitoring target at the current moment can be obtained through formula (4), formula (5) and formula (1) based on the estimation of the state and profile of the monitoring target at the previous moment).

在本实施例中,预测量测值是监测目标在当前时刻下的预测状态所对应的量测值(即监测目标在t时刻的预测状态对应的量测值)。In this embodiment, the predicted measurement value is the measurement value corresponding to the predicted state of the monitoring target at the current moment (ie, the measurement value corresponding to the predicted state of the monitoring target at moment t).

步骤S102,对监测目标的状态估计值以及轮廓估计值进行粒子采样,得到在当前时刻下监测目标的预测量测值的带权粒子。Step S102, performing particle sampling on the state estimation value and the profile estimation value of the monitoring target to obtain weighted particles of the predicted measurement value of the monitoring target at the current moment.

可选地,对监测目标的状态估计值以及轮廓估计值进行粒子采样,得到在当前时刻下监测目标的预测量测值的带权粒子的步骤,包括:对状态估计值以及轮廓估计值进行粒子采样,生成多个粒子;基于预设参数,计算每个粒子的初始粒子权重;基于所有初始粒子权重,对各初始粒子权重进行归一化处理,得到每个粒子的目标粒子权重,其中,将携带有目标粒子权重的粒子表征为带权粒子;更新带权粒子,得到监测目标在当前时刻下关联预测量测值的带权粒子。Optionally, the step of performing particle sampling on the state estimation value and the profile estimation value of the monitoring target to obtain weighted particles of the predicted measurement value of the monitoring target at the current moment includes: performing particle sampling on the state estimation value and the profile estimation value to generate multiple particles; calculating the initial particle weight of each particle based on preset parameters; normalizing each initial particle weight based on all initial particle weights to obtain the target particle weight of each particle, wherein particles carrying the target particle weight are characterized as weighted particles; updating the weighted particles to obtain weighted particles associated with the predicted measurement value of the monitoring target at the current moment.

在本发明实施例中,可以对各监测目标的状态估计值以及轮廓估计值进行粒子采样,以得到在当前时刻下监测目标的预测量测值的带权粒子。具体为:In an embodiment of the present invention, particle sampling may be performed on the state estimation value and the profile estimation value of each monitoring target to obtain weighted particles of the predicted measurement value of the monitoring target at the current moment. Specifically:

对于每一组目标状态以及目标轮廓(xk(t),Xk(t))进行粒子采样,可以生成J个粒子(即对状态估计值以及轮廓估计值进行粒子采样,生成多个粒子),得到一组带权粒子

Figure BDA0004011310020000106
其中,wk,j表示粒子j的粒子权重。For each set of target states and target profiles (x k (t), X k (t)), particle sampling can generate J particles (i.e., particle sampling is performed on the state estimation value and the profile estimation value to generate multiple particles), and a set of weighted particles is obtained.
Figure BDA0004011310020000106
Where w k,j represents the particle weight of particle j.

粒子权重的计算方式如下:Particle weights are calculated as follows:

wk,j=f(xk,j|xk(t),cXk(t))(6);w k,j =f(x k,j |x k (t),cX k (t))(6);

其中,c>1为放大因子(即预设参数),所有权重计算完成后进行归一化处理(即基于预设参数,计算每个粒子的初始粒子权重,再基于所有初始粒子权重,对各初始粒子权重进行归一化处理,得到每个粒子的目标粒子权重):Among them, c>1 is the amplification factor (i.e., the preset parameter), and all weights are normalized after calculation (i.e., based on the preset parameters, the initial particle weight of each particle is calculated, and then based on all the initial particle weights, each initial particle weight is normalized to obtain the target particle weight of each particle):

Figure BDA0004011310020000111
Figure BDA0004011310020000111

在本实施例中,可以将携带有目标粒子权重的粒子表征为带权粒子。之后更新带权粒子,以得到监测目标在当前时刻下关联预测量测值的带权粒子(即对采样得到的带权粒子进行更新,得到各目标在t时刻的预测状态对应量测值的带权粒子)。对带权粒子中的状态进行更新的公式为:In this embodiment, the particles carrying the weight of the target particles can be characterized as weighted particles. The weighted particles are then updated to obtain weighted particles with associated predicted measurement values of the monitoring target at the current moment (i.e., the weighted particles obtained by sampling are updated to obtain weighted particles with corresponding measurement values of the predicted state of each target at time t). The formula for updating the state in the weighted particles is:

xk,j=Fxk,j(8);x k,j =Fx k,j (8);

其中,F表示状态转移矩阵。Where F represents the state transfer matrix.

可选地,基于在当前时刻下监测目标的预测量测值的带权粒子,对每个传感器的量测数据进行处理之前,还包括:获取每个传感器的量测数据,其中,量测数据至少包括:量测向量集合、量测总数;基于量测数据,构建面向监测目标的目标关联变量,其中,目标关联变量用于表示预设监测目标在当前时刻下产生的第一预设位数的量测指示的量测索引;基于量测数据,构建面向量测的量测关联变量,其中,量测关联变量用于表示在当前时刻下由量测索引所映射的监测目标产生的第二预设位数的量测。Optionally, based on the weighted particles of the predicted measurement values of the monitored target at the current moment, before processing the measurement data of each sensor, it also includes: obtaining the measurement data of each sensor, wherein the measurement data at least includes: a measurement vector set and a total number of measurements; based on the measurement data, constructing a target-associated variable for the monitored target, wherein the target-associated variable is used to represent a measurement index indicating a first preset number of bits of measurement generated by a preset monitored target at the current moment; based on the measurement data, constructing a measurement-associated variable for vector measurement, wherein the measurement-associated variable is used to represent a second preset number of bits of measurement generated by the monitored target mapped by the measurement index at the current moment.

在本发明实施例中,对任一传感器s,监测目标k最多产生

Figure BDA0004011310020000112
个量测,为了区分监测目标k产生哪些量测,传感器s使用面向监测目标的关联变量进行描述:In the embodiment of the present invention, for any sensor s, the monitoring target k generates at most
Figure BDA0004011310020000112
In order to distinguish which measurements are generated by the monitoring target k, the sensor s is described using associated variables oriented to the monitoring target:

Figure BDA0004011310020000113
Figure BDA0004011310020000113

为了表示出量测是源于目标还是源于杂波,传感器s定义如下面向量测的关联变量:In order to indicate whether the measurement is from the target or from clutter, the sensor s defines the following variables associated with the measurement:

Figure BDA0004011310020000121
Figure BDA0004011310020000121

在本实施例中,可以先定义传感器索引为s∈{1,2,…,S},S为传感器总个数;目标索引为k∈{1,2,…,K},K为目标总数量,并假定已知。采用zs表示传感器s量测的集合,

Figure BDA0004011310020000122
表示传感器s产生的量测的总个数。假设各目标产生的最大量测数为
Figure BDA0004011310020000123
面向目标k的关联变量αkq表示目标k产生的第q个量测对应的量测索引,其值范围为
Figure BDA0004011310020000124
其中,0表示未产生对应的量测;面向量测的关联变量βm=(k,q),采用二元组表示索引为m的量测对应为目标k产生的第q个量测。In this embodiment, the sensor index can be defined as s∈{1,2,…,S}, where S is the total number of sensors; the target index can be defined as k∈{1,2,…,K}, where K is the total number of targets, and is assumed to be known. z s is used to represent the set of measurements of sensor s,
Figure BDA0004011310020000122
represents the total number of measurements generated by sensor s. Assume that the maximum number of measurements generated by each target is
Figure BDA0004011310020000123
The associated variable α kq for target k represents the measurement index corresponding to the qth measurement generated by target k, and its value range is
Figure BDA0004011310020000124
Wherein, 0 indicates that no corresponding measurement is generated; the associated variable β m =(k,q) of the face vector measurement uses a two-tuple to indicate that the measurement with index m corresponds to the qth measurement generated for target k.

在本实施例中,可以先获取每个传感器的量测数据,该量测数据至少包括:量测向量集合zs、量测总数

Figure BDA0004011310020000125
然后根据量测数据,构建面向监测目标的目标关联变量αkq,该目标关联变量用于表示预设监测目标在当前时刻下产生的第一预设位数(例如,第q个)的量测指示的量测索引,并根据量测数据,构建面向量测的量测关联变量βm=(k,q),该量测关联变量用于表示在当前时刻下由量测索引所映射的监测目标产生的第二预设位数(例如,第q个)的量测。In this embodiment, the measurement data of each sensor can be obtained first, and the measurement data at least includes: a measurement vector set zs, a measurement total number
Figure BDA0004011310020000125
Then, based on the measurement data, a target-associated variable α kq for the monitoring target is constructed, and the target-associated variable is used to represent the measurement index of the measurement indication of the first preset number of bits (for example, the qth) generated by the preset monitoring target at the current moment. And based on the measurement data, a measurement-associated variable β m = (k, q) for vector measurement is constructed, and the measurement-associated variable is used to represent the second preset number of bits (for example, the qth) of measurement generated by the monitoring target mapped by the measurement index at the current moment.

步骤S103,基于在当前时刻下监测目标的预测量测值的带权粒子,对每个传感器的量测数据进行处理,得到在当前时刻下与传感器对应的所有监测目标的当前状态估计值以及当前轮廓估计值。Step S103, based on the weighted particles of the predicted measurement values of the monitoring targets at the current moment, the measurement data of each sensor is processed to obtain the current state estimation values and the current contour estimation values of all monitoring targets corresponding to the sensor at the current moment.

可选地,基于在当前时刻下监测目标的预测量测值的带权粒子,对每个传感器的量测数据进行处理,得到在当前时刻下与传感器对应的所有监测目标的当前状态估计值以及当前轮廓估计值的步骤,包括:基于在当前时刻下监测目标的预测量测值的带权粒子以及量测向量集合,对目标关联变量进行量测评估,得到似然概率;在似然概率属于预设概率阈值范围的情况轮廓估计值下,基于目标关联变量以及量测关联变量,构建关联因子图;基于关联因子图,进行迭代计算,得到第一参数值以及第二参数值;基于第一参数值以及第二参数值,更新传感器的监测目标的量测数据,得到更新结果;基于更新结果,得到与传感器对应的在当前时刻下监测目标的当前轮廓估计值,并基于更新结果,对带权粒子进行加权计算,得到与传感器对应的监测目标的当前状态估计值。Optionally, based on the weighted particles of the predicted measurement values of the monitored targets at the current moment, the measurement data of each sensor is processed to obtain the current state estimation values and the current profile estimation values of all the monitored targets corresponding to the sensor at the current moment, including: based on the weighted particles of the predicted measurement values of the monitored targets at the current moment and the set of measurement vectors, the target-associated variables are measured and evaluated to obtain the likelihood probability; in the case where the likelihood probability belongs to the profile estimation value within the preset probability threshold range, an association factor graph is constructed based on the target-associated variables and the measurement-associated variables; based on the association factor graph, an iterative calculation is performed to obtain a first parameter value and a second parameter value; based on the first parameter value and the second parameter value, the measurement data of the monitored targets of the sensor are updated to obtain an updated result; based on the updated result, the current profile estimation value of the monitored targets corresponding to the sensor at the current moment is obtained, and based on the updated result, the weighted particles are weighted calculated to obtain the current state estimation value of the monitored targets corresponding to the sensor.

在本发明实施例中,所有传感器s可以并行进行内部的数据关联迭代,从而各自得到t时刻各目标的状态以及轮廓参数的估计(即可以基于在当前时刻下监测目标的预测量测值的带权粒子,对每个传感器的量测数据进行处理,得到在当前时刻下与传感器对应的所有监测目标的当前状态估计值以及当前轮廓估计值)。具体为:In the embodiment of the present invention, all sensors s can perform internal data association iterations in parallel, so that each can obtain the state of each target at time t and the estimation of the profile parameters (that is, the measurement data of each sensor can be processed based on the weighted particles of the predicted measurement values of the monitoring targets at the current moment to obtain the current state estimation values and current profile estimation values of all monitoring targets corresponding to the sensor at the current moment). Specifically:

(1)量测评估,使用目标在t时刻的预测状态对应量测值的带权粒子通过计算对应的似然概率对测量数据进行评估(即基于在当前时刻下监测目标的预测量测值的带权粒子以及量测向量集合,对目标关联变量进行量测评估,得到似然概率):(1) Measurement evaluation: Use the weighted particles corresponding to the measurement value of the target's predicted state at time t to evaluate the measurement data by calculating the corresponding likelihood probability (i.e., based on the weighted particles and measurement vector set of the predicted measurement value of the monitoring target at the current time, the target-associated variable is measured and evaluated to obtain the likelihood probability):

使用

Figure BDA0004011310020000131
对于各目标k的关联变量αkq进行量测评估,具体评估公式如下:use
Figure BDA0004011310020000131
The associated variable α kq of each target k is measured and evaluated. The specific evaluation formula is as follows:

Figure BDA0004011310020000132
Figure BDA0004011310020000132

其中,

Figure BDA0004011310020000133
表示带权粒子,J为粒子总数;jkq(xk,j,Xkkq;zs)计算公式如下:in,
Figure BDA0004011310020000133
represents weighted particles, J is the total number of particles; the calculation formula of j kq (x k,j ,X kkq ;z s ) is as follows:

Figure BDA0004011310020000134
Figure BDA0004011310020000134

其中,Xk表示目标k的目标轮廓,本实施例中采用的为随机矩阵;

Figure BDA0004011310020000135
Figure BDA0004011310020000136
在以xk,j为均值,Xk为协方差矩阵的高斯分布下的概率;μfa为在测量场景中虚警量测出现个数的均值,本实施例中假设均值的数量服从泊松分布,在量测范围内服从均匀分布;
Figure BDA0004011310020000137
表示传感器s产生的量测的总个数,
Figure BDA0004011310020000138
表示传感器s产生的第m个量测,
Figure BDA0004011310020000139
表示对应量测的虚警概率。Wherein, X k represents the target profile of target k, and a random matrix is used in this embodiment;
Figure BDA0004011310020000135
for
Figure BDA0004011310020000136
The probability under Gaussian distribution with x k,j as the mean and X k as the covariance matrix; μfa is the mean of the number of false alarm measurements in the measurement scenario. In this embodiment, it is assumed that the number of means follows a Poisson distribution and follows a uniform distribution within the measurement range;
Figure BDA0004011310020000137
represents the total number of measurements produced by sensor s,
Figure BDA0004011310020000138
represents the mth measurement produced by sensor s,
Figure BDA0004011310020000139
Represents the false alarm probability of the corresponding measurement.

(2)迭代数据关联,通过对面向监测目标以及面向量测的关联变量构建数据关联因子图,通过迭代计算得到各关联情况对应的概率(即在似然概率属于预设概率阈值范围的情况轮廓估计值下,基于目标关联变量以及量测关联变量,构建关联因子图,并基于关联因子图,进行迭代计算,得到第一参数值以及第二参数值):(2) Iterative data association, by constructing a data association factor graph for the association variables facing the monitoring target and the measurement, and obtaining the probability corresponding to each association situation through iterative calculation (that is, under the situation profile estimation value where the likelihood probability belongs to the preset probability threshold range, the association factor graph is constructed based on the target association variable and the measurement association variable, and iterative calculation is performed based on the association factor graph to obtain the first parameter value and the second parameter value):

所有传感器可以并行进行执行,分别使用各自的量测数据进行迭代数据关联。定义迭代计算索引为

Figure BDA0004011310020000141
其中,nit表示循环计算的总次数,每个迭代计算的过程中,对于所有可行的二元组
Figure BDA0004011310020000142
以及所有量测
Figure BDA0004011310020000143
进行如下计算:All sensors can be executed in parallel, using their own measurement data for iterative data association. Define the iterative calculation index as
Figure BDA0004011310020000141
Among them, nit represents the total number of loop calculations. During each iterative calculation, for all feasible binary pairs
Figure BDA0004011310020000142
And all measurements
Figure BDA0004011310020000143
Perform the following calculations:

Figure BDA0004011310020000144
Figure BDA0004011310020000144

Figure BDA0004011310020000145
Figure BDA0004011310020000145

其中,

Figure BDA0004011310020000146
表示对除选定的二元组(k,q)进行和积算法的运算;in,
Figure BDA0004011310020000146
It means to perform the sum-product algorithm operation on the selected two-tuple (k,q);

Figure BDA0004011310020000147
Figure BDA0004011310020000147

Figure BDA0004011310020000148
Figure BDA0004011310020000149
表示αk包含的有效量测个数,
Figure BDA00040113100200001410
表示有效量测个数对应的概率。
Figure BDA0004011310020000148
Figure BDA0004011310020000149
represents the number of valid measurements contained in α k ,
Figure BDA00040113100200001410
Represents the probability corresponding to the number of valid measurements.

对于

Figure BDA00040113100200001411
进行如下计算:for
Figure BDA00040113100200001411
Perform the following calculations:

Figure BDA00040113100200001412
Figure BDA00040113100200001412

并且对于迭代计算进行如下初始计算:And the following initial calculation is performed for the iterative calculation:

Figure BDA00040113100200001413
Figure BDA00040113100200001413

利用因子图的特性对该迭代数据关联计算进行优化,可以将上述计算公式替换如下:By using the characteristics of factor graphs to optimize the iterative data association calculation, the above calculation formula can be replaced as follows:

Figure BDA0004011310020000151
Figure BDA0004011310020000151

Figure BDA0004011310020000152
Figure BDA0004011310020000152

Figure BDA0004011310020000153
Figure BDA0004011310020000153

其中,

Figure BDA0004011310020000154
Figure BDA0004011310020000155
in,
Figure BDA0004011310020000154
Figure BDA0004011310020000155

(3)量测更新,利用迭代数据关联得到的结果(即第一参数值

Figure BDA0004011310020000156
以及第二参数值
Figure BDA0004011310020000157
),对加权粒子进行测量更新,同时根据关联概率选择合适的一部分量测值对目标的轮廓进行估计(即基于第一参数值以及第二参数值,更新传感器的监测目标的量测数据,得到更新结果,并基于更新结果,得到与传感器对应的在当前时刻下监测目标的当前轮廓估计值):(3) Measurement update, using the result obtained by iterative data association (i.e., the first parameter value)
Figure BDA0004011310020000156
And the second parameter value
Figure BDA0004011310020000157
), measure and update the weighted particles, and select a suitable part of the measured values according to the associated probability to estimate the contour of the target (that is, based on the first parameter value and the second parameter value, update the measurement data of the monitoring target of the sensor to obtain the updated result, and based on the updated result, obtain the current contour estimation value of the monitoring target corresponding to the sensor at the current moment):

经过迭代数据关联计算后,得到

Figure BDA0004011310020000158
以及
Figure BDA0004011310020000159
进行如下量测更新:After iterative data association calculation, we get
Figure BDA0004011310020000158
as well as
Figure BDA0004011310020000159
Perform the following measurement updates:

Figure BDA0004011310020000161
Figure BDA0004011310020000161

Figure BDA0004011310020000162
Figure BDA0004011310020000162

Figure BDA0004011310020000163
Figure BDA0004011310020000163

对于各目标选取用于更新轮廓的量测

Figure BDA0004011310020000164
并进行如下更新:For each target, select the measurement used to update the contour
Figure BDA0004011310020000164
And update as follows:

Figure BDA0004011310020000165
Figure BDA0004011310020000165

Figure BDA0004011310020000166
Figure BDA0004011310020000166

其中,γ表示关联概率,Zk表示所选量测对应的协方差矩阵,τk以及νk为随机矩阵更新中的参数。Among them, γ represents the association probability, Zk represents the covariance matrix corresponding to the selected measurement, τk and νk are parameters in the random matrix update.

(4)目标估计,对带权粒子进行加权计算,得到对目标状态的估计(即基于更新结果,对带权粒子进行加权计算,得到与传感器对应的监测目标的当前状态估计值):(4) Target estimation: weighted calculation is performed on the weighted particles to obtain an estimate of the target state (i.e., based on the update result, weighted calculation is performed on the weighted particles to obtain an estimate of the current state of the monitoring target corresponding to the sensor):

使用得到的量测更新结果,结合各组带权粒子对估计各目标的状态:Use the obtained measurement update results and combine each group of weighted particle pairs to estimate the state of each target:

Figure BDA0004011310020000167
Figure BDA0004011310020000167

步骤S104,基于当前状态估计值以及当前轮廓估计值,筛选与每个监测目标对应的目标传感器集合。Step S104: based on the current state estimation value and the current profile estimation value, filter the target sensor set corresponding to each monitoring target.

可选地,基于当前状态估计值以及当前轮廓估计值,筛选与每个监测目标对应的目标传感器集合的步骤,包括:将当前状态估计值表征为预设高斯分布的均值,并将当前轮廓估计值表征为预设高斯分布的协方差矩阵;基于均值以及协方差矩阵,计算每个传感器对监测目标的估计概率;对所有估计概率进行归一化处理,得到目标估计概率;在目标估计概率大于预设概率阈值的情况下,将与目标估计概率对应的传感器加入至监测目标对应的目标传感器集合。Optionally, based on the current state estimation value and the current profile estimation value, the step of screening the target sensor set corresponding to each monitoring target includes: characterizing the current state estimation value as the mean of a preset Gaussian distribution, and characterizing the current profile estimation value as the covariance matrix of a preset Gaussian distribution; based on the mean and the covariance matrix, calculating the estimated probability of each sensor for the monitoring target; normalizing all estimated probabilities to obtain the target estimated probability; and when the target estimated probability is greater than a preset probability threshold, adding the sensor corresponding to the target estimated probability to the target sensor set corresponding to the monitoring target.

在本发明实施例中,融合中心可以根据各传感器得到的估计,基于投票机制选择正确率较高的估计(即基于当前状态估计值以及当前轮廓估计值,筛选与每个监测目标对应的目标传感器集合)。具体为:In the embodiment of the present invention, the fusion center can select the estimate with higher accuracy based on the voting mechanism according to the estimate obtained by each sensor (that is, based on the current state estimate value and the current profile estimate value, filter the target sensor set corresponding to each monitoring target). Specifically:

对于各传感器得到的目标状态以及目标轮廓的估计

Figure BDA0004011310020000171
由融合中心对于各估计进行选择,过程如下:Estimation of target state and target contour obtained by each sensor
Figure BDA0004011310020000171
The fusion center makes a selection for each estimate, and the process is as follows:

对于各目标k,传感器s,将其余传感器对状态以及轮廓的估计看作是高斯分布的均值和协方差矩阵,计算对应的概率,进行交叉验证投票(即将当前状态估计值表征为预设高斯分布的均值,并将当前轮廓估计值表征为预设高斯分布的协方差矩阵;基于均值以及协方差矩阵,计算每个传感器对监测目标的估计概率

Figure BDA0004011310020000172
):For each target k, sensor s, the remaining sensors' estimates of the state and profile are regarded as the mean and covariance matrix of the Gaussian distribution, the corresponding probabilities are calculated, and cross-validation voting is performed (i.e., the current state estimate is represented as the mean of the preset Gaussian distribution, and the current profile estimate is represented as the covariance matrix of the preset Gaussian distribution; based on the mean and covariance matrix, the estimated probability of each sensor for the monitored target is calculated
Figure BDA0004011310020000172
):

Figure BDA0004011310020000173
Figure BDA0004011310020000173

计算完成后对同一目标k所有传感器的概率进行归一化(即对所有估计概率进行归一化处理,得到目标估计概率),并选择大于阈值ε的传感器估计进行后续融合,目标k选择的传感器索引集合表示为Sk(即在目标估计概率大于预设概率阈值ε的情况下,将与目标估计概率对应的传感器加入至监测目标对应的目标传感器集合)。After the calculation is completed, the probabilities of all sensors of the same target k are normalized (that is, all estimated probabilities are normalized to obtain the target estimated probability), and the sensor estimates greater than the threshold ε are selected for subsequent fusion. The sensor index set selected by target k is represented as S k (that is, when the target estimated probability is greater than the preset probability threshold ε, the sensor corresponding to the target estimated probability is added to the target sensor set corresponding to the monitoring target).

步骤S105,融合目标传感器集合中与每个目标传感器对应的当前状态估计值以及当前轮廓估计值,得到在当前时刻下监测目标的目标状态值以及目标轮廓值。Step S105 , fusing the current state estimation value and the current contour estimation value corresponding to each target sensor in the target sensor set to obtain the target state value and the target contour value of the monitored target at the current moment.

可选地,融合目标传感器集合中与每个目标传感器对应的当前状态估计值以及当前轮廓估计值,得到在当前时刻下监测目标的目标状态值以及目标轮廓值的步骤,包括:计算所有与每个目标传感器对应的当前状态估计值的平均值,得到目标状态值;基于当前轮廓估计值,确定监测目标的当前目标轮廓的第一半轴长值、第二半轴长值以及偏转角度值;基于所有目标传感器的协方差矩阵,确定目标协方差矩阵;基于目标协方差矩阵,确定目标偏转角度值;基于目标偏转角度值,计算每个目标传感器的目标第一半轴长值以及目标第二半轴长值;基于所有目标第一半轴长值以及目标第二半轴长值,得到目标半轴长值;基于目标半轴长值以及目标偏转角度值,得到目标轮廓值。Optionally, the step of fusing the current state estimation value and the current contour estimation value corresponding to each target sensor in the target sensor set to obtain the target state value and the target contour value of the monitored target at the current moment includes: calculating the average of all the current state estimation values corresponding to each target sensor to obtain the target state value; determining the first semi-axis length value, the second semi-axis length value and the deflection angle value of the current target contour of the monitored target based on the current contour estimation value; determining the target covariance matrix based on the covariance matrix of all target sensors; determining the target deflection angle value based on the target covariance matrix; calculating the target first semi-axis length value and the target second semi-axis length value of each target sensor based on the target deflection angle value; obtaining the target semi-axis length value based on all target first semi-axis length values and target second semi-axis length values; obtaining the target contour value based on the target semi-axis length value and the target deflection angle value.

在本发明实施例中,对优选的目标传感器集合中的估计进行融合,最终得到t时刻各目标的状态以及轮廓参数的估计(即融合目标传感器集合中与每个目标传感器对应的当前状态估计值以及当前轮廓估计值,得到在当前时刻下监测目标的目标状态值以及目标轮廓值)。In an embodiment of the present invention, the estimates in the preferred target sensor set are fused to ultimately obtain estimates of the state and profile parameters of each target at time t (i.e., the current state estimate value and the current profile estimate value corresponding to each target sensor in the target sensor set are fused to obtain the target state value and the target profile value of the monitored target at the current moment).

对于目标k,根据确定的传感器索引集合Sk进行状态和轮廓的估计。For target k, the state and contour are estimated according to the determined sensor index set Sk .

对于目标状态,采用算数平均的方法进行计算(即计算所有与每个目标传感器对应的当前状态估计值的平均值,得到目标状态值):For the target state, the arithmetic average method is used for calculation (that is, the average value of all current state estimates corresponding to each target sensor is calculated to obtain the target state value):

Figure BDA0004011310020000181
Figure BDA0004011310020000181

对于目标的轮廓,可以先将其转换为半轴长加偏转角度的形式(即基于当前轮廓估计值,确定监测目标的当前目标轮廓的第一半轴长值、第二半轴长值以及偏转角度值,其中,第一半轴长值即为长半轴长,第二半轴长值即为短半轴长):For the contour of the target, it can be first converted into the form of semi-axis length plus deflection angle (that is, based on the current contour estimation value, the first semi-axis length value, the second semi-axis length value and the deflection angle value of the current target contour of the monitored target are determined, wherein the first semi-axis length value is the major semi-axis length, and the second semi-axis length value is the minor semi-axis length):

Figure BDA0004011310020000182
Figure BDA0004011310020000182

其中,ls1,ls2s分别为长半轴长、短半轴长以及偏转角度。Among them, l s1 ,l s2s are the major semi-axis length, minor semi-axis length and deflection angle respectively.

对于各传感器对轮廓的估计,将其看作是高斯分布的协方差矩阵,将这些高斯分布

Figure BDA0004011310020000183
进行相乘得到最终的协方差矩阵σ2(即基于所有目标传感器的协方差矩阵,确定目标协方差矩阵),高斯分布相乘的协方差矩阵公式如下:For the estimation of the contour by each sensor, it is regarded as the covariance matrix of the Gaussian distribution.
Figure BDA0004011310020000183
The final covariance matrix σ 2 is obtained by multiplication (i.e., the target covariance matrix is determined based on the covariance matrices of all target sensors). The formula of the covariance matrix of Gaussian distribution multiplication is as follows:

Figure BDA0004011310020000184
Figure BDA0004011310020000184

通过得到的协方差矩阵σ2确定目标偏转角度的估计

Figure BDA0004011310020000185
(即基于目标协方差矩阵,确定目标偏转角度值),由此可以计算各传感器对应目标长短半轴的估计(即基于目标偏转角度值,计算每个目标传感器的目标第一半轴长值
Figure BDA0004011310020000186
以及目标第二半轴长值
Figure BDA0004011310020000187
):The estimated target deflection angle is determined by the obtained covariance matrix σ2
Figure BDA0004011310020000185
(i.e., based on the target covariance matrix, determine the target deflection angle value), from which the estimation of the target long and short semi-axis corresponding to each sensor can be calculated (i.e., based on the target deflection angle value, calculate the target first semi-axis length value of each target sensor
Figure BDA0004011310020000186
And the target second semi-axis length value
Figure BDA0004011310020000187
):

Figure BDA0004011310020000188
Figure BDA0004011310020000188

Figure BDA0004011310020000191
Figure BDA0004011310020000191

最终通过算术平均的方法得到长短半轴长的估计(即基于所有目标第一半轴长值以及目标第二半轴长值,得到目标半轴长值

Figure BDA0004011310020000192
再基于目标半轴长值以及目标偏转角度值,得到目标轮廓值):Finally, the estimation of the major and minor semi-axis lengths is obtained by the arithmetic mean method (i.e., the target semi-axis length value is obtained based on all the target first semi-axis length values and the target second semi-axis length values).
Figure BDA0004011310020000192
Then, based on the target semi-axis length value and the target deflection angle value, the target contour value is obtained):

Figure BDA0004011310020000193
Figure BDA0004011310020000193

图2是根据本发明实施例的一种可选的面向生命线安全的多传感器多扩展目标融合算法原理的流程图,如图2所示,包括如下流程:FIG2 is a flow chart of an optional multi-sensor multi-extended target fusion algorithm principle for lifeline safety according to an embodiment of the present invention. As shown in FIG2 , the flow chart includes the following steps:

步骤(1):t-1时刻目标的估计,即得到t-1时刻各目标的状态以及轮廓参数的估计;Step (1): Estimation of the target at time t-1, that is, obtaining the state of each target at time t-1 and the estimation of the profile parameters;

步骤(2):粒子采样,即对各目标的估计进行粒子采样;Step (2): Particle sampling, i.e., performing particle sampling on the estimation of each target;

步骤(3):状态更新,即对步骤(2)中的采样得到的带权粒子进行更新,得到各目标在t时刻的预测状态对应的量测值的带权粒子;Step (3): state update, that is, updating the weighted particles obtained by sampling in step (2) to obtain the weighted particles of the measurement values corresponding to the predicted state of each target at time t;

步骤(4):传感器并行数据关联迭代,即所有传感器1,…,传感器s等并行进行内部的数据关联迭代,各自得到t时刻各目标的状态以及轮廓参数的估计;Step (4): sensor parallel data association iteration, that is, all sensors 1, ..., sensor s, etc. perform internal data association iteration in parallel, and each obtains the state of each target at time t and the estimation of the profile parameters;

步骤(5):优选估计,即融合中心根据各传感器得到的估计,基于投票机制选择正确率较高的估计;Step (5): Optimizing the estimation, i.e., the fusion center selects the estimation with the higher accuracy based on the estimation obtained by each sensor based on the voting mechanism;

步骤(6):估计融合,即对步骤(5)中优选的估计进行融合,最终得到t时刻目标的估计(即t时刻各目标的状态以及轮廓参数的估计)。Step (6): estimation fusion, that is, fusing the preferred estimates in step (5) to finally obtain an estimate of the target at time t (that is, an estimate of the state of each target and the profile parameters at time t).

本发明实施例中,能够基于因子图优化实现传感器多扩展目标跟踪,并通过使用投票机制优选传感器估计结果之后,再实现多传感器的扩展目标估计融合,达到了如下有益效果:(1)基于因子图优化实现多扩展目标跟踪,计算速度快,并且不需要进行聚类等预处理操作;(2)投票机制优选以及多传感器扩展目标估计融合,使用交叉验证投票能够确保优选效果,并且有效提高了利用多传感器的量测数据得到的估计结果的精准性。In the embodiments of the present invention, sensor multi-extended target tracking can be realized based on factor graph optimization, and after optimizing sensor estimation results by using a voting mechanism, multi-sensor extended target estimation fusion can be realized, thereby achieving the following beneficial effects: (1) multi-extended target tracking is realized based on factor graph optimization, with fast calculation speed and no need for preprocessing operations such as clustering; (2) the voting mechanism optimization and multi-sensor extended target estimation fusion can ensure the optimization effect by using cross-validation voting, and effectively improve the accuracy of the estimation results obtained using the measurement data of multiple sensors.

下面结合另一实施例进行详细说明。The following is a detailed description in conjunction with another embodiment.

实施例二Embodiment 2

本实施例中提供的一种基于多传感器的监测目标融合装置包含了多个实施单元,每个实施单元对应于上述实施例一中的各个实施步骤。A multi-sensor based monitoring target fusion device provided in this embodiment includes multiple implementation units, each implementation unit corresponds to each implementation step in the above-mentioned embodiment 1.

图3是根据本发明实施例的一种可选的基于多传感器的监测目标融合装置的示意图,如图3所示,该融合装置可以包括:获取单元30,采样单元31,处理单元32,筛选单元33,融合单元34,其中,FIG3 is a schematic diagram of an optional multi-sensor based monitoring target fusion device according to an embodiment of the present invention. As shown in FIG3 , the fusion device may include: an acquisition unit 30, a sampling unit 31, a processing unit 32, a screening unit 33, and a fusion unit 34, wherein:

获取单元30,用于获取在前一时刻下各监测目标的目标状态的状态估计值以及目标轮廓的轮廓估计值;An acquisition unit 30 is used to acquire a state estimation value of a target state of each monitoring target and a contour estimation value of a target contour at a previous moment;

采样单元31,用于对监测目标的状态估计值以及轮廓估计值进行粒子采样,得到在当前时刻下监测目标的预测量测值的带权粒子;The sampling unit 31 is used to perform particle sampling on the state estimation value and the profile estimation value of the monitoring target to obtain weighted particles of the predicted measurement value of the monitoring target at the current moment;

处理单元32,用于基于在当前时刻下监测目标的预测量测值的带权粒子,对每个传感器的量测数据进行处理,得到在当前时刻下与传感器对应的所有监测目标的当前状态估计值以及当前轮廓估计值;A processing unit 32 is used to process the measurement data of each sensor based on the weighted particles of the predicted measurement value of the monitoring target at the current moment, and obtain the current state estimation value and the current profile estimation value of all monitoring targets corresponding to the sensor at the current moment;

筛选单元33,用于基于当前状态估计值以及当前轮廓估计值,筛选与每个监测目标对应的目标传感器集合;A screening unit 33, configured to screen a target sensor set corresponding to each monitoring target based on a current state estimation value and a current profile estimation value;

融合单元34,用于融合目标传感器集合中与每个目标传感器对应的当前状态估计值以及当前轮廓估计值,得到在当前时刻下监测目标的目标状态值以及目标轮廓值。The fusion unit 34 is used to fuse the current state estimation value and the current contour estimation value corresponding to each target sensor in the target sensor set to obtain the target state value and the target contour value of the monitored target at the current moment.

上述融合装置,可以通过获取单元30获取在前一时刻下各监测目标的目标状态的状态估计值以及目标轮廓的轮廓估计值,通过采样单元31对监测目标的状态估计值以及轮廓估计值进行粒子采样,得到在当前时刻下监测目标的预测量测值的带权粒子,通过处理单元32基于在当前时刻下监测目标的预测量测值的带权粒子,对每个传感器的量测数据进行处理,得到在当前时刻下与传感器对应的所有监测目标的当前状态估计值以及当前轮廓估计值,通过筛选单元33基于当前状态估计值以及当前轮廓估计值,筛选与每个监测目标对应的目标传感器集合,通过融合单元34融合目标传感器集合中与每个目标传感器对应的当前状态估计值以及当前轮廓估计值,得到在当前时刻下监测目标的目标状态值以及目标轮廓值。在本发明实施例中,可以先对前一时刻下各监测目标的目标状态的状态估计值以及目标轮廓的轮廓估计值进行粒子采样,以得到当前时刻下监测目标的预测量测值的带权粒子,然后根据带权粒子对每个传感器的量测数据进行处理,以得到当前时刻下每个传感器对应的所有监测目标的当前状态估计值以及当前轮廓估计值,再筛选每个监测目标的目标传感器集合,之后融合每个目标传感器对应的当前状态估计值以及当前轮廓估计值,以得到当前时刻下监测目标的目标状态值以及目标轮廓值,能够实时计算监测目标的状态值以及轮廓值,并且利用筛选结果能够高效得到各监测目标的状态值以及轮廓值的准确估计,进而解决了相关技术中计算监测目标的状态值以及轮廓值的时效性较差以及准确率较低的技术问题。The above-mentioned fusion device can obtain the state estimation value of the target state and the contour estimation value of the target contour of each monitored target at the previous moment through the acquisition unit 30, perform particle sampling on the state estimation value and the contour estimation value of the monitored target through the sampling unit 31, and obtain the weighted particles of the predicted measurement value of the monitored target at the current moment, and process the measurement data of each sensor based on the weighted particles of the predicted measurement value of the monitored target at the current moment through the processing unit 32 to obtain the current state estimation value and the current contour estimation value of all monitored targets corresponding to the sensor at the current moment, and screen the target sensor set corresponding to each monitored target based on the current state estimation value and the current contour estimation value through the screening unit 33, and fuse the current state estimation value and the current contour estimation value corresponding to each target sensor in the target sensor set through the fusion unit 34 to obtain the target state value and the target contour value of the monitored target at the current moment. In an embodiment of the present invention, particle sampling can be first performed on the state estimation value of the target state and the contour estimation value of the target contour of each monitored target at the previous moment to obtain weighted particles of the predicted measurement value of the monitored target at the current moment, and then the measurement data of each sensor is processed according to the weighted particles to obtain the current state estimation value and current contour estimation value of all monitored targets corresponding to each sensor at the current moment, and then the target sensor set of each monitored target is screened, and then the current state estimation value and current contour estimation value corresponding to each target sensor are fused to obtain the target state value and target contour value of the monitored target at the current moment, so that the state value and contour value of the monitored target can be calculated in real time, and the screening result can be used to efficiently obtain accurate estimates of the state value and contour value of each monitored target, thereby solving the technical problems of poor timeliness and low accuracy in calculating the state value and contour value of the monitored target in the related art.

可选地,融合装置还包括:第一确定模块,用于在获取在前一时刻下各监测目标的目标状态的状态估计值以及目标轮廓的轮廓估计值之后,基于在前一时刻下监测目标的状态估计值,确定监测目标在当前时刻下的状态预测值;第二确定模块,用于基于在前一时刻下监测目标的轮廓估计值,确定监测目标在当前时刻下的轮廓预测值;第一输出模块,用于基于状态预测值以及轮廓预测值,得到在当前时刻下监测目标的预测量测值,其中,预测量测值是监测目标在当前时刻下的预测状态所对应的量测值。Optionally, the fusion device also includes: a first determination module, which is used to determine the state prediction value of the monitoring target at the current moment based on the state estimation value of the monitoring target at the previous moment after obtaining the state estimation value of the target state of each monitoring target and the contour estimation value of the target contour at the previous moment; a second determination module, which is used to determine the contour prediction value of the monitoring target at the current moment based on the contour estimation value of the monitoring target at the previous moment; and a first output module, which is used to obtain the predicted measurement value of the monitoring target at the current moment based on the state prediction value and the contour prediction value, wherein the predicted measurement value is the measurement value corresponding to the predicted state of the monitoring target at the current moment.

可选地,采样单元包括:第一生成模块,用于对状态估计值以及轮廓估计值进行粒子采样,生成多个粒子;第一计算模块,用于基于预设参数,计算每个粒子的初始粒子权重;第一处理模块,用于基于所有初始粒子权重,对各初始粒子权重进行归一化处理,得到每个粒子的目标粒子权重,其中,将携带有目标粒子权重的粒子表征为带权粒子;第一更新模块,用于更新带权粒子,得到监测目标在当前时刻下关联预测量测值的带权粒子。Optionally, the sampling unit includes: a first generation module, used to perform particle sampling on state estimation values and profile estimation values to generate multiple particles; a first calculation module, used to calculate the initial particle weight of each particle based on preset parameters; a first processing module, used to normalize each initial particle weight based on all initial particle weights to obtain a target particle weight for each particle, wherein particles carrying target particle weights are characterized as weighted particles; a first update module, used to update weighted particles to obtain weighted particles associated with predicted measurement values of the monitoring target at the current moment.

可选地,融合装置还包括:第一获取模块,用于基于在当前时刻下监测目标的预测量测值的带权粒子,对每个传感器的量测数据进行处理之前,获取每个传感器的量测数据,其中,量测数据至少包括:量测向量集合、量测总数;第一构建模块,用于基于量测数据,构建面向监测目标的目标关联变量,其中,目标关联变量用于表示预设监测目标在当前时刻下产生的第一预设位数的量测指示的量测索引;第二构建模块,用于基于量测数据,构建面向量测的量测关联变量,其中,量测关联变量用于表示在当前时刻下由量测索引所映射的监测目标产生的第二预设位数的量测。Optionally, the fusion device also includes: a first acquisition module, which is used to acquire the measurement data of each sensor before processing the measurement data of each sensor based on the weighted particles of the predicted measurement value of the monitoring target at the current moment, wherein the measurement data at least includes: a measurement vector set and a total number of measurements; a first construction module, which is used to construct a target-associated variable for the monitoring target based on the measurement data, wherein the target-associated variable is used to represent a measurement index indicating a measurement indication of a first preset number of bits generated by a preset monitoring target at the current moment; and a second construction module, which is used to construct a measurement-associated variable for vector measurement based on the measurement data, wherein the measurement-associated variable is used to represent a second preset number of bits of measurement generated by the monitoring target mapped by the measurement index at the current moment.

可选地,处理单元包括:第一评估模块,用于基于在当前时刻下监测目标的预测量测值的带权粒子以及量测向量集合,对目标关联变量进行量测评估,得到似然概率;第三构建模块,用于在似然概率属于预设概率阈值范围的情况轮廓估计值下,基于目标关联变量以及量测关联变量,构建关联因子图;第二计算模块,用于基于关联因子图,进行迭代计算,得到第一参数值以及第二参数值;第二更新模块,用于基于第一参数值以及第二参数值,更新传感器的监测目标的量测数据,得到更新结果;第三计算模块,用于基于更新结果,得到与传感器对应的在当前时刻下监测目标的当前轮廓估计值,并基于更新结果,对带权粒子进行加权计算,得到与传感器对应的监测目标的当前状态估计值。Optionally, the processing unit includes: a first evaluation module, which is used to measure and evaluate the target associated variables based on the weighted particles and measurement vector set of the predicted measurement values of the monitored target at the current moment to obtain the likelihood probability; a third construction module, which is used to construct a correlation factor graph based on the target associated variables and the measurement associated variables under the profile estimation value when the likelihood probability belongs to the preset probability threshold range; a second calculation module, which is used to perform iterative calculation based on the correlation factor graph to obtain the first parameter value and the second parameter value; a second update module, which is used to update the measurement data of the monitored target of the sensor based on the first parameter value and the second parameter value to obtain an updated result; a third calculation module, which is used to obtain the current profile estimation value of the monitored target corresponding to the sensor at the current moment based on the updated result, and to perform weighted calculation on the weighted particles based on the updated result to obtain the current state estimation value of the monitored target corresponding to the sensor.

可选地,筛选单元包括:第一表征模块,用于将当前状态估计值表征为预设高斯分布的均值,并将当前轮廓估计值表征为预设高斯分布的协方差矩阵;第四计算模块,用于基于均值以及协方差矩阵,计算每个传感器对监测目标的估计概率;第二处理模块,用于对所有估计概率进行归一化处理,得到目标估计概率;第一加入模块,用于在目标估计概率大于预设概率阈值的情况下,将与目标估计概率对应的传感器加入至监测目标对应的目标传感器集合。Optionally, the screening unit includes: a first characterization module, used to characterize the current state estimation value as the mean of a preset Gaussian distribution, and to characterize the current profile estimation value as the covariance matrix of a preset Gaussian distribution; a fourth calculation module, used to calculate the estimated probability of each sensor for the monitored target based on the mean and the covariance matrix; a second processing module, used to normalize all estimated probabilities to obtain the target estimated probability; and a first adding module, used to add the sensor corresponding to the target estimated probability to the target sensor set corresponding to the monitored target when the target estimated probability is greater than a preset probability threshold.

可选地,融合单元包括:第五计算模块,用于计算所有与每个目标传感器对应的当前状态估计值的平均值,得到目标状态值;第三确定模块,用于基于当前轮廓估计值,确定监测目标的当前目标轮廓的第一半轴长值、第二半轴长值以及偏转角度值;第四确定模块,用于基于所有目标传感器的协方差矩阵,确定目标协方差矩阵;第五确定模块,用于基于目标协方差矩阵,确定目标偏转角度值;第六计算模块,用于基于目标偏转角度值,计算每个目标传感器的目标第一半轴长值以及目标第二半轴长值;第二输出模块,用于基于所有目标第一半轴长值以及目标第二半轴长值,得到目标半轴长值;第三输出模块,用于基于目标半轴长值以及目标偏转角度值,得到目标轮廓值。Optionally, the fusion unit includes: a fifth calculation module, used to calculate the average value of all current state estimation values corresponding to each target sensor to obtain a target state value; a third determination module, used to determine the first semi-axis length value, the second semi-axis length value and the deflection angle value of the current target contour of the monitored target based on the current contour estimation value; a fourth determination module, used to determine the target covariance matrix based on the covariance matrix of all target sensors; a fifth determination module, used to determine the target deflection angle value based on the target covariance matrix; a sixth calculation module, used to calculate the target first semi-axis length value and the target second semi-axis length value of each target sensor based on the target deflection angle value; a second output module, used to obtain the target semi-axis length value based on all target first semi-axis length values and target second semi-axis length values; a third output module, used to obtain the target contour value based on the target semi-axis length value and the target deflection angle value.

上述的融合装置还可以包括处理器和存储器,上述获取单元30,采样单元31,处理单元32,筛选单元33,融合单元34等均作为程序单元存储在存储器中,由处理器执行存储在存储器中的上述程序单元来实现相应的功能。The above-mentioned fusion device may also include a processor and a memory. The above-mentioned acquisition unit 30, sampling unit 31, processing unit 32, screening unit 33, fusion unit 34, etc. are all stored in the memory as program units, and the processor executes the above-mentioned program units stored in the memory to realize corresponding functions.

上述处理器中包含内核,由内核去存储器中调取相应的程序单元。内核可以设置一个或以上,通过调整内核参数来融合目标传感器集合中与每个目标传感器对应的当前状态估计值以及当前轮廓估计值,得到在当前时刻下监测目标的目标状态值以及目标轮廓值。The processor includes a kernel, which retrieves the corresponding program unit from the memory. One or more kernels can be set, and the current state estimation value and the current contour estimation value corresponding to each target sensor in the target sensor set are fused by adjusting the kernel parameters to obtain the target state value and the target contour value of the monitored target at the current moment.

上述存储器可能包括计算机可读介质中的非永久性存储器,随机存取存储器(RAM)和/或非易失性内存等形式,如只读存储器(ROM)或闪存(flash RAM),存储器包括至少一个存储芯片。The above-mentioned memory may include non-permanent memory in a computer-readable medium, random access memory (RAM) and/or non-volatile memory, such as read-only memory (ROM) or flash RAM, and the memory includes at least one storage chip.

本申请还提供了一种计算机程序产品,当在数据处理设备上执行时,适于执行初始化有如下方法步骤的程序:获取在前一时刻下各监测目标的目标状态的状态估计值以及目标轮廓的轮廓估计值,对监测目标的状态估计值以及轮廓估计值进行粒子采样,得到在当前时刻下监测目标的预测量测值的带权粒子,基于在当前时刻下监测目标的预测量测值的带权粒子,对每个传感器的量测数据进行处理,得到在当前时刻下与传感器对应的所有监测目标的当前状态估计值以及当前轮廓估计值,基于当前状态估计值以及当前轮廓估计值,筛选与每个监测目标对应的目标传感器集合,融合目标传感器集合中与每个目标传感器对应的当前状态估计值以及当前轮廓估计值,得到在当前时刻下监测目标的目标状态值以及目标轮廓值。The present application also provides a computer program product, which, when executed on a data processing device, is suitable for executing a program that is initialized with the following method steps: obtaining a state estimation value of a target state and a contour estimation value of a target contour of each monitored target at a previous moment, performing particle sampling on the state estimation value and the contour estimation value of the monitored target to obtain weighted particles of a predicted measurement value of the monitored target at a current moment, processing the measurement data of each sensor based on the weighted particles of the predicted measurement value of the monitored target at the current moment to obtain a current state estimation value and a current contour estimation value of all monitored targets corresponding to the sensor at the current moment, screening a set of target sensors corresponding to each monitored target based on the current state estimation value and the current contour estimation value, fusing the current state estimation value and the current contour estimation value corresponding to each target sensor in the target sensor set to obtain a target state value and a target contour value of the monitored target at the current moment.

根据本发明实施例的另一方面,还提供了一种计算机可读存储介质,计算机可读存储介质包括存储的计算机程序,其中,在计算机程序运行时控制计算机可读存储介质所在设备执行上述的基于多传感器的监测目标融合方法。According to another aspect of an embodiment of the present invention, a computer-readable storage medium is also provided, the computer-readable storage medium including a stored computer program, wherein when the computer program is running, the device where the computer-readable storage medium is located is controlled to execute the above-mentioned multi-sensor-based monitoring target fusion method.

根据本发明实施例的另一方面,还提供了一种电子设备,包括一个或多个处理器和存储器,存储器用于存储一个或多个程序,其中,当一个或多个程序被一个或多个处理器执行时,使得一个或多个处理器实现上述的基于多传感器的监测目标融合方法。According to another aspect of an embodiment of the present invention, there is also provided an electronic device, comprising one or more processors and a memory, wherein the memory is used to store one or more programs, wherein when the one or more programs are executed by one or more processors, the one or more processors implement the above-mentioned multi-sensor based monitoring target fusion method.

图4是根据本发明实施例的一种用于基于多传感器的监测目标融合方法的电子设备(或移动设备)的硬件结构框图。如图4所示,电子设备可以包括一个或多个(图4中采用402a、402b,……,402n来示出)处理器402(处理器402可以包括但不限于微处理器MCU或可编程逻辑器件FPGA等的处理装置)、用于存储数据的存储器404。除此以外,还可以包括:显示器、输入/输出接口(I/O接口)、通用串行总线(USB)端口(可以作为I/O接口的端口中的一个端口被包括)、网络接口、键盘、电源和/或相机。本领域普通技术人员可以理解,图4所示的结构仅为示意,其并不对上述电子装置的结构造成限定。例如,电子设备还可包括比图4中所示更多或者更少的组件,或者具有与图4所示不同的配置。FIG. 4 is a hardware structure block diagram of an electronic device (or mobile device) for a multi-sensor-based monitoring target fusion method according to an embodiment of the present invention. As shown in FIG. 4 , the electronic device may include one or more (402a, 402b, ..., 402n are used in FIG. 4 to illustrate) processors 402 (the processor 402 may include but is not limited to a processing device such as a microprocessor MCU or a programmable logic device FPGA), and a memory 404 for storing data. In addition, it may also include: a display, an input/output interface (I/O interface), a universal serial bus (USB) port (which may be included as one of the ports of the I/O interface), a network interface, a keyboard, a power supply and/or a camera. It can be understood by those skilled in the art that the structure shown in FIG. 4 is only for illustration and does not limit the structure of the above-mentioned electronic device. For example, the electronic device may also include more or fewer components than those shown in FIG. 4 , or have a configuration different from that shown in FIG. 4 .

上述本发明实施例序号仅仅为了描述,不代表实施例的优劣。The serial numbers of the above embodiments of the present invention are only for description and do not represent the advantages or disadvantages of the embodiments.

在本发明的上述实施例中,对各个实施例的描述都各有侧重,某个实施例中没有详述的部分,可以参见其他实施例的相关描述。In the above embodiments of the present invention, the description of each embodiment has its own emphasis. For parts that are not described in detail in a certain embodiment, reference can be made to the relevant descriptions of other embodiments.

在本申请所提供的几个实施例中,应该理解到,所揭露的技术内容,可通过其它的方式实现。其中,以上所描述的装置实施例仅仅是示意性的,例如所述单元的划分,可以为一种逻辑功能划分,实际实现时可以有另外的划分方式,例如多个单元或组件可以结合或者可以集成到另一个系统,或一些特征可以忽略,或不执行。另一点,所显示或讨论的相互之间的耦合或直接耦合或通信连接可以是通过一些接口,单元或模块的间接耦合或通信连接,可以是电性或其它的形式。In the several embodiments provided in this application, it should be understood that the disclosed technical content can be implemented in other ways. Among them, the device embodiments described above are only schematic. For example, the division of the units can be a logical function division. There may be other division methods in actual implementation. For example, multiple units or components can be combined or integrated into another system, or some features can be ignored or not executed. Another point is that the mutual coupling or direct coupling or communication connection shown or discussed can be through some interfaces, indirect coupling or communication connection of units or modules, which can be electrical or other forms.

所述作为分离部件说明的单元可以是或者也可以不是物理上分开的,作为单元显示的部件可以是或者也可以不是物理单元,即可以位于一个地方,或者也可以分布到多个单元上。可以根据实际的需要选择其中的部分或者全部单元来实现本实施例方案的目的。The units described as separate components may or may not be physically separated, and the components shown as units may or may not be physical units, that is, they may be located in one place or distributed on multiple units. Some or all of the units may be selected according to actual needs to achieve the purpose of the present embodiment.

另外,在本发明各个实施例中的各功能单元可以集成在一个处理单元中,也可以是各个单元单独物理存在,也可以两个或两个以上单元集成在一个单元中。上述集成的单元既可以采用硬件的形式实现,也可以采用软件功能单元的形式实现。In addition, each functional unit in each embodiment of the present invention may be integrated into one processing unit, or each unit may exist physically separately, or two or more units may be integrated into one unit. The above-mentioned integrated unit may be implemented in the form of hardware or in the form of software functional units.

所述集成的单元如果以软件功能单元的形式实现并作为独立的产品销售或使用时,可以存储在一个计算机可读取存储介质中。基于这样的理解,本发明的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的全部或部分可以以软件产品的形式体现出来,该计算机软件产品存储在一个存储介质中,包括若干指令用以使得一台计算机设备(可为个人计算机、服务器或者网络设备等)执行本发明各个实施例所述方法的全部或部分步骤。而前述的存储介质包括:U盘、只读存储器(ROM,Read-Only Memory)、随机存取存储器(RAM,Random Access Memory)、移动硬盘、磁碟或者光盘等各种可以存储程序代码的介质。If the integrated unit is implemented in the form of a software functional unit and sold or used as an independent product, it can be stored in a computer-readable storage medium. Based on this understanding, the technical solution of the present invention, in essence, or the part that contributes to the prior art, or all or part of the technical solution can be embodied in the form of a software product, and the computer software product is stored in a storage medium, including a number of instructions for a computer device (which can be a personal computer, a server or a network device, etc.) to perform all or part of the steps of the method described in each embodiment of the present invention. The aforementioned storage medium includes: U disk, read-only memory (ROM, Read-Only Memory), random access memory (RAM, Random Access Memory), mobile hard disk, magnetic disk or optical disk and other media that can store program codes.

以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也应视为本发明的保护范围。The above is only a preferred embodiment of the present invention. It should be pointed out that for ordinary technicians in this technical field, several improvements and modifications can be made without departing from the principle of the present invention. These improvements and modifications should also be regarded as the scope of protection of the present invention.

Claims (10)

1. A multi-sensor-based monitoring target fusion method, comprising:
acquiring a state estimation value of a target state of each monitoring target and a contour estimation value of a target contour at the previous moment;
sampling particles of the state estimation value and the contour estimation value of the monitoring target to obtain weighted particles of the predicted measurement value of the monitoring target at the current moment;
processing the measurement data of each sensor based on weighted particles of the predicted measurement values of the monitoring targets at the current time to obtain current state estimation values and current contour estimation values of all the monitoring targets corresponding to the sensors at the current time;
screening a target sensor set corresponding to each monitoring target based on the current state estimated value and the current contour estimated value;
And fusing the current state estimated value and the current contour estimated value corresponding to each target sensor in the target sensor set to obtain a target state value and a target contour value of the monitoring target at the current moment.
2. The fusion method according to claim 1, further comprising, after acquiring the state estimation value of the target state and the contour estimation value of the target contour of each monitoring target at the previous time, the steps of:
determining a state predicted value of the monitoring target at the current time based on the state estimated value of the monitoring target at the previous time;
determining a contour predicted value of the monitoring target at the current time based on the contour estimated value of the monitoring target at the previous time;
and obtaining a predicted measurement value of the monitoring target at the current time based on the state predicted value and the contour predicted value, wherein the predicted measurement value is a measurement value corresponding to the predicted state of the monitoring target at the current time.
3. The fusion method according to claim 2, wherein the step of sampling the state estimation value and the contour estimation value of the monitoring target to obtain weighted particles of the predicted measurement value of the monitoring target at the current time, includes:
Performing particle sampling on the state estimation value and the contour estimation value to generate a plurality of particles;
calculating an initial particle weight of each particle based on preset parameters;
normalizing each initial particle weight based on all the initial particle weights to obtain a target particle weight of each particle, wherein the particles carrying the target particle weights are characterized as weighted particles;
and updating the weighted particles to obtain weighted particles of the monitoring target associated with the predicted measurement value at the current moment.
4. The fusion method of claim 1, wherein prior to processing the measurement data for each sensor based on weighted particles of the predicted measurement value of the monitored target at the current time, further comprising:
acquiring the measurement data of each sensor, wherein the measurement data at least comprises: a set of measurement vectors, a total number of measurements;
based on the measurement data, constructing a target associated variable facing a monitoring target, wherein the target associated variable is used for representing a measurement index of a measurement instruction of a first preset bit number generated by a preset monitoring target at the current moment;
And constructing a measurement-oriented measurement related variable based on the measurement data, wherein the measurement related variable is used for representing measurement of a second preset bit number generated by a monitoring target mapped by the measurement index at the current time.
5. The fusion method according to claim 4, wherein the step of processing the measurement data of each sensor based on weighted particles of the predicted measurement values of the monitoring targets at the current time to obtain the current state estimation values and the current contour estimation values of all the monitoring targets corresponding to the sensors at the current time includes:
based on weighted particles of the predicted measurement value of the monitoring target at the current time and the measurement vector set, carrying out measurement evaluation on the target associated variable to obtain likelihood probability;
under the condition that the likelihood probability belongs to a preset probability threshold range, constructing a correlation factor graph based on the target correlation variable and the measurement correlation variable;
performing iterative computation based on the associated factor graph to obtain a first parameter value and a second parameter value;
updating the measurement data of the monitoring target of the sensor based on the first parameter value and the second parameter value to obtain an updating result;
And based on the updating result, obtaining a current contour estimated value of the monitoring target corresponding to the sensor at the current time, and based on the updating result, carrying out weighted calculation on the weighted particles to obtain a current state estimated value of the monitoring target corresponding to the sensor.
6. The fusion method according to claim 1, wherein the step of screening the target sensor set corresponding to each of the monitoring targets based on the current state estimation value and the current contour estimation value comprises:
characterizing the current state estimation value as a mean value of a preset Gaussian distribution, and characterizing the current contour estimation value as a covariance matrix of the preset Gaussian distribution;
calculating the estimated probability of each sensor to the monitoring target based on the mean value and the covariance matrix;
normalizing all the estimated probabilities to obtain target estimated probabilities;
and adding the sensor corresponding to the target estimated probability to the target sensor set corresponding to the monitoring target under the condition that the target estimated probability is larger than a preset probability threshold.
7. The fusion method according to claim 6, wherein the step of fusing the current state estimation value and the current contour estimation value corresponding to each target sensor in the set of target sensors to obtain a target state value and a target contour value of the monitoring target at the current time, includes:
calculating the average value of all the current state estimation values corresponding to each target sensor to obtain the target state value;
determining a first half-axis length value, a second half-axis length value and a deflection angle value of a current target profile of the monitoring target based on the current profile estimation value;
determining a target covariance matrix based on the covariance matrices of all the target sensors;
determining a target deflection angle value based on the target covariance matrix;
calculating a target first half-axis length value and a target second half-axis length value of each target sensor based on the target deflection angle values;
obtaining a target half-axis length value based on all the target first half-axis length values and the target second half-axis length values;
and obtaining the target profile value based on the target half-axis length value and the target deflection angle value.
8. A multi-sensor based monitoring target fusion device, comprising:
the acquisition unit is used for acquiring a state estimation value of a target state of each monitoring target and a contour estimation value of a target contour at the previous moment;
the sampling unit is used for sampling particles of the state estimation value and the contour estimation value of the monitoring target to obtain weighted particles of the predicted measurement value of the monitoring target at the current moment;
the processing unit is used for processing the measurement data of each sensor based on weighted particles of the predicted measurement value of the monitoring target at the current time to obtain current state estimation values and current contour estimation values of all the monitoring targets corresponding to the sensors at the current time;
a screening unit, configured to screen a target sensor set corresponding to each monitoring target based on the current state estimation value and the current contour estimation value;
and the fusion unit is used for fusing the current state estimated value and the current contour estimated value corresponding to each target sensor in the target sensor set to obtain a target state value and a target contour value of the monitoring target at the current time.
9. A computer readable storage medium, characterized in that the computer readable storage medium comprises a stored computer program, wherein the computer program, when run, controls a device in which the computer readable storage medium is located to perform the multi-sensor based monitoring target fusion method according to any one of claims 1 to 7.
10. An electronic device comprising one or more processors and a memory for storing one or more programs, wherein the one or more programs, when executed by the one or more processors, cause the one or more processors to implement the multi-sensor based monitoring target fusion method of any of claims 1-7.
CN202211649534.8A 2022-12-21 2022-12-21 Multi-sensor based monitoring target fusion method, device, and electronic equipment Pending CN116049759A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202211649534.8A CN116049759A (en) 2022-12-21 2022-12-21 Multi-sensor based monitoring target fusion method, device, and electronic equipment

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202211649534.8A CN116049759A (en) 2022-12-21 2022-12-21 Multi-sensor based monitoring target fusion method, device, and electronic equipment

Publications (1)

Publication Number Publication Date
CN116049759A true CN116049759A (en) 2023-05-02

Family

ID=86117310

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211649534.8A Pending CN116049759A (en) 2022-12-21 2022-12-21 Multi-sensor based monitoring target fusion method, device, and electronic equipment

Country Status (1)

Country Link
CN (1) CN116049759A (en)

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20120194873A1 (en) * 2011-01-27 2012-08-02 Guoyi Fu Image registration parameters and confidence estimation from sensor data
CN112146648A (en) * 2020-09-23 2020-12-29 河北工业大学 A multi-target tracking method based on multi-sensor data fusion
CN112991397A (en) * 2021-04-19 2021-06-18 深圳佑驾创新科技有限公司 Traffic sign tracking method, apparatus, device and storage medium
DE102020116026A1 (en) * 2020-06-17 2021-12-23 Bayerische Motoren Werke Aktiengesellschaft Method and device for determining and / or tracking the contour of an object
CN115310526A (en) * 2022-07-28 2022-11-08 杭州电子科技大学 Sensor data rapid fusion method based on Sigma point belief propagation

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20120194873A1 (en) * 2011-01-27 2012-08-02 Guoyi Fu Image registration parameters and confidence estimation from sensor data
DE102020116026A1 (en) * 2020-06-17 2021-12-23 Bayerische Motoren Werke Aktiengesellschaft Method and device for determining and / or tracking the contour of an object
CN112146648A (en) * 2020-09-23 2020-12-29 河北工业大学 A multi-target tracking method based on multi-sensor data fusion
CN112991397A (en) * 2021-04-19 2021-06-18 深圳佑驾创新科技有限公司 Traffic sign tracking method, apparatus, device and storage medium
CN115310526A (en) * 2022-07-28 2022-11-08 杭州电子科技大学 Sensor data rapid fusion method based on Sigma point belief propagation

Similar Documents

Publication Publication Date Title
CN111562571B (en) A Maneuvering Multi-target Tracking and Track Maintenance Method with Unknown New Strength
CN113688287A (en) Enterprise risk infection path analysis method, device, equipment and storage medium
CN112348011B (en) Vehicle damage assessment method and device and storage medium
CN108197795B (en) Malicious group account identification method, device, terminal and storage medium
CN118761520A (en) A method and device for predicting the progress of engineering projects based on neural network
Sheng et al. Network traffic anomaly detection method based on chaotic neural network
Sun et al. Air Quality Prediction and Multi-Task Offloading based on Deep Learning Methods in Edge Computing
CN114462925B (en) Abnormal inventory asset identification method, device and terminal equipment
CN116049759A (en) Multi-sensor based monitoring target fusion method, device, and electronic equipment
CN116070149A (en) Data analysis method, system and cloud platform based on artificial intelligence
CN119377803B (en) Construction method of application efficiency evaluation model, application efficiency evaluation method and device
CN119917796B (en) A Spatiotemporal Graph Neural Network-Based Water Quality Prediction Method Based on DEM Mapping
CN118432952B (en) Anomaly detection method, electronic device and storage medium in zero-trust environment
CN108898067A (en) Determine the method, apparatus and computer readable storage medium of people and the object degree of association
Casaer et al. Analysing space use patterns by Thiessen polygon and triangulated irregular network interpolation: a non-parametric method for processing telemetric animal fixes
Cai CAMs as Shapley Value-based Explainers
CN113256422B (en) Method and device for identifying bin account, computer equipment and storage medium
Alguliyev et al. Weighted clustering for anomaly detection in big data
CN106156470B (en) Time series abnormity detection and labeling method and system
CN115327503A (en) Fixed single-station passive positioning method based on Gaussian particle filtering and related device
CN116823003A (en) Blockchain digital resource processing method, device, computer equipment and storage medium
US20230300046A1 (en) Packet information analysis method and network traffic monitoring device
CN118761557B (en) Water-Energy-Food System Nexus Analysis Method
Leite e Lopes et al. Applying complex network metrics to individual-tree diameter growth modeling.
CN120579706B (en) A digital collectibles grading and virtual-real interactive trading system based on multi-technology integration

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination