Abstract

In the current scenario of high-range resolution radar and noncooperative target, the rotational motion parameters of the target are unknown and migration through resolution cells (MTRC) is apparent in the obtained inverse synthetic aperture radar (ISAR)images, in both slant-range and cross-range directions. In the case of the high-speed maneuvering target with a small value of rotation, the phase retrieval algorithm can be applied to compensate for the translational motion to form an autofocusing image. However, when the target has a relatively large rotation angle during the coherent integration time, phase retrieval method cannot get an acceptable image for viewing and analysis as the location of the scatterer will not be true due to the Doppler shift imposed by the target’s rotational motion. In this paper, a novel ISAR imaging method for maneuvering targets based on phase retrieval and keystone transform is proposed, which can effectively solve the above problems. First, the keystone transform is used to solve the MTRC effects caused by the rotation component. Next, phase errors caused by the remaining translational motion will be removed by employing phase retrieval algorithm, allowing the scatterers are always kept in their range cells. Finally, the Doppler frequency shifts of scatterers will be time invariant in the phase of the received signal. Furthermore, this approach does not need to estimate the motion parameters of the target, which simplifies the processing steps. The simulated results demonstrate the validity of this method.

1. Introduction

Inverse synthetic aperture radar imaging is a radar technique to obtain a high-resolution image of moving targets [16]. In general, the target is assumed to have both translational and rotational motion, where the target’s rotational motion is used to attain cross-range resolution. The conventional ISAR motion compensation algorithms [7, 8] are dedicated to the imaging of the maneuvering target with simple motion. Once the target contains complex motions such as large angle rotation, high speed, and vibration, these algorithms cannot achieve high-resolution images.

Various techniques are performed to mitigate or eliminate these unwanted motion effects in ISAR imaging [911], among which some of the popular algorithms include cross-correlation method [12, 13] and minimum entropy method [1416]. However, the computational complexity of the minimum entropy algorithm is too large, and the cross-correlation method can only measure the target’s velocity within a small interval, and the accuracy of the algorithm is low. Phase retrieval is a method developed for optical imaging in recent years, and it is gradually applied in various fields [17, 18]. The principle of phase retrieval is to calculate the phase distribution of the light field by measuring the intensity distribution of the optical field. Shi et al. [19, 20] believe that the ISAR echo module of the maneuvering target is not affected by its radial motion, that is, phase retrieval can be applied to compensate for translational motion (radial velocity and radial acceleration). Furthermore, the method does not need to estimate the target motion parameters, and the autofocus images can be obtained under the condition that the rotation of the target is small for the coherent integration time of the radar. However, when the rotation angle becomes larger, phase retrieval algorithm has difficulties in eliminating the MTRC caused by the rotational motion.

Keystone transform was originally proposed in Ref. [21] for SAR imaging of moving targets, and it was called keystone remapping. To obtain high-resolution SAR images, the keystone transform is used to rescale the time axis for each frequency, and it corrects range migration in SAR images. The keystone transform has been widely used in ISAR autofocus to correct range cell migration correction (RCMC). For example, keystone transform (KT) are utilized to transform the energy of all the scatterers into one range cell [22]. Modified keystone transforms can be employed to simultaneously transform all multicomponent linear frequency modulation subechoes into multicomponent single-frequency signals [23]. Moreover, Wang and Kasilingam [24, 25] apply the first-order generalized keystone transform to remove the linear range cell migration induced by the target’s uniform rotation. The effects of higher-order terms can be solved by performing second-order or higher-order generalized keystone transforms [26, 27]. Usually, the unfavorable effect induced by the rotational motion is less severe than that induced by translational motion. In this way, keystone transform can be first applied for removing the range cell migration due to the target’s rotation. However, keystone transform cannot overcome the issue of severe scatters’ walk with respect to slow time imposed by translational acceleration and large translational velocity, which makes the target’s resultant image defocused in range and cross-range directions. To make the Doppler frequency shifts constant among the range cells in radar signal processing, it is necessary to combine the keystone transform and phase retrieval.

Based on the above analysis, a novel ISAR imaging algorithm for maneuvering targets based on phase retrieval and keystone transform is proposed. First, keystone transform is performed to correct the MTRC due to the rotation motion. Then, we use phase retrieval to compensate for the target’s translational motion component and solve the range walk effect. Finally, the returned signals from the same scatterer are always kept in the same range cell, and the time-invariant Doppler shift can be obtained.

The rest of the paper is organized as follows: maneuvering target echo signal analysis is shown in Section 2. In Section 3, we focus on the theory of keystone transform and phase retrieval. Experimental simulation analysis is presented in Section 4, and conclusions are given in Section 5.

2. Maneuvering Target Echo Signal Analysis

The space-fixed or global system and the body-fixed or local system , which is rigidly fixed in the body, are the two coordinate systems commonly used to describe a maneuvering target as shown in Figure 1. To describe the rotation of the target, another reference coordinate system is introduced, which is parallel to the global system with its origin at the origin of the target body-fixed system. The uniformly rotating target contains scatterers, and the random scatterer on the target has translational motion and rotational motion.

The baseband echo signal received by the radar can be expressed as where is fast time; is slow time, where represents pulse repetition interval; and are the reflectivity density function from any point scatterer at and the distance from the radar at instantaneous time , respectively; denotes the normalized echo envelope; and is the carrier frequency.

Transforming from fast time domain to baseband frequency domain, it can be represented to yield where, represents the Fourier transform of .

The target’s translational range distance, , can be modelled in terms of motion parameters:

Here, and are the target’s translational velocity and the acceleration, respectively. is the initial range of the target. Therefore, the distance from the scatterer to the radar can be written as where represents the rotational angle of the target with respect to the radar line of sight axis, and , where denotes the angular velocity of the target. Assuming that the initial angle of the target is equal to zero and the CPI (coherent processing interval) is short, we can get

Therefore, the backscattered signal can be rewritten as where, , which is a constant.

The first term in the phase represents the range distribution of scatters and can be ignored for the imaging process. The second and third terms denote echo envelope shift caused by the translational motion and the rotational motion, respectively. The last term indicates the carrier phase changes caused by the Doppler effect. To have a motion-free range-Doppler image of the target, the envelope shift in the second and third terms cannot expand to more than one range resolution cell; otherwise, the range of the scatterer, , will migrate from its original position (i.e., ) to a wrong position (i.e., ). It is worth noting that there is a coupling between signal frequency and slow time in equation (6), which can be regarded as the MTRC effect or range drift shift. Therefore, the RCMC procedure should be applied to the second and third phase errors terms to suppress them.

According to phase retrieval theory, we know that autofocus imaging can be achieved when the raw data magnitude of the maneuvering target is approximate to the stationary one in the ISAR system [17]. Obviously, when the “” is not small, this rotational motion has to be compensated for by applying an initial compensating procedure, that is, using phase retrieval alone does not eliminate the MTRC caused by the third term in equation (6). To effectively remove the rotational motion effects, we will perform the keystone transform.

3. Algorithm Theory and Application

3.1. Keystone Transform

The keystone transform was used to rescale the slow time axis for each frequency [21] by where is the rescaled slow time axis.

Then, the frequency domain signal becomes

As is much larger than , we have

We can see that frequency and slow time have been decoupled in the second term of the phase function by applying the keystone transform, and the envelope shift position of the echo signal stays at the range . In fact, the quadratic term can certainly be neglected compared with the normal size of the range resolution cell (e.g., 1 m). For example, if , , and , the quadratic term is less than 0.1 m. The Doppler frequency shift due to carrier phase changes caused by motion can be calculated by taking the time derivative of the phase in equation (9) as

The Doppler frequency shift changes linearly due to the acceleration during the CPI, resulting in frequency dispersion of the echo signal and negatively affecting the range resolution. Additionally, when the target’s translational velocity becomes larger, Doppler blurring occurs in the resulting ISAR image, and more importantly, the image smearing can be drastic for fast moving targets in range direction. Therefore, an algorithm needs to be applied to overcome such the range walk issue and make the Doppler frequency shifts constant.

3.2. Phase Retrieval Algorithm

Phase retrieval is the process of finding the optimal solution to the phase problem of the recovered signal. Given a complex signal with the amplitude and phase : where is an M-dimensional spatial frequency coordinate and is an M-dimensional spatial coordinate. The principle of phase retrieval is to find a set of constrained phases that satisfy the measured amplitude, that is, the known Fourier transform information is applied to recover the Fourier phase information , then the distribution function is obtained by inverse Fourier transform.

However, in signal processing, the correctness of the recovery result cannot be guaranteed due to the lack of a priori information, that is, it is likely to lead to nonuniqueness of the phase retrieval result. There are two reasons for this: the first is caused by the spatial displacement, conjugate inversion, and global phase shift and the second is related to the dimension of the unknown signal. In order to eliminate the phase ambiguity and ensure the uniqueness of the recovery result when performing the phase retrieval, we use the error phase of the blurred image obtained by the traditional ISAR imaging algorithm (keystone transform) as a priori phase information. In this paper, the classic phase retrieval algorithm, namely, the oversampling smoothness (OSS) phase retrieval [19], is applied and a priori information is added to it. Besides, the support domain size of the improved OSS algorithm is set with respect to the blurred target image.

Based on the support domain constraints of the traditional hybrid input-output (HIO) algorithm [28], the OSS algorithm adds iterative steps of frequency domain filtering. To search for the phase information of the recovered signal in solution space, it iterates back and forth between the Fourier and image domains. The following steps are the process of the OSS algorithm from the th to the th iteration at each run. (1) is the input signal with the random phase. Obtain a Fourier pattern by applying the Fourier transform(2)Retain the phase information of and replace the magnitude of it with a known Fourier magnitude , then a new complex—valued function can be formed. It should be pointed out that is the magnitude of the ISAR echo signal (3)Perform inverse Fourier transform to and get a new image (4)Revise based on HIO equation and generate the new where represents a finite support and the value of the parameter is between 0.5 and 1(5)Calculate the image for the th iteration where is the Fourier transform of and is a normalized Gaussian function in the Fourier domain. It is defined as

The smoothing filter is only applied to the density outside the support domain. By changing parameter , the width of the Gaussian filter can be adjusted to overcome the issue of high-frequency information outside the support, while the density inside the support domain is not disturbed.

By applying OSS phase retrieval algorithm, the maneuvering targets module in (9) becomes

From (16), ISAR echo module is not affected by the target’s translational velocity and the acceleration with respect to the rescaled slow time . As an autofocus processing procedure, OSS phase retrieval effectively removes the effects of Doppler blurring and range walk. Specifically, the first term in the absolute value can be ignored for the imaging process as mentioned in Section 2. Noting that the Doppler frequency shift is equal to , which is constant between the scatterers. Therefore, the module in (16) is approximate to the ideal module of the target autofocus imaging, so the RCMC procedure has been performed by OSS phase retrieval. In other words, phase retrieval performs a fine phase compensation, and the target’s resultant image will be completely focused in range and cross-range directions.

The overall flow diagram of the proposed method is shown in Figure 2.

4. Simulation Analysis

To verify the feasibility of this method, we conduct several experiments. A finite number of perfect point scatterers are assumed to be located as shown in Figure 3(a), and the module of the raw data is plotted in Figure 3(b). Table 1 shows the radar parameters that are used in the simulation.

4.1. Target Imaging Results with Different Translational Motion Parameters

To obtain a more intuitive view, the different parameters from each experiment are listed in Table 2.

When the angular velocity of the target is large, Figure 4(a) shows the ISAR image with applying the RD algorithm, where the blur occurs due to the range cell migration caused by the rotational motion. As can be seen from Figure 4(b), because of the small translational velocity of the target, keystone transform can correct the MTRC caused by rotation and translation and obtain a focused ISAR image. Keeping the translation parameters unchanged and increasing the rotation velocity, then the resultant rotation motion-compensated ISAR image after applying keystone transform are still relatively focused in Figure 5(b), which proves that keystone transform can perform the RCMC procedure in this case. When the radial velocity becomes larger, by utilizing keystone transform, the ISAR image is still smearing in Figure 6(b), which is induced by larger carrier phase changes with respect to Figure 4(b) imposed by the Doppler effect. In addition, due to the existence of radial acceleration, Figure 7(b) shows the keystone transform cannot handle the difficulties of the nonlinear Doppler shift. The proposed method can effectively eliminate carrier phase errors caused by the target translational motion while compensating for the rotational motion to obtain a high-resolution image.

In Figure 8, the range and azimuth profiles of the scattering point on the aircraft head marked by the red circle related to Figures 4(a), 4(b), and 4(c) are plotted, which intuitively demonstrates profile comparison of the target imaging through RD, keystone transform, and the proposed method. Due to the target’s translational motion and rotational motion, the peak sidelobe ratio produced by RD algorithm is the largest, indicating that the point scatterer has the most serious range and azimuth defocusing. Obviously, the profiles with respect to keystone transform show the range broadening are alleviated, which indicates that keystone transform removes the envelope shift related to fast time in equation (6) under the condition of small velocity. As can be clearly seen from the figures, the proposed method has the minimum peak sidelobe ratio and therefore possesses good focusing performance in the process of reconstructing images.

4.2. Target Imaging Results with Different Rotational Motion Parameters

Figures 9(a) and 9(b) show RD algorithm produces degradations such as blurring and smearing. In the case of high speed and small rotation angle, the scatters’ range walk caused by the target’s rotation motion does not exceed a range cell. Therefore, as shown in Figures 9(b) and 9(c), high-quality images can be obtained by using both the minimum entropy method and the OSS phase retrieval algorithm. However, by comparing the two resultant motion-compensated ISAR images, the one with applying the OSS phase retrieval algorithm is more focused. When the rotation velocity is much larger, the minimum entropy method has difficulty in eliminating the MTRC caused by the rotation motion. Figure 10(b) clearly demonstrates the rotational motion-based defocusing is noted in the ISAR image. In this case, the application of phase retrieval alone does not remove the echo envelope shift induced by the rotational component in equation (6), so some scattering points of the target in Figure 10(c) are missing. In Figure 10(d), the resultant motion-free ISAR image is well focused, which confirms that the proposed method can perform autofocus imaging of the maneuvering target with high-speed and large angle.

Figures 11(a) and 11(b) show profile comparison of the scatterer point on the airplane head marked by the red circle in Figures 10(a), 10(b), and 10(d). Seen from the peak sidelobe ratio of the two figures above, the imaging result of the minimum entropy method clearly demonstrates the broadening of the range direction and azimuth direction, which can also indicate the poor performance of the minimum entropy method under the condition of the high-speed maneuvering target with large angle. Compared with the RD algorithm and the minimum entropy method, the proposed method has a smaller peak sidelobe ratio, that is, it obtains a more focused image in both directions.

5. Conclusions

In this paper, we introduce a new method for maneuvering target ISAR imaging based on phase retrieval and keystone transform. By applying the keystone transform, the range cell migration caused by the target’s rotation motion can be corrected to obtain an initial phase compensation. Furthermore, phase retrieval is dedicated to performing the remaining translational motion compensation with respect to slow time, so that the time-invariant Doppler shifts can be obtained. Finally, the method can obtain resultant motion-free ISAR image after completely removing the phase error of the received signal, which clearly demonstrates that the unwanted effects due to the target’s motion are eliminated.

In the experimental simulation, we compare the imaging results of the keystone transform, phase retrieval, and minimum entropy method with the imaging results of the proposed method. In the case of the high-speed maneuvering target with a large value of rotation, the former three traditional methods cannot successfully eliminate the motion effects. Consequently, the motion-compensated ISAR image is obtained while the effect of motion in the scattered field can be mitigated by applying the proposed method, wherein the target’s scattering centers are very well displayed with good resolution.

In summary, the results of this study provide a new way of thinking for ISAR imaging problem of the maneuvering target with complex motion. At the same time, it should be pointed out that the approach has a high accuracy, and it does not need to estimate the motion parameters of the target.

Data Availability

No data were used to support this study.

Disclosure

The funding sponsors had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; and in the decision to publish the results.

Conflicts of Interest

The authors declare no conflict of interest.

Authors’ Contributions

All the authors made significant contributions to this work. Hongyin Shi proposed the novel ISAR imaging method. Ting Yang performed the experiments and wrote the paper. Yue Liu performed part of the experiments. Jingjing Si revised the manuscript.

Acknowledgments

This work was supported by the National Natural Science Foundation of China under Grants 61571388 and 61701429 and Natural Science Foundation of Hebei Province under Grant F2016203251.