CN110849372B - Underwater multi-target track association method based on EM clustering - Google Patents

Underwater multi-target track association method based on EM clustering Download PDF

Info

Publication number
CN110849372B
CN110849372B CN201911188834.9A CN201911188834A CN110849372B CN 110849372 B CN110849372 B CN 110849372B CN 201911188834 A CN201911188834 A CN 201911188834A CN 110849372 B CN110849372 B CN 110849372B
Authority
CN
China
Prior art keywords
track
target
clustering
model
underwater
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201911188834.9A
Other languages
Chinese (zh)
Other versions
CN110849372A (en
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.)
Harbin Engineering University
Original Assignee
Harbin Engineering University
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 Harbin Engineering University filed Critical Harbin Engineering University
Priority to CN201911188834.9A priority Critical patent/CN110849372B/en
Publication of CN110849372A publication Critical patent/CN110849372A/en
Application granted granted Critical
Publication of CN110849372B publication Critical patent/CN110849372B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01CMEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C21/00Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
    • G01C21/20Instruments for performing navigational calculations
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V9/00Prospecting or detecting by methods not provided for in groups G01V1/00 - G01V8/00
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00Pattern recognition
    • G06F18/20Analysing
    • G06F18/23Clustering techniques

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Remote Sensing (AREA)
  • General Physics & Mathematics (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Data Mining & Analysis (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Theoretical Computer Science (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Evolutionary Biology (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Artificial Intelligence (AREA)
  • Geophysics (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Automation & Control Theory (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

The invention provides an underwater multi-target track association method based on EM clustering. Preprocessing the data; constructing and constructing a track quality grading model by introducing an information entropy, and preferentially performing track association matching by taking a point with good quality as a center when a track is associated; establishing a Gaussian mixture model through the topological information among the tracks, and simultaneously taking the points with good quality as the mass center to obtain a Gaussian probability density function; establishing a mixed integer nonlinear programming model, and reducing association deviation by using a recursive idea; and solving an extreme value of the unknown parameters through EM clustering to set a correlation judgment threshold, simultaneously carrying out uniqueness processing on the result, solving to obtain a corresponding relation of the underwater target track, and finally matching the track of the underwater target. Under the conditions of different target numbers, different sensor angle and distance measurement errors, different sensor detection probabilities and the like, the method has good positive correlation rate and has certain superiority and robustness.

Description

Underwater multi-target track association method based on EM clustering
Technical Field
The invention relates to an underwater multi-target track association method.
Background
The ocean resources of the earth are very abundant, and the ocean area accounts for about seventy percent of the total surface area of the earth. However, as the development of the method is slower in the exploration and utilization link of the ocean, a large development space is left for ocean resources. The development of the ocean is undoubtedly the most important in the development of each country in the world nowadays, but before the ocean resources are developed and utilized, the ocean needs to be relatively comprehensively understood through a large amount of ocean data, and the process needs a series of ocean development support technologies including track correlation technology. Meanwhile, the detection and tracking technology of underwater targets is greatly improved and the space is improved, the targets can be detected in time and accurately tracked, the development and utilization of the ocean are very important, meanwhile, in the field of multi-sensor information fusion, track association is a premise and a basis of multi-sensor information fusion, the purpose is to judge whether tracks reported by different sensors originate from the same target, when the sensors have monitoring areas, missed detection and random false alarms which are not identical, the reported targets of the sensors are not identical, so that a corresponding track does not exist in a track set reported by another sensor, and the difficulty of the original complex track association problem is increased. In addition, system errors generally exist in the detection process of the sensor and are influenced by various factors, so that deviation exists between the position state estimation of the target and the real target position, and the performance of the traditional track correlation algorithm is seriously deteriorated. If the target tracks can be accurately correlated under the complex environment with multiple noise points under water, the method is a premise and a basis for target tracking under high level water. If the high-resolution track association can be carried out in the multi-target track association, the subsequent link influence caused by noise and false alarm can be reduced, a new target can be perceived by finding a new track, and the influence caused by track omission is reduced.
Common algorithms adopted by the multi-sensor track association include Bayes (Bayes) estimation, least square estimation and evidence combination theory (Dempster-Shafer, D-S) in estimation algorithm aspect; BP neural network series in the aspect of artificial intelligence algorithm, GA genetic algorithm optimization algorithm, weighted fusion algorithm and the like. The Bayesian estimation is the earliest algorithm, and has great limitation because the prior measurement variance of the sensor is needed, so that the algorithm is difficult to be widely applied. Although the D-S evidence theory does not need prior information and needs fewer conditions than Bayes reasoning, the combination condition of evidence is very strict and cannot be applied to the premise background of evidence conflict and evidence basic probability determination, and a plurality of scholars have researched the D-S evidence theory problem based on evidence conflict at present. Although the information fusion algorithm based on the neural network type has good parallel processing capability and stability, when the number of layers of the network is large, the convergence cannot be fast, and the global optimal solution is difficult to obtain. The information fusion method for optimizing by using the genetic algorithm has the advantages of implicit parallel capability, strong robustness and the like, but the time complexity of the algorithm is high. Meanwhile, for the conditions of high noise and more false alarms of the underwater environment, an algorithm which can still accurately and effectively solve the problem of track association of underwater multiple targets when a certain system error exists in the sensor is needed.
Disclosure of Invention
The invention aims to provide an underwater multi-target track association method based on EM clustering, which has better precision and tolerance.
The purpose of the invention is realized as follows:
the method comprises the following steps: carrying out classification pretreatment on the track, and simultaneously carrying out classification according to the track quality entropy;
step two: the GMM model is established by the method,
2.1 grading the Pre-treatment uncertainty level α j (k) Level 1 track set X B The locus set X of other levels is regarded as the centroid of Gaussian model in GMM A A set of sample points considered to be GMM;
2.2 mixing different data points by Gaussian distribution, calculating the mean vector, covariance matrix and mixing weight in GMM model, obtaining model probability density expression, carrying out vector estimation on sensor deviation, and carrying out vector estimation according to bias vector eta k Predicting the covariance of the deviation vector by the followed dynamic model to obtain the optimal corresponding deviation estimation;
step three: evaluating the maximum likelihood and establishing an MINLP model for recursion;
step four: and performing EM clustering, and finally realizing accurate target classification and identification.
The present invention may further comprise:
1. in step one, the set of quality exclusion events for the tracked trace is: and (4) according to the track scanning state, whether the collection is in the extrapolation time, the track state estimation stable condition and the track stable condition, processing the idea of introducing weights into the four event sets to obtain track quality entropy for grading.
2. In the second step, a GMM model is established, corresponding movement is carried out on the center of mass of the GMM to the sample point set according to the neighborhood topological structure, if the distance between the final center of mass and the sample point is smaller, the higher the association degree between the tracks is, and after the optimal matching relation is obtained in a certain mode, the matching association relation between the track sets is obtained by utilizing posterior probability.
3. In the second step, a Gaussian radial basis function is obtained:
Figure GDA0003937185210000021
X A a track sample point set is obtained, and K is a dimension; delta 2 Introducing a uniform distribution for the covariance in the Gaussian model
Figure GDA0003937185210000022
Then, a probability density function is obtained:
Figure GDA0003937185210000023
where ω is a uniformly distributed weight coefficient.
4. In the second step, the sensor deviation is estimated by vector, and the offset vector of the sensor is represented by eta k Is expressed as eta k Different for each sensor, the offset vector η is for the time-varying case k Following the following dynamic model η k =F k-1,η η k-1 +w k-1,η In the formula F k-1,η Is a transition matrix, and w k-1,η Is a zero mean; an initial bias estimate and corresponding covariance are
Figure GDA0003937185210000031
For
Figure GDA0003937185210000032
Obtaining an optimal corresponding bias estimate using maximum likelihood rules
Figure GDA0003937185210000033
WhereinU is a correspondence matrix.
5. The maximum likelihood is evaluated in the third step, the solution is realized by a linear binary distribution problem and a least square continuous optimization problem, and the estimation is based on the current deviation
Figure GDA0003937185210000034
Determining a correspondence matrix
Figure GDA0003937185210000035
i, j are track points and are based on passing current
Figure GDA0003937185210000036
Calculating a deviation estimate
Figure GDA0003937185210000037
Recursion is continued until the matrix and bias estimates no longer change.
6. The correlation method in the fourth step is that the EM clustering algorithm completes the correlation matching of the underwater target track, wherein the EM algorithm comprises the step E and the step M,
e, calculating the posterior probability of the feature vector to the GMM model on the basis of the initial parameters,
Figure GDA0003937185210000038
m is that the posterior probability obtained by E is used to respectively calculate the extreme value of the unknown parameter set, if the target track i monitored on the sensor A and the target track j monitored on the sensor B come from the same underwater target, then there is a target
Figure GDA0003937185210000039
In the formula:
Figure GDA00039371852100000310
is to carry out the judgment threshold of the track association and simultaneously to carry out
Figure GDA00039371852100000311
The largest corresponding trajectory matches.
The invention relates to the technical field of underwater multi-target track association, in particular to a multi-sensor underwater multi-target track association method based on EM clustering application recursion idea modeling. Compared with the prior art, the invention has the advantages that: a. the underwater environment is complex and changeable, the noise interference is serious, the accuracy of a single sensor is low, and in order to improve the probability and the accuracy of detection and tracking, the invention adopts multiple sensors for correlation; b. according to the method, by a preprocessing method of the track grading condition, the characteristic information representing the track quality is used for establishing the entropy grade, so that the entropy weight is obtained for track association, the track fusion accuracy is effectively improved, the risk that the track quality is reduced after fusion, the track fusion robustness is enhanced, and the fusion precision is improved. d. The method has the advantages that the MINLP is used for recursion, deviation estimation is reduced, and the tolerance can be effectively improved in the environment with multiple underwater noise points and high false alarm. The method can solve the problem of poor anti-noise capability of the traditional underwater target correlation, can effectively improve the underwater multi-target track correlation accuracy rate, and has certain applicability.
Drawings
FIG. 1 is a flow chart of the method of the present invention;
FIG. 2 is a simulation of a simulation experiment of the present method;
FIG. 3 is a graph comparing performance of an upper limit algorithm for sensor angle measurement error;
FIG. 4 is a graph comparing the performance of the sensor range error ceiling algorithm;
FIG. 5 is a graph comparing performance of different detection probability algorithms of a sensor;
FIG. 6 is a comparison graph of performance of an algorithm for detecting the number of different targets by a sensor;
FIG. 7 is a plot of root mean square error for different sensor goniometric system errors;
FIG. 8 is a plot of the root mean square error under different sensor ranging system errors.
Detailed Description
The invention is described in more detail below by way of example.
With reference to fig. 1, the specific steps of the present invention are as follows:
(1) Track grading pretreatment for underwater target track
The set of quality exclusion events for the traced trajectory is: h is H A {H A1 The track status is scanned and changed, H A2 The scanning condition of the track is unchanged }, H B {H B1 Track is within extrapolation time, H B2 Track outside of maximum extrapolation time }, H C {H C1 Trajectory state estimation is stable, H C2 Trajectory state estimation jitter, H D {H D1 Stable track, H D2 Track hunting }.
The event set H can reflect the track condition, i.e. the track quality, in a period of time, the above four event sets are all repulsive, and whether the track state detected by the sensor is stable and the deviation degree is known through the four event sets. The situation of the quality entropy is judged by the frequency of the occurrence of the above H event set. While preprocessing the data for trajectory correlation.
Entropy I of the quality classification of the trace points is
Figure GDA0003937185210000041
In the formula C nn Normalized weight for the nth event; m nn The total number of sub-events for the nth event; p (m) represents the probability of the mth sub-event of the nth event. Assuming that the k time has O tracks in common, the mass entropy E (k) { E) is obtained 1 (k),E 2 (k),...,E U (k) And normalizing the obtained product, namely:
Figure GDA0003937185210000042
in the formula E j (k) For the tracking quality entropy of the trace j at time k, j ∈ {1,2, \8230;,O}。
k time trajectory j uncertainty level alpha j (k) The method is divided into the following steps: e is not less than 0 j (k) 1 when the content is less than or equal to 0.2; e is more than 0.2 j (k) 2 when the content is less than or equal to 0.5; e is more than 0.5 j (k) When the content is less than or equal to 1, the content is 3. Therefore, data preprocessing is graded, the lower the grade number is, the better the tracking quality is, and the higher the grade number is, the worse the tracking quality is, so that in the process of track association, a point with good quality can be taken as a center to carry out track association matching.
(2) Construction of GMM model
Ranking the pre-processing to a level of uncertainty α j (k) Level 1 track set X B The locus set X of other levels is regarded as the centroid of Gaussian model in GMM A Considering the sample point set of the GMM, their gaussian radial basis functions are given by:
Figure GDA0003937185210000051
in the formula: k is dimension; delta 2 Is the covariance in the gaussian model. Meanwhile, when the sensor is in a noisy state, a system error state and the like, the target of the two sensors is inconsistent, namely a non-homologous track appears in the track set, so that a uniform distribution is brought:
Figure GDA0003937185210000052
thus obtaining the track
Figure GDA0003937185210000053
The GMM probability density function of (a) is:
Figure GDA0003937185210000054
in the formula: ω is a uniformly distributed weight coefficient. Since the smaller the difference of the neighborhood topology of the trajectory, the higher the probability that they are from the same target, the higher the corresponding weight ratio in the GMM, and the probability of the proportional weight of each gaussian component is described as:
Figure GDA0003937185210000055
in the formula:
Figure GDA0003937185210000056
the resulting GMM probability density function is therefore:
Figure GDA0003937185210000057
(3) Evaluating maximum likelihood function and performing MINLP modeling recursion
(1) Vector estimation of sensor bias
The offset vector of the sensor can be represented by eta k Is expressed as eta k May be different for each sensor. For time-varying cases, the offset vector η k Following the following dynamic model
η k =F k-1,η η k-1 +w k-1,η
In the formula F k-1,η Is a transition matrix, and w k-1,η Is a zero mean; the initial bias estimate and corresponding covariance are
Figure GDA0003937185210000061
(2) Evaluating a maximum likelihood function
If one wants to optimize the model, the first solution is the evaluation of the likelihood function, for which the two parts P (z) of the likelihood function are considered separately k,η |U kk ) And P (Z) k ,z k,η |U kk ) And (4) evaluating.
1) Evaluation of P (z) k,η |U kk )
For
Figure GDA0003937185210000062
Under the Gaussian assumption, the likelihood function z k,η Can be represented by the following formula
Figure GDA0003937185210000063
In the formula n η Denotes a deviation parameter, P (z) k,η |U kk ) And matrix U k Are independent of each other.
2) Evaluation of P (Z) k ,z k,η |U kk ),
Figure GDA0003937185210000064
Figure GDA0003937185210000065
Figure GDA0003937185210000066
Is a single point distance from the point at which,
Figure GDA0003937185210000067
as deviation of the track
A target density over the monitored volume defined by the target number per unit volume is beta t The monitoring probability of the sensor is
Figure GDA0003937185210000068
β t Can be regarded as beta t =1/4D 2 Where D represents the average target distance from its nearest target point. The likelihood function relating the two local orbits can be evaluated, and the local orbit estimate can be obtained as
Figure GDA0003937185210000069
Maximizing the likelihood function P (Z) k ,z k,η |U kk ) Corresponding to minimizing its negative logarithm. Bonding with
Figure GDA00039371852100000610
Comprises the following steps:
Figure GDA00039371852100000611
while removing constants not related to parameter estimation to obtain the following optimization model
Figure GDA0003937185210000071
This problem then becomes the MINLP problem. It can be solved by a linear binary allocation problem and a least-squares continuous optimization problem.
(3) Recursive MINLP
Estimation of bias based on current
Figure GDA0003937185210000072
Determining a correspondence matrix
Figure GDA0003937185210000073
And (3) eliminating the influence of invariant parameters to obtain an optimization model as follows:
Figure GDA0003937185210000074
by the current time
Figure GDA0003937185210000075
Calculating a deviation estimate
Figure GDA0003937185210000076
Figure GDA0003937185210000077
Obtaining a corresponding covariance:
Figure GDA0003937185210000078
the MINLP model can be recursive, the pseudo-code is as follows:
inputting the deviation estimated value of the current time k
Figure GDA0003937185210000079
1) By the present deviation estimation
Figure GDA00039371852100000710
Calculating a correspondence matrix
Figure GDA00039371852100000711
2) In view of the current matrix
Figure GDA00039371852100000712
Updating bias estimates
Figure GDA00039371852100000713
3) Repeating (1) and (2) until the corresponding matrix is not changed any more.
4) Outputting a final deviation estimate
Figure GDA00039371852100000714
And a corresponding matrix
Figure GDA00039371852100000715
The transformation relationship between sets of point sets may be expressed as follows:
X A =X B +v(X B )
in the formula v (X) n ) For the offset function, a regularization term is added to allow the trajectory set to move as a whole. Therefore, a regularization function is added to the regenerated kernel Hilbert space
Figure GDA0003937185210000081
Obtain an expectation of a set of parameters and a set of sample points as
Figure GDA0003937185210000082
In the formula:
Figure GDA0003937185210000083
{ω,δ 2 v is the unknown parameter set;
Figure GDA0003937185210000084
is the posterior probability; λ is the weight coefficient of the regularization term. Describing the form of the v function by maximizing it according to the variational method
Figure GDA0003937185210000085
Wherein: omega j Is a weight coefficient matrix W = [ omega ] 12 ,...,ω m ] T An element of (1);
Figure GDA0003937185210000086
is a Gaussian kernel matrix; beta is a smoothness coefficient. Thus, the matching incidence relation among the target track point set is obtained and is X A =X B + GW, while substituting into the above formula
Figure GDA0003937185210000087
In the formula: g j,n Is a row vector of a gaussian kernel matrix.
(4) Performing EM clustering to perform underwater target track association
At the expected step, the posterior probability is known as follows according to Bayesian theorem:
Figure GDA0003937185210000088
in the formula: j =1,2. If a posterior probability
Figure GDA0003937185210000089
Higher, this means that the probability that the object trajectories i and j are from the same object is higher, while the matrix
Figure GDA0003937185210000091
Is the associated probability matrix of the set of trajectories. Therefore, the probability that non-homologous trajectories can be obtained is:
Figure GDA0003937185210000092
in the maximization step, in order to maximize Q, the unknown parameter sets are respectively extremized. Solution (II)
Figure GDA0003937185210000093
Figure GDA0003937185210000094
So as to obtain the composite material,
W=[diag(P1)G+λδ 2 I] -1 (PX A -diag(P1)X B )
Figure GDA0003937185210000095
Figure GDA0003937185210000096
in the formula: 1 is a full 1 column vector; and I is an identity matrix.
From the above equation, if the target trajectory i monitored on sensor a and the target trajectory j monitored on sensor B are from the same underwater target, then there are:
Figure GDA0003937185210000097
in the formula:
Figure GDA0003937185210000098
is the judgment threshold for performing the track association. Meanwhile, because one-to-many error association conditions exist, in order to ensure the uniqueness of the association result, the method can be used for solving the problem that the error association between the two or more error association conditions exists
Figure GDA0003937185210000099
The largest corresponding trajectory matches.
In order to verify the effectiveness of the EM clustered underwater multi-target track association method, N underwater target track association problems are detected by using two sensors in a global rectangular coordinate system, and a simulation diagram is shown in FIG. 2.
The effectiveness of each algorithm under different angle measurement errors of a sensor is detected, and fig. 3 is an algorithm effectiveness comparison graph under different angle measurement system errors.
FIG. 4 is a comparison graph of performance of the proposed algorithm under different ranging system errors, and when the ranging system errors become large, the accuracy rate of the FFT and REP algorithms decreases significantly, and the non-homologous trajectory is wrongly associated. The algorithm reduces the correlation deviation through MINLP modeling recursion, can better perform track correlation and simultaneously has excellent tolerance and can keep a relatively stable correlation rate under the condition of increasing system errors in an underwater environment with multiple noise points and a high false alarm environment.
FIG. 5 is a comparison graph of performance of the proposed algorithm under different detection probabilities, and we only change the detection probability of the sensor under the condition that other conditions are not changed, so that the detection probability is changed at intervals of 0.5-1, and other parameters are not changed. As can be seen from the figure, when the detection probability of the sensor is small, the algorithm has obvious advantages in positive correlation rate and strong tolerance.
FIG. 6 is a positive correlation efficiency graph of an algorithm for the number of different underwater targets, and since the algorithm is adaptive to a threshold after passing through neighborhood topology information and an optimization model, when non-homologous tracks increase, the algorithm can still identify the correlation matching relationship of each target track, and the influence of the non-homologous tracks on a matching result is reduced.
It can be known from fig. 7 and 8 that the estimation accuracy of the algorithm of the present invention is more accurate than that of the other three algorithms in terms of the system error of the sensor, the average estimation deviation between the angle measurement and the distance measurement is small, and when the system error between the angle measurement and the distance measurement is fixed, the experimental result shows that the algorithm of the present invention has higher accuracy, because the algorithm of the present invention continuously performs model optimization on the deviation in the MINLP recursion, and reduces the influence caused by the sensor deviation, the accuracy of the algorithm of the present invention is higher than that of the other three algorithms in the continuous iteration process.
TABLE 1 Algorithm average run time Table(s)
Figure GDA0003937185210000101
As can be seen from the running time of the algorithms in table 1, the average running time of each algorithm gradually becomes higher as the number of targets increases. The FFT algorithm continuously and alternately iterates the target track position, and the time consumption is shortest. The algorithm needs to be updated iteratively all the time in the clustering process of the EM algorithm when the track association is carried out, so the running time of the algorithm is relatively long. However, under the condition that the time is similar to the REP algorithm FCM algorithm, the tolerance of the algorithm is superior to the three algorithms, although the running time of the algorithm is greater than that of the FFT algorithm, under the condition that the underwater track correlation noise is high, certain requirements are required on the accuracy and the tolerance of the track correlation, and the tolerance of the FFT algorithm obviously does not meet the actual requirements, so that the algorithm has certain superiority in the track correlation for underwater targets.

Claims (5)

1. An underwater multi-target track association method based on EM clustering is characterized by comprising the following steps:
the method comprises the following steps: carrying out classification pretreatment on the track: grading according to the track quality entropy;
step two: establishing a GMM model;
2.1 grading the preprocessing into a trajectory set X with uncertainty level 1 B The locus set X of other levels is regarded as the centroid of Gaussian model in GMM A A set of sample points considered to be GMM; obtaining a gaussian radial basis function:
Figure FDA0003937185200000011
X A a track sample point set is obtained, and K is a dimension; delta. For the preparation of a coating 2 Introducing a uniform distribution for covariance in Gaussian model
Figure FDA0003937185200000012
Then, a probability density function is obtained:
Figure FDA0003937185200000013
wherein ω is a uniformly distributed weight coefficient;
2.2 mixing different data points by Gaussian distribution, calculating the mean vector, covariance matrix and mixing weight in GMM model, obtaining model probability density expression, carrying out vector estimation on sensor deviation, and carrying out vector estimation according to bias vector eta k Predicting the covariance of the deviation vector by the followed dynamic model to obtain the optimal corresponding deviation estimation;
step three: evaluating the maximum likelihood and establishing an MINLP model for recursion;
for two parts of maximum likelihood P (z) k,η |U kk ) And P (Z) k ,z k,η |U kk ) Evaluating, converting the optimized dynamic model into an MINLP problem, and recursion is carried out on the MINLP: estimation of bias based on current
Figure FDA0003937185200000014
Determining a correspondence matrix
Figure FDA0003937185200000015
To trace points, based on passing current
Figure FDA0003937185200000016
Calculating a deviation estimate
Figure FDA0003937185200000017
Recursion is continued until the matrix and the bias estimate no longer change; adding regularization function in regeneration kernel Hilbert space
Figure FDA0003937185200000018
Enabling the track set to move integrally;
step four: and performing EM clustering, and finally realizing accurate target classification and identification.
2. The underwater multi-target track association method based on EM clustering as claimed in claim 1, wherein: in step one, the set of quality exclusion events for the tracked trace is: and processing the idea of introducing weights into the four event sets according to the track scanning state, the fact that whether the collection is within extrapolation time, the track state estimation stability condition and the track stability condition, and grading the obtained track quality entropy.
3. The underwater multi-target track association method based on EM clustering as claimed in claim 2, wherein: in the second step, a GMM model is established, corresponding movement is carried out on the center of mass of the GMM to the sample point set according to the neighborhood topological structure, if the distance between the final center of mass and the sample point is smaller, the higher the association degree between the tracks is, and after the optimal matching relation is obtained, the matching association relation between the track sets is obtained by utilizing posterior probability.
4. The underwater multi-target track association method based on EM clustering as claimed in claim 3, wherein: in step two, the sensor deviation is subjected to vector estimation, and the offset vector of the sensor uses eta k Is expressed as eta k Different for each sensor, the offset vector eta for the time-varying case k Following the following dynamic model η k =F k-1,η η k-1 +w k-1,η In the formula F k-1,η Is a transition matrix, w k-1,η Is a zero mean; an initial bias estimate and corresponding covariance are
Figure FDA0003937185200000021
For
Figure FDA0003937185200000022
Obtaining an optimal corresponding bias estimate using maximum likelihood rules
Figure FDA0003937185200000023
Wherein U is k Is a correspondence matrix.
5. The underwater multi-target track association method based on EM clustering as claimed in claim 4, wherein: the correlation method in the fourth step is that the EM clustering algorithm completes correlation matching of underwater target tracks, wherein the EM algorithm comprises the step E and the step M;
e, calculating the posterior probability of the feature vector to the GMM model on the basis of the initial parameters,
Figure FDA0003937185200000024
m is that the posterior probability obtained by E is used to respectively calculate the extreme value of the unknown parameter set, if the target track i monitored on the sensor A and the target track j monitored on the sensor B come from the same underwater target, then there is a target
Figure FDA0003937185200000025
In the formula:
Figure FDA0003937185200000026
is to carry out the judgment threshold of the track association and simultaneously to carry out the judgment of the track association
Figure FDA0003937185200000027
The largest corresponding trajectory matches.
CN201911188834.9A 2019-11-28 2019-11-28 Underwater multi-target track association method based on EM clustering Active CN110849372B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911188834.9A CN110849372B (en) 2019-11-28 2019-11-28 Underwater multi-target track association method based on EM clustering

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911188834.9A CN110849372B (en) 2019-11-28 2019-11-28 Underwater multi-target track association method based on EM clustering

Publications (2)

Publication Number Publication Date
CN110849372A CN110849372A (en) 2020-02-28
CN110849372B true CN110849372B (en) 2023-02-14

Family

ID=69605830

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911188834.9A Active CN110849372B (en) 2019-11-28 2019-11-28 Underwater multi-target track association method based on EM clustering

Country Status (1)

Country Link
CN (1) CN110849372B (en)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112598061B (en) * 2020-12-23 2023-05-26 中铁工程装备集团有限公司 Clustering and grading method for tunnel surrounding rock
CN112906746B (en) * 2021-01-25 2024-04-02 北京工业大学 Multi-source track fusion evaluation method based on structural equation model
CN113793327B (en) * 2021-09-18 2023-12-26 北京中科智眼科技有限公司 Token-based high-speed rail foreign matter detection method
CN114219022A (en) * 2021-12-13 2022-03-22 南京理工大学 Multi-sensor multi-target tracking method combining cluster analysis and particle swarm optimization algorithm
CN116381607B (en) * 2023-04-11 2023-10-27 哈尔滨工程大学 Multi-target water-striking sound characteristic association method

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102724631A (en) * 2012-06-28 2012-10-10 上海交通大学 Position information updating method for position-based routing of vehicular ad hoc network
CN106959618A (en) * 2017-05-05 2017-07-18 国网山东省电力公司电力科学研究院 A kind of voltage control method for coordinating for optimizing weight based on ladder
CN107066806A (en) * 2017-02-15 2017-08-18 中国人民解放军海军航空工程学院 Data Association and device
CN108286971A (en) * 2017-10-18 2018-07-17 北京航空航天大学 A kind of forecast Control Algorithm that the Inspector satellite based on the optimization of MIXED INTEGER second order cone is evaded
CN109858526A (en) * 2019-01-08 2019-06-07 沈阳理工大学 Sensor-based multi-target track fusion method in a kind of target following
CN110188951A (en) * 2019-05-30 2019-08-30 重庆大学 A kind of method for building up of the optimizing scheduling of the brick field ferry bus based on least energy consumption
CN110361744A (en) * 2019-07-09 2019-10-22 哈尔滨工程大学 RBMCDA underwater multi-target tracking based on Density Clustering

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9977123B2 (en) * 2014-05-20 2018-05-22 Bae Systems Information And Electronic Systems Integration Inc. Automated track projection bias removal using frechet distance and road networks

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102724631A (en) * 2012-06-28 2012-10-10 上海交通大学 Position information updating method for position-based routing of vehicular ad hoc network
CN107066806A (en) * 2017-02-15 2017-08-18 中国人民解放军海军航空工程学院 Data Association and device
CN106959618A (en) * 2017-05-05 2017-07-18 国网山东省电力公司电力科学研究院 A kind of voltage control method for coordinating for optimizing weight based on ladder
CN108286971A (en) * 2017-10-18 2018-07-17 北京航空航天大学 A kind of forecast Control Algorithm that the Inspector satellite based on the optimization of MIXED INTEGER second order cone is evaded
CN109858526A (en) * 2019-01-08 2019-06-07 沈阳理工大学 Sensor-based multi-target track fusion method in a kind of target following
CN110188951A (en) * 2019-05-30 2019-08-30 重庆大学 A kind of method for building up of the optimizing scheduling of the brick field ferry bus based on least energy consumption
CN110361744A (en) * 2019-07-09 2019-10-22 哈尔滨工程大学 RBMCDA underwater multi-target tracking based on Density Clustering

Non-Patent Citations (16)

* Cited by examiner, † Cited by third party
Title
3D Multi-Drone-Cell Trajectory Design for Efficient IoT Data Collection;Weisen Shi等;《ICC 2019 - 2019 IEEE International Conference on Communications (ICC)》;20190715;全文 *
An outer-approximation algorithm for generalized maximum entropy sampling;Han-Lim Choi等;《2008 American Control Conference》;20080805;全文 *
Hybridization of Nonlinear and Mixed-Integer Linear Programming for Aircraft Separation With Trajectory Recovery;Jérémy Omer等;《IEEE Transactions on Intelligent Transportation Systems》;20131231;全文 *
joint track-to-track association and sensor registration at the track level;HongyanZhu等;《Digital SignalProcessing41(2015)》;20150323;48-59 *
Mixed Integer Nonlinear Programming Framework for Fixed Path Coordination of Multiple Underwater Vehicles Under Acoustic Communication Constraints;Pramod Abichandani等;《IEEE Journal of Oceanic Engineering》;20150112;第40卷(第4期);864-873 *
地空协同防空目标抗差跟踪算法;崔亚奇等;《航空学报》;20140430;第35卷(第04期);1079-1090 *
基于m-best数据关联和小轨迹关联多目标跟踪算法;谷晓琳等;《***工程与电子技术》;20170731;第39卷(第07期);1640-1646 *
基于松弛标号算法的多传感抗差航迹关联;朱洪艳等;《控制与决策》;20150430;第30卷(第04期);593-598 *
基于航迹矢量分级聚类的雷达与电子支援措施抗差关联算法;李保珠等;《电子与信息学报》;20190630;第41卷(第06期);1310-1316 *
基于高斯混合模型的航迹抗差关联算法;李保珠等;《航空学报》;20190630;第40卷(第6期);226-234 *
复杂场景中运动目标的检测与跟踪;丁婕;《中国优秀硕士学位论文全文数据库 (信息科技辑)》;20120615;I138-1947 *
复杂背景下多目标跟踪技术研究;李伟;《中国优秀硕士学位论文全文数据库 (信息科技辑)》;20190115;I138-3322 *
多目标跟踪中联合概率数据关联优化算法;廖辉荣等;《计算机仿真》;20111130;第28卷(第11期);14-18 *
差分进化算法及应用研究;吴亮红;《中国优秀硕士学位论文全文数据库 (信息科技辑)》;20070415;I138-39 *
混合整数线性规划形式的抗差状态估计方法;陈艳波等;《电力自动化设备》;20150731;第35卷(第07期);26-31+49 *
过失误差识别和数据校正的MILP模型的新视角;孙少超等;《华东理工大学学报(自然科学版)》;20110131;第37卷(第01期);95-99 *

Also Published As

Publication number Publication date
CN110849372A (en) 2020-02-28

Similar Documents

Publication Publication Date Title
CN110849372B (en) Underwater multi-target track association method based on EM clustering
CN109188423B (en) Distributed multi-target tracking method based on multi-source clustering
CN103729859B (en) A kind of probability nearest neighbor domain multi-object tracking method based on fuzzy clustering
Clark et al. Multi-target state estimation and track continuity for the particle PHD filter
CN109508000A (en) Isomery multi-sensor multi-target tracking method
CN104156984A (en) PHD (Probability Hypothesis Density) method for multi-target tracking in uneven clutter environment
CN105137418A (en) Multi-object tracking and data interconnection method based on whole neighborhood fuzzy clustering
CN112098992B (en) Multi-hypothesis multi-target track starting method based on grid clustering
CN111929653B (en) Target detection and tracking method and system based on unmanned ship navigation radar
CN114509750A (en) Water target tracking method based on multi-navigation radar
CN113076686A (en) Aircraft trajectory prediction method based on social long-term and short-term memory network
CN117075631A (en) Detection tracking identification method based on unmanned aerial vehicle cluster target
CN113378998B (en) Stratum lithology while-drilling identification method based on machine learning
CN109671096A (en) A kind of space-time neighbour target detection and Grid Clustering measure more extension method for tracking target under dividing
CN113836818B (en) Ocean current motion prediction algorithm based on BP neural network prediction model
Dai et al. Clustering of DOA data in radar pulse based on SOFM and CDbw
CN115130523A (en) Flight target behavior intention prediction method based on hidden Markov model
CN115345245A (en) Tropical cyclone rapid reinforcement forecasting method based on random forest and transfer learning
CN115015908A (en) Radar target data association method based on graph neural network
CN109581305B (en) Multi-radar error correction method based on historical data
CN112308229A (en) Dynamic multi-objective evolution optimization method based on self-organizing mapping
Aoki et al. Suboptimal JPDA for tracking in the presence of clutter and missed detections
CN111811515A (en) Multi-target track extraction method based on Gaussian mixture probability hypothesis density filter
Xie et al. Moving Object Detection Algorithm Based on Adaptive Clustering
Chu et al. A Multi-Target Track Association Algorithm in Underwater Multi-Sensors Environments

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
GR01 Patent grant
GR01 Patent grant