Next Article in Journal
3D Airborne EM Forward Modeling Based on Time-Domain Spectral Element Method
Previous Article in Journal
Multi-Source EO for Dynamic Wetland Mapping and Monitoring in the Great Lakes Basin
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Time-Frequency Reversion-Based Spectrum Analysis Method and Its Applications in Radar Imaging

1
National Laboratory of Radar Signal Processing, Xidian University, Xi’an 710071, China
2
Center for Advanced Communications, Villanova University, Villanova, PA 19085, USA
*
Author to whom correspondence should be addressed.
Submission received: 30 December 2020 / Revised: 31 January 2021 / Accepted: 1 February 2021 / Published: 8 February 2021
(This article belongs to the Section Remote Sensing Image Processing)

Abstract

:
Spectrum analysis (SA) plays an important role in radar signal processing, especially in radar imaging algorithm design. Because it is usually hard to obtain the analytical expression of spectrum by the Fourier integral directly, principle of stationary phase (POSP)-based SA is applied to approximate this integral. However, POSP requires the phase of the signal to vary rapidly, which is not the case in circular synthetic aperture radar (SAR) and turntable inverse SAR (ISAR). To solve this problem, a new SA method based on time-frequency reversion (TFRSA) is proposed, which utilizes the relationship of the Fourier transform pairs and their corresponding signal phases. In addition, the connection between the imaging geometry and time-frequency relationship is also analyzed and utilized to help solve the time-frequency reversion. When the TFRSA is applied to the linear trajectory SAR, the obtained spectrum expression is the same as the result of POSP. When it is applied to ISAR, the spectrum expressions of near-field and far-field are derived and their difference is found to be position-independent. Based on this finding, an extended polar format algorithm (EPFA) for near-field ISAR imaging is proposed, which can solve the distortion and defocusing problems caused by traditional ISAR imaging algorithms. When it is applied to the circular SAR (CSAR), a new and efficient imaging method based on EPFA is proposed, which can solve the low efficiency problem of conventional BP-based CSAR imaging algorithms. The simulations and real data processing results are provided to validate the effectiveness of proposed method.

Graphical Abstract

1. Introduction

Radar can realize long distance target detection, localization, imaging, classification and recognition in all-day and all-weather conditions, and has become a significant tool in many applications. Among these applications, Synthetic aperture radar (SAR) [1,2,3] and inverse SAR (ISAR) [4,5,6] imaging attracted much attention, because they can provide two-dimension high-resolution images, which contain many features, for example, the size and shape. Different from optical imaging, radar imaging produces images by complex computation. The resolution of radar image in the range dimension is accomplished by transforming the large signal bandwidth, whereas the fine resolution in azimuth stems from the difference in observation angles. With the development of advanced hardware, the available radar bandwidth is becoming larger, which in turn produces ultrahigh-resolution radar images. Recently, the microwave-photonic technique [7,8,9], which combines the advantages of electromagnetic wave and photonics techniques, is applied to the radar system to produce an ultra-large 10 GHz bandwidth radar signal. This microwave-photonic radar can generate radar images with 1.5 cm resolutions [10,11]. In addition to the large bandwidth, wide-angle radar imaging is also a worthwhile study direction, because it can not only produce ultrahigh azimuth resolution but also provide anisotropy characteristic, which is very useful for target classification and recognition.
The design of radar imaging algorithm is usually based on the spectrum of the radar received signal [3]. Different constructions of virtual synthetic apertures always lead to different imaging algorithms. Thus, the spectrum analysis (SA) is very important to radar imaging. The Fourier integral is a direct way to perform SA. However, the integral is usually difficult to evaluate, and approximation is often applied to obtain the Fourier integral. Principle of stationary phase (POSP) is one of the popular methods used for asymptotic analysis [12]. Because the phase of the radar signal changes rapidly, the cancellation of sinusoids would cause the integral to be zero except near the stationary point. Due to its simplicity, POSP is widely used in radar signal processing and imaging algorithm design [13,14,15], especially for the linear trajectory SAR (LTSAR). Because the phase is assumed to be changing rapidly, POSP is considered applicable when the signal’s time and bandwidth product (TBP) is larger than 100, which is satisfied in most of the LTSAR applications.
For curved trajectory SAR (CTSAR), like the circular SAR (CSAR), the TBP condition cannot be satisfied for returns of scatterers near the scene center, and thus POSP cannot be applied to CSAR. This leads to the difficulty of obtaining the spectrum expression and designing efficient Doppler domain imaging algorithm for CSAR. Back-projection algorithm (BPA) is a well-known and accurate time-domain imaging algorithm, which is applicable to any trajectories of SAR. However in BPA, the value of every image pixel is calculated by integrating using the accurate slant range, this incurs huge computational costs. In order to improve the efficiency of imaging, many accelerated BPAs are proposed [16,17]. However, their efficiencies are still not comparable to that of frequency-domain algorithms, for example, chirp scaling algorithm (CSA), frequency scaling algorithm (FSA), and range migration algorithm (RMA) [3].
Different from SAR, which usually observes the stationary scene by moving the radar along the predesigned trajectory, ISAR usually observes non-cooperative moving targets in the static platform, which means its corresponding virtual synthetic aperture is more complex and this leads to complicated Doppler spectrum. In traditional ISAR imaging, the distance between the radar and target is usually considered much larger than the size of the target. Thus, the plane wave assumption is applied to decompose the slant range into the two uncoupled parts: translational motion (TM) and rotational motion (RM). TM is the same for all scatterers of the target and has no contribution to ISAR imaging. There are many methods proposed to compensate the TM to produce an equivalent turntable model for ISAR imaging [18,19,20]. Moreover, many imaging algorithms were proposed for the turntable ISAR imaging under the plane wave assumption, for example, the range-Doppler algorithm (RDA), keystone transform (KT)-based imaging algorithm, and polar format algorithm (PFA) [5].
However, with the increasing demand for finer image resolution, the accuracy requirement of error compensation for radar imaging has also increased. Besides, the distance between radar and target for some applications are becoming smaller, for example, target detection and ranging for autonomous vehicles, harbor surveillances, and airport security checks [21,22,23]. In these situations, termed near-field ISAR imaging, the plane wave assumption may not be accurate enough for high-resolution ISAR imaging and the corresponding image quality would be degraded. In the near-filed situations, traditional ISAR imaging algorithms, like KT-based method and PFA, would produce images with geometric distortion and defocusing problems. Besides, for the near-field turntable ISAR, where a target is placed on a turntable with only rotational motion, the POSP is also inapplicable because the returns of the scatterers near the turntable center are narrow-band signals.
In order to overcome the limitations of POSP, a new SA method based on time-frequency reversion (TFR) is proposed in this paper. By utilizing the relationship between the Fourier transform (FT) pairs and their corresponding signal phases, the frequency-time distributed lines (FTDL) [24,25] can be obtained and the spectrum can be evaluated by integrating the FTDL. Besides, the connection between the time-frequency relationship and the imaging geometry is also analyzed, leading to a simplification of the time-frequency reversion process, and at the same time, providing a physical insight of the problem. When the proposed TFR-based spectrum analysis (TFRSA) method is applied to LTSAR, the derived spectrum expression is the same as the result of POSP, which validates the accuracy of TFRSA. When TFRSA is applied to the near-field ISAR, its spectrum expression is also derived. It is worth noting that the derived difference between near-field and far-field spectrums is a position-independent expression. With this finding, an extended polar format algorithm (EPFA) is proposed for the near-field ISAR imaging. When the TFRSA is applied to the CSAR, its spectrum is approximated to a near-field ISAR spectrum with position-dependent constant error, and an efficient CSAR imaging algorithm is proposed based on the EPFA.
The remainder of the paper is organized as follows. In Section 2.1, the signal model and adopted imaging geometry are given. In Section 2.2, the proposed TFRSA method is introduced. In Section 2.3, the TFRSA is applied to the LTSAR to derive its spectrum expression and this expression is compared with the result of POSP. In Section 2.4, the TFRSA is applied to the near-field ISAR to derive its spectrum expression and an EPFA is also proposed for near-field ISAR imaging. In Section 2.5, the TFRSA is applied to the CSAR, an efficient imaging method based on EPFA is proposed for CSAR. In Section 3, simulations and real data processing results are provided to validate the effectiveness of the proposed TFRSA and EPFA imaging methods. In Section 4, an explanation on the performances of different algorithms is given and the future research directions are also discussed. Finally, a conclusion is given in Section 5.

2. Materials and Methods

2.1. Radar Imaging Geometry and Signal Model

Different from optical imaging, which is based on the accumulation of photons, radar imaging produces images by complex computation. Radar transmits electromagnetic wave and receives the returns and samples them into digital data. Then, imaging processing is applied to the digital data to generate a radar image.
For SAR, its imaging geometry is illustrated in Figure 1a, in which radar is observing the static scene and moving along a predesigned trajectory. A Cartesian coordinate system X O Y is built on the imaging plane, which is defined by the radar velocity vector and the line of sight (LOS) vector. The origin is set at the center of the scene. The observation angle θ of radar is defined by the angle between the radar position vector and the Y-axis, which is demonstrated in Figure 1a. The i-th scatterer on the ground, whose coordinates are x i g , y i g on the ground coordinate system X g O Y g , is projected onto the imaging plane with the coordinates x i , y i . Thus, the slant range between radar and i-th scatterer in the scene can be expressed as
R i S A R θ = r θ cos θ y i 2 + r θ sin θ x i 2 = r 2 θ 2 r θ x i sin θ + y i cos θ + x i 2 + y i 2 ,
where r θ denotes the distance between the radar and the origin.
For ISAR, its imaging geometry is illustrated in Figure 1b, in which the target is moving along arbitrary trajectory whereas the radar is static. In this case, the origin of the Cartesian coordinate system is set at the target center and the Y-axis is set parallel to the LOS. The motion of the target can be decomposed into two parts of rotation motion and radial translation, which corresponding to RM and TM, respectively. In order to express the rotational angle between two different positions of the target, the coordinate system X 1 O Y 1 is first rotated to be parallel with X O Y , which is demonstrated by step 1 in Figure 1b. Then, X 1 O Y 1 is moved along the LOS to overlap with X O Y , which is demonstrated by step 2. The rotational angle is defined by the ϑ demonstrated in Figure 1b. For the i-th scatterer in the target, its slant range can be expressed by
R i I S A R ϑ = r ϑ x i sin ϑ y i cos ϑ 2 + x i cos ϑ y i sin ϑ 2 = r 2 ϑ 2 r ϑ x i sin ϑ + y i cos ϑ + x i 2 + y i 2 ,
where r ϑ denotes the distance between radar and the target center and x i , y i depicts the coordinates of the i-th scatterer on the target when the rotational angle is 0. Comparing Equation (1) with Equation (2), we see that they have the same form, which means they are equivalent and they can be transformed into each other. Thus, for convenience, in the followings, we only take SAR for an example to introduce the signal model.
Assuming that radar transmits linear frequency modulation (LFM) signal and there are K scatterers, after down conversion and range compression, the received signal in range wavenumber and azimuth angle domain can be expressed as [26]
s k r , θ = i = 1 K A i rect k r Δ k r exp j k r R i θ = i = 1 K A i rect k r Δ k r exp j k r r 2 θ 2 r θ x i sin θ + y i cos θ + x i 2 + y i 2 ,
where A i signifies the back scatter coefficient of the i-th scatterer, rect   represents the rectangular window, and k r and Δ k r denote the range wavenumber and range wavenumber bandwidth, respectively. The angle θ is an azimuth variable and it can also be replaced by other azimuth variables like azimuth time t or azimuth radar position X r in the specific situations. The signal spectrum of the radar returns in this paper refers to applying azimuth FT along the azimuth variable of the signal in Equation (3).

2.2. Time-Frequency Reversion-Based Spectrum Analysis Method

Assuming θ is a function of azimuth time t and it is expressed as θ t , we can rewrite Equation (3) using the azimuth variable t as
s k r , t = i = 1 K A i rect k r Δ k r exp j k r r 2 t 2 r t x i sin θ t + y i cos θ t + x i 2 + y i 2 ,
By performing azimuth FT for the variable t, the Doppler spectrum of Equation (4) can be obtained by
S k r , f D = s k r , t exp j 2 π f D t d t ,
where f D signifies the Doppler frequency. Because the Fourier integral in Equation (5) is usually difficult to evaluate, an approximation is often applied to Equation (5) to obtain a simple result.
POSP is one of the common techniques to make this approximation. The main idea of POSP relies on the cancellation of sinusoids with rapidly varying phases. The phase of the radar signal usually changes rapidly due to the high carrier frequency, so the POSP is suitable for radar signal processing and it is widely used in SAR imaging applications, for instance, the FT of the LFM signal and LTSAR Doppler spectrum analysis.
Taking the FT for an example, assuming the analyzed signal is expressed as
g t = A t exp j φ t ,
where A t and φ t stand for the amplitude and phase of the signal g t , respectively, and the FT of g t in Equation (6) can be evaluated by
S ω = A t exp j φ t exp j ω t d t = A t exp j φ t ω t d t = A t exp j ϕ t d t ,
where ω signifies the frequency variable. Assuming the stationary phase point of Equation (7) is t = t * and approximating the phase term ϕ t by its second-order Taylor series, which yields
ϕ t φ t * ω t * + 1 2 φ t * t t * 2 ,
where ϕ denotes the second derivative of ϕ t about t. By applying the POSP to the integral of Equation (7), the result can be expressed by
S ω A t * exp j φ t * ω t * + j 1 2 φ t * t t * 2 d t = 2 π ϕ t * A t * exp j φ t * ω t * ± π 4 ,
where 2 π ϕ t * is a constant number and it is often ignored. From Equation (9) we find that POSP can simplify the FT integral and produce good approximated result, so it is widely used in radar signal processing.
However, the phase of the radar signal does not always change rapidly, for example, the returns of the CSAR and turntable ISAR. In these cases, the POSP would cause large error and the analyzed results are unusable. In order to solve this problem, we propose a new TFRSA method, which utilizes the relationship of the Fourier pairs and their corresponding signal phase. Here, we take the Fourier pairs of time and frequency as an example. The frequency can be expressed as the derivative of the signal phase of Equation (6), which is given by
ω = d φ t d t = φ t ,
where φ depicts the derivative of the signal phase. Equation (10) is also called the time-frequency distribution lines (TFDL) [24,25], which is used to describe the characteristic of the signal. Moreover, for the spectrum of Equation (7), the time can also be expressed by the negative derivative of its phase, which is given by
t ω = d φ S ω d ω ,
where φ S ω stands for the phase of the signal spectrum. Equation (11) is also called the FTDL. TFDL can be regarded as a 90 degrees rotation of the FTDL [27,28], which means that we can obtain the FTDL from the TFDL. Then, the spectrum phase can be calculated by integrating the FTDL along the frequency variable ω , which yields
φ S ω = t ω d ω = φ 1 ω d ω ,
where φ 1 ω denotes the inverse function of φ t . Because the amplitude A t is a low-frequency narrow band signal, and its FT can be regarded as an impulse function near the zero frequency, its effect on the signal spectrum can be regarded as a constant number and it is often ignored. This spectrum analysis method is called TFRSA. Furthermore, the FTDL can also be obtained from the imaging geometry, which utilizes the relationship between the instantaneous squint angle (ISA) and the spectrum variable. The ISA of the i-th scatterer is defined as the angle between the normal vector v of radar velocity and the vector from radar pointing to the scatterer, which is illustrated in Figure 1a. The ISA can also be expressed by the spectrum variable as follows
f D i = 2 v i n s λ sin θ i i n s sin θ i i n s = λ f D i 2 v i n s ,
where f D i and θ i i n s depict the Doppler frequency and the ISA of the i-th scatterer, respectively, λ and v i n s signify the wavelength of the transmitted signal and the instantaneous velocity of the radar, respectively. Equation (13) reveals that there is a one-to-one mapping relationship between frequency and ISA. Combing the geometry relationship and Equation (13) can help us simplify the spectrum derivation and have a better understanding of the TFRT. In the following, we will apply the TFRSA to three specific applications to validate its effectiveness.

2.3. Linear Trajectory SAR Spectrum Analysis Based on TFRSA

Linear trajectory sweep mode is one of the most common modes in the SAR, and its spectrum result has been well studied by the POSP and validated in many applications. In this subsection, we apply the TFRSA to the LTSAR and compare the derived spectrum with the result of POSP. The imaging geometry of the LTSAR is illustrated in Figure 2. Here, we use the radar azimuth position X r as the azimuth variable and the corresponding spectrum azimuth variable is the azimuth wavenumber k x . From the geometry relationship in Figure 2, the ISA of the i-th scatterer can be expressed as
tan θ i i n s = X r x i R b y i ,
where R b signifies the range between the scene center and the linear trajectory. Using the spectrum azimuth variable k x to express the ISA and it yields
tan θ i i n s = k x k r 2 k x 2 ,
where the radial wavenumber vector is parallel to the line from the radar pointing to the scatterer and the azimuth wavenumber vector is parallel to the radar velocity. Combining Equations (14) and (15), we can obtain
X r k x = R b y i k x k r 2 k x 2 + x i ,
The phase of the spectrum can be calculated by
φ S k x = X r k x d k x = R b y i k x k r 2 k x 2 x i d k x = k r 2 k x 2 R b y i k x x i ,
So, the spectrum derived by TFRSA is
S k r , k x = i = 1 K A i rect k r Δ k r exp j k r 2 k x 2 R b y i + k x x i ,
Equation (18) is the same as the result derived by the POSP [3], which means the TFRSA has the same accuracy as POSP for LTSAR, but it is easier to understand physically. In the following, the TFRTAR is applied to the situations where POSP cannot be used, which are near-field ISAR and CSAR applications.

2.4. Near-Field ISAR Spectrum Analysis Based on TFRSA and EPFA for Near-Field ISAR Imaging

Because the plane wave assumption is not valid in the near-field ISAR, the traditional ISAR imaging algorithm would cause distortion and defocusing for near-field ISAR imaging. To solve these problems, a new imaging algorithm should be designed based on the spectrum of the near-field ISAR. For simplicity, here we take the turntable near-field ISAR as an example to analyze its spectrum. The geometry of the near-field ISAR is demonstrated in Figure 3 and the slant range of i-th scatterer in the target can be expressed by
R i N F ϑ = R 0 x i sin ϑ y i cos ϑ 2 + x i cos y i sin ϑ 2 ,
where R 0 stands for the distance between the radar and the turntable center. From the geometry given in Figure 3, we can calculate the ISA of the i-th scatterer by
tan θ i i n s = x i cos ϑ y i sin ϑ R 0 x i sin ϑ y i cos ϑ ,
Because the relationship of Equation (15) still holds here, by combining Equations (15) and (20), we can obtain
x i cos ϑ y i sin ϑ R 0 x i sin ϑ y i cos ϑ = k x k r 2 k x 2 ,
Thus, the rotational angle ϑ can be expressed by the azimuth wavenumber k x as
ϑ N F k x = arccos R 0 k x k r x i 2 + y i 2 + arctan x i k x y i k r 2 k x 2 y i k x + x i k r 2 k x 2 ,
where arccos   and arctan   signify the arccosine and arctangent functions.
We perform FT on the rotation angle ϑ with the angular wavenumber ζ as a variable. Transforming the imaging geometry of near-field ISAR into the CSAR imaging geometry means that the radar is moving along a circle centered at the turntable. Then, the relationship between radar aperture position X r and rotational angle ϑ is X r = R 0 ϑ . According to the property of FT, the relationship between ζ and k x can be expressed by
k x = ζ / R 0     ζ = k x R 0 ,
Substituting Equation (23) into Equation (22), the rotational angle can be rewritten as
ϑ N F ζ = arccos ζ k r x i 2 + y i 2 + arctan x i ζ y i R 0 2 k r 2 ζ 2 y i ζ + x i R 0 2 k r 2 + ζ 2 ,
The spectrum phase of the near-field ISAR can be calculated by
φ N F S i ζ = ϑ N F ζ d ζ = arccos ζ k r x i 2 + y i 2 + arctan x i ζ y i R 0 2 k r 2 ζ 2 y i ζ + x i R 0 2 k r 2 ζ 2 d ζ = k r 2 x i 2 + y i 2 ζ 2 R 0 2 k r 2 ζ 2 ζ arccos ζ k r x i 2 + y i 2 ζ arctan x i ζ y i R 0 2 k r 2 ζ 2 y i ζ + x i R 0 2 k r 2 ζ 2 ,
where φ N F S i depicts the spectrum phase of the i-th scatterer in near-field ISAR situation when x i , y i 0 , 0 . When the scatterer is located at the turntable center, according to the imaging geometry, the range between radar and the scatterer is unchanged. Thus, the phase of the scatterer is constant and the spectrum of the scatterer is an impulse function. The spectrum of the i-th scatterer in the near-field can be expressed by
S N F k r , ζ = i = 1 K S N F i k r , ζ ,
with
S N F i k r , ζ = A i rect k r Δ k r δ ζ , x i , y i = 0 , 0 A i rect k r Δ k r rect ζ ζ N F i c Δ ζ N F i exp j k r 2 x i 2 + y i 2 ζ 2 R 0 2 k r 2 ζ 2 × exp j ζ arccos ζ k r x i 2 + y i 2 + arctan x i ζ y i R 0 2 k r 2 ζ 2 y i ζ + x i R 0 2 k r 2 ζ 2 ,   otherwise ,
ζ N F i = k r d R i ϑ d ϑ = k r R 0 x i cos ϑ y i sin ϑ R i N F ϑ ,
Δ ζ N F i = max ζ N F i min ζ N F i ,
ζ N F i c = k r R 0 x i R 0 y i 2 + x i 2 ,
where S N F i denotes the near-field spectrum of the i-th scatterer, δ   stands for the impulse function, ζ N F i signifies the angle wavenumber of the i-th scatterer, and Δ ζ N F i and ζ N F i c depict the bandwidth and the center angle wavenumber of the i-th scatterer in the near-field case. Observing Equation (27) reveals that different scatterers have different spectrums, namely, position-dependent, and the spectrum is more sensitive to the X-coordinate of the scatterer.
From electromagnetic theory, we know that a spherical wave can be decomposed into a summation of plane waves [29]. Thus, instead of trying to design an imaging algorithm directly from the near-field spectrum, we decided to design an algorithm from these plane wave components. According to the plane wave assumption, the far-field slant range of i-th scatterer can be expressed by
R i F F ϑ = R 0 x i sin ϑ y i cos ϑ ,
By applying the TFRSA to the far-field situation, we can obtain its TFDL by
ζ F F i = k r d R i F F ϑ d ϑ = k r x i cos ϑ y i sin ϑ ,
Rotating the TFDL by 90 degrees, the FTDL can be obtained by
ϑ F F ζ = arccos ζ k r x i 2 + y i 2 arctan y i x i ,
According to the TFRSA, the phase of the far-field spectrum can be calculated by
φ F F S i ζ = ϑ F F ζ d ζ = arccos ζ k r x i 2 + y i 2 + arctan y i x i d ζ = k r 2 x i 2 + y i 2 ζ 2 ζ arccos ζ k r x i 2 + y i 2 + ζ arctan y i x i ,
where φ F F S i depicts the spectrum phase of the i-th scatterer in far-field ISAR satiation. Like near-field ISAR, the spectrum of the i-th scatterer in the far-field can be expressed by
S F F k r , ζ = i = 1 K S F F i k r , ζ ,
with
S F F i k r , ζ = A i rect k r Δ k r δ ζ , x i , y i = 0 , 0 A i rect k r Δ k r rect ζ ζ F F i c Δ ζ F F i exp j k r 2 x i 2 + y i 2 ζ 2 × exp j ζ arccos ζ k r x i 2 + y i 2 + arctan x i y i ,   otherwise ,
Δ ζ F F i = max ζ F F i min ζ F F i ,
ζ F F i c = k r x i ,
where S F F i denotes the far-field spectrum of the i-th scatterer, ζ F F i signifies the angular wavenumber of the i-th scatterer, and Δ ζ F F i and ζ F F i depict the bandwidth and the center angle wavenumber of the i-th scatterer in the far-field case, respectively. Comparing Equations (27) and (36), the spectrums in the near-field and far-field situations are very similar except for a phase difference. This phase difference can be obtained by subtracting Equation (25) from Equation (34), and it yields
φ D i f S i ζ = φ F F S i ζ φ N F S i ζ = ζ arctan y i x i + R 0 2 k r 2 ζ 2 + ζ arctan x i ζ y i R 0 2 k r 2 ζ 2 y i ζ + x i R 0 2 k r 2 ζ 2 = R 0 2 k r 2 ζ 2 + ζ arcsin ζ k r R 0 ,
The detailed derivation of the above result is provided in Appendix A. Observing Equation (39) we can find that the phase difference is position-independent, which means we can transform the near-field spectrum into the far-field spectrum just by multiplying a compensation function in the spectrum domain. The constructed compensation function is given by
H N 2 F k r , ζ = exp j φ D i f S i ζ j k r R 0 = exp j R 0 2 k r 2 ζ 2 + ζ arcsin ζ k r R 0 k r R 0 ,
where the term k r R 0 is added to avoid changing the position of range profile.
After spectrum compensation, the near-field ISAR signal can then be modeled as an equivalent far-field ISAR signal. Transforming the compensated signal into the range wavenumber and azimuth angle domains via the azimuth inverse FT and it yields
s C k r , ϑ = i = 1 K A i rect k r Δ k r rect ϑ ϑ i c Δ ϑ i exp j k r x i sin ϑ + y i cos ϑ ,
with
ϑ i c = x i R 0 2 + x i 2 ,
Δ ϑ i = max ϑ i min ϑ i ,
where ϑ i c denotes the angle bias of the i-th scatterer caused by the compensation function H N 2 F and Δ ϑ i depicts the radar observation angle range of the i-th scatterer. Observing Equation (41) reveals that after compensating, the signals of different scatterers would shift along the azimuth dimension.
By applying the polar reformatting to Equation (41), namely, k α = k r sin ϑ ,   k β = k r cos ϑ , it can be rewritten as
s C k β , k α = i = 1 K A i rect k β k β i c Δ k β i rect k α k α i c Δ k α i exp j k α x i + k β y i ,
where k β and k α denote the range wavenumber and cross-range wavenumber, respectively; k β i c and k α i c signify the range and cross-range wavenumber centers of the i-th scatterer, respectively; and Δ k β i and Δ k α i depict the range and cross-range wavenumber bandwidths of the i-th scatterer, respectively. Finally, by performing two-dimensional FFT (2-D FFT) to Equation (44), the focused near-field ISAR image can be obtain by
I x , y = FFT 2 s k β , k α = i = 1 K A k sinc x x i 2 π / Δ k α i   sinc y y i 2 π / Δ k β i   ,
where FFT 2   depicts the 2-D FFT. Observing Equation (45) reveals that the resolutions of the image in range and cross-range dimensions are ρ r i = 2 π / Δ k β i and ρ a i = 2 π / Δ k α i , which means different scatterers have different resolutions. Besides, the image pixel size in cross-range and range dimensions are
p i x α = 2 π / Δ k α p i x β = 2 π / Δ k β ,
with
Δ k α = max k α i min k α i ,   for   all   i = 1 , 2 , , K Δ k β = max k β i min k β i ,     for   all   i = 1 , 2 , , K ,
where p i x α and p i x β are used to perform 2-D image scaling for the obtained image. This new imaging algorithm for near-field ISAR imaging is called extended polar format algorithm (EPFA), and its flowchart is given in Figure 4. The whole algorithm includes three steps: spectrum correction, polar reformatting, and 2-D image scaling.

2.5. CSAR Spectrum Analysis Based on TFRSA and EPFA-Based CSAR Imaging Method

In the CSAR imaging geometry, which is illustrated in Figure 5, the radar moves along a circular trajectory and continuously observes the same scene. This leads to a much large observation angle, which could produce ultrahigh resolution radar image of the scene. The CASR imaging geometry is equivalent to the geometry of the near-field turntable ISAR without consideration of the height of the trajectory. This means we can apply the EPFA to the CSAR imaging with some additional compensation processing. According to the geometry in Figure 5, the slant range of the i-th scatterer can be expressed by
R i C S A R θ , H = r 0 cos γ cos θ y i 2 + r 0 cos γ sin θ x i 2 + H 2 r 0 x i sin θ cos γ y i cos θ cos γ 2 + x i cos θ cos γ y i sin θ cos γ 2 + x i 2 + y i 2 sin 2 γ r 0 y i cos γ 2 + x i 2 cos 2 γ ,
where θ and γ denote the rotational angle of the radar along the circular orbit and the elevation angle of the orbit, respectively, which are illustrated in Figure 5. r 0 stands for the distance between the radar and the scene center and H signifies the height of the orbit, which can be expressed by H = r 0 sin γ . Observing the approximation result of Equation (48) and comparing it with Equation (19), the CSAR slant range has the same form as the near-field ISAR slant range without consideration of the constant error. Because the constant error does not influence the signal spectrum, we can still apply the EPFA to CSAR to obtain the equivalent far-field signal. However, the error would influence the polar reformatting result, which would cause distortion and defocusing to the final image.
Substituting Equation (48) into the signal model of Equation (3), the received CSAR signal after down conversion and range compression can be expressed as
s k r , θ = i = 1 K A i rect k r Δ k r exp j k r r 0 y i cos θ cos γ 2 + r 0 x i sin θ cos γ 2 × exp j k r x i 2 + y i 2 sin 2 γ r 0 y i cos γ 2 + x i 2 cos 2 γ ,
Performing EPFA to Equation (49), we can obtain
s k β , k α = i = 1 K A i rect k β k β i c Δ k β i rect k α k α i c Δ k α i exp j k α x i cos γ + k β y i cos γ × exp j k α 2 + k β 2 x i 2 + y i 2 sin 2 γ r 0 y i cos γ 2 + x i 2 cos 2 γ i = 1 K A i rect k β k β i c Δ k β i rect k α k α i c Δ k α i exp j k α x i cos γ + k β y i cos γ exp j φ e 1 + φ e 2 ,
with
φ e 1 = k β x i 2 sin 2 γ r 0 2 + x i 2 cos 2 γ ,
φ e 2 = k α 2 + k β c 2 y i 2 sin 2 γ r 0 y i cos γ ,
where φ e 1 and φ e 2 signify the decomposed range profile error and phase error, which would cause distortion and defocusing, respectively. k β c stands for the center of the range wavenumber k β .
In order to generate the focused image without distortion, we first compensate the range profile error. Transforming Equation (50) to range wavenumber and azimuth position domains by applying azimuth FT, we obtain
s k β , x = i = 1 K A i rect k β k β i c Δ k β i sinc x x i cos γ η 2 π / Δ k α i   exp j k β y i cos γ exp j φ e 1 ,
where η stands for the azimuth main-lobe expansion factor caused by the phase error φ e 2 . Observing the sinc function term of Equation (53) reveals that the azimuth positions of scatterers are scaled by the factor cos γ . We construct the range profile compensation functions in the range wavenumber and azimuth position domain by
H E C 1 k β , x n = exp j k β x n 2 sin 2 γ / cos 2 γ   r 0 2 + x n 2 ,
with
x n = 1 cos γ n N a 2 2 π Δ k α , n = 1 , 2 , , N a ,
where x n denotes the real X-coordinate of the n-th azimuth pixel, N a signifies the number of the pixels in azimuth dimension, and the term 2 π / Δ k α stands for the size of the image pixel in azimuth dimension.
After compensating the range profile error, the signal is transformed into the range position and azimuth wavenumber domains, and it yields
s y , k α = i = 1 K A i sinc y y i cos γ + y i 2 sin 2 γ / R 0 y i cos γ 2 π / Δ k β i rect k α k α i c Δ k α i × exp j k α x i cos γ exp j φ e 2 ,
From Equation (56), we can find that the focus positions of scatterers in range dimension are modulated by a complex function, which is given by
y ^ i = M y i = y i cos γ + y i 2 sin 2 γ / R 0 y i cos γ   ,
where y ^ i stands for the real focus position and M   signifies the focus position modulation function. We can solve the real Y-coordinates of image pixels by using the inverse function of M   , and it yields
y m = M 1 y ^ m ,
with
y ^ m = m N r 2 2 π Δ k β , m = 1 , 2 , , N r ,
where y ^ m denotes the focus position of the m-th range pixel and y m depicts the real Y-coordinate of the m-th range pixel. Thus, the phase error compensation function can be constructed by
H E C 2 k α , y m = exp j k α 2 + k β c 2 y m 2 sin 2 γ r 0 y m cos γ ,
After compensating the phase error, the signal can be expressed as
s y , k α = i = 1 K A i sinc y y i cos γ + y i 2 sin 2 γ / R 0 y i cos γ   2 π / Δ k β i   rect k α k α i c Δ k α i exp j k α x i cos γ ,
Finally, the focused CSAR imaging result can be obtained by applying the azimuth FT to Equation (61), and it yields
I x , y = i = 1 K A i sinc y y i cos γ + y i 2 sin 2 γ / R 0 y i cos γ   2 π / Δ k β i   sinc x x i cos γ 2 π / Δ k α i ,
Besides, the vectors x n and y m can be used to scale the image in azimuth and range dimensions, respectively. Finally, the flowchart of the EPFA-based CSAR imaging method is illustrated in Figure 6.

3. Results

To evaluate the performance of the proposed TFRSA and EPFA, we use both simulated and real data. Results achieved by the proposed EPFA are compared with those obtained by traditional PFA and BPA. Apart from comparing the visual images, the entropy, impulse response width (IRW), and peak side lobe ratio (PSLR) of the final focused images are measured to assess image qualities and the running times are measured to assess their efficiencies. Because the proposed method is not revolved in parameters estimation, the performance comparison between the proposed methods and other methods under various SNR conditions will not have a much difference.

3.1. Simulations

3.1.1. Simulation for Near-Field ISAR Imaging

The simulated radar parameters for near-field ISAR are listed in Table 1. The target point model of the airplane is shown in Figure 7a and the distance between radar and target is 10 m. The imaging results of different algorithms are presented in Figure 7b–d, which correspond to the PFA, BPA, and EPFA, respectively. Observing Figure 7b reveals that the imaging result of PFA is severely distorted and defocused when it is compared with the target point model in Figure 7a, whereas the BPA and EPFA in Figure 7c,d both show perfect airplane shape without obvious defocusing. It is worth mentioning that because the BPA is an accurate time-domain imaging algorithm, it is used as reference in this paper. In this simulation, the running times of BPA and PFA are 141.89 s and 11.29 s, respectively. For EPFA, we use the nonuniform FFT (NUFFT) [30] to speed up the polar reformatting operation and its running time is 2.22 s. It is clear that EPFA is highly efficient and is very suitable for real-time applications.
In order to have a more detailed comparison, the contour maps of three edge points, which are located at the head, right-wing, and left tail of the airplane, are plotted in Figure 8. Results from different algorithms are presented in the different columns, beginning with the first column for PFA, then the second column for BPA, and then the third column for EPFA. Different points are demonstrated in the different rows. It can be observed that, the contour maps of PFA are defocused severely, whereas the contour maps of EPFA are similar to those of BPA, that all of them have elliptical main-lobes and they are well separated from their side-lobes. The azimuth slices of the contour maps are also plotted in Figure 9, from which we can find the main-lobes of EPFA are narrower than those of PFA and similar to those of BPA. This signifies that EPFA has high accuracy in near-filed ISAR imaging.
In order to make a quantitative comparison, some parameters, which can evaluate the image quality, are calculated and listed in Table 2. The entropy [31,32,33] is a statistical measure of randomness that can be used to characterize the focus performance of the radar image and it is defined by
E n = i = 1 M N p i log p i ,
with
p i = I i 2 / i = 1 M N I i 2 ,
where I i and p i denote the value and the normalized power of the i-th element of the image. The smaller the entropy, the higher the image quality. The impulse response width (IRW) and the peak side lobe ratio (PLSR) of the image slice are used to measure the focusing performance. From Table 2 we can find that the entropies, IRWs and PLSRs of EPFA are smaller than PFA, whereas they are similar to those of BPA. Thus, the simulation result validate the effectiveness of the proposed EPFA.

3.1.2. Simulation for CSAR Imaging

The simulated radar parameters for CSAR are listed in Table 3. The scene of a 11 × 11 points array is shown in Figure 10a and the distance from radar to the center of the scene is 1100 m. The imaging results of different algorithms are shown in Figure 10b–d, which corresponding to the PFA, BPA, and the proposed method, respectively. Observing Figure 10b reveals that the imaging result of PFA is obviously distorted with a trapezoid shape, whereas the imaging results of BPA and the proposed method both show the standard square shape in Figure 10c,d. Besides, the running time for the proposed method is 8.78 s by using the NUFFT acceleration, whereas they are 3242.70 s and 81.96 s for BPA and PFA, respectively. We can find there is a large improvement in the efficiency using the proposed method.
The contour maps of four edge points are plotted to make a detailed comparison, which are shown in Figure 11. Different methods are presented in the different columns, beginning with the first column for PFA, then the second column for BPA, and then the third column for the proposed method. Different points are demonstrated in the different rows. Observing these contour maps, for PFA, they are defocused severely and their main lobes are connected with the side-lobes, whereas for the proposed method, they are similar to those of BPA, all of them have elliptical main lobes and separated side-lobes. The azimuth slices of the contour maps are also illustrated in Figure 12, which shows the proposed method has similar focus ability to the BPA.
In order to make a quantitative comparison, the IRWs and PLSRs of the azimuth slices are calculated and listed in Table 4. Besides, the image entropies are also calculated. From Table 4, we can find that the entropy, IRW, and PLSR of EPFA are similar to those of BPA while they are smaller than those of PFA. The simulation results validate the effectiveness of the proposed EPFA-based CSAR imaging method.

3.2. Real Data Processing Results

Real near-field ISAR data are collected in a microwave anechoic chamber, in which an aircraft model target is placed on the turntable model as shown in Figure 13a. The radar is fixed, and the height of radar is set a little bit higher than the turntable surface. The horizontal distance between the radar and the turntable surface is 7 m. The rotational angle of turntable can be precisely controlled, and it is changed discretely with step of 0.2 degree. The radar transmits and receives signal at every discrete rotational angle. The parameters of the radar are listed in the Table 5.
The imaging results based on the real data processed by different algorithms are illustrated in Figure 13b–d, which corresponding to the PFA, BPA, and EPFA, respectively. Observing Figure 13b, we can find that the head of the aircraft, which is marked with a rectangular, shows an expanded main-lobe for PFA, whereas in Figure 13c,d, their aircraft head parts show more focused main lobes. The contour maps of the aircraft’s head are plotted in Figure 14, in which the main lobe of PFA is of an oval shape whereas they are round for BPA and EPFA. The azimuth slices of the contour maps are plotted in Figure 15 for further comparison, in which the main lobe of the PFA is obviously wider than that of BPA and PFA. The performance parameters of image entropy, IRW, PLSR, and running time are calculated and listed in Table 6 for a quantitative comparison. We can find that the results of EPFA are similar to the results of BPA, while they are smaller than those of PFA. In addition, the running time for EFPA is much smaller than BPA. The real data processing result validate the effectiveness of the proposed EPFA.

4. Discussion

PFA is one of the common imaging algorithms for SAR and ISAR because of its simplicity and efficiency. However, due to the plane wave assumption of PFA, it would produce large slant range error for near-field ISAR and wide scene SAR, which would generate range profiles migration and azimuth phase error, and thus cause the image to be distorted and defocused. Because the BPA constructs user-defined grids, the accurate slant range of every grid can be calculated by the known trajectory information. Then, radar returns can be accumulated coherently for every grid by using this accurate slant range. BPA is an accurate imaging algorithm for all modes of SAR and ISAR, so it would produce focused image without distortion. The drawback of BPA is its large computation cost. On the other hand, according to the introduction of the proposed TFRSA in the Section 2.2, the influence of the envelopes of the scatterer’s returns can be ignored. This approximation is applicable and often used in radar signal processing because the envelope is usually a low-frequency and narrow bandwidth signal. Thus, the proposed TFRSA-based EPFA has high accuracy in radar imaging and generates comparably high image quality to BPA.
For the future research, we will focus on the possibility of applying the TFRSA to more applications and consider the influence of the motion error in the realistic environments. Now, we are trying to apply the proposed method to the helical SAR(HSAR)2D imaging, HSAR and tomographic CSAR 3D imaging, and the primary results is prove that the proposed method is valid to these applications.

5. Conclusions

In this paper, a time-frequency reversion-based spectrum analysis (TFRSA) method is proposed, which utilizes the relationship of the Fourier pairs and their corresponding signal phase. Based on this method, the spectrum of the linear trajectory SAR is analyzed. The derived spectrum expression is shown to be the same as the result of the POSP. In addition, the TFRSA is also applied to the near-field ISAR, which the POSP cannot be applied to. The spectrum expression of near-field ISAR is derived and it is compared with the far-field ISAR spectrum expression. The spectrum difference is found to be position-independent. Based on this finding, an extended PFA is proposed to solve the distortion and defocusing problems caused by traditional ISAR imaging algorithms. Because the CSAR has similar imaging geometry to the near-field ISAR, its spectrum is also derived by TFRSA and an EPFA-based imaging method is proposed to improve the imaging efficiency. Simulations and real data processing results validate the effectiveness of the proposed TFRSA, EPFA, and EPFA-based CSAR imaging method.

Author Contributions

Conceptualization, J.F., G.S., and M.X.; methodology, J.F., G.S., and M.X.; software, J.F. and G.S.; validation, J.F., G.S., and M.X.; formal analysis, J.F., G.S., and M.X.; investigation, M.X.; resources, J.F. and G.S.; writing—original draft preparation, G.S. and M.X.; writing—review and editing, J.F. and G.S.; visualization, J.F. and G.S.; supervision, M.X.; project administration, M.X.; funding acquisition, M.X. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported in part by the National Science Fund for Distinguished Young Scholars (Grant number 61825105) and in part by the Foundation for Innovative Research Groups of the National Natural Science Foundation of China (61621005).

Acknowledgments

The authors would like to acknowledge the support of the China Scholarship Council. The authors would also like to acknowledge the help of Tat Soon Yeo for his reviewing and improving the paper. Fu was a visiting scholar at the Center for Advanced Communications (CAC), Villanova University, in 2019–2020 and would like to acknowledge the support of CAC and the director Moeness G. Amin.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
BPABack-projection Algorithm
CSAChirp Scaling Algorithm
CSARCircular SAR
EPFAExtended Polar Format Algorithm
FSAFrequency Scaling Algorithm
FTDLFrequency-time Distribution Line
IRWImpulse Response Width
ISAInstantaneous Squint Angle
ISARInverse Synthetic Aperture Radar
KTkeystone Transform
LFMlinear frequency modulation
LOSLight of Sight
LTSARLinear Trajectory SAR
PFAPolar Format Algorithm
PLSRPeak Side Lobe Ratio
POSPPrinciple of Stationary Phase
RDRange-Doppler
RMRotational Motion
RMARange Migration Algorithm
SASpectrum Analysis
SARSynthetic Aperture Radar
TFDLTime-frequency Distribution Line
TFRTime-frequency reversion
TFRSATFR-based Spectrum Analysis
TMTranslational Motion

Appendix A

In order to simplify Equation (39), we apply the angle sum tan trigonometric identity to the two a r c t a n   terms, which is
tan A + B = tan A + tan B 1 tan A tan B ,
Performing the t a n   operation to the two a r c t a n   terms and applying Equation (A1) yields
tan arctan y i x i + arctan x i ζ y i R 0 2 k r 2 ζ 2 y i ζ + x i R 0 2 k r 2 ζ 2 = y i x i + x i ζ y i R 0 2 k r 2 ζ 2 y i ζ + x i R 0 2 k r 2 ζ 2 / 1 y i x i x i ζ y i R 0 2 k r 2 ζ 2 y i ζ + x i R 0 2 k r 2 ζ 2   = y i y i ζ + x i R 0 2 k r 2 ζ 2 + x i x i ζ y i R 0 2 k r 2 ζ 2 x i y i ζ + x i R 0 2 k r 2 ζ 2 y i x i ζ y i R 0 2 k r 2 ζ 2 = x i 2 + y i 2 ζ x i 2 + y i 2 R 0 2 k r 2 ζ 2 = ζ R 0 2 k r 2 ζ 2 ,
Applying the a r c t a n   operation to two sides of Equation (A2), we can obtain
arctan y i x i arctan x i ζ + y i R 0 2 k r 2 ζ 2 y i ζ + x i R 0 2 k r 2 ζ 2 = arctan ζ R 0 2 k r 2 ζ 2 = arcsin ζ R 0 k r ,

References

  1. Wehner, D.R. High Resolution Radar; Artech House: Norwood, MA, USA, 1987. [Google Scholar]
  2. Mensa, D.L. High Resolution Radar Imaging; Artech House: Norwood, MA, USA, 1981. [Google Scholar]
  3. Cumming, I.; Wong, F. Digital Processing of Synthetic Aperture Radar Data: Algorithms and Implementation; Artech House: Norwood, MA, USA, 2005. [Google Scholar]
  4. Chen, C.C.; Andrews, H.C. Target-motion-induced radar imaging. IEEE Trans. Aerosp. Electron. Syst. 1980, 2–14. [Google Scholar] [CrossRef]
  5. Chen, V.C. Inverse Synthetic Aperture Radar Imaging: Principles, Algorithms and Applications; Institution of Engineering and Technology: London, UK, 2014. [Google Scholar]
  6. Chen, V.C.; Ling, H. Time-Frequency Transforms for Radar Imaging and Signal Analysis; Artech House: Norwood, MA, USA, 2002. [Google Scholar]
  7. Yao, J. Microwave photonics. J. Light. Technol. 2009, 27, 314–335. [Google Scholar] [CrossRef]
  8. McKinney, J.D. Photonics illuminates the future of radar. Nature 2014, 507, 310–312. [Google Scholar] [CrossRef] [PubMed]
  9. Capmany, J.; Novak, D. Microwave photonics combines two worlds. Nat. Photonics 2007, 1, 319–330. [Google Scholar] [CrossRef]
  10. Gao, Y.; Xing, M.; Zhang, Z.; Guo, L. Ultrahigh Range Resolution ISAR Processing by Using KT–TCS Algorithm. IEEE Sens. J. 2018, 18, 8311–8317. [Google Scholar] [CrossRef]
  11. Yang, L.; Xing, M.; Zhang, L.; Sun, G.; Gao, Y.; Zhang, Z.; Bao, Z. Integration of Rotation Estimation and High-Order Compensation for Ultrahigh-Resolution Microwave Photonic ISAR Imagery. IEEE Trans. Geosci. Remote. Sens. 2020, 1–21. [Google Scholar] [CrossRef]
  12. Bleistein, N.; Handelsman, R.A. Asymptotic Expansions of Integrals; Courier Corporation: North Chelmsford, MA, USA, 1986. [Google Scholar]
  13. Li, J.; Zhang, S.; Chang, J. Bistatic forward-looking SAR imaging based on two-dimensional principle of stationary phase. In Proceedings of the 2012 International Workshop on Microwave and Millimeter Wave Circuits and System Technology, Chengdu, China, 19–20 April 2012; pp. 1–4. [Google Scholar]
  14. Andre, S.; Ruben, R. Analysis of scattering arrays using the principle of stationary phase. In Proceedings of the 1964 Antennas and Propagation Society International Symposium, Long Island, NY, USA, 21–24 September 1964; Volume 2, pp. 19–23. [Google Scholar]
  15. Chassande-Mottin, E.; Flandrin, P. On the stationary phase approximation of chirp spectra. In Proceedings of the IEEE-SP International Symposium on Time-Frequency and Time-Scale Analysis (Cat. No. 98TH8380), Pittsburgh, PA, USA, 9 October 1998; pp. 117–120. [Google Scholar]
  16. Yang, Z.; Xing, M.; Zhang, Y.D.; Sun, G.; Bao, Z. Factorised polar-format back-projection algorithm. IET Radar Sonar Navig. 2015, 9, 875–880. [Google Scholar] [CrossRef]
  17. Li, H.; Zhang, L.; Xing, M.; Bao, Z. Expediting back-projection algorithm for circular SAR imaging. Electron. Lett. 2015, 51, 785–787. [Google Scholar] [CrossRef]
  18. Mengdao, X.; Renbiao, W.; Jinqiao, L.; Zheng, B. Migration through resolution cell compensation in ISAR imaging. IEEE Geosci. Remote. Sens. Lett. 2004, 1, 141–144. [Google Scholar] [CrossRef]
  19. Zhang, L.; Sheng, J.-l.; Duan, J.; Xing, M.-d.; Qiao, Z.-j.; Bao, Z. Translational motion compensation for ISAR imaging under low SNR by minimum entropy. EURASIP J. Adv. Signal Process. 2013, 2013, 33. [Google Scholar] [CrossRef] [Green Version]
  20. Li, D.; Zhan, M.; Liu, H.; Liao, Y.; Liao, G. A robust translational motion compensation method for ISAR imaging based on keystone transform and fractional Fourier transform under low SNR environment. IEEE Trans. Aerosp. Electron. Syst. 2017, 53, 2140–2156. [Google Scholar] [CrossRef]
  21. Ji, Z.; Zhu, X. Application of DBF in 77 GHz Automotive Millimeter-wave Radar. In IOP Conference Series: Materials Science and Engineering; IOP Publishing: Bristol, UK, 2019; Volume 490, p. 072066. [Google Scholar] [CrossRef]
  22. Keziah, J.; Muthukumaran, N. Design of K Band Transmitting Antenna for Harbor Surveillance Radar Application. Int. J. Appl. Electr. Electron. Eng. 2016, 2, 16–20. [Google Scholar]
  23. Turin, F.; Pastina, D.; Lombardo, P.; Corucci, L. ISAR imaging of ground targets with an X-band FMCW Radar system for airport surveillance. In Proceedings of the IEEE 2014 15th International Radar Symposium (IRS), Gdansk, Poland, 16–18 June 2014; pp. 1–5. [Google Scholar] [CrossRef]
  24. Amin, M.G. Spectral decomposition of time-frequency distribution kernels. IEEE Trans. Signal Process. 1994, 42, 1156–1165. [Google Scholar] [CrossRef]
  25. Cohen, L. Time-frequency distributions-a review. Proc. IEEE 1989, 77, 941–981. [Google Scholar] [CrossRef] [Green Version]
  26. Xu, G.; Xing, M.; Zhang, L.; Duan, J.; Chen, Q.; Bao, Z. Sparse Apertures ISAR Imaging and Scaling for Maneuvering Targets. IEEE J. Sel. Top. Appl. Earth Obs. Remote. Sens. 2014, 7, 2942–2956. [Google Scholar] [CrossRef]
  27. Bailey, D.H.; Swarztrauber, P.N. The fractional Fourier transform and applications. SIAM Rev. 1991, 33, 389–404. [Google Scholar] [CrossRef]
  28. Ozaktas, H.M.; Kutay, M.A. The fractional fourier transform. In Proceedings of the 2001 European Control Conference (ECC), Porto, Portugal, 4–7 September 2001; pp. 1477–1483. [Google Scholar]
  29. Rafaely, B. Plane-wave decomposition of the sound field on a sphere by spherical convolution. J. Acoust. Soc. Am. 2004, 116, 2149–2157. [Google Scholar] [CrossRef] [Green Version]
  30. Greengard, L.; Lee, J.Y. Accelerating the nonuniform fast Fourier transform. SIAM Rev. 2004, 46, 443–454. [Google Scholar] [CrossRef] [Green Version]
  31. Flores, B.C.; Ugarte, A.; Kreinovich, V. Choice of an entropy-like function for range-Doppler processing. In Proceedings of the Automatic Object Recognition III, International Society for Optics and Photonics, Orlando, FL, USA, 11–16 April 1993; Volume 1960, pp. 47–56. [Google Scholar]
  32. Xi, L.; Guosui, L.; Ni, J. Autofocusing of ISAR images based on entropy minimization. IEEE Trans. Aerosp. Electron. Syst. 1999, 35, 1240–1252. [Google Scholar] [CrossRef]
  33. Kragh, T.J.; Kharbouch, A.A. Monotonic iterative algorithm for minimum-entropy autofocus. In Proceedings of the Adaptive Sensor Array Processing (ASAP) Workshop, Lexington, MA, USA, 6–7 June 2006; Volume 40, pp. 1147–1159. [Google Scholar]
Figure 1. Radar imaging geometry. (a) Synthetic aperture radar (SAR) imaging geometry. (b) Inverse synthetic aperture radar (ISAR) imaging geometry.
Figure 1. Radar imaging geometry. (a) Synthetic aperture radar (SAR) imaging geometry. (b) Inverse synthetic aperture radar (ISAR) imaging geometry.
Remotesensing 13 00600 g001
Figure 2. Imaging geometry of linear trajectory SAR.
Figure 2. Imaging geometry of linear trajectory SAR.
Remotesensing 13 00600 g002
Figure 3. Imaging geometry of near-field ISAR.
Figure 3. Imaging geometry of near-field ISAR.
Remotesensing 13 00600 g003
Figure 4. The flowchart of extended polar format algorithm (EPFA) for near-field ISAR imaging.
Figure 4. The flowchart of extended polar format algorithm (EPFA) for near-field ISAR imaging.
Remotesensing 13 00600 g004
Figure 5. Imaging geometry of CSAR.
Figure 5. Imaging geometry of CSAR.
Remotesensing 13 00600 g005
Figure 6. Flowchart of EPFA-based CSAR imaging method.
Figure 6. Flowchart of EPFA-based CSAR imaging method.
Remotesensing 13 00600 g006
Figure 7. Target point model and imaging results. (a) Airplane point model. (b) ISAR Imaging result of polar format algorithm (PFA). (c) ISAR Imaging result of back-projection algorithm (BPA). (d) ISAR Imaging result of EPFA.
Figure 7. Target point model and imaging results. (a) Airplane point model. (b) ISAR Imaging result of polar format algorithm (PFA). (c) ISAR Imaging result of back-projection algorithm (BPA). (d) ISAR Imaging result of EPFA.
Remotesensing 13 00600 g007
Figure 8. Contour maps of the selected points in the imaging results. (ac) Contour map of point 1 in the PFA result. (df) Contour map of point 1 in the BPA result. (gi) Contour map of point 1 in the EPFA result.
Figure 8. Contour maps of the selected points in the imaging results. (ac) Contour map of point 1 in the PFA result. (df) Contour map of point 1 in the BPA result. (gi) Contour map of point 1 in the EPFA result.
Remotesensing 13 00600 g008
Figure 9. Azimuth slices of the selected points by different imaging algorithms. (a) Azimuth slices of point 1. (b) Azimuth slices of point 2. (c) Azimuth slices of point 3.
Figure 9. Azimuth slices of the selected points by different imaging algorithms. (a) Azimuth slices of point 1. (b) Azimuth slices of point 2. (c) Azimuth slices of point 3.
Remotesensing 13 00600 g009
Figure 10. Target point model and imaging results. (a) Point array Scene. (b) CSAR Imaging result of PFA. (c) CSAR Imaging result of BPA. (d) CSAR Imaging result of proposed method.
Figure 10. Target point model and imaging results. (a) Point array Scene. (b) CSAR Imaging result of PFA. (c) CSAR Imaging result of BPA. (d) CSAR Imaging result of proposed method.
Remotesensing 13 00600 g010
Figure 11. Contour maps of the selected points. (ac) Contour map of point 1 by different imaging algorithms. (df) Contour map of point 2 by different imaging algorithms. (gi) Contour map of point 3 by different imaging algorithms. (jl) Contour map of point 4 by different imaging algorithms.
Figure 11. Contour maps of the selected points. (ac) Contour map of point 1 by different imaging algorithms. (df) Contour map of point 2 by different imaging algorithms. (gi) Contour map of point 3 by different imaging algorithms. (jl) Contour map of point 4 by different imaging algorithms.
Remotesensing 13 00600 g011
Figure 12. Azimuth slices of the selected points by different imaging algorithms. (a) Azimuth slices of point 1. (b) Azimuth slices of point 2. (c) Azimuth slices of point 3. (d) Azimuth slices of point 4.
Figure 12. Azimuth slices of the selected points by different imaging algorithms. (a) Azimuth slices of point 1. (b) Azimuth slices of point 2. (c) Azimuth slices of point 3. (d) Azimuth slices of point 4.
Remotesensing 13 00600 g012
Figure 13. Target point model and imaging results. (a) Real experiment scene. (b) Real data ISAR Imaging result of PFA. (c) Real data ISAR Imaging result of BPA. (d) Real data ISAR Imaging result of EPFA.
Figure 13. Target point model and imaging results. (a) Real experiment scene. (b) Real data ISAR Imaging result of PFA. (c) Real data ISAR Imaging result of BPA. (d) Real data ISAR Imaging result of EPFA.
Remotesensing 13 00600 g013
Figure 14. Contour maps of the aircraft’s head. (a) Contour map of the aircraft’s head by PFA. (b) Contour map of the aircraft’s head by BPA. (c) Contour map of the aircraft’s head by EPFA.
Figure 14. Contour maps of the aircraft’s head. (a) Contour map of the aircraft’s head by PFA. (b) Contour map of the aircraft’s head by BPA. (c) Contour map of the aircraft’s head by EPFA.
Remotesensing 13 00600 g014
Figure 15. Azimuth slices of the aircraft’s head by different algorithms.
Figure 15. Azimuth slices of the aircraft’s head by different algorithms.
Remotesensing 13 00600 g015
Table 1. Near-field ISAR imaging simulation parameters.
Table 1. Near-field ISAR imaging simulation parameters.
Near-Field ISAR Simulation ParametersValuesNear-Field ISAR Simulation ParametersValues
Bandwidth B1 GHzPulse Width T p 1 μs
Sampling Frequency F s 1.1 GHzCarrier Frequency f c 10 GHz
Pulse Repeat Frequency f P R F 500 HzRotational Angle Δ θ 47 °
Table 2. Image performance parameters comparison of near-field ISAR simulation.
Table 2. Image performance parameters comparison of near-field ISAR simulation.
AlgorithmsEntropyIRW (m)PLSR (dB)Time (s)
P1P2P3P1P2P3
PFA11.29Null0.190.058 0.54 2.97 3.61 11.29
BPA9.290.0270.0450.013 13.05 14.58 7.68 141.88
EPFA9.410.0270.0470.013 12.94 14.51 10.72 2.22
Table 3. CSAR simulation parameters.
Table 3. CSAR simulation parameters.
CSAR Simulation ParametersValuesCSAR Simulation ParametersValues
Bandwidth B500 MHzPulse Width T p 600 μs
Sampling Frequency F s 600 MHzCarrier Frequency f c 10 GHz
Pulse Repeat Frequency f P R F 8000 HzRotational Angle Δ θ 14.67 °
Reference Range r 0 1100 mTrajectory Height H214.60 m
Table 4. Image performance parameters comparison of CSAR simulation.
Table 4. Image performance parameters comparison of CSAR simulation.
MethodsEntropyIRW (m)PLSR (dB)Time (s)
P1P2P3P4P1P2P3P4
PFA9.720.110.400.350.30 4.93 0.61 0.23 0.52 81.96
BPA8.210.0460.0460.0590.059 13.21 13.16 13.27 13.21 3242.70
Proposed8.220.0460.0460.0590.059 12.17 13.23 13.27 12.59 3.24
Table 5. Near-field ISAR radar parameters.
Table 5. Near-field ISAR radar parameters.
Radar ParametersValuesRadar ParametersValues
Bandwidth B4 GHzCarrier Frequency f c 10 GHz
Reference Range R 0 7 mRotational Angle Δ θ 40 °
Table 6. Image performance parameters comparison of near-field ISAR real data processing.
Table 6. Image performance parameters comparison of near-field ISAR real data processing.
MethodsEntropyIRW (m)PLSR (dB)Time (s)
PFA9.750.040 12.25 0.45
BPA9.700.019 12.63 1.12
EPFA9.690.019 13.02 0.04
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Fu, J.; Xing, M.; Sun, G. Time-Frequency Reversion-Based Spectrum Analysis Method and Its Applications in Radar Imaging. Remote Sens. 2021, 13, 600. https://0-doi-org.brum.beds.ac.uk/10.3390/rs13040600

AMA Style

Fu J, Xing M, Sun G. Time-Frequency Reversion-Based Spectrum Analysis Method and Its Applications in Radar Imaging. Remote Sensing. 2021; 13(4):600. https://0-doi-org.brum.beds.ac.uk/10.3390/rs13040600

Chicago/Turabian Style

Fu, Jixiang, Mengdao Xing, and Guangcai Sun. 2021. "Time-Frequency Reversion-Based Spectrum Analysis Method and Its Applications in Radar Imaging" Remote Sensing 13, no. 4: 600. https://0-doi-org.brum.beds.ac.uk/10.3390/rs13040600

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