Next Article in Journal
Displacement Sensor Based on a Small U-Shaped Single-Mode Fiber
Next Article in Special Issue
Convolutional Recurrent Neural Network-Based Event Detection in Tunnels Using Multiple Microphones
Previous Article in Journal
Application of an Improved Correlation Method in Electrostatic Gait Recognition of Hemiparetic Patients
Previous Article in Special Issue
Energy-Efficient Spatial Query-Centric Geographic Routing Protocol in Wireless Sensor Networks
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Target Identification Method for the Millimeter Wave Seeker via Correlation Matching and Beam Pointing

1
No. 203 Research Institute of China Ordnance Industries, Xi’an 710065, China
2
School of Computer Science, Shaanxi Normal University, Xi’an 710119, China
3
National Laboratory of Radar Signal Processing, Xidian University, Xi’an 710071, China
*
Author to whom correspondence should be addressed.
Sensors 2019, 19(11), 2530; https://doi.org/10.3390/s19112530
Submission received: 5 May 2019 / Revised: 29 May 2019 / Accepted: 31 May 2019 / Published: 3 June 2019
(This article belongs to the Special Issue Sensors In Target Detection)

Abstract

:
Target identification is a challenging task under land backgrounds for the millimeter wave (MMW) seeker, especially under complex backgrounds. Focusing on the problem, an effective method combining correlation matching and beam pointing is proposed in this paper. In the beginning, seeker scanning for target detection is conducted in two rounds, and target information of the detected targets is stored for correlation matching. Point or body feature judgment is implemented by using high resolution range profile (HRRP). Then, the error distribution zone is constructed with the beam pointing as the origin. In the end, we identify the target by searching the one which lies in the closest error distribution from the beam pointing center. The effectiveness of the proposed method is verified by using mooring test-fly and real flight data.

1. Introduction

The millimeter wave (MMW) seeker [1,2,3] is capable of working under all weather conditions, day and night, and has attracted increasing popularity. However, target identification for the MMW seeker is very difficult under land backgrounds, especially under complex backgrounds. Moreover, for the MMW seeker with strict limitations of the non-deviation flight trajectory, it is unable to image the target in two dimensionalities as synthetic aperture radar (SAR) [4,5,6]. That is to say, for small range tactical missiles, one has to realize target identification just by using the high-resolution range profile (HRRP) in the range direction [7,8,9]. As a result, the performance of target identification for the MMW seeker is not satisfying in plenty of application situations.
The key point for target identification under complex backgrounds is the false alarms caused by the clutters. How to eliminate or ease the negative influences by the clutter distractions is the precondition of successful target identification. Study of the land clutters is very complex, since it is influenced by various factors, such as geographical factors, climate factors, and radio waves, etc. Focusing on the problem, plenty of effective models have been proposed to model the distribution probability of the clutter. The essential point of these models is to generate random sequences under certain principles [10,11,12].
Fortunately, through massive data of mooring test-fly experiments and real flight experiments, we have found that the most explicit characteristics of the land clutters are unstable. Said another way, for any clutter, its point feature or body feature for the seeker detection is not certain. It may appear as a point under a certain incident angle, but the same clutter distraction may appear as a body target under another incident angle. Besides, for the same target, the target or body feature may change with different time or different weather conditions. To summarize, the clutter appears a fluctuating characteristic, or its property is quite sensitive to external conditions.
This unstable body or point feature of the clutters is helpful for target identification of the MMW seeker, since the point or body feature of the interested target is usually quite stable. Or better, the target feature is much more robust under various conditions including the backgrounds and the seeker status.
Making use of the different body/point features of the clutters and the target, a target identification method using correlation matching and beam pointing of the target is presented in this paper. Since body or point features of the clutters are not stable, we conduct two round scanning of the seeker to eliminate clutter distractions. Then, point or body feature identification is further evaluated by HRRP. In the end, error distribution zone is calculated on the basis of beam pointing. The beam pointing center will be viewed as the origin. Final target identification is realized by searching the target who lies in the closest error distribution zone of the beam pointing of the seeker. The performance of the target identification can be dramatically improved by using the proposed algorithm. The proposed method can realize satisfying target identification under complex land backgrounds, which provides further support for target tracking and precise attack for the MMW guidance missile.

2. Target Identification via Correlation Matching and Beam Pointing

The proposed method consists of two main steps. In the first step, we conduct scanning correlation for the targets detected in both scanning rounds. Additionally, in the second step, we compare the radial and azimuth range errors with the baseline parameters. The target that lies in the closest error zone of beam pointing of the target will be regarded as the final target identification result.

2.1. Correlation Matching

As is known, the seeker starts to work when it reaches the effective working range, it will scan and capture the target, and output the information needed by the guidance law, such as the pitching and yawing line-of-sight (LOS) rates. The diagram of the working process of the seeker is shown in Figure 1. In the following, we will discuss the scanning correlation matching strategy of the MMW seeker in detail, which will suppress distractions effectively. Firstly, the seeker starts scanning from the left to the right by moving the gimbal angles by the servo mechanism (we name it the first round scanning in this paper), and the seeker will implement point/body feature identification based on the radar echoes, and store the information of the detected target. Then, the seeker will conduct the second round scanning, i.e., scanning from the right to the left by moving the gimbal. The same as the first round, point or body identification will also be implemented in the second round, and related target information will also be stored in the meantime for the upcoming correlation matching.
Assuming that, targets ( T 1 1 , T 2 1 , , T i 1 , , T N 1 1 ) are detected at different times ( t 1 1 , t 2 1 , , t i 1 , , t N 1 1 ) in the first round scanning by using the ordered statistic constant false alarm rate (OS-CFAR) detection [13], where i = 1 , 2 , , N 1 , N 1 represents the total number of all the detected targets. The information of seeker to target distances ( R 1 1 , R 2 1 , , R i 1 , , R N 1 1 ) and the gimbal angles [ ( ϑ g 1 1 , φ g 1 1 ) , ( ϑ g 2 1 , φ g 2 1 ) , , ( ϑ g i 1 , φ g i 1 ) , , ( ϑ g N 1 1 , φ g N 1 1 ) ] in the missile body coordinate system are stored. Similarly, interested targets ( T 1 2 , T 2 2 , , T j 2 , , T N 2 2 ) are detected at time ( t 1 2 , t 2 2 , , t j 2 , , t N 2 2 ) in the second round scanning, where j = 1 , 2 , , N 2 , N 2 is the total number of detected targets in the second round scanning. The targets are also detected by using the OS-CFAR detection just as what is implemented in the first round scanning. The corresponding information ( R 1 2 , R 2 2 , , R j 2 , , R N 2 2 ) and [ ( ϑ g 1 2 , φ g 1 2 ) , ( ϑ g 2 2 , φ g 2 2 ) , , ( ϑ g j 2 , φ g j 2 ) , , ( ϑ g N 2 2 , φ g N 2 2 ) ] associated with the seeker and the target for the second round scanning are stored for correlation matching.
Having the target information from the seeker by the two round scanning process in hand, we will carry out correlation matching in the following. However, since the target information is obtained at different times, we have to realize time synchronization in the beginning. In this paper, we will compensate the parameters obtained from the first round scanning to the second round scanning. In other words, we will update the parameters of the detected targets ( T 1 1 , T 2 1 , , T i 1 , , T N 1 1 ) in the first round scanning to expected values at moments ( t 1 2 , t 2 2 , , t i 2 , , t N 1 2 ) . Then, we will compare these parameters to the ones obtained by the seeker itself in the second round scanning. If the correlation matching for the targets has been successful, we will view the successfully matched targets as the interested ones.
Here, we display the whole scanning matching process again to improve clarity. Firstly, target positions obtained in the first round scanning by the seeker in the inertial coordinate system are calculated by using coordinate system transformation. Then, target positions are updated by using the velocity information obtained from the inertial navigation system (INS), calculating target positions at the moments of target detection in the second round scanning in real time. In the following, update parameters of the seeker to target distances, the pitching gimbal angles and the yawing gimbal angles are calculated using trigonometric functions. In the end, correlation matching is implemented by making comparisons between the updated parameters and the ones outputted by the seeker itself in the second round scanning under certain restricted conditions. To further clarify the process, we will give a more detailed discussion on the assumption of only one detected target in the first round scanning.
The target position obtained in the first round scanning in the tracker coordinate system can be expressed as ( x m , y m , z m )
{ x m = R t 1 1 cos ε m cos β m y m = R t 1 1 sin ε m z m = R t 1 1 cos ε m sin β m ,
where R t 1 1 denotes the seeker to target distance at time t 1 , ε m denotes the misalignment angle in the pitching direction, and β m represents the misalignment angle in the yawing direction.
After coordinate transformation, target position ( x m 1 , y m 1 , z m 1 ) in the body coordinate system can be given by
[ x m 1 y m 1 z m 1 ] = A 1 [ x m y m z m ] ,
where A 1 is the direction cosine matrix.
A 1 = [ cos ϑ g cos φ g sin ϑ g cos ϑ g sin φ g sin ϑ g cos φ g cos ϑ g sin ϑ g sin φ g sin φ g 1 cos ϑ g ] ,
where φ g represents the yawing gimbal angular, and ϑ g represents the pitching gimbal angular.
As previously discussed, so as to achieve scanning correlation matching, we need to realize time synchronization firstly. Said another way, we need to update the seeker parameters associated with the same target (seeker to target distance and the gimbal angles in the body coordinate system) by the time detected in the second round scanning to realize further correlation matching.
Target position ( x m 2 , y m 2 , z m 2 ) in the navigation coordinate system can be expressed as [14]
[ x m 2 y m 2 z m 2 ] = A 2 [ x m 1 y m 1 z m 1 ] ,
where
A 2 = [ cos ϑ cos φ sin ϑ cos ϑ sin φ sin ϑ cos φ cos γ + sin φ sin γ cos ϑ cos γ sin ϑ sin φ cos γ + cos φ sin γ sin ϑ cos φ sin γ + sin φ cos γ cos ϑ sin γ sin ϑ sin φ sin γ + cos φ cos γ ] ,
where ( ϑ , φ , γ ) denotes the attitude angle set, which represents the pitching angle, the yawing angle, and the rolling angle, respectively. These parameters represent the attitude of the seeker, which can be obtained by the INS.
Target position ( x m 3 , y m 3 , z m 3 ) can be updated in real time during the flight by
{ x m 3 = x m 2 t 1 t v x ( t ) d t y m 3 = y m 2 t 1 t v y ( t ) d t z m 3 = z m 2 t 1 t v z ( t ) d t ,
where v x ( t ) , v y ( t ) , and v z ( t ) represents the north velocity, vertical velocity, and east velocity of the seeker, respectively. t is the time when the seeker detects the target in the second round scanning. The scanning speed is very fast, and the scanning angular scope is limited. Therefore, the time interval Δ t = t t 1 is usually quite small. The variance of the speed vector of the seeker can be neglected, and Equation (6) can be simplified by
{ x m 3 = x m 2 v x Δ t y m 3 = y m 2 v y Δ t z m 3 = z m 2 v z Δ t .
So far, we have obtained the updated target position at time t . Or better, ( x m 3 , y m 3 , z m 3 ) can be viewed as the compensated result for the first round scanning. What follows is to calculate the pitching and yawing gimbal angles. The target position ( x m 4 , y m 4 , z m 4 ) in the body coordinate system can be expressed as
[ x m 4 y m 4 z m 4 ] = A 2 T [ x m 3 y m 3 z m 3 ] .
Seeker parameters associated with the target can be calculated by using ( x m 4 , y m 4 , z m 4 ) . The seeker to target distance can be given by
R = x m 4 2 + y m 4 2 + z m 4 2 .
The beam pointing angle of the target along the pitching direction ϑ g can be given by
ϑ g = a r c sin ( y m 4 / x m 4 2 + y m 4 2 + z m 4 2 ) .
Additionally, the beam pointing angle of the target along the yawing direction φ g can be given by
φ g = a r c tan ( z m 4 / x m 4 ) .
Correlation matching is implemented by comparing the parameters of R , ϑ g and φ g at time t with the corresponding parameters outputted by the seeker at this time.
{ | R R | Δ R | ϑ g ϑ g | Δ ϑ g | φ g φ g | Δ φ g ,
where R is the seeker to target distance, ϑ g is the gimbal angle in the pitching direction, and φ g is the gimbal angle in the yawing direction outputted by the seeker itself at this time, Δ R , Δ ϑ g , and Δ φ g are the thresholds for correlation matching.
We have to note that these values are quite different under different circumstances. We have to determine them according to the real system in practice. For different INS equipment, the measurement precision of v x , v y , and v z are different, leading to different beam pointing directions. Moreover, the error of the gimbal angle output and the angle measurement error are also inevitable. Briefly, the determination of the parameters of Δ R , Δ ϑ g , and Δ φ g greatly relies on real application conditions, which have to be adjusted in practice.

2.2. Error Calculation of Radial Range and Azimuth Range

For complex land backgrounds, after correlation matching between two round scanning, multiple targets may still meet the restriction of correlation matching. We cannot realize target identification just by using the proposed method presented in Section 2.1. Assuming that, there are K targets { T c 1 , T c 2 , , T c K } which satisfy the condition restricted by Equation (12), the corresponding parameters of the seeker are [ ( ϑ g c 1 , φ g c 1 ) , ( ϑ g c 2 , φ g c 2 ) , , ( ϑ g c K , φ g c K ) ] and ( R c 1 , R c 2 , , R c K ) , respectively. In the following, we calculate the radial range error and the azimuth range error. Target position ( x t , y t , z t ) calculated in the navigation coordinate system can be expressed as
x t = R · cos ( θ n 0 ) · cos ( φ n 0 ) y t = R · sin ( θ n 0 ) z t = R · cos ( θ n 0 ) · sin ( φ n 0 ) } ,
where R represents the seeker to target position calculated by the fire-control system, θ n 0 is the pitching angle, and φ n 0 is the azimuth angle.
Additionally, target position ( x t 0 , y t 0 , z t 0 ) in the launching system can be given by
[ x t 0 y t 0 z t 0 ] = A ( φ b x ) [ x t y t z t ] ,
where φ b x represents the north angle, the direction cosine matrix A can be given by
A ( φ b x ) = [ cos φ b x 0 sin φ b x 0 1 0 sin φ b x 0 cos φ b x ] .
Since the position ( x m , y m , z m ) can be obtained by the INS mounted on the seeker, the seeker to target distance R can be calculated in real time
R = ( x t 0 x m ) 2 + ( y t 0 y m ) 2 + ( z t 0 z m ) 2 .
For the k th ( k = 1 , 2 , , K ) target, the radial range error Δ R c k r can be given by
Δ R c k r = R c k R .
According to the geometry relationship between the seeker and the target, the beam pointing direction of the yaw angle φ p t can be expressed as
φ p t = a r c tan ( z t 0 z m x t 0 x m ) .
Additionally, the azimuth range error can be given by
Δ R c k a = R c k · ( φ g c K φ p t ) .
Hereto, we obtained the range errors of the successfully matched K targets { ( Δ R c k r , Δ R c k a ) , ( Δ R c k r , Δ R c k a ) , , ( Δ R c k r , Δ R c k a ) , , ( Δ R c K r , Δ R c K a ) } . What follows is to calculate the location of the targets at the error distribution zone. The area was divided into 1 σ , 2 σ , and 3 σ error distribution zones empirically according to massive mooring test-fly data and real flight data exploration with the consideration of various error introduction terms. Final target identification is conducted by locating the target, which lies in the closest error distribution zone with respect to the beam pointing of the seeker. The flow diagram of the proposed method is displayed in Figure 2.

3. Experimental Results and Analysis

Mooring test-fly data and real flight data were used to verify the effectiveness of the proposed method.

3.1. Mooring Test-Fly Experiment

Firstly, we tested the proposed method on a mooring test-fly experiment. Here, we will give an explanation of the mooring test-fly experiment. As is known, a seeker is very expensive, especially for the MMW seeker. As a result, before a seeker is mounted on the missile for launching, we have to evaluate its performance carefully. To fulfill the task, we usually conduct the mooring test-fly experiment. In the mooring test-fly experiment, we mounted the seeker on a platform, such as a helicopter with pilots or an unmanned aerial vehicle (UAV) with remote control on the ground. The platform with the seeker mounted on will move according to the missile trajectory. When the distance from the seeker to the target is within the effective range of the seeker, the seeker will start to illuminate an electromagnetic wave, and implement all the processing just as the real missile fight case. The track will continue until the platform approaches the target, i.e., the distance from the seeker to the target is zero. Then, the platform will fly high and far from the target, getting ready for the next round trip for seeker data collection and performance evaluation. Since the seeker still exists, we can conduct similar processes again and again to search satisfying parameters. For the MMW seeker that we focus on, the platform with the MMW seeker mounted on flies a simulated trajectory. The thresholds for correlation matching are set to be 15 m for range, and 2° for angles, respectively.
The seeker starts to work when reaching the effective working range. Two round scanning of the seeker was conducted subsequently. A total of 11 targets were detected in the first round scanning. The corresponding target information is tabulated in Table 1, target 1- c ( c = 1 , 2 , 11 ) denotes the c th detected target. We can know that target 1-3 is the accurate target by the prior information calculation from target location and seeker position obtained by the high-precision global position system (GPS) mounted on the target and the INS mounted on the seeker [15,16].
As previously discussed, point/body feature of the target can be determined by using the HRRP [7,8,9]. The HRRP of the MMW echo is given in Figure 3. Additionally, the profile after target detection is shown in Figure 4. As can be seen, there are two objects that exceed the threshold. One appears as a body target, which coincides with the size of a tank, and the other one appears as a point target.
As can be seen in Table 1, targets that appear as point features will be discarded (targets 1-5 and 1-7 in Table 1), and the five strongest reflection targets with body features will be chosen for correlation matching, i.e., targets 1-1, 1-2, 1-3, 1-4, and 1-6.
Then, the seeker implemented the second round scanning, and 10 targets were detected in this round scanning. The corresponding target information is shown in Table 2. Similar to Table 1, 2- n ( n = 1 , 2 , 10 ) denotes the n th target in the second round scanning.
Likewise, discarding the targets appearing as point features (targets 2-7, 2-8, and 2-10 in Table 2), and choosing the five strongest reflection targets with body features, we can get targets 2-1, 2-2, 2-3, 2-4 and 2-5 for correlation matching.
Correlation matching results for targets 1-1, 1-2, 1-3, 1-4, and 1-6 are given in Table 3, Table 4, Table 5, Table 6 and Table 7, respectively.
From Table 3, we can see that only target 2-1 can realize satisfying matching with target 1-1. The range difference is only 6 m, and the value is the smallest of all the candidate targets. We can tell that the range difference may be viewed as the top choice for correlation matching. In addition, the second effective factor is the pitching angle difference. The azimuth angle can provide further assistance.
From Table 4, we can see that only target 2-2 realizes satisfying matching with target 1-2, and we can further find that the range difference of target 2-4 is also quite small, which is 24 m. Although the range difference is the top choice, we have to set the threshold appropriately to avoid misjudging. Combining the pitching and azimuth angles provides great discriminating power.
From Table 5, we can see that target 2-3 realizes satisfying matching with target 1-3. In addition, from Table 6, we can see that there is no target that matches target 1-4. All the targets selected in the second round scanning cannot match target 1-4. In other words, target 1-4 will not arrive at the final stage of target identification. From Table 7, we can see that target 1-6 realizes satisfying correlation with target 2-4.
To summarize, after two round target correlation matching, the following target pairs realize satisfying matching: targets 1-1 and 2-1, targets 1-2 and 2-2, targets 1-3 and 2-3, targets 1-6 and 2-4. Here, we label the successfully matched targets as target 1, target 2, target 3, and target 4, respectively. The radial range and azimuth range errors with respect to the beam pointing center are given in Figure 5. In Figure 5, the blue line and the red line represent the distance differences in radial and azimuth directions with respect to the origin, i.e., the beam pointing center, respectively.
Then, we calculated error distances associated with the error distribution zone. Just as the determination of the thresholds set for correlation matching, the calculation of the error distribution zone is even more complex. The determination of the parameters along the range and azimuth directions involves plenty of error terms, such as target orientation error by the fire control system, floating error of the INS, and angle determination error of the seeker, etc. Here, we set the parameters empirically through massive data exploration again. The 1 σ error distribution rectangular zone is formed by 35 and 50 m in the radial and azimuth directions, respectively. Said another way, for a given seeker, the determination of these parameters is achieved by excessive times of mooring test-fly experiments to support the real missile flight.
According to the geometry between the interested targets and the beam pointing center, the error distribution zone with the four successfully matched targets are shown in Figure 6. The abscissa denotes the azimuth range error, and the ordinate denotes the radial range error. From Figure 6, we can see that target 1 and target 2 lie in the 2 σ error zone of the beam pointing center. Target 4 lies in the 3 σ error zone of the beam pointing center. Only target 3 lies in the 1 σ error zone of the beam pointing center. From Figure 6, we can see that target 3 lies in the closest error distribution zone of the beam pointing of the seeker. Target 3 will be viewed as the interested target for the seeker. We can see that target identification can be obtained, and the proposed method is capable of providing powerful assistance and support for precise attack of the missile.
Here, we will give more explanations of the parameters, which will affect the final identification result. Although satisfying results can be obtained by using the proposed method, parameter setting is of great importance, and appropriate parameter setting is the precondition of accurate identification. When the parameters are not well-selected, wrong target identification will emerge. In this part, we will display examples in situations with inappropriate parameter setting. Firstly, we take the azimuth angle difference threshold as an example. From Table 3, we can see that if we set the azimuth angle difference threshold to be less than 1.48°, the correlation matching for the interested target 3 will fail. In other words, target 3 will not be selected for further processing. Wrong identification will inevitably occur.
In the following, we take the error distribution zone threshold as another example of inappropriate parameter setting for illustration. In the experiment, the 1 σ error distribution zone is given by the inside rectangular with radial range threshold of 35 m and azimuth range threshold of 50 m, as shown in Figure 6. If we set the radial range threshold to be 50 m, the error distribution zone for the targets will change, as shown in Figure 7. Comparing Figure 6 with Figure 7, we can see that for the four interested targets, although target 1 and target 3 still lie in the same error distribution zone with different radial range thresholds, target 2 and target 4 change their error distribution zones. Target 2 lies in the 2 σ error distribution zone in Figure 6, whereas it lies in the 1 σ error distribution zone in Figure 7. As for target 4, it lies in the 3 σ error distribution zone in Figure 6, but it approaches to the 2 σ error distribution zone in Figure 7.
As can be seen, both target 2 and target 3 lie in the 1 σ error distribution zone in Figure 7. We will choose the target with larger amplitude to be the final identification result, when more than one target lies in the same error distribution zone. That is to say, wrong target identification will emerge in this situation. Target 2 corresponds to target 1-2 in the first round and target 2-2 in the second round, whereas target 3 corresponds to target 1-3 in the first round and target 2-3 in the second round. From Table 1, we can see that the amplitude of target 1-2 is 120, and the amplitude of target 1-3 is 100. From Table 2, we can see that the amplitude of target 2-2 is 120, whereas the amplitude of target 2-3 is 111. As can be seen, the amplitude of target 2 is larger than that of target 3 in both rounds. In other words, target 2 will be viewed to be more threatening than target 3. This is definitely what we would like to avoid.
Above all, appropriate parameter setting is of great importance to target identification by using the proposed method. Inappropriate parameter setting will lead to undesirable failure. That is to say, one has to adjust the parameters according to the real seeker system in practice.

3.2. Real Data

In the following, we verify the effectiveness of the proposed method by using real flight data of the seeker. The background of the targets is not as complex as the one in mooring test-fly experiment. A total of 14 targets were detected in the first round scanning, and the target information is shown in Table 8. In the experiment, target 1-1 is the interested target. Fifteen targets were detected in the second round scanning, and the corresponding target information is shown in Table 9. Firstly, we will select the five largest targets with body features as correlation matching candidates just as the case in Section 3.1. That is to say, we would like to implement correlation matching between target sets of 1-1, 1-2, 1-3, 1-4 1-5 detected in the first round scanning and target sets of 2-1, 2-2, 2-3, 2-4, and 2-5 in the second round scanning. The correlation matching results for targets 1-1, 1-2, 1-3, 1-4, and 1-5 are shown in Table 10, Table 11, Table 12, Table 13 and Table 14, respectively.
From the experimental results we can tell that all the targets fail to realize successful correlation matching except target 1-1. We can see that from Table 10, target 1-1 realizes satisfying matching with target 2-2. The seeker will view target 1-1 as the interested target to attack. The radial and azimuth range errors are 50 and 116.7 m, respectively. The thresholds for error distribution zone are set to be the same as the mooring test-fly experiment (actually, as previously discussed, the parameters are determined by excessive mooring test-fly experiments). In other words, the target lies in the 3 σ error zone of the beam pointing center. Real flight data further verifies the effectiveness of the proposed method.

4. Conclusions

So as to overcome the great obstacle of target identification in land backgrounds for the MMW seeker, especially under complex conditions, an effective method taking the unstable body or point feature characteristic of the clutters into account is presented in this paper. Clutter distractions are discarded effectively through information correlation matching by using two round scanning of the seeker. A satisfying target identification result is obtained by combining the prior beam pointing information.
Moreover, appropriate parameter setting of the thresholds as well as the appropriate construction of the error distribution zone is the precondition of the proposed method. However, the determination of the parameters is restricted by various factors. The parameters of the proposed method presented in this paper is determined by the exploration and summation of plenty of real data including mooring test-fly and real flight experimental data. How to construct a precise model for the determination of these parameters deserves further study.

Author Contributions

S.C. proposed the idea of the algorithm; M.L. and S.C. conceived and performed the experiments; S.C. wrote the paper; F.L. and M.X. supervised the work and provided valuable suggestions.

Funding

This work was funded by the national natural science foundation of China (No. 61701289), the natural science foundation of Shaanxi province (No. 2018JQ6083, No. 2018JQ6087), and the fundamental research funds for the central universities (No. GK201903084).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Lu, F.; Chen, S.; Liu, M.; Wang, J.; Ma, F.; Yang, T. A Target Recapturing Method for the Millimeter Wave Seeker with Narrow Beamwidth. In Proceedings of the IGARSS 2018-2018 IEEE International Geoscience and Remote Sensing Symposium, Valencia, Spain, 22–27 July 2018; pp. 2841–2844. [Google Scholar]
  2. Cheng, Y.; Hu, F.; Gui, L.; Su, J.; Qi, B.; Liu, S.; Huang, M. Linear polarisation property and fusion method for target recognition in passive millimetre-wave polarimetric imaging. Electron. Lett. 2016, 52, 1221–1223. [Google Scholar] [CrossRef]
  3. Cecchini, L.; Pizzingrilli, E.; Russo, S.; D’Elia, U. MMW Active Phased Array Seeker Project for Hit to Kill Engagement. In Proceedings of the 2008 IEEE Radar Conference, Rome, Italy, 26–30 May 2008; pp. 1–6. [Google Scholar]
  4. Bi, H.; Zhang, B.; Zhu, X.; Hong, W.; Sun, J.; Wu, Y. L1-regularization-based SAR imaging and CFAR detection via complex approximated message passing. IEEE Trans. Geosci. Remote Sens. 2017, 55, 3426–3440. [Google Scholar] [CrossRef]
  5. Chang, W.; Tao, H.; Sun, G.; Wang, Y.; Bao, Z. A novel multi-angle SAR imaging system and method based on an ultrahigh speed platform. Sensors 2019, 19, 1701. [Google Scholar] [CrossRef] [PubMed]
  6. Liu, A.; Wang, F.; Xu, H.; Li, L. N-SAR: A new multichannel multimode polarimetric airborne SAR. IEEE J. Sel. Topics Appl. Earth Observ. Remote Sens. 2018, 11, 3155–3166. [Google Scholar] [CrossRef]
  7. Xing, Y.; You, P.; Yong, S. Parameter estimation of micro-motion targets for high-resolution-range radar using online measured reference. Sensors 2018, 18, 2773. [Google Scholar] [CrossRef] [PubMed]
  8. Zhu, L.; Xiong, G.; Yu, W. Radar HRRP Group-Target Recognition Based on Combined Methods in the Background of Sea Clutter. In Proceedings of the 2018 International Conference on Radar (RADAR), Brisbane, QLD, Australia, 27–31 August 2018; pp. 1–6. [Google Scholar]
  9. Zhang, H.; Ding, D.; Fan, Z.; Chen, R. Adaptive neighborhood-preserving discriminant projection method for HRRP-based radar target recognition. IEEE Antennas Wirel. Propag. Lett. 2014, 14, 650–653. [Google Scholar] [CrossRef]
  10. Zhu, G.; Chen, Y.; Yin, H. Analysis of Typical Ground Clutter Statistical Characteristics. In Proceedings of the 2017 International Applied Computational Electromagnetics Society Symposium (ACES), Suzhou, China, 1–4 August 2017; pp. 1–2. [Google Scholar]
  11. Jiang, H.; Liao, G.; Qu, Y. Compensation of Clutter Spectrum for Airborne Forward-Looking Radar Based on the Clutter Subspace Transformation. In Proceedings of the 2009 International Radar Conference “Surveillance for a Safer World” (RADAR 2009), Bordeaux, France, 12–16 October 2009; pp. 1–4. [Google Scholar]
  12. Liu, Z.; Chen, Y.; Xie, R.; Zhang, L. Clutter Property Analysis and Clutter Suppression for Hypersonic Airborne Radar. In Proceedings of the 2013 International Conference on Radar, Adelaide, SA, Australia, 9–12 September 2013; pp. 374–379. [Google Scholar]
  13. Hyun, E.; Lee, J. A New OS-CFAR Detector Design. In Proceedings of the 2011 First ACIS/JNU International Conference on Computers, Networks, Systems and Industrial Engineering, Jeju Island, South Korea, 23–25 May 2011; pp. 133–136. [Google Scholar]
  14. Yu, Y.; Zhang, Y.; Liu, Z. Coordinate systems for use on guide bomb’s dynamics of flight. Acta Armamentarii 2002, 23, 201–204. [Google Scholar]
  15. Lu, F.; Chen, S.; Wang, J.; Liu, M.; Yang, T. A MMW Seeker Performance Evaluation Method for Moving Targets via RTK Technology. In Proceedings of the IGARSS 2018-2018 IEEE International Geoscience and Remote Sensing Symposium, Valencia, Spain, 22–27 July 2018; pp. 1078–1081. [Google Scholar]
  16. Malleswaran, M.; Deborah, S.; Manjula, S.; Vaidehi, V. Integration of INS and GPS Using Radial Basis Function Neural Networks for Vehicular Navigation. In Proceedings of the 2010 11th International Conference on Control Automation Robotics & Vision, Singapore, 7–10 December 2010; pp. 2427–2430. [Google Scholar]
Figure 1. Diagram of the seeker working process.
Figure 1. Diagram of the seeker working process.
Sensors 19 02530 g001
Figure 2. The flow diagram of the proposed algorithm.
Figure 2. The flow diagram of the proposed algorithm.
Sensors 19 02530 g002
Figure 3. High resolution range profile (HRRP) of a detected target.
Figure 3. High resolution range profile (HRRP) of a detected target.
Sensors 19 02530 g003
Figure 4. HRRP after target detection.
Figure 4. HRRP after target detection.
Sensors 19 02530 g004
Figure 5. Radial and azimuth range differences.
Figure 5. Radial and azimuth range differences.
Sensors 19 02530 g005
Figure 6. Error distribution zone of the interested targets.
Figure 6. Error distribution zone of the interested targets.
Sensors 19 02530 g006
Figure 7. Error distribution zone of the interested targets with radial range threshold of 50 m.
Figure 7. Error distribution zone of the interested targets with radial range threshold of 50 m.
Sensors 19 02530 g007
Table 1. Detected target information in the first round scanning.
Table 1. Detected target information in the first round scanning.
Target1-11-21-31-41-51-61-71-81-91-101-11
Point/Bodybodybodybodybodypointbodypointbodybodybodybody
Amplitude121120111106105102101100949292
Pitching angle (°)1.171.111.141.251.231.181.181.151.141.291.31
Azimuth angle (°)−3.24−3.71−4.71−0.570.77−7.25−7.26−0.93−6.661.491.01
Seeker to target distance (m)16951806173317541695183718941652164218961839
Table 2. Detected target information in the second round scanning.
Table 2. Detected target information in the second round scanning.
Target2-12-22-32-42-52-62-72-82-92-10
Point/Bodybodybodybodybodybodypointpointpointbodypoint
Amplitude12112011210310210299999378
Pitching angle (°)1.251.141.21.191.181.431.431.281.171.19
Azimuth angle (°)−2.14−2.97−3.23−5.9−0.76−6.68−6.32−5.711.92.85
Seeker to target distance (m)1689179817351830164418871689163518861797
Table 3. Correlation matching results for target 1-1.
Table 3. Correlation matching results for target 1-1.
Target2-12-22-32-42-52-9
Pitching angle difference (°)0.080.030.030.020.010
Azimuth angle difference (°)1.10.270.012.662.485.14
Range difference (m)61034013551191
Table 4. Correlation matching results for target 1-2.
Table 4. Correlation matching results for target 1-2.
Target2-12-22-32-42-52-9
Pitching angle difference (°)0.140.030.090.080.070.06
Azimuth angle difference (°)1.570.740.482.192.955.61
Range difference (m)1178712416280
Table 5. Correlation matching results for target 1-3.
Table 5. Correlation matching results for target 1-3.
Target2-12-22-32-42-52-9
Pitching angle difference (°)0.1100.060.050.040.03
Azimuth angle difference (°)2.571.741.481.193.956.61
Range difference (m)446529789153
Table 6. Correlation matching results for target 1-4.
Table 6. Correlation matching results for target 1-4.
Target2-12-22-32-42-52-9
Pitching angle difference (°)00.110.050.060.070.08
Azimuth angle difference (°)1.572.42.665.330.192.47
Range difference (m)65441976110132
Table 7. Correlation matching results for target 1-6.
Table 7. Correlation matching results for target 1-6.
Target2-12-22-32-42-52-9
Pitching angle difference (°)0.070.040.020.010.000.01
Azimuth angle difference (°)5.114.284.021.356.499.15
Range difference (m)14839102719349
Table 8. Detected target information in the first round scanning for real flight data.
Table 8. Detected target information in the first round scanning for real flight data.
Target1-11-21-31-41-51-61-71-81-91-101-111-121-131-14
Point/Bodybodybodybodybodybodybodybodybodypointpointbodybodypointpoint
Amplitude11096959595949494939292919090
Pitching angle (°)−6.98−7.01−7.05−7.04−7.04−7.1−7.09−7.01−7.05−7.04−6.96−6.97−7.1−7.1
Azimuth angle (°)−1.011.172.623.913.914.12−5.121.172.62−4.47−0.77−2.31−5.21−5.21
Seeker to target distance (m)22982163220623372304240622052262243024372363225722352364
Table 9. Detected target information in the second round scanning for real flight data.
Table 9. Detected target information in the second round scanning for real flight data.
Target2-12-22-32-42-52-62-72-82-92-102-112-122-132-142-15
Point/Bodybodybodybodybodybodybodybodybodybodybodypointbodybodypointbody
Amplitude1079897969695959595959494949494
Pitching angle (°)−6.99−7−7−7.1−7.09−7−7.09−7.04−7.01−7.01−6.88−7.03−7−6.97−6.99
Azimuth angle (°)1.32−2.73−0.264.555.222.335.223.50.391.11−3.81−2.170.9−3.61−3.38
Seeker to target distance (m)229422952185244122612150230323412230239523432153226124112380
Table 10. Correlation matching results for target 1-1 for real flight data.
Table 10. Correlation matching results for target 1-1 for real flight data.
Target2-12-22-32-42-52-9
Pitching angle difference (°)0.010.020.020.020.120.11
Azimuth angle difference (°)2.331.720.752.665.566.23
Range difference (m)43113143374
Table 11. Correlation matching results for target 1-2 for real flight data.
Table 11. Correlation matching results for target 1-2 for real flight data.
Target2-12-22-32-42-5
Pitching angle difference (°)0.020.010.010.090.08
Azimuth angle difference (°)0.153.91.433.384.05
Range difference (m)1311322227898
Table 12. Correlation matching results for target 1-3 for real flight data.
Table 12. Correlation matching results for target 1-3 for real flight data.
Target2-12-22-32-42-5
Pitching angle difference (°)0.060.050.050.050.04
Azimuth angle difference (°)1.35.352.881.932.6
Range difference (m)888911023555
Table 13. Correlation matching results for target 1-4 for real flight data.
Table 13. Correlation matching results for target 1-4 for real flight data.
Target2-12-22-32-42-5
Pitching angle difference (°)0.050.040.040.060.05
Azimuth angle difference (°)2.596.644.170.641.31
Range difference (m)434215210476
Table 14. Correlation matching results for target 1-5 for real flight data.
Table 14. Correlation matching results for target 1-5 for real flight data.
Target2-12-22-32-42-5
Pitching angle difference (°)0.050.040.040.060.05
Azimuth angle difference (°)2.596.644.170.641.31
Range difference (m)10911913743

Share and Cite

MDPI and ACS Style

Chen, S.; Liu, M.; Lu, F.; Xing, M. A Target Identification Method for the Millimeter Wave Seeker via Correlation Matching and Beam Pointing. Sensors 2019, 19, 2530. https://doi.org/10.3390/s19112530

AMA Style

Chen S, Liu M, Lu F, Xing M. A Target Identification Method for the Millimeter Wave Seeker via Correlation Matching and Beam Pointing. Sensors. 2019; 19(11):2530. https://doi.org/10.3390/s19112530

Chicago/Turabian Style

Chen, Shichao, Ming Liu, Fugang Lu, and Mengdao Xing. 2019. "A Target Identification Method for the Millimeter Wave Seeker via Correlation Matching and Beam Pointing" Sensors 19, no. 11: 2530. https://doi.org/10.3390/s19112530

APA Style

Chen, S., Liu, M., Lu, F., & Xing, M. (2019). A Target Identification Method for the Millimeter Wave Seeker via Correlation Matching and Beam Pointing. Sensors, 19(11), 2530. https://doi.org/10.3390/s19112530

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop