Design of Insect Target Tracking Algorithm in Clutter Based on the Multidimensional Feature Fusion Strategy

: Entomological radar is an effective means of monitoring insect migration, and can realize long-distance and large-scale rapid monitoring. The stable tracking of individual insect targets is the basic premise underlying the identiﬁcation of insect species and the study of insect migration mechanisms. However, the complex motion trajectory and large number of false measurements decrease the performance of insect target tracking. In this paper, an insect target tracking algorithm in clutter was designed based on the multidimensional feature fusion strategy (ITT-MFF). Firstly, multiple feature parameters of measurements were fused to calculate the membership of measurements and target, thereby improving the data association accuracy in the presence of clutter. Secondly, a distance-correction factor was introduced to the probabilistic data association (PDA) algorithm to accomplish multi-target data association with a low computational cost. Finally, simulation scenarios with different target numbers and clutter densities were constructed to verify the effectiveness of the proposed method. The tracking result comparisons of the experimental data acquired from a Ku-band entomological radar also indicate that the proposed method can effectively reduce computational cost while maintaining high tracking precision, and is suitable for engineering implementation.


Introduction
Insects are the most abundant and economically important group of terrestrial migrants [1]. They transfer enormous quantities of parasites, propagules, nutrients, and energy between regions, with substantial effects on the essential ecosystem [2]. However, some migratory pests can cause serious plant diseases [3] and virus transmission [4], creating great losses for human society. According to statistics, in China, pests cause an average of 17.5 billion kilograms of lost crop production annually [5]. Observations of migratory flight have utility for crop protection that aim to predict damaging pests and reduce their impacts [6]. Radar, having the advantages of a wide detection range all-time, is an effective technology to realize the large-scale and rapid monitoring of flying insects [7]. Meanwhile, the tracking of insect targets via high-resolution radar is of great significance for the study of insect migration mechanisms [8]. However, for the insect target tracking in cluttered environments [9], complex flight trajectories (such as parallel and crossing trajectories) and large numbers of false measurements not only make it difficult to accurately determine the sources of the measurements but also increase the calculation burden.
The Bayesian multi-target probability data association algorithm is widely used in engineering due to its superior tracking performance [10]. Multiple hypothesis tracking (MHT) is recognized theoretically as the optimal approach in Bayesian tracking, in which the hypotheses are propagated into the future in anticipation that subsequent data will resolve the uncertainty [11]; however, this process requires high computational complexity. Furthermore, the joint probabilistic data association algorithm (JPDA), where a track is updated by a weighted sum of all measurements in its gate, can maintain excellent performance in the case of dense clutter and missed detections [12]. Nevertheless, the intensive computations limit its practical applications.
In terms of reducing computational complexity, in [13], a cheap JPDA algorithm was designed using an ad hoc formula to simplify association probability calculations. On this basis, a suboptimal JPDA algorithm was developed [14], in which joint event probability is simplified by assuming that the probability of detection is one, or very nearly one, using partial joint events. Although cheap JPDA and suboptimal JPDA can greatly reduce computational complexity, the tracking performance degrades in the presence of dense clutter. Furthermore, from the perspective of fuzzy mathematics theory, a joint probabilistic data association algorithm based on all-neighbor fuzzy clustering (ANFCJPDA) was developed [15]. In this algorithm, the predicted position of the target is set as a cluster center, and the association probabilities are calculated according to fuzzy clustering. In [16], a rough set probabilistic data association (RS-PDA) algorithm was presented, in which rough set theory is introduced to judge the origination of measurements in the validation regions, and the measurements in the intersection region can then be processed discriminatingly. The above two algorithms can maintain a good tracking performance in the presence of clutter, but their time consumption is only slightly less than that of JPDA. Therefore, it is difficult to achieve a balance between computational cost and tracking performance in the existing optimized JPDA algorithm, which makes it difficult to ensure excellent performance in insect target tracking scenarios.
Another challenging problem for multitarget tracking is differentiating between measurements arising from the target of interest and measurements originating from other target returns or clutter. As a result, accurate judgments of measurement origin would help enhance tracking accuracy. In addition to target position measurements, the Doppler component [17] is widely used in combination with kinetic information as an additional discriminant of measurement origin. However, when target Doppler is included in the measurement vector, a nonlinear filter process is required to deal with the nonlinearity between the target's kinematic state and Doppler measurement, which brings additional computation. Moreover, in many studies, several non-kinematic features were utilized to assist in determining the measurement source. In [18,19], four different types of sensor data were introduced for target tracking: kinematic measurements (target position and velocity), features (such as radar cross section (RCS), signal strength, and wing span), attributes (such as the number of aircraft engines and type of radar), and categorical features (such as wing span, with only a few types of targets). In [20], target class information was integrated into the data association process, and the association results was improved when the kinematic likelihoods were similar for different targets. However, for migrating insect groups, the individuals in the group are usually of the same category, potentially invalidating the class information. Moreover, amplitude is also a vital information to discriminate against false measurements and further improve tracking performance. In [21], the target RCS was used to track multiple Rayleigh targets. Further, in [22], based on the probabilistic data association filter with amplitude information [23], appropriate amplitude likelihoods were specified to cope with K-distributed clutter. Nevertheless, in the above algorithm, only a single-dimensional feature is utilized, making it difficult to achieve excellent performance in complex insect target tracking scenarios.
For a rotating-polarization radar [24], which is generally used to observe insect migration, many feature parameters related to the insect target characteristics can be measured, besides the RCS parameters. Therefore, in this paper, an insect target tracking algorithm in clutter was designed based on the multidimensional feature fusion strategy, which is mainly aimed at data association processing, one of the most important aspects of tracking. For the echo data collected by high resolution and fully polarimetric entomological radar, multidimensional feature parameters of the measurements are extracted, and then these features are fused based on the proposed multidimensional feature fusion strategy. The obtained membership of the multidimensional feature fusion can effectively improve the Remote Sens. 2021, 13, 3744 3 of 20 accuracy of measurement-to-track associations. Further, in this work, a distance-correction factor was introduced to PDA to modify the association probability of the measurements falling into the crossing gate, thereby avoiding the split of association probability matrix in JPDA and ensuring the multitarget tracking performance under the premise of low computation.
The remainder of this paper is organized as follows. In Section 2, the multiple target tracking problem is described and the problems in the insect target tracking scenario are analyzed. Section 3 presents the principles of the proposed insect target tracking algorithm based on the multidimensional feature fusion strategy. Section 4 introduces comparisons of simulations and experimental data results between the proposed method and traditional methods. The discussion and future work are presented in Section 5, and the conclusions are given in Section 6.

Multiple Target Tracking Problem Formation
Assume that the number of targets in the surveillance is n t , and the dynamic and measurement models for target t, t= 1, 2, · · · , n t , are defined as where x t (k) denotes the state vector of target t at time k, and z t (k) denotes the measurement vector. F t (k) is the state transition matrix and H t (k) is the measurement transition matrix. G t (k) is the process noise distribution matrix. v t (k) and w t (k) are zero-mean mutually dependent white Gaussian noise vectors with covariance matrices Q t (k) and R t (k), respectively. Suppose the number of validated measurements at time k is m k . The validated measurement set at time k is denoted as . Generally, JPDA [12] is the preferred method to handle multitarget tracking in clutter. The following is a brief introduction to the process of multi-target tracking based on JPDA.
The state prediction and the measurement prediction at time k are defined aŝ The predicted covariance is defined as The innovation covariance is The filter gain is A validation gate is designed to select the validated measurements of the target t, and is represented as where the value of parameter g can be obtained from tables of chi-square distribution based on the probability P G that the true measurements will fall in the elliptical gate. On this basis, the following m k × (n t + 1) validation matrix Ω can be defined as follows: where ω jt refers to the binary elements that indicate whether measurement j lies in the validation gate for target t. Index t = 0 indicates "no target" and the corresponding column of Ω has all units, which represents each measurement may originate from clutter or a false alarm. The state estimation for target t iŝ e t (k) is combined innovation and defined as β jt e jt (k) (11) where e jt (k) = z j (k) −ẑ t (k|k − 1). β jt is the association probability that the validated measurement j originates from target t at time k and is calculated by where θ i (k) is the ith feasible event, and θ i jt (k) is the event that the measurement j originates from target t in the ith feasible event. Z k {Z(1), Z(2), · · · , Z(k)} denotes the cumulative measurement set up to time k, and n k is the number of the feasible events.
For multi-target tracking based on JPDA, the probabilities of all feasible joint events should be calculated. As the number of targets and clutter density increase, the joint events also increase exponentially, which result in huge amounts of calculations. Moreover, only the position distribution is considered in Pr θ i (k) Z k . Consequently, it is difficult to achieve excellent tracking performance for insect target tracking in a cluttered background.

Insect Target Tracking Algorithm Based on Multidimensional Feature Fusion
To improve the tracking accuracy of insect targets in cluttered environments with low computational cost, we designed an insect target tracking algorithm based on the multidimensional feature fusion strategy. The flow diagram is shown in Figure 1. Based on the framework of the probabilistic data association algorithm, two key modules are added to better deal with insect target tracking scenarios, as shown in the grey boxes in Figure 1. On the one hand, a multidimensional feature fusion strategy based on fuzzy logic synthesis is proposed to calculate the membership of multi-dimensional feature fusion between the measurements and the existing target track, thereby improving the accuracy of measurement-to-track associations. On the other hand, a distance correction factor is introduced to correct the association probability of the measurements falling into the crossing gate, enabling the PDA [23] to handle multitarget tracking with a low computational load. Finally, the association probability of PDA is modified based on the membership of multidimensional feature fusion and the distance correction factor. Remote Sens. 2021, 13, x FOR PEER REVIEW

Multidimensional Feature Fusion Strategy
A key problem of multitarget tracking in cluttered environments is the measurement-to-track data associations, which means it is difficult to corre the sources of measurements by only relying on kinematic information (targ velocity). According to the previous literature [24,25], it can be concluded polarization pattern and scattering matrix (SM) are closely related to the ins and composition, which can be measured using a fully polarimetric radar. Th can measure all four linear-polarization receive-transmit (HH, HV, VH, VV) combinations. The corresponding feature parameters reflect the inherent ch the insect target and are beneficial to distinguish the target and clutter meas Therefore, a multidimensional feature fusion strategy based on fuzzy theory is proposed to obtain the membership of the multidimensional featu flow diagram of this strategy is shown in Figure 2. It mainly includes fou as shown in the grey boxes. Firstly, the feature parameters of measurement in real time. Secondly, a new membership function, the fitting probabili function (PDF) ratio of the target and clutter feature parameters, is design the characteristics of insect targets. Meanwhile, the analytic hierarchy pr employed to calculate the weight of each feature. In the end, multidimen fusion is realized by fuzzy logic synthesis.

Multidimensional Feature Fusion Strategy
A key problem of multitarget tracking in cluttered environments is the uncertainty of measurement-to-track data associations, which means it is difficult to correctly determine the sources of measurements by only relying on kinematic information (target position and velocity). According to the previous literature [24,25], it can be concluded that an insect's polarization pattern and scattering matrix (SM) are closely related to the insect's geometry and composition, which can be measured using a fully polarimetric radar. This type of radar can measure all four linear-polarization receive-transmit (HH, HV, VH, VV) backscattering combinations. The corresponding feature parameters reflect the inherent characteristics of the insect target and are beneficial to distinguish the target and clutter measurements.
Therefore, a multidimensional feature fusion strategy based on fuzzy logic synthesis theory is proposed to obtain the membership of the multidimensional feature fusion. The flow diagram of this strategy is shown in Figure 2. It mainly includes four key modules, as shown in the grey boxes. Firstly, the feature parameters of measurements are extracted in real time. Secondly, a new membership function, the fitting probability distribution function (PDF) ratio of the target and clutter feature parameters, is designed to adapt to the characteristics of insect targets. Meanwhile, the analytic hierarchy process (AHP) is employed to calculate the weight of each feature. In the end, multidimensional feature fusion is realized by fuzzy logic synthesis. as shown in the grey boxes. Firstly, the feature parameters of measuremen in real time. Secondly, a new membership function, the fitting probabili function (PDF) ratio of the target and clutter feature parameters, is design the characteristics of insect targets. Meanwhile, the analytic hierarchy pr employed to calculate the weight of each feature. In the end, multidime fusion is realized by fuzzy logic synthesis.

Feature Parameter Exaction
The feature parameters extracted by the polarization pattern and SM are represented as a 0 , a 1 , d and v, and their definitions and physical meanings are as follows.
According to [26], the insect's polarization pattern σ(α) can be expressed as where α is the direction of linear polarization, and a 0 represents the average of the RCS over all polarization angles (360 • ). θ 1 and θ 2 are parameters of theoretical model in natural coordinates. a 1 and a 2 are magnitudes of harmonic modulations representing elongated component and cruciform component respectively. The backscatter of a target is defined by the scattering matrix. For monostatic radar, the general form of the insect SM is modeled as follows [26,27]: s 12 e jβ s 12 e jβ s 22 e jγ (14) where s 11 , s 12 , and s 22 represent the square roots of the RCSs, σ xx , σ xy , and σ yy (unit: m 2 ), respectively. σ xx and σ yy are the maximum of the copolar polarization pattern and the value in the orthogonal direction, respectively, and σ xy is the cross polarisation element magnitudes. They can be calculated from a 0 , a 1 , and a 2 [26]. β and γ are the phases related to target properties. The eigenvalue of S is defined as λ s1 and λ s2 . For fully polarimetric radar, by transmitting orthogonally polarized signals, the inset SM elements can be measured from the echo signals of different polarization directions. According to [26], d and v are the invariant target parameters of SM, which can be calculated by the Graves matrix. The Graves matrix [28] is defined as where superscript H represents combined conjugate and transpose operations. λ 1 and λ 2 are the two eigen values of the Graves matrixDet G [29,30]. Assuming that λ 1 ≥ λ 2 , then where The determinant of the Graves matrix can be expressed as the product of its two eigenvalues The determinant is represented by variable d: An "eigen parameter" v is defined as where, φ represents the phase difference between λ s1 and λ s2 .

Membership Function Definition
The membership function [31] is the core element when using fuzzy synthesis theory to solve practical problems. Commonly used membership functions include normal distribution and Γ distribution. Considering the characteristics of insect targets, a new membership function is proposed. For each feature, the ratio of the target and the clutter feature parameter's fitting PDF is defined as the membership function. Correspondingly, the membership µ p of the pth feature is expressed as where, pd f T p is the target PDF of the pth feature, pd f Cp is the clutter PDF of the pth feature, and u p is the feature parameter value calculated by the measurement information. As a result, the statistical models for the target and clutter measurement feature parameters need to be known a priori to achieve multi-dimensional feature fusion. For this reason, the PDF of the feature parameters defined in Section 3.1.1 was modeled based on the experimental data collected by the fully polarimetric entomological radar in Ku-band in August 2019 in Lancang, Yunnan province, China. The experimental scenario and parameters are detailed in Section 4.3.1. The measurements, after undertaking constant false alarm rate (CFAR) detection process, included not only kinetic information but also the feature parameters extracted in Section 3.1.1. Here, the number of target trajectory samples is 5000, and the number of measurements of each trajectory is greater than 100. The ground clutter data were collected by the radar at different elevation angles and distances.
Next, based on the least squares fitting criterion and Kolmogorov-Smirnov (K-S) fitting goodness test methods [32], four typical distribution models (Lognormal, Weibull, Gamma, and Normal distribution) were employed to fit the statistical characteristics of the four feature parameters. The least squares fitting error and K-S test parameters are shown in Figures 3 and 4, respectively. Then, on the basis of the fitting results, the fitting model of the four feature parameters was selected, as shown in Table 1. Furthermore, Figure 5 shows the PDF fitting curves for the different feature parameters of insect targets and clutter. It can be seen that the extracted features present obvious distinction between the target and clutter measurements.          (c) (d) Figure 5. PDF fitting curves for the different feature parameters of insect targets and clutter: (a) d ,

Feature Weight Assignment
As shown in Figure 5, the overlapping areas of target and clutter PDF are different for different features. The smaller the overlap area is, the better the distinction. Therefore, it is necessary to define a scheme to achieve feature weight assignment.
The analytic hierarchy process is an effective way to determine the weight assignment [33,34] and is widely employed to quantitatively describe the weight value of each feature. The first step of the AHP method is to evaluate the relative importance of each pair of factors and to build a pairwise comparison matrix. In this section, the overlapping area of the PDF fitting model of the target and clutter feature parameters is used to measure the discrimination of those features, expressed as d  Table 2.

Feature Weight Assignment
As shown in Figure 5, the overlapping areas of target and clutter PDF are different for different features. The smaller the overlap area is, the better the distinction. Therefore, it is necessary to define a scheme to achieve feature weight assignment.
The analytic hierarchy process is an effective way to determine the weight assignment [33,34] and is widely employed to quantitatively describe the weight value of each feature. The first step of the AHP method is to evaluate the relative importance of each pair of factors and to build a pairwise comparison matrix. In this section, the overlapping area of the PDF fitting model of the target and clutter feature parameters is used to measure the discrimination of those features, expressed as A d , A v , A a 0 , and A a 1 . Then the ratios of the PDF overlapping areas between multiple features are taken as the corresponding elements of the comparison matrix, as shown in Table 2.
Based on the comparison matrix, an arithmetic mean method is adopted to calculate the assigned weight b p of each feature. Specifically, we use the arithmetic mean value of all column vectors to calculate the weight vector. The expression is as follows: where A pj is the element in the comparison matrix, and n is the dimension of matrix, i.e., the number of feature parameters.

Fuzzy Logic Synthesis
Fuzzy logic synthesis [35] is implemented using the weighted average method. Consequently, the multidimensional feature fusion membership of each measurement j that falls into the target t validation gate at time k can be obtained, as follows: where m k,t is the number of measurements in the target t validation gate, and C 0 is the normalization factor. µ jp is the pth feature membership of measurement j and target t, which can be calculated by (22).

Distance-Correction Factor Calculation
For the scenario of insect target tracking, the JPDA, due to its combinatorial nature, comprises the bulk of the computational load in multitarget tracking processing. The PDA uses all of the validated measurements with different probabilities to update the target states, but it does not take into account the measurements in the crossing gate and is also not suitable for multitarget tracking. Therefore, here, we introduce a distance-correction factor into PDA to modify the association probability of the measurement in the crossing gate, enabling the algorithm to better handle multiple target tracking and achieve a compromise between the quantity of calculations and the tracking performance. The calculation steps for the distance-correction factor are described as follows: Firstly, the target set is divided based on the validation matrix defined in (9). Assume that n j is the number of validation gates that measurement j lies in, i.e., n j = n t ∑ t=1 ω jt . n j > 1 indicates that measurement j might have originated from multiple targets. In this case, the corresponding target index set is defined as T_index = t ii t satis f y ω jt = 1 ii = 1, 2, · · · , n j , n j ≤ n t .
Then, the Mahalanobis distance [36] is employed to measure the distance between measurement j and the center of the validation gate t, the size of which measures the competition degree between multiple crossing gates and the involved measurements. The mathematical expression is Here, the smaller the distance, the higher the relevance degree between measurement j and target t. Therefore, the distance-correction factor between measurement j and target t is defined as , t ∈ T_index

Association Probability Matrix Modification
The number of candidate measurements of target t is denoted as m k,t , and the validation measurement set originating from target t at time k can be obtained based on the validation matrix denoted as Z t (k): The association probability that measurement j originated from target t, as calculated by the PDA algorithm [23], is defined as where, θ jt (k) is the feasible event that the measurement j originates from target t, and Z k−1 t denotes the measurement set up to frame k − 1. However, in a cluttered environment, the number of measurements is usually greater than the number of targets. Thus, the source of measurements within the crossing gate needs to be considered. Moreover, only the position distribution is considered in (28). It is, therefore, difficult to deal with complex insect target tracking.
Thus, in the proposed ITT-MFF, the association probability β jt , defined in (28), is modified by the multidimensional feature fusion membership in (24) and the distancecorrection factor in (26). The utilization of more dimensional features is beneficial to improving the accuracy of measurement-to-track associations, and the introduction of the distance-correction factor can realize multi-target tracking under the premise of a low computational cost. The modified expression is defined as

Simulations and Experimental Results
In this section, the tracking performance of the proposed insect target tracking algorithm is first evaluated in the simulation scenarios with different target numbers and clutter densities. Then, based on the real observational data acquired from the Ku-band entomological radar, the tracking results in a cluttered environment also verify the effectiveness of the proposed method.

Tracking Performance Evaluation Criteria
Three evaluation criteria were considered to evaluate the tracking performance of the algorithms.

• Success Tracking Ratio (STR)
If a track can maintain more than N T frames and ensure that the mean of the tracking error is less than the standard deviation of measurement noise σ R , the trajectory is considered to be tracked successfully, thereby satisfying: wherex(k) is the state estimation at time k, and x(k) is the real value. N T is the number of scan frames.
• Optimal Subpattern Assignment (OSPA) Statistic The OSPA statistic [37] gives an evaluation of the performance of a tracking algorithm. The metric includes two parameters, c and p, where c limits the penalty assigned for false or lost tracks, and p determines the penalty assigned for estimates that are not close to any true target.
The OSPA statistic between the set of estimated state vectorsx(k|k) and the set of true target states x(k) can be expressed as where represents the distance betweenx i (k|k) and x π(i) (k);x i (k|k) is the i-th state ofx(k|k); Π n represents the set of all permutations of the true target states; n is the number of measurements; and m is the number of targets. In the following simulation, we set the values as p = 1 and c 2 = 10, 000 m 2 .
After M times Monte Carlo simulation runs, the mean OSPA (MOSPA) is defined by taking the average of the errors, which is expressed as • Running Time To assess the computational requirements of the proposed method, the running time (denoted by T R ) is obtained by computing the average CPU time in MATLAB 2017 on an Intel Core 2.40 GHz two-CPU computer operating under Windows 10.

Simulation Scenarios
In this section, a crossing target scenario is simulated. The number of targets is set to two, four, and six, as shown in Figure 6. The targets move towards each other, cross at a certain moment, and then move away from each other. In each simulation scenario, the three clutter density values are set to 10 −5 , 10 −4 , and 10 −3 , respectively. The tracking filter parameters are shown in Table 3. The PDF of the target and the clutter feature parameters satisfy the statistical model in Section 3.1.2.
Remote Sens. 2021, 13, x FOR PEER REVIEW 13 of 21 assigned for false or lost tracks, and p determines the penalty assigned for estimates that are not close to any true target.
The OSPA statistic between the set of estimated state vectors ( ) | x k k and the set of true target states ( ) xk can be expressed as

Simulation Scenarios
In this section, a crossing target scenario is simulated. The number of targets is set to two, four, and six, as shown in Figure 6. The targets move towards each other, cross at a certain moment, and then move away from each other. In each simulation scenario, the three clutter density values are set to 10 −5 , 10 −4 , and 10 −3 , respectively. The tracking filter parameters are shown in Table 3. The PDF of the target and the clutter feature parameters satisfy the statistical model in Section 3.1.2.

Result Comparisons
Based on the performance evaluation criteria listed in Section 4.1, the proposed ITT-MFF algorithm is compared with several probability data association algorithms (JPDA [12], cheap JPDA [13] and suboptimal JPDA [14]) in different simulation scenarios. At the same time, for fairness, a comparison with the ITT-NoMFF (only the distancecorrection factor is introduced) algorithm is added to verify the performance in the absence of multidimensional feature fusion. In addition, when the clutter density is 10 −3 and the number of targets is six, the result of the JPDA algorithm is not shown due to the high memory requirements and time-consuming.

• STR
The comparison results of the algorithm's STR in different simulation scenarios are shown in Table 4. For JPDA, suboptimal JPDA, cheap JPDA, and ITT-NoMFF, without the aid of multidimensional features, the STR decreases with an increase of clutter density. With an increase in the number of targets, the magnitude of this decline increases. Notably, when the clutter density is 10 −3 and the number of targets is greater than four, the STR drops to less than 50%. However, the ITT-NoMFF algorithm is superior to cheap JPDA and suboptimal JPDA, and is basically the same as the JPDA algorithm, which verifies the validity of the distance-correction factor. Besides, the proposed ITT-MFF algorithm is robust in cluttered background due to the use of multidimensional feature fusion. Moreover, when the clutter density is 10 −3 , the successful tracking ratio of ITT-MFF can still be guaranteed as greater than 90%. Therefore, compared to other algorithms, ITT-MFF offers more stable tracking performance. • MOSPA The MOPSA metric is an effective means to evaluate the algorithm's tracking precision. To avoid the impact of unsuccessful tracking results on the evaluation in the process of multiple Monte Carlo simulations, only the tracking accuracy of successful tracking is evaluated. The MOPSA distance comparison results in different simulation scenarios are shown in Figures 7-9. It can be seen that with an increase in clutter density, the MOSPA distances of JPDA, suboptimal JPDA, Cheap JPDA, and ITT-NoMFF gradually increase. However, ITT-NoMFF outperforms the other three methods.
Remote Sens. 2021, 13, x FOR PEER REVIEW 15 of 21 the MOSPA distances of JPDA, suboptimal JPDA, Cheap JPDA, and ITT-NoMFF gradually increase. However, ITT-NoMFF outperforms the other three methods. When the multi-dimensional feature fusion strategy is introduced, the proposed ITT-MFF algorithm still guarantees high tracking accuracy with an increase in clutter density, which indicates that the extracted features can effectively improve the distinction between the target and clutter measurements.  Table 5. For the JPDA algorithm, as the number of targets increases, the running time increases rapidly. Moreover, the greater the clutter density is, the larger the growth rate. Such increases are unacceptable in practical applications. Since cheap JPDA and suboptimal JPDA simplify the calculations of association probability, the running time is reduced greatly. However, the tracking precision in dense clutter is difficult to meet the demand. Although the running time of the proposed ITT-MFF algorithm increased slightly compared to the two When the multi-dimensional feature fusion strategy is introduced, the proposed ITT-MFF algorithm still guarantees high tracking accuracy with an increase in clutter density, which indicates that the extracted features can effectively improve the distinction between the target and clutter measurements. •

Running Time
Running time comparisons of different algorithms are shown in Table 5. For the JPDA algorithm, as the number of targets increases, the running time increases rapidly. Moreover, the greater the clutter density is, the larger the growth rate. Such increases are unacceptable in practical applications. Since cheap JPDA and suboptimal JPDA simplify the calculations of association probability, the running time is reduced greatly. However, the tracking precision in dense clutter is difficult to meet the demand. Although the running time of the proposed ITT-MFF algorithm increased slightly compared to the two When the multi-dimensional feature fusion strategy is introduced, the proposed ITT-MFF algorithm still guarantees high tracking accuracy with an increase in clutter density, which indicates that the extracted features can effectively improve the distinction between the target and clutter measurements. •

Running Time
Running time comparisons of different algorithms are shown in Table 5. For the JPDA algorithm, as the number of targets increases, the running time increases rapidly. Moreover, the greater the clutter density is, the larger the growth rate. Such increases are unacceptable in practical applications. Since cheap JPDA and suboptimal JPDA simplify the calculations of association probability, the running time is reduced greatly. However, the tracking precision in dense clutter is difficult to meet the demand. Although the running time of the proposed ITT-MFF algorithm increased slightly compared to the two When the multi-dimensional feature fusion strategy is introduced, the proposed ITT-MFF algorithm still guarantees high tracking accuracy with an increase in clutter density, which indicates that the extracted features can effectively improve the distinction between the target and clutter measurements.

Running Time
Running time comparisons of different algorithms are shown in Table 5. For the JPDA algorithm, as the number of targets increases, the running time increases rapidly. Moreover, the greater the clutter density is, the larger the growth rate. Such increases are unacceptable in practical applications. Since cheap JPDA and suboptimal JPDA simplify the calculations of association probability, the running time is reduced greatly. However, the tracking precision in dense clutter is difficult to meet the demand. Although the running time of the proposed ITT-MFF algorithm increased slightly compared to the two simplified algorithms, it will not, like JPDA, increase rapidly with an increase in clutter density or number of targets. Moreover, compared to ITT-NoMFF, the multidimensional feature fusion processing does not entail much extra computation. The engineering practicability of the proposed ITT-MFF algorithm is verified based on experimental field data. In August 2019, the echo data were collected using a Kuband high resolution and fully polarimetric entomological radar (shown in Figure 10) in Lancang, Yunnan province, China. The experimental scene is shown in Figure 11. The basic parameters of the radar system are shown in Table 6. simplified algorithms, it will not, like JPDA, increase rapid density or number of targets. Moreover, compared to ITT-N feature fusion processing does not entail much extra compu

Experiment Scenarios
The engineering practicability of the proposed ITT-MF on experimental field data. In August 2019, the echo data we high resolution and fully polarimetric entomological rad Lancang, Yunnan province, China. The experimental scene basic parameters of the radar system are shown in Table 6.    After CFAR detection processing, the target measurement sets were acquired, as shown in Figure 12. For this study, we selected six insect target trajectories, including parallel and crossing situations. To evaluate the algorithm's tracking performance in different clutter densities, a semi-physical simulation experiment [38] was adopted. We assumed that the positions of the clutter measurements were uniformly distributed around the target measurements, and set the clutter density to 10 −3 . The target feature parameters were obtained in real time based on the theoretical formulas in Section 3.1.1, and the clutter feature parameters were subjected to the PDF distribution model obtained in Section 3.  After CFAR detection processing, the target measurement sets were acquired, as shown in Figure 12. For this study, we selected six insect target trajectories, including parallel and crossing situations. To evaluate the algorithm's tracking performance in different clutter densities, a semi-physical simulation experiment [38] was adopted. We assumed that the positions of the clutter measurements were uniformly distributed around the target measurements, and set the clutter density to 10 −3 . The target feature parameters were obtained in real time based on the theoretical formulas in Section 3.1.1, and the clutter feature parameters were subjected to the PDF distribution model obtained in Section 3.1.2.

Result Comparisons
Based on the practical measurements of insect targets, a comparative analysis of various algorithms was conducted based on the filter result diagram and running time.
The filtering results are shown in Figure 13. It can be seen that because of the interference of the clutter points, the cheap JPDA, suboptimal JPDA, JPDA, and ITT-NoMFF algorithms all have false associations. Especially around the 70th frame, when the target trajectory crosses and the distance is small, the performance significantly degrades. However, the proposed ITT-MFF algorithm can effectively distinguish target and clutter measurements based on multidimensional feature fusion, thereby reducing the number of false associations and achieving stable tracking in a cluttered environment.

Discussion
For insect target tracking in a cluttered environment, complex flight trajectories (such Correspondingly, the running time comparison results are shown in Table 7. It can be seen that, although the running time of ITT-MFF is slightly higher than that of cheap JPDA and suboptimal JPDA, it is still much lower than that of JPDA.

Discussion
For insect target tracking in a cluttered environment, complex flight trajectories (such as parallel and crossing trajectories) and large numbers of false measurements not only make it difficult to accurately determine the origins of the measurements but also increase the calculation burden. The JPDA algorithm can maintain an excellent performance in the case of dense clutter, but the expensive computational cost is a critical problem in applying the algorithm. Therefore, determining how to achieve a tradeoff between computation and tracking accuracy is a primary problem in achieving insect target tracking. To this end, we designed an insect target tracking algorithm based on multidimensional feature fusion strategy.
The simulation experiments show that the time consumption of JPDA algorithm increases exponentially with an increase in the number of targets. Although the suboptimal JPDA and cheap JPDA offer high computing efficiency, their performance degrades seriously in crossing trajectory and densely cluttered environments. The main reason for this decrease in tracking performance is the difficulty in achieving the correct discrimination of measurement sources in clutter. The method proposed in this paper fully considers the scattering characteristics of insect targets and improves the data association accuracy in cluttered environments based on multidimensional feature fusion strategy. A distancecorrection factor was also introduced into the PDA algorithm to ensure the effectiveness of multi-target tracking processing. The experimental data processing results further show that, compared to other algorithms, this method can realize the stable tracking of insect targets with a low computational cost.
The accurate extraction of insect target trajectories is of great significance to the study of insect migration mechanisms. However, the research in this paper mainly focuses on the tracking algorithm for individual insect targets. In future research, we will consider using the movement characteristics of group targets to track larger numbers of insects. In addition, multi-sensor tracking systems will be employed to obtain more target information and deal with more complex tracking scenarios.

Conclusions
In this paper, an insect target tracking algorithm in clutter based on the multidimensional feature fusion strategy was proposed. Based on the simulation results under different scenarios, this method achieved a higher successful tracking rate and tracking accuracy, which means that multidimensional feature fusion can effectively improve the accuracy of data association. In addition, compared to the JPDA algorithm, the running time was greatly reduced and there was no exponential explosion in computations, which indicates that introducing the distance-correction factor avoided splitting the confirmation matrix, and also provided excellent performance in the crossing track. To summarize, the proposed method can guarantee an excellent tracking performance in cluttered background with a low computing cost and is suitable for insect target tracking scenarios.