Next Article in Journal
A Novel Method for Designing General Window Functions with Flexible Spectral Characteristics
Previous Article in Journal
Unobtrusive Vital Sign Monitoring in Automotive Environments—A Review
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Bistatic ISAR Imaging with a V-FM Waveform Based on a Dual-Channel-Coupled 2D-CS Algorithm

1
State Key Laboratory of Complex Electromagnetic Environment Effects on Electronics and Information System, National University of Defense Technology, Changsha 410073, China
2
Naval Research Academy, PLA (NVRA), Beijing 100000, China
*
Author to whom correspondence should be addressed.
Sensors 2018, 18(9), 3082; https://doi.org/10.3390/s18093082
Submission received: 27 June 2018 / Revised: 21 August 2018 / Accepted: 5 September 2018 / Published: 13 September 2018
(This article belongs to the Section Remote Sensors)

Abstract

:
Due to the sparsity of the space distribution of point scatterers and radar echo data, the theory of Compressed Sensing (CS) has been successfully applied in Inverse Synthetic Aperture Radar (ISAR) imaging, which can recover an unknown sparse signal from a limited number of measurements by solving a sparsity-constrained optimization problem. In this paper, since the V style modulation(V-FM) signal can mitigate the ambiguity apparent in range and velocity, the dual-channel, two-dimension, compressed-sensing (2D-CS) algorithm is proposed for Bistatic ISAR (Bi-ISAR) imaging, which directly deals with the 2D signal model for image reconstruction based on solving a nonconvex optimization problem. The coupled 2D super-resolution model of the target’s echoes is firstly established; then, the 2D-SL0 algorithm is applied in each channel with different dictionaries, and the final image is obtained by synthesizing the two channels. Experiments are used to test the robustness of the Bi-ISAR imaging framework with the two-dimensional CS method. The results show that the framework is capable accurately reconstructing the Bi-ISAR image within the conditions of low SNR and low measured data.

1. Introduction

The inverse synthetic aperture radar (ISAR) has many important detecting and recognizing applications in many civilian and military fields such as air traffic control, harbor and river traffic surveillance, and the remote sensing of satellites under the conditions all-weather and all-day [1,2,3]. Compared to monostatic ISAR, Bistatic ISAR (Bi-ISAR) has many advantages, such as platform flexibility, and immunity to electronic countermeasures and physical attack. The configurations in which the transmitter and the receiver are installed on separate platforms can offer more freedom for acquiring complementary target information and avoiding blind velocities [4,5,6]. In conventional Range-Doppler imaging, by means of intra-pulse frequency modulation for radar waveforms, the range resolution and the detection range can be improved via the pulse compression technique and large time-bandwidth product waveforms [7]. These large time and band product signals, which include chirp signal, frequency-stepped signal, and frequency-stepped chirp signal, are commonly utilized in modern radar systems [8,9,10,11]. However, the drawback of the chirp signal is that the ambiguity apparent in range and velocity, which can be seen from the “ridge” ambiguity function, and frequency-stepped waveform need a relatively long coherent processing interval (CPI) to complete the transmitting and receiving process of each burst, and the wide aspect angle with long CPI may lead to range cell migration (RCM). These disadvantages greatly degrade the ISAR imaging quality. Thus, ISAR imaging based on the V-style modulation (V-FM) waveform has been proposed because of the “thumbtack” ambiguity function, which can mitigate the ambiguity apparent in range and velocity [12].
Compressed Sensing (CS) theory states that sparse or compressible signals are able to be reconstructed exactly from limited measurements by solving a sparsity-constrained optimization problem with a high probability, which has been used in many fields such as medical imaging and radar [13]. CS has been adopted to achieve high resolution and well-focused ISAR images using sparse sampling data [14,15,16]. ISAR imaging with wideband V-FM waveforms has been proposed in [17]. Gu and Pan et al. in [12] adopted dual-channel CS-based dechirping (CS-D) to achieve high resolution range profiles (HRRPs) of V-FM waveforms. In the framework, the HRRPs reconstruct themselves by solving the a l 1 norm convex optimization problem and using Fourier transform (FT) in azimuth direction to obtain the final 2D images. An improved CS-based high azimuth imaging method is developed to overcome strong noise and clutter in [18]. However, the above methods are the sparsity-based high range or azimuth resolution imaging methods, which are very time consuming because of the high computational complexity and storage system requirements. To avoid the above disadvantages, two-dimension (2D) coupled imaging methods have been proposed [19,20,21,22]. The traditional 2D Fourier Transform method (2D-FFT) needs long data collection, which in reality, is unachievable. In [19], sensing matrices are designed in the cross range and range direction to implement sparse sampling in both directions and 2D-SL0 algorithm is used to achieve high-resolution, fully polarimetric ISAR images. A novel two-dimensional (2D) group primal dual active set with the continuation (2DGPDASC) algorithm is used to recover an ISAR image in [10]. In [22], the sparsity-based imaging algorithms used to achieve range and azimuth resolutions simultaneously are proposed, but the ISAR imaging matrix vectoring to the one-dimensional (1D) model results in extensive memory usage.
In summary, existing research on ISAR imaging with the V-FM signal can achieve high-resolution images based on 1D sparsity-driven algorithms, but there is little research on V-FM signal Bi-ISAR 2D-CS imaging. In this framework, we establish Bi-ISAR echoes as a 2D-coupled, super-resolution model in the range and azimuth domains by analyzing the characteristics of echoes. Then, we adopt a dual-channel, two dimensional, smoothed L0 (2D-SL0) model to achieve Bi-ISAR imaging, which belongs to a nonconvex optimization problem. Using the steepest descent method and gradient projection principle, the nonconvex approach is advantageous in that it usually yields sparser solutions for a given size of samples, as well as short reconstruction time and low calculation load.
The remainder of this paper is organized as follows. In Section 2, the signal model and Bi-ISAR echoes are introduced. The CS theory is presented in Section 3, the 2D-coupled SR models of V-FM echoes in the dual channel are constructed, and the improved dual-channel 2D-SL0 algorithm for Bi-ISAR imaging is proposed. Section 4 presents the experimental results to validate the effectiveness of the proposed method, and, finally, conclusions are drawn and presented in Section 5.

2. Signal Model and Imaging

The relative motion between radar and target in ISAR imaging consists of translation motion and rotation motion. In this paper, we suppose that the target’s translation motion, which is defined as the movement of the target along the range axis of the radar, has been well completed previously [2]. The target has been converted into a turntable model. The geometry of bistatic ISAR imaging is illustrated in Figure 1, the local coordinate system. x o y is on the moving target, and the point o of the moving target is chosen as the origin. The moving target has the circular motion with a constant rotation rate of w rad/s; r i represents the distance from point i to reference point o , and R R and R T are the distances from point o to the receiver and transmitter, respectively. The dash line represents the initial position of the target, the solid line represents the position after a rotation angle Δ θ i , and β is the bistatic angle, which is kept constant during short CPI; the instantaneous distance sum from the scattering point i to the transmitter and receiver is given as (see Appendix A.1)
R R T + R R + 2 cos β 2 [ y i cos Δ θ i + x i sin Δ θ i ]
Generally, the observation time of Bi-ISAR imaging is very short, and the rotation angle can be approximated by the first-order series
R R T + R R + 2 cos β 2 [ y i + x i w t m ]
Assume that the V-FM waveforms are composed of two chirp signals with opposite slopes as depicted in Figure 2, which follows
υ ( t ^ ) = rect h a l f ( t ^ T p ) υ 1 ( t ^ ) + rect h a l f + ( t ^ T p ) υ 2 ( t ^ )
in which
rect h a l f = { 1 t ^ [ T p 2 , 0 ] 0 e l s e
rect h a l f + = { 1 t ^ [ 0 , T p 2 ] 0 e l s e
υ 1 ( t ^ ) = { exp ( j π γ t ^ 2 ) t ^ [ , 0 ] 0 e l s e
υ 2 ( t ^ ) = { exp ( j π γ t ^ 2 ) t ^ [ 0 , ] 0 e l s e
and γ is the chirp rate, T p is the pulse duration and bandwidth B = γ T P , and t ^ represents the fast time. Figure 2 shows the time-frequency of V-FM waveforms.
The transmitted V-FM waveform can be expressed as follows:
s ( t ^ ) = υ ( t ^ ) exp ( j 2 π f 0 t ) = rect h a l f ( t ^ T p ) exp ( j 2 π ( f 0 t 1 2 γ t ^ 2 ) ) + rect h a l f + ( t ^ T p ) exp ( j 2 π ( f 0 t + 1 2 γ t ^ 2 ) )
in which f 0 is the frequency of carrier wave, t m is the slow time, t is the full time, and t ^ = t t m . If the pulse repetition interval (PRI) of radar is T P R I , then t m = m T P R I .
To compare the difference between linear frequency modulation (LFM) signal and V-FM signal, the ambiguity function is calculated and simulated, respectively. The definition of the ambiguity function of μ ( t ) is as follows:
χ ( τ , ξ ) = μ ( t ) μ ( t + τ ) e j 2 τ ξ t d t
in which τ is the delay, which represents range, and ξ is the Doppler frequency, which represents velocity. The results of both signals are shown in Figure 3; the horizontal axis represents range and velocity, and the vertical axis represents the normalized value of | χ ( τ , ξ ) | , which shows that the “thumbtack” ambiguity function of the V-FM signal mitigates the ambiguity apparent in range and velocity compared with the “ridge” ambiguity function of the LFM signal.
Supposing that the target consists of I strong scatterers, the echo of the scattering center i can be written as
s r ( t ^ , t m ) = i = 1 K σ i s ( t ^ R c , t m ) = i = 1 K σ i rect h a l f ( t ^ R c T p ) exp ( j 2 π ( f 0 ( t R c ) 1 2 γ ( t ^ R c ) 2 ) ) + i = 1 K σ i rect h a l f + ( t ^ R c T p ) exp ( j 2 π ( f 0 ( t R c ) + 1 2 γ ( t ^ R c ) 2 ) )
in which σ i is the scattering coefficient of point i , R represents the distance sum from the i-th scattering center to the transmitter and the receiver, and c is the speed of the electromagnetic wave.
The reference signals of dual-channel dechirping can be expressed as
s r e f 1 ( t ^ , t m ) = rect h a l f ( t ^ R r e f c T r e f ) × exp ( j 2 π ( f 0 ( t T r e f ) 1 2 γ ( t ^ T r e f ) 2 ) )
s r e f 2 ( t ^ , t m ) = rect h a l f + ( t ^ R r e f c T r e f ) × exp ( j 2 π ( f 0 ( t T r e f ) + 1 2 γ ( t ^ T r e f ) 2 ) )
in which R r e f is the reference range R r e f = R R + R T , T r e f is the reference pulse duration, and T r e f = R r e f c ; then, the echoed signal in dual channels after dechirping can be expressed as follows:
s i f 1 ( t ^ , t m ) = i = 1 K σ i rect h a l f ( t ^ R c T p ) × exp ( j 2 π γ c ( t ^ R r e f c ) R Δ ) × exp ( j π γ c 2 R Δ 2 ) × exp ( j 2 π f 0 R Δ c )
s i f 2 ( t ^ , t m ) = i = 1 K σ i rect h a l f + ( t ^ R c T p ) × exp ( j 2 π γ c ( R r e f c t ^ ) R Δ ) × exp ( j π γ c 2 R Δ 2 ) exp ( j 2 π f 0 R Δ c )
in which R Δ = R R r e f , and we set R r e f = R R + R T . The first exponentials of s i f 1 ( t ^ , t m ) and s i f 2 ( t ^ , t m ) are the range items that produce beat frequencies f 1 i = 2 γ R Δ c and f 2 i = 2 γ R Δ c , respectively. The second exponentials of s i f 1 ( t ^ , t m ) and s i f 2 ( t ^ , t m ) are residual video phases (RVPs), and the third exponentials account for Doppler shift. The RVP can be removed by the translation motion compensation (TMC) algorithm.
Taking the time of reference point as a basis, applying Fourier transform (FT) and multiplying by a constant conversion coefficient 2 γ / c , we have
s i f 1 ( r , t m ) = i = 1 K T p σ i 2 sin c ( γ T p c ( r R Δ ) ) × exp ( j 2 π f 0 R Δ c j π γ R 2 Δ c 2 )
s i f 2 ( r , t m ) = i = 1 K T p σ i 2 sin c ( γ T p c ( r + R Δ ) ) × exp ( j 2 π f 0 R Δ c + j π γ R 2 Δ c 2 )
in which r represents the down-range domain. It is noted that two channels have the same amplitude but are symmetrical regarding zero in the down range domain. After envelope alignment and removing the residual video phase (RVP), the final synthesized signal can be obtained as
s ( r , t m ) = s i f 1 ( r , t m ) + f l i p l r ( s i f 2 ( r , t m ) ) = i = 1 K T p σ i sin c ( γ T p c ( r R Δ ) ) × exp ( j 2 π f 0 R Δ c )
in which f l i p l r ( ) represents flipping the left part of the signal in a channel to the right part.
Substituting (2) into (13), we have
s ( r , t m ) = i = 1 K σ i T p sin c ( γ T p c ( r R Δ ) ) × exp ( j 2 π f d t m ) × exp ( j 4 π λ y cos β 2 )
in which f d = 2 x w cos ( β / 2 ) λ represents the Doppler frequency and λ is the wavelength. The signal in the range domain, ignoring the constant phase term, becomes
s ( r , t m ) = i = 1 K σ i T p sin c ( γ T p c ( r R Δ ) ) × exp ( j 2 π f d t m )
The cross range compression is obtained by applying FT to (15) in terms of t m as follows:
s ( r , f ) = T p T a σ i sin c ( γ T p c ( r R Δ ) ) × sin c ( T a ( f f d ) )
in which f [ f r / 2 , f r / 2 ] is the Doppler extent, f r is the pulse repetition frequency (PRF), and T a is the observation time. From (16), it can be seen that the target image is sparse in the Range-Doppler (RD) plane, and there are a few strong scatter points in the whole plane. It can be inferred that RD imaging may fail when measured data is too limited to provide sufficient resolution. Long data-collect time is needed for high resolution in ISAR imaging, but because the targets are noncooperative and usually maneuvering, it is difficult to collect sufficient data and have long observation times, which motivates the study of high-resolution ISAR imaging algorithms based on V-FM waveforms with limited data.

3. 2D-CS-Based ISAR Imaging

3.1. Review of CS

CS theory provides a new framework for sampling and compressing signals and is able to recover a sparse signal from limited measurements with a high probability. It can overcome the limitation of Nyquist sampling theorem. A signal of interest is expressed as a vector α C N which is sparse in some orthonormal basis, and can be represented as α = Ψ s , where Ψ C L × L is a linear and orthonormal basis and s C L is a sparse vector with less K ( K < L ) non-zero elements. In other words, the inverse of Ψ represents a sparse transformation that transforms α into a sparse vector s . Suppose that the signal is collected through an inner product of signal with a different sensing matrix, which can be expressed as y = Φ α = Φ Ψ s , in which Φ C Q × L ( Q < L ) is known as the measurement matrix, and y C Q is a vector of measurement samples. It is worth noting that the number of observations Q cannot be infinitely small, satisfying Q O ( K lgL ) . Because Q is much smaller than L , the recovery of the signal s from y is ill-posed. This requires that the dictionary matrix Φ Ψ obeys the restricted isometry principle (RIP)
( 1 δ ) s 2 2 Φ Ψ s 2 2 ( 1 + δ ) s 2 2
in which δ ( 0 , 1 ) is a constant, and 2 represents the l 2 norm. The coefficient α can be exactly recovered from limited measurements by solving the l 0 norm minimization, which can be described as
{ min ( s 0 ) s . t y Φ Ψ s 2 < ε α = Ψ s
in which 0 represents the l 0 norm, min ( ) denotes the minimization, and ε is the noise level. Herein, we can conclude that the CS theory combines three key points to guarantee its performance; the signal to reconstruct is sparse or compressible, the sensing matrix makes Φ Ψ so as to satisfy the RIP, and the reconstruction algorithm is efficient and has low computational complexity.

3.2. 2D-CS Based ISAR Imaging

According to the practical scattering responses of a target, the final 2D range-azimuth image is always sparse or compressible, since only a few prominent scattering centers occupy the whole imaging plane. We establish the 2D-coupled SR model of the target’s echoes (9) and (10) based on this fact. After the RVP and constant terms are removed, we get dual channels target’s echoes as follows:
s i f 1 ( t ^ , t m ) = ( x , y ) I σ ˜ 1 ( x , y ) e j 4 π w t m λ x e j 4 π γ t ^ c y
s i f 2 ( t ^ , t m ) = ( x , y ) I σ ˜ 2 ( x , y ) e j 4 π w t m λ x e j 4 π γ t ^ c y
in which σ ˜ 1 ( x , y ) and σ ˜ 2 ( x , y ) denotes the complex amplitude of the point i ( x i , y i ) in dual channel (see Appendix A.2). In the case of limited bandwidth and small angle, the scattering points do not have Migration Through Resolution Cell (MTRC), so (19) and (20) ignore the influence of w t m on the range image.
The discrete expressions of s i f 1 ( t ^ , t m ) and s i f 2 ( t ^ , t m ) are S 1 = [ s i f 1 n m ] N × M and S 2 = [ s i f 2 n m ] N × M ( n = 0 , 1 , , N 1 ;   m = 0 , 1 , , M 1 ) , in which N is the number of range samples (range bins) and M is the number of emitted pulses (azimuth samples). Scattering distribution function σ ˜ ( x , y ) is discretely represented as a two-dimensional scattering rate distribution matrix X 1 = [ σ ˜ 1 ( x h , y k ) ] N r × N a and X 2 = [ σ ˜ 2 ( x h , y k ) ] N r × N a ( k = 0 , 1 , , N r 1 ;   h = 0 , 1 , , N a 1 ) , in which N r > N , N a > M . Then, we let ρ x denote the azimuth resolution, ρ y represent range resolution, Δ θ represent the imaging accumulation angle, B represent the bandwidth, and Δ f represent the frequency sampling interval. Thus, we get the following equation:
x h = h ρ x , y k = k ρ y ρ x = λ / ( 2 Δ θ ) , ρ y = c / ( 2 B ) B = N r Δ f , γ t ^ = n Δ f Δ θ = w N a T p }
Taking (21) into (19) and (20), we get the discrete form of the ISAR observation signal model
S 1 = [ s i f 1 n m ] N × M = h , k σ ˜ 1 k h e i 2 π m h N a e j 2 π n k N r
S 2 = [ s i f 2 n m ] N × M = h , k σ ˜ 2 k h e i 2 π m h N a e j 2 π n k N r
Taking the noise into consideration, the 2D discrete expression ISAR acquisition signals are described as follows:
S 1 = Φ X 1 Ψ T a 1 + E
S 2 = Φ X 2 Ψ T a 2 + E
in which Φ = [ ϕ 0 , ϕ 1 , , ϕ N 1 ] T C N × N r ( N < N r ) is the over-complete Fourier dictionaries in the range direction, in which ϕ n = e j 2 π n k N r . Ψ a 1 = [ φ 1 0 , φ 1 1 , , φ 1 M 1 ] T C N a × M , in which φ m = e j 2 π m h N a and Ψ a 2 = [ φ 0 , φ 1 , , φ M 1 ] T C N a × M denote the over-complete Fourier dictionaries of dual channel, respectively, in the cross-range dictionaries. Depicted in the range-Doppler plane, X 1 and X 2 C N r × N a represent the interested complex ISAR image, which contains N r range bins and N a Doppler frequency cells. E C N × M denotes the synthetic additive noise matrix. It is worth noting that the popular method to solve the 2D reconstruction problem is to use the Kronecker product and vectorization operation convert 2D-coupled SR model into 1D reconstruction problem. Taking Channel 1 as an example, then (24) is rewritten as s 1 = Θ 1 x 1 , in which s 1 = v e c ( S 1 ) , x 1 = v e c ( X 1 ) , Θ 1 = Ψ a 1 Φ , and v e c ( ) denotes the vectorization of a matrix, which converts the matrix into a column vector; stands for Kronecker product. We note that Θ 1 has N M × N r N a elements, while Φ and Ψ a 1 is determined by N N r + N a M elements. Thus, the conventional 1D signal model has a large computing memory. This paper focuses on directly solving the 2D-coupled SR model based on a 2D nonconvex optimization problem to avoid this disadvantage.
We solve the ISAR imaging problem as
min X 1 X 1 0 s . t .   S 1 Φ X 1 Ψ a 1 T F 2 ε
min X 2 X 2 0 s . t .   S 2 Φ X 2 Ψ a 2 T F 2 ε
X ^ = X ^ 1 + X ^ 2 = [ σ 1 k h + σ 2 k h ] N × M
in which 0 denotes the number of non-zero components in X 1 and X 2 , F represents the Frobenius norm of a matrix, and then ε is a small constant.
The optimization problem in (26) and (27) can be solved via 2D-SL0 algorithm described in [23]. SL0 is not based on minimizing the l 1 norm, but it directly minimizes the l 0 norm of the solution. It introduces a continuous Gaussian function to approximate the discontinuous l 0 norm, and then we can get the image X in (28) using a projective steepest ascent approach; the details of this algorithm can be found in [23]. As for V-FM waveform, we give the Two-dimensional SL0 (2D-SL0) Algorithm 1 in channel 1 shown as follows:
Algorithm 1. Two-dimensional SL0 (2D-SL0) algorithm
Initialization
Let X ^ 1 0 = Φ S 1 ( Ψ a 1 ) T .
Choose a suitable decreasing sequence for σ , [ σ 1 σ J ] .
For j = 1 , , J :
 (a).
Let σ = σ j .
 (b).
Maximize (approximately) the function F σ ( X 1 ) = i , j exp ( x 1 i j 2 / 2 σ 2 ) on the feasible set { X 1 | S 1 = Φ X 1 Ψ a 1 T } using L iterations of the steepest ascent algorithm (followed by projection onto the feasible set):
-Initialization: X 1 = X ^ 1 j 1 .
-For l = 1 , , L (Loop L times):
(1) Let Δ = [ δ i j ] , in which δ i j = exp ( x 1 i j / 2 σ 2 ) .
(2) Let X 1 X 1 μ Δ ( μ is a small positive constant).
(3) Project X 1 back onto the feasible set: X 1 X 1 Φ ( Φ X 1 Ψ T X 1 ) ( B ) T
  (c). Set X ^ 1 j = S .
Final answer is X ^ 1 = X ^ J .
in which Φ = Φ T ( Φ Φ T ) 1 and Ψ a 1 = Ψ T a 1 ( Ψ a 1 Ψ T a 1 ) 1 is the pseudoinverse of Φ and Ψ a 1 , respectively. Finally, we can get the ISAR image by (28).

4. Experimental Results

In this section, we will implement numerical simulations to verify the performance of dual-channel 2D sparse signal reconstruction in V-FM bistatic ISAR imaging; targets modeled as point scatterers are utilized in the following, as shown in Figure 4a. The V-FM signal ISAR works at the X-band with a carrier frequency of 10 GHz and bandwidth of 300 MHz. The pulse-width is 100 us, the pulse frequency is 1 kHz, and the rotating velocity is 0.05 Hz; we set N = M = 256 and N a = N r = 300 . The main simulated parameters are listed in Table 1. The bistatic ISAR imaging results via the traditional RD method and 2D-SL0 are shown in Figure 4b,c; it can be noted that the 2D-SL0 based image shown in Figure 4c has a finer resolution than by another method.
To provide a quantitative evaluation for imaging performance, we take three indexes related to radar target recognition. The first evaluation index is the entropy of the image; the larger the entropy is, the more complex the image information is and the worse the imaging result. The second evaluation index is the normalized root mean square error (RMSE) between the reconstructed image and the target model is computed, which is
R M S E = 1 N M n = 1 N m = 1 M ( X ( n , m ) X ^ ( n , m ) ) 2
In addition, the correlation coefficient (CC) is the third index, which can be expressed as follows:
C C = C o v ( X , X ^ ) D ( X ) D ( X ^ )
in which C o v ( ) and D ( ) represent the covariance and variance coefficient between two matrices, respectively.
In the following imaging experiments, these indexes are adopted as the quantitative indexes to evaluate the Bi-ISAR image quality with different reconstruction methods, including 2DFFT, FFT in down-range domain, SL0 in cross-range domain (FFTSL0), and 2DSL0.

4.1. Imaging Results Versus Bi-Angle

First, we will examine the performance of different reconstruction methods versus bi-angle (60°, 90°, 120°), and the imaging results are shown in Figure 5; the first, second, third columns represent 2DFFT, FFTSL0, 2DSL0 images, respectively, and the first, second, third rows represent images with bi-angle 60°, 90°, and 120°, respectively. We can note that the imaging results become defocused as the bi-angle increases, and 2DSL0 method can provide better focal quality with high resolution than other methods, and, when the bi-angle is 120°, the quality of the image becomes worse. From Figure 5, it implies that the traditional Fourier transform or 2D-CS method cannot eliminate the influence of the bistatic angle and another method is required to overcome this problem.
To characterize the performance quantitatively, the corresponding curves of entropy, RMSE, and CC are shown in Figure 6. In Figure 6a, it is obvious that the entropy obtained by using 2DSL0 method is lower than that of both 2DFFT and FFTSL0 method. Lower entropy suggests that the quality of imaging result is better. In Figure 6b, the 2DSL0 method performs best among three methods, since it retains the smallest reconstruction errors. From Figure 6c, we find that when the bi-angle is lower than 30°, the reconstructed Bi-ISAR image with the 2DSL0 method is very close to the original image, and CC values increase when bi-angle decreases.

4.2. Performance Versus SNRs

To prove that the robustness of different reconstruction methods degrades with increasing noise, we analyzed the effect of different SNRs on Bi-ISAR imaging; simulations with Gaussian distributed complex noise with three SNRs of 5 dB, 0 dB, and −5 dB added into the simulated VFM signals were performed, and the imaging results are shown in Figure 7. The first, second, and third columns represent 2DFFT, FFTSL0, and 2DSL0 images, respectively, and first, second, and third rows represent images with SNR of 5 dB, 0 dB, and −5 dB, respectively. With the increase in SNR in experiments, the focusing performance of the 2DSL0 method degrades slightly; the 2DSL0 can eliminate much more noise than other methods.
Entropy, RMSE, and CC are also used to analyze their performance quantitatively, as shown in Figure 8. In Figure 8a, the entropy values of 2DSL0 are much lower than those of both 2DFFT and FFTSL0 method, which means 2DSL0-based imaging results have lower chaos than others methods do, and have the best imaging as a result. In Figure 8b, the 2DSL0 method performs best among three methods, since it retains the smallest reconstruction errors, and the RMSE values decrease under an increase of SNR. From Figure 6c, it notes that when the SNR is higher than 15 dB, the reconstructed Bi-ISAR image with 2DSL0 method is very close to the original image, and the results have a positive correlation with SNR.

4.3. Experiments with Different Echo Data

In the following, the Bi-ISAR images are generated with different measurements and constant SNR (20 dB), and the measured data used changes from 40% to 80% of echo data. Figure 9 demonstrates the Bi-ISAR imaging results with different reconstruction methods under different measurements (50%, 65%, and 78% of echo data). The first, second, and third columns represent 2DFFT, FFTSL0, 2DSL0 images, respectively, and the first, second, and third rows represent images with measurements (50%, 65%, and 78% of echo data, respectively). It is noticeable that among these Bi-ISAR images, the 2D resolution and focusing energy obtained by the 2DSL0 method is the best, and the 2DFFT method the worst. There is aliasing in cross-range domain under incomplete measurements when using traditional Fourier transform.
Taking three indexes to analyze the quality of imaging results with different reconstruction versus measurements, as shown in Figure 10. In Figure 10a, the entropy values of 2DSL0 are much lower than those of both 2DFFT and FFTSL0 method, especially when the measurements of echo data are low. In Figure 8b, the 2DSL0 method performs best among three methods, since it retains the smallest reconstruction errors, and more echo data provides more information for imaging. In Figure 10c, the CC values of the 2DSL0 method are larger than those both methods. The larger CC is higher than the similarity between the reconstructed Bi-ISAR image and the referenced image. This implies that the 2DSL0 method can perform more accurate reconstruction than the other two methods.

5. Conclusions

In this paper, based on the characteristics of the V-FM signal, we have proposed a dual-channel 2D bistatic ISAR imaging method, which has different dictionaries in each channel and solves a nonconvex optimization problem that can be handled efficiently and reliably. By analyzing the signal and image model of VFM signal bistatic ISAR via traditional Fourier transform and then establishing the 2D SR model of the target’s echo, one can reconstruct the images using the 2D-CS method. Using the simulated data of an aircraft to validate the 2DSL0 based imaging processing, we show that the proposed method maintains better performance than other methods in a low SNR environment and with low measurement ratios of echo data.

Author Contributions

Conceptualization, J.C., W.W. and X.P.; Methodology, J.C. and L.X.; Software, J.C. and X.P.; Validation, J.C. and X.P.; Formal Analysis, J.C.; Investigation, L.X.; Resources, W.W.; Data Curation, J.C.; Writing-Original Draft Preparation, J.C.; Writing-Review & Editing, J.C. and X.P.; Visualization, J.C.; Funding Acquisition, X.P. and W.W.

Funding

This research was funded by National Natural Science Foundation of China under Grant 61701507 and Grant 61571451.

Acknowledgments

The authors would like to thank the editors and receivers for their valuable comments and suggestions to improve the quality of this paper.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Appendix A.1. Derivation of Instantaneous Range in (1)

During short CPI, assume bistatic angle remains constant, the instantaneous distance in Figure 1 can be represented as
R = R T 2 + r i 2 2 R T r i cos ( π 2 + β 2 + θ i ) + R R 2 + r i 2 2 R R r i cos ( π 2 β 2 + θ i )
in which θ i = θ 0 + Δ θ i , in far field condition. We have R T > r i , R R > r i , then (A1) can be approximately expressed by
R R T r i cos ( π 2 + β 2 + θ i ) + R R r i cos ( π 2 β 2 + θ i )
in which x i = r i cos θ 0 , y i = r i sin θ 0 , then (A2) can be rewritten as
R R T + R R + 2 cos β 2 [ y i cos Δ θ i + x i sin Δ θ i ]

Appendix A.2. Derivation of the Complex Amplitude in (19) and (20)

Take (2) into (9) and (10), we will get the amplitude shown as follows:
σ ˜ 1 ( x i , y i ) = σ i rect ( t ^ R c T p ) exp ( j 2 π γ ( t ^ T r e f ) c x i w t m )
σ ˜ 2 ( x i , y i ) = σ i rect ( t ^ R c T p ) exp ( j 2 π γ ( T r e f t ^ ) c x l i w t m )

References

  1. Chen, V.C.; Martorella, M. Inverse Synthetic Aperture Radar Imaging: Principles, Algorithms and Applications; SciTech Publishing: Edison, NJ, USA, 2014. [Google Scholar]
  2. 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]
  3. Chen, V.C.; Miceli, W.J. Simulation of ISAR imaging of moving targets. IEEE Proc. Radar Sonar Navig. 2001, 148, 160–166. [Google Scholar] [CrossRef]
  4. Chen, V.C.; Rosiers, A.D.; Lipps, R. Bi-static ISAR Range-Doppler Imaging and Resolution Analysis. In Proceedings of the 2009 IEEE Radar Conference, Pasadena, CA, USA, 4–8 May 2009; pp. 1–5. [Google Scholar]
  5. Martorella, M.; Palmer, J.; Littleton, B.; Longstaff, I.D. On bistatic inverse synthetic aperture. IEEE Trans. Aerosp. Electron. Syst. 2007, 43, 1125–1134. [Google Scholar] [CrossRef]
  6. Ozdemir, C.; Demirci, S.; Yilmaz, B.; Ak, C.; Yigit, E. A new and practical formulation of bistatic inverse synthetic aperture radar imaging and verification of the formulation using numerical examples. In Proceedings of the 2009 International Conference on Electrical and Electronics Engineering-ELECO 2009, Bursa, Turkey, 5–8 November 2009; pp. II-161–II-164. [Google Scholar]
  7. Skolink, M.I. Introduction to Radar Systems; McGraw-Hill: New York, NY, USA, 2001. [Google Scholar]
  8. Sitler, K.L.; Temple, M.A.; Novack, R.C.; Hughes, J.A. High range resolution profiling using phase-coded, stepped-frequency waveforms. Electron. Lett. 2002, 38, 46–48. [Google Scholar] [CrossRef]
  9. Kim, K.T. Focusing of high range resolution profiles of moving targets using stepped frequency waveforms. IET Radar Sonar Navig. 2010, 4, 564–575. [Google Scholar] [CrossRef]
  10. Pan, X.; Wang, W.; Feng, D.; Liu, Y.; Wang, G. Signatures extraction of rotating targets based on a fraction of HRRPs. IEEE Trans. Antennas Propag. 2015, 63, 585–592. [Google Scholar] [CrossRef]
  11. Pan, X.; Wang, W.; Wang, G. Sub-Nyquist Sampling Jamming against ISAR with CS-Based HRRP Reconstruction. IEEE Sens. J. 2016, 16, 1597–1602. [Google Scholar] [CrossRef]
  12. Gu, Z.; Liu, J.; Pan, X. ISAR Imaging with Wideband V-FM Waveforms via Dual-channel CS-D. IEEE Access 2017, 5, 18557–18563. [Google Scholar] [CrossRef]
  13. Li, S.; Zhao, G.; Zhang, W.; Qiu, Q.; Sun, H. ISAR Imaging by Two-dimensional Convex Optimization-Based Compressive Sensing. IEEE Sens. J. 2016, 16, 7088–7093. [Google Scholar]
  14. Patel, V.M.; Easley, G.R.; Healy, D.M., Jr.; Chellappa, R. Compressed synthetic aperture radar. IEEE J. Sel. Top. Signal Process. 2010, 4, 244–254. [Google Scholar] [CrossRef]
  15. Zhang, L.; Qiao, Z.J.; Xing, M.D.; Sheng, J.L.; Guo, R.; Bao, Z. High resolution ISAR imaging by exploiting sparse apertures. IEEE Trans. Antennas Propag. 2012, 60, 997–1008. [Google Scholar] [CrossRef]
  16. Hu, J.M.; Zhou, W.; Fu, Y.W.; Li, X.; Jing, N. Uniform rotational motion compensation for ISAR based on phase cancellation. IEEE Geosci. Remote Sens. Lett. 2011, 8, 636–640. [Google Scholar] [CrossRef]
  17. Feng, D.; Pan, X.; Gu, Z. High-Resolution ISAR imaging with Wideband V-FM Waveforms. Int. J. Antennas Propag. 2017, 2017, 1908204. [Google Scholar] [CrossRef]
  18. Zhang, L.; Xing, M.D.; Qiu, C.W.; Li, J.; Sheng, J.L.; Li, Y.C.; Bao, Z. Resolution enhancement for inversed synthetic aperture radar imaging under low SNR via improved compressive sensing. IEEE Trans. Geosci. Remote Sens. 2010, 48, 3824–3838. [Google Scholar] [CrossRef]
  19. Qiu, W.; Zhao, H.; Zhou, J.; Fu, Q. High-resolution fully polarimetric ISAR imaging based on compressive sensing. IEEE Trans. Geosci. Remote Sens. 2014, 52, 6119–6131. [Google Scholar] [CrossRef]
  20. He, X.; Tong, N.; Hu, X.; Feng, W. High-Resolution ISAR imaging based on two-dimensional group sparse recovery. IET Radar Sonar Naving. 2017, 12, 82–86. [Google Scholar] [CrossRef]
  21. Zhang, S.; Zhang, W.; Zong, Z.; Tian, Z.; Yeo, T.S. High-resolution bistatic ISAR imaging based on two-dimensional compressed sensing. IEEE Trans. Antennas Propag. 2015, 63, 2098–2111. [Google Scholar] [CrossRef]
  22. Li, G.; Zhang, H.; Wang, X.; Xia, X.G. ISAR 2-D imaging of uniformly rotating targets via matching pursuit. IEEE Trans. Aerosp. Electron. Syst. 2012, 48, 1838–1846. [Google Scholar] [CrossRef]
  23. Ghaffari, A.; Babaie-Zadeh, M.; Jutten, C. Sparse decomposition of two dimensional signals. In Proceedings of the 20th European Signal Processing Conference, Taipei, Taiwan, 19–24 April 2009; pp. 3157–3160. [Google Scholar]
Figure 1. ISAR (Inverse Synthetic Aperture Radar) imaging geometric model.
Figure 1. ISAR (Inverse Synthetic Aperture Radar) imaging geometric model.
Sensors 18 03082 g001
Figure 2. Sketch map of V-FM (V style modulation) waveforms.
Figure 2. Sketch map of V-FM (V style modulation) waveforms.
Sensors 18 03082 g002
Figure 3. (a) Ambiguity function of LFM signal; (b) ambiguity function of VFM signal.
Figure 3. (a) Ambiguity function of LFM signal; (b) ambiguity function of VFM signal.
Sensors 18 03082 g003
Figure 4. (a) Aircraft model, (b) RD-based image, and (c) 2D-SL0-based image.
Figure 4. (a) Aircraft model, (b) RD-based image, and (c) 2D-SL0-based image.
Sensors 18 03082 g004
Figure 5. Bi-ISAR (Bistatic ISAR) imaging results with different reconstruction methods versus bi-angle. Column 1 represents 2DFFT images, Column 2 represents FFTSL0 images, and Column 3 represents 2DSL0 images. Row 1 represents bi-angle 60°, Row 2 represents bi-angle 90°, and Row 3 represents bi-angle 120°.
Figure 5. Bi-ISAR (Bistatic ISAR) imaging results with different reconstruction methods versus bi-angle. Column 1 represents 2DFFT images, Column 2 represents FFTSL0 images, and Column 3 represents 2DSL0 images. Row 1 represents bi-angle 60°, Row 2 represents bi-angle 90°, and Row 3 represents bi-angle 120°.
Sensors 18 03082 g005
Figure 6. Curves of imaging performance under different bi-angles. (a) Entropy versus bi-angle; (b) RMSE (root mean square error) versus bi-angle; (c) CC (correlation coefficient) versus bi-angle.
Figure 6. Curves of imaging performance under different bi-angles. (a) Entropy versus bi-angle; (b) RMSE (root mean square error) versus bi-angle; (c) CC (correlation coefficient) versus bi-angle.
Sensors 18 03082 g006
Figure 7. Bi-ISAR imaging results with different reconstruction methods under different SNRs. Column 1 represents 2DFFT images, Column 2 represents FFTSL0 images, and Column 3 represents 2DSL0 images. Row 1 represents SNR is 5 dB, Row 2 represents SNR is 0 dB, and Row 3 represents SNR is −5 dB.
Figure 7. Bi-ISAR imaging results with different reconstruction methods under different SNRs. Column 1 represents 2DFFT images, Column 2 represents FFTSL0 images, and Column 3 represents 2DSL0 images. Row 1 represents SNR is 5 dB, Row 2 represents SNR is 0 dB, and Row 3 represents SNR is −5 dB.
Sensors 18 03082 g007
Figure 8. Curves of imaging performance under different SNRs. (a) Entropy versus SNRs, (b) RMSE versus SNRs, and (c) CC versus SNRs.
Figure 8. Curves of imaging performance under different SNRs. (a) Entropy versus SNRs, (b) RMSE versus SNRs, and (c) CC versus SNRs.
Sensors 18 03082 g008
Figure 9. Bi-ISAR imaging results with different reconstruction methods versus measurements. Column 1 represents 2DFFT images, Column 2 represents FFTSL0 images, and Column 3 represents 2DSL0 images. Row 1 represents 50% of echo data, Row 2 represents 65% of echo data, and Row 3 represents 78% of echo data.
Figure 9. Bi-ISAR imaging results with different reconstruction methods versus measurements. Column 1 represents 2DFFT images, Column 2 represents FFTSL0 images, and Column 3 represents 2DSL0 images. Row 1 represents 50% of echo data, Row 2 represents 65% of echo data, and Row 3 represents 78% of echo data.
Sensors 18 03082 g009
Figure 10. Curves of imaging performance under different measurements. (a) Entropy versus measurements, (b) RMSE versus measurements, and (c) CC versus measurements.
Figure 10. Curves of imaging performance under different measurements. (a) Entropy versus measurements, (b) RMSE versus measurements, and (c) CC versus measurements.
Sensors 18 03082 g010
Table 1. Simulated parameters.
Table 1. Simulated parameters.
f 0 10 GHz R T 500 km
B 300 MHz R R 600 km
T p 100   μ s β 30°
P R F 1 kHz w 0.05 Hz

Share and Cite

MDPI and ACS Style

Chen, J.; Pan, X.; Xu, L.; Wang, W. Bistatic ISAR Imaging with a V-FM Waveform Based on a Dual-Channel-Coupled 2D-CS Algorithm. Sensors 2018, 18, 3082. https://doi.org/10.3390/s18093082

AMA Style

Chen J, Pan X, Xu L, Wang W. Bistatic ISAR Imaging with a V-FM Waveform Based on a Dual-Channel-Coupled 2D-CS Algorithm. Sensors. 2018; 18(9):3082. https://doi.org/10.3390/s18093082

Chicago/Turabian Style

Chen, Jiyuan, Xiaoyi Pan, Letao Xu, and Wei Wang. 2018. "Bistatic ISAR Imaging with a V-FM Waveform Based on a Dual-Channel-Coupled 2D-CS Algorithm" Sensors 18, no. 9: 3082. https://doi.org/10.3390/s18093082

APA Style

Chen, J., Pan, X., Xu, L., & Wang, W. (2018). Bistatic ISAR Imaging with a V-FM Waveform Based on a Dual-Channel-Coupled 2D-CS Algorithm. Sensors, 18(9), 3082. https://doi.org/10.3390/s18093082

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