Next Article in Journal
A Novel Pre-Processing Technique for Original Feature Matrix of Electronic Nose Based on Supervised Locality Preserving Projections
Previous Article in Journal
Joint Parameter Estimation for the Two-Wave with Diffuse Power Fading Model
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Parameter Estimation of a Ground Moving Target Using Image Sharpness Optimization

School of Marine Science and Technology, Northwestern Polytechnical University, Xi’an 710072, China
*
Author to whom correspondence should be addressed.
Sensors 2016, 16(7), 1017; https://doi.org/10.3390/s16071017
Submission received: 13 April 2016 / Revised: 16 June 2016 / Accepted: 25 June 2016 / Published: 30 June 2016
(This article belongs to the Section Remote Sensors)

Abstract

:
Motion parameter estimation of a ground moving target is an important issue in synthetic aperture radar ground moving target indication (SAR-GMTI) which has significant applications for civilian and military. The SAR image of a moving target may be displaced and defocused due to the radial and along-track velocity components, respectively. The sharpness cost function presents a measure of the degree of focus of the image. In this work, a new ground moving target parameter estimation algorithm based on the sharpness optimization criterion is proposed. The relationships between the quadratic phase errors and the target’s velocity components are derived. Using two-dimensional searching of the sharpness cost function, we can obtain the velocity components of the target and the focused target image simultaneously. The proposed moving target parameter estimation method and image sharpness metrics are analyzed in detail. Finally, numerical results illustrate the effective and superior velocity estimation performance of the proposed method when compared to existing algorithms.

Graphical Abstract

1. Introduction

Ground moving target indication (GMTI) combined with synthetic aperture radar (SAR) has become a well-established technique for civil and military applications. The challenge of GMTI includes both the detection of targets and the estimation of their velocity components. Recently, the topic of moving target detection has been extensively studied [1,2,3,4,5]. In conventional SAR images, the target will be displaced and defocused because of the target’s radial (cross-track) and along-track velocities, respectively [6]. The target’s moving parameters, such as radial velocity, along-track velocity, and original azimuth position, play an important role in monitoring ground vehicles. Using these parameters, we can reposition them to the true azimuth location and extrapolate a target’s future position [7,8]. Generally, the interferometric phase is employed to determine the radial velocity while the along-track velocity is acquired by the target’s Doppler parameters. For a space-borne SAR system, the speed of ground slow moving targets are rather slow compared with the high speed of the satellite platform (typical platform velocities are in the range of 7000 m/s to 7500 m/s), therefore, the influence of the target’s along-track velocity can be ignored. In contrast, for airborne SAR systems, the target’s along-track velocity component will inevitably cause the image to be defocused, consequently reducing the detection capability.
In this paper, the main effort lays in the estimation of the motion parameters of a point moving target whose along-track velocity must be considered. Assuming that the moving targets have been detected using the displaced phase center antenna (DPCA) method, the along-track interferometric (ATI) method, or some other method, a new moving target parameter estimation algorithm based on the image sharpness optimization criterion is proposed. The target’s velocity components can be obtained by optimizing the sharpness cost function. The algorithm shows effective and superior velocity estimation performance in a number of simulations.
The paper is arranged as follows: Section 2 provides the echoes model of a typical SAR-GMTI system. In Section 3, we provide a brief introduction to the image sharpness metrics. In Section 4, the new moving target parameter estimation algorithm is described in detail. In Section 5, the algorithm is applied to the simulated SAR data; the performance of the algorithm is proved to be effective. In the closing section, a summary of the findings is presented.

2. Echoes Model of a Ground Moving Target

Recently, the majority of SAR-GMTI algorithms utilize multiple apertures to provide an additional degree of freedom with which unwanted clutters may be suppressed. For simplicity and without loss of generality, we consider a single aperture SAR system model for the echoes backscattered from a ground point moving target [9,10]. Assume that the SAR system functions in side-looking mode with a fixed pointing angle orthogonal to the flight path. The moving target and radar geometry relation is illustrated in Figure 1, where the projection of the aircraft flying direction is defined as the X-axis, and the Y-axis is the radial (cross-track) direction, and the Z-axis represents the altitude. The platform moves at a constant speed va and a fixed altitude H. An arbitrary point ground moving target P is assumed to be at position (X0,Y0,0) at t = 0 and moving with along-track and cross-track velocity components vx and vy, respectively. Suppose that when t = 0, the slant range of radar to target P is R 0 = H 2 + X 0 2 + Y 0 2 . For a small antenna azimuth beamwidth, it is possible to approximate the instantaneous range between radar and the point target by a second-order Taylor series:
R ( t ) R 0 + v y Y 0 R 0 t ( v a v x ) X 0 R 0 t + v y 2 + ( v a v x ) 2 2 R 0 t 2
Assuming the radar transmits a linear frequency modulated signal, after range compression, the signal can be expressed as:
S 1 ( t ) = A T r e c t ( t T ) exp ( j 4 π λ R ( t ) )
where T is the signal pulse period and λ is the wavelength.
After Fourier transformation, the signal can be expressed as follows:
S 1 ( f ) = A T 1 | f d r | r e c t ( f + f d c B ) exp ( j 2 π f d c f f d r ) exp ( j π f 2 f d r )
where AT is the magnitude, B is the transmission bandwidth, fdc is the Doppler center frequency of the target, and fdr is the Doppler rate of the target. Derived from Equation (1), fdc and fdr can be given as:
f d c = 2 λ × d R ( t ) d t = 2 λ R 0 [ X 0 ( v a v x ) v y Y 0 ]
f d r = 2 λ × d 2 R ( t ) d t 2 = 2 λ R 0 [ ( v a v x ) 2 + v y 2 ]
For a stationary scene, the Doppler center frequency and Doppler rate are:
f d c 0 = f d c | v x = 0 , v y = 0 = 2 X 0 v a λ R 0
f d r 0 = f d r | v x = 0 , v y = 0 = 2 v a 2 λ R 0
In SAR imaging, the reference function used in azimuth processing is constructed using the stationary target’s Doppler parameters. The reference function can be formed as follows:
s r e f a ( f ) = exp ( j π f 2 / f d r 0 )
Synthetic aperture radar can image a stationary scene with fine resolution. However, for a moving target, smearing, defocusing, and displacement are inevitable because the Doppler parameters of a moving target are different from those of a stationary scene.
Consequently, the signal after azimuth processing in the Doppler domain can be expressed as:
S 1 ( f ) = A T 1 | f d r | r e c t ( f + f d c B ) exp ( j 2 π f d c f f d r ) exp ( j π f 2 f d r + j π f 2 f d r 0 )
From Equation (9), for a moving target, the quadratic phase term is not equal to zero because fdrfdr0. In other words, the Doppler rate of the moving target is not completely compensated and there still exists quadratic phase errors. For SAR images, the linear phase errors may cause image displacement, whereas the quadratic phase errors may introduce image defocus and reduce the signal amplitude. If the quadratic phase errors can be estimated and fully compensated for the signal, then the image will be focused and the signal energy will be converged in the center, see Equation (11).
Let γ represent the quadratic phase error in Equation (9):
γ = j π f d r j π f d r 0 = j π { λ R 0 2 X 0 v a λ R 0 2 [ ( v a v x ) 2 + v y 2 ] }
After phase error compensation, the focused image in the Doppler domain can be expressed as:
S c o m p e n ( f ) = S 1 ( f ) exp ( γ f 2 )
After compensation, the quadratic phase is equal to zero, see Equation (12):
j π f 2 f d r + j π f 2 f d r 0 + γ = 0
One of the measurements of the degree of image focus is named image sharpness. In the next section, we will provide a brief introduction to image sharpness metrics.

3. Image Sharpness Metrics

Sharpness is considered as the ratio of the difference in luminance of an object and its immediate surroundings. In SAR images, sharpness can be considered as a measurement of the degree of the focus of the image because the sharpness allows one to emphasize the difference in the intensity of the scene. The amplitude of the focused image is concentrated in the center pixels. As a result, the focused image has several pronounced peaks corresponding to the scatters of the target that produce high fluctuations of the image intensity. Conversely, the amplitude of the unfocused image exhibits small fluctuations around its mean value. Thus, we expect a high sharpness value of the focused image because there are great differences in the intensity. Instead, the amplitude of an unfocused image is concentrated around its mean value and the sharpness value is low.
Because of unknown platform movement and ground target motion, SAR suffers from image degradation due to the presence of phase errors in the received signal. These phase errors act as a blurring filter, resulting in loss of resolution, spurious targets and a decrease in image sharpness. If the quadratic phase errors can be estimated and compensated for the signal, then the image can be focused and the sharpness value will reach the highest value. Determining how to estimate the phase errors is the key problem for image focusing.
Figure 2 shows a small part of real SAR image; it contains 5 × 5 pixels. The amplitude of each pixel is presented by different colors.
Let Ii,j(γ) represent the complex value of pixel (i,j) in the phase compensated image S c o m p e n , which is a function of γ.
Image sharpness is used to evaluate the degree of focus [11,12,13,14,15]. Because of the point-like nature of the SAR image model, maximizing sharpness is found to increase the image focus. Because the complex value of the compensated image pixel is a function of γ, the sharpness can also be expressed as a function of γ.
In this paper, we study the standard deviation intensity sharpness function, as proposed in [9]:
C ( γ ) = A { [ i = 1 M j = 1 N I i j 2 ( γ ) A { i = 1 M j = 1 N I i j 2 ( γ ) } ] 2 } A { i = 1 M j = 1 N I i j 2 ( γ ) }
where M and N are the number of pixels in range and azimuth directions, respectively. A{•} is the spatial mean operator. As an example, application of A{•} to a real sequence x(p) with p = 1,…,p is shown in Equation (14):
A [ x ( p ) ] = 1 p p = 1 P x ( p )
Sharpness optimization is a method by which a quadratic phase error estimate is chosen to maximize a particular cost function, see Equation (15):
γ = arg max γ C ( γ )
The properties of the sharpness function C(γ) have been analyzed in detail in [15]. When C(γ) achieves the maximum value, the quadratic phase error is completely compensated and the image will be focused. The variable γ corresponds to the estimated quadratic phase errors, which is presented by γ .
For example, Figure 3 illustrates the relationship between the sharpness value and the quadratic phase. For a point moving target, when the quadratic phase error is fully compensated, the quadratic phase coefficient is equal to zero, and the sharpness value is maximized.
In SAR imaging, sharpness metrics are always used for image autofocus. Through the process of estimating the phase errors and compensating them for the unfocused image, a higher definition image can be acquired. For a moving target, the velocities introduce a phase error, thus causing the image to be unfocused. The phase error of a moving target can also be estimated using sharpness optimization metrics. Thus, inspired by the concept of sharpness, we proposed a moving target parameter estimate method in the next section.

4. Parameter Estimation Using Sharpness Optimization

4.1. Coarse Parameter Estimation

Traditionally, after clutter suppression and moving target detection in the image domain, the motion parameter of the target can be estimated using ATI [16,17] and DPCA [18,19] techniques.
For ATI techniques, the velocities of the ground moving target can be obtained via the interferometric phase, using Equations (16) and (17):
X 0 = λ R 0 ϕ ( f ) 2 π d + d 4
v y = v a X 0 Y 0 λ R 0 f 2 Y 0
In which ϕ ( f ) is the interferometric phase, d is the distance between phase centers, and the other parameters are the same as mentioned before.
Because the image of moving target is unfocused, the interferometric phase ϕ ( f ) is also inaccurate. The motion parameters obtained by Equations (16) and (17) may have large errors. To improve the estimation accuracy, we propose an accurate target motion parameter estimation method that estimates the parameters and obtains the focused image simultaneously.

4.2. Accurate Parameter Estimation

Derived from Equation (10), the quadratic phase error γ is a function of the along-track velocity vx and the radial velocity vy. From Figure 3, if the quadratic phase error is completely compensated, then the sharpness function reaches peak amplitude and the image is fully focused. Thus, the estimation of the values of vx and vy involve an optimization problem of the function of Equation (15).
The properties of the sharpness function C(γ) have been analyzed in detail in [13]. C(γ) has been proven to be a unimodal function, which has only one maximum value. The maximum value corresponds to a focused image.
The unknown variables in C(γ) are vx and vy, which are independent variables. The problem of Equation (15) can be simplified as two independent one-dimensional unimodal function optimization problems. Searching vx and vy, when γ equals to the real quadratic phase errors, C(γ) reaches the peak value. The corresponding vx and vy are the true velocity components of the target.
From Equations (16) and (17), a coarse estimation of vy can be acquired. As a result, the searching intervals of vy can be set centered on v y , which can substantially reduce the computational load. The searching intervals of v x can be set adequately large. As the searching process is constrained to a small image area containing the moving target, the algorithm is not time consuming.
Because C(γ) contains only one maximum, unimodal function optimization methods can be applied to effectively solve the optimization problem, such as the advance-and-retreat search method [20] and the golden-section search method [21]. Taking as an example the golden-section search method, whose flowchart is illustrated in Figure 4.
The searching steps can be summarized below:
(1)
Choose the search interval of vy: [vya,vyb].
(2)
Choose the search interval of vx: [vxa,vxb].
(3)
Let vy = vya; search vx in [vxa,vxb], set ε = 0.001. According to the process shown in Figure 4, obtain the maximum value C a * and the corresponding v a x * .
(4)
Let vy = vyi = vya + Δy, where Δy is the step size acquired by the flowchart in Figure 4. According to the process shown in Figure 4, obtain the maximum value C i * and the corresponding v i x * .
(5)
Repeat Step (4); until vy = vyb, and then obtain C b * and v b x * .
(6)
Let C max * = max { C a * , , C b * , the corresponding vx and vy are the estimates of the target’s velocity components.

4.3. Flow Chart of Parameter Estimation of a Ground Moving Target

Based on the principle presented above, the flowchart of the proposed approach is illustrated in Figure 5. The main steps of the ground moving target parameter estimation algorithm using sharpness optimization are summarized as follows:
(1)
Imaging of the raw SAR data.
(2)
Clutter suppression and detection of the moving target in the image domain.
(3)
Extraction of the small part of the image containing the moving target.
(4)
Corse estimation of the motion parameters.
(5)
Joint searching of vx and vy to obtain the quadratic phase error.
(6)
Calculation of the sharpness cost function C(γ).
(7)
Repeat Steps (5) and (6); when C(γ) achieves the maximum value, stop. The corresponding vx and vy are the accurate along-track and radial velocity of the target, respectively.

5. Experimental Results

In this section, we present and discuss a number of results obtained by applying the proposed moving target parameter estimation method to simulated SAR data. The simulation conditions are: aircraft velocity = 150 m/s, transmission wavelength = 0.03 m, pulse repetition frequency = 500 Hz, and distance between phase centers = 0.96 m. The parameters of the moving target are: radial velocity = 2 m/s, along-track velocity = 4.5 m/s and original azimuth position = 130 m. Following the procedure shown in Figure 5, following imaging, detection, and extraction, the image of the moving target is depicted in Figure 6a.
Because of the along-track velocity, the image is unfocused and spread along the azimuth direction. Both the advance-and-retreat search method and the golden-section search method are applied to solve the optimization problem. The iteration process is shown in Figure 7. From Figure 7, we can infer that the golden-section search method has a superior convergence rate than the advance-and-retreat method. For the moving target in Figure 6a, under the same computer conditions, the golden-section search method takes 1.292 s, whereas the advance-and-retreat search method takes 1.873 s to obtain the optimal estimation of the target velocities. The results show that the moving target parameter estimation proposed in this paper has a very low computation load. The estimated moving target parameters are: along-track velocity = 4.503 m/s, radial velocity = 2.032 m/s and original azimuth location = 124.39 m. The estimated parameters are in agreement with the simulated parameters. Figure 6b shows the focused image obtained by compensating the quadratic phase errors. We also compared the performance of the proposed parameter estimation method with the traditional ATI method. Figure 8a shows the radial velocity estimation error versus the signal-to-clutter ratio (SCR), and Figure 8b shows the along-track velocity estimation error versus the SCR. The proposed method is found to exhibit higher velocity estimation accuracy than the traditional ATI method.

6. Conclusions

This paper introduced a novel motion parameter estimation method using contrast optimization for a ground moving target. The echo model of a target in a SAR image was analyzed and the relationships between the quadratic phase errors and the target’s velocity components are derived first. Next, the concepts of image sharpness metrics were discussed in detail. Based on the sharpness optimization criterion, the peak magnitude of the sharpness cost function corresponds to the real phase errors, which are functions of the target’s velocity components. Using two-dimensional searching, we can obtain the along-track and radial velocities of the target. The proposed method can obtain the motion parameters and the focused image simultaneously, and obviously enhance the estimation accuracy. In the simulation section, the method was applied to simulated SAR data, and the results illustrated that the proposed method can effectively estimate the moving target parameters with better accuracy than the conventional ATI method.

Acknowledgments

The work described in this paper is supported by the National Natural Science Foundation of China (No. 51409214) and Natural Science Foundation research project of Shaanxi Province (No. 2016JQ6071). The authors would appreciate the anonymous reviewers for their valuable comments and suggestions for improving this paper.

Author Contributions

Jing Yu designed the approach and performed the computer simulation; Yaan Li analyzed the simulation data. Jing Yu wrote the paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Zhang, S.; Xing, M.; Xia, X.; Guo, R. Robust clutter suppression and moving target imaging approach for multichannel in azimuth high-resolution and wide-swath synthetic aperture radar. IEEE Trans. Geosci. Remote Sens. 2015, 53, 687–709. [Google Scholar] [CrossRef]
  2. Melvin, W.L. A STAP Overview. IEEE AES Mag. 2004, 19, 19–35. [Google Scholar] [CrossRef]
  3. Shu, Y.; Liao, G.; Yang, Z. Robust radial velocity estimation of moving targets based on adaptive data reconstruction and subspace projection algorithm. IEEE Geosci. Remote Sens. Lett. 2014, 11, 1101–1105. [Google Scholar] [CrossRef]
  4. Feng, Z.; Li, Z.; Zheng, B. A new approach to ground moving target detection and location based on two-look processing for the single channel SAR system. J. Xidian Univ. 2006, 33, 673–681. [Google Scholar]
  5. Zheng, M.; Yang, R. SAR moving targets detection based of DPCA and interferometric processing. J. Electron. Inform. Technol. 2003, 25, 1525–1530. [Google Scholar]
  6. Raney, R.K. Synthetic aperture imaging radar and moving targets. IEEE Trans. Aerosp. Electron. Syst. 1971, 7, 499–505. [Google Scholar] [CrossRef]
  7. Yang, S.; Li, Y.; Zhang, K.; Tang, W. Multiple-parameter estimation method based on spation-temporal 2-D processing for bistatic MIMO radar. Sensors 2015, 15, 31442–31452. [Google Scholar] [CrossRef] [PubMed]
  8. Wang, Y. Radar imaging of non-uniformly rotating targets via a novel approach for multi-component AM-FM signal parameter estimation. Sensors 2015, 15, 6905–6923. [Google Scholar] [CrossRef] [PubMed]
  9. Chen, C.; Andrews, H.C. Target-motion-induced radar imaging. IEEE Trans. Aerosp. Electron. Syst. 1980, 16, 2–14. [Google Scholar] [CrossRef]
  10. Elachi, C.; Bicknell, T.; Rolando, L.J.; Wu, C. Spaceborne synthetic-aperture imaging radars: Applications, techniques, and technology. Proc. IEEE 1982, 70, 1174–1209. [Google Scholar] [CrossRef]
  11. Berizzi, F.; Corsini, G. Autofocusing of inverse synthetic aperture radar images using contrast optimization. IEEE Trans. Aerosp. Electron. Syst. 1996, 32, 1185–1191. [Google Scholar] [CrossRef]
  12. Gao, Y.; Yu, W.; Liu, Y.; Wang, R. Autofocus algorithm for SAR imagery based on sharpness optimization. Electron. Lett. 2014, 50, 830–832. [Google Scholar] [CrossRef]
  13. Barbarossa, S.; di Lorenzo, P.; Vecchiarelli, P. Polynomial phase signals: An application to SAR imaging of moving targets. IEEE Trans. Signal Process. 2014, 62, 4375–4389. [Google Scholar] [CrossRef]
  14. Schulz, T.J. Optimal sharpness function for SAR autofocus. IEEE Signal Process. Lett. 2007, 14, 27–30. [Google Scholar] [CrossRef]
  15. Morrison, R.L.; Do, M.N.; Munson, D.C. SAR image autofocus by sharpness optimization: A theoretical study. IEEE Trans. Image Process. 2007, 16, 2309–2321. [Google Scholar] [CrossRef] [PubMed]
  16. Kim, D.; Moon, W.M.; Moller, D.; Imel, D.A. Measurements of ocean surface waves and currents using L- and C-band along-track interferometric SAR. IEEE Trans. Geosci. Remote Sens. 2003, 41, 2821–2832. [Google Scholar]
  17. Wang, W.; Shao, H. Two-antenna SAR with waveform diversity for ground moving target indication. IEEE Trans. Geosci. Remote Sens. 2014, 11, 2154–2158. [Google Scholar] [CrossRef]
  18. Cerutti-Maori, D.; Sikaneta, I. A generalization of DPCA processing for multichannel SAR/GMTI radars. IEEE Trans. Geosci. Remote Sens. 2013, 51, 560–572. [Google Scholar] [CrossRef]
  19. Dawidowicz, B.; Kulpa, K.S.; Malanowski, M.; Misiurewicz, J.; Kulpa, K.S. DPCA detection of moving targets in airborne passive radar. IEEE Trans. Aerosp. Electron. Syst. 2012, 48, 1347–1357. [Google Scholar] [CrossRef]
  20. Zhang, J.; Xiao, Z.; Tan, Y.; He, X. Hybrid particle swarm optimizer with advance and retreat strategy and clonal mechanism for global numerical optimization. IEEE Congr. Evol. Comput. 2008, 28, 2059–2066. [Google Scholar]
  21. Guler, O. Foundations of Optimization; Springer: New York, NY, USA; Dordrecht, The Netherland; Heidelberg, Germany; London, UK, 2010. [Google Scholar]
Figure 1. Geometry of a SAR system with a moving target.
Figure 1. Geometry of a SAR system with a moving target.
Sensors 16 01017 g001
Figure 2. A part of a SAR image.
Figure 2. A part of a SAR image.
Sensors 16 01017 g002
Figure 3. Relationship between the sharpness value and the quadratic phase.
Figure 3. Relationship between the sharpness value and the quadratic phase.
Sensors 16 01017 g003
Figure 4. Flowchart of the golden-section search method.
Figure 4. Flowchart of the golden-section search method.
Sensors 16 01017 g004
Figure 5. Flowchart of the moving target parameter estimation process.
Figure 5. Flowchart of the moving target parameter estimation process.
Sensors 16 01017 g005
Figure 6. (a) Unfocused image of the moving target; (b) Focused image of the moving target.
Figure 6. (a) Unfocused image of the moving target; (b) Focused image of the moving target.
Sensors 16 01017 g006
Figure 7. Estimation errors of the quadratic phase coefficient versus the iteration times.
Figure 7. Estimation errors of the quadratic phase coefficient versus the iteration times.
Sensors 16 01017 g007
Figure 8. (a) Radial velocity estimation error versus the SCR; (b) Along-track velocity estimation error versus the SCR.
Figure 8. (a) Radial velocity estimation error versus the SCR; (b) Along-track velocity estimation error versus the SCR.
Sensors 16 01017 g008

Share and Cite

MDPI and ACS Style

Yu, J.; Li, Y. Parameter Estimation of a Ground Moving Target Using Image Sharpness Optimization. Sensors 2016, 16, 1017. https://doi.org/10.3390/s16071017

AMA Style

Yu J, Li Y. Parameter Estimation of a Ground Moving Target Using Image Sharpness Optimization. Sensors. 2016; 16(7):1017. https://doi.org/10.3390/s16071017

Chicago/Turabian Style

Yu, Jing, and Yaan Li. 2016. "Parameter Estimation of a Ground Moving Target Using Image Sharpness Optimization" Sensors 16, no. 7: 1017. https://doi.org/10.3390/s16071017

APA Style

Yu, J., & Li, Y. (2016). Parameter Estimation of a Ground Moving Target Using Image Sharpness Optimization. Sensors, 16(7), 1017. https://doi.org/10.3390/s16071017

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