Xidian University, National Laboratory of Radar Signal Processing, 2 South Taibai Road, Xi’an 710071, China
In inverse synthetic aperture radar (ISAR) imaging of targets with complex motion, the azimuth echoes have to be modeled as multicomponent cubic phase signals (CPSs) after motion compensation. For the CPS model, the chirp rate and the quadratic chirp rate deteriorate the ISAR image quality due to the Doppler frequency shift; thus, an effective parameter estimation algorithm is required. This paper focuses on a parameter estimation algorithm for multicomponent CPSs based on the local polynomial ambiguity function (LPAF), which is simple and can be easily implemented via the complex multiplication and fast Fourier transform. Compared with the existing parameter estimation algorithm for CPS, the proposed algorithm can achieve a better compromise between performance and computational complexity. Then, the high-quality ISAR image can be obtained by the proposed LPAF-based ISAR imaging algorithm. The results of the simulated data demonstrate the effectiveness of the proposed algorithm.
The inverse synthetic aperture radar (ISAR) imaging technique for moving targets has attracted the attention of many radar researchers in the past three decades due to its significance in both civil and military applications.1–6 Generally, in order to obtain a well-focused ISAR image, the first procedure is to implement motion compensation, which includes range alignment (the translational and rotational migrations) and phase adjustment (the Doppler phase caused by the translation).6 Then the conventional range-Doppler (RD) algorithm can be used to generate a focused ISAR image, which indicates that, during the coherent integration interval, the Doppler frequency shift must remain constant. However, instances of targets that are usually uncooperative and engaged in complicated maneuvers often occurs.7–14 In this case, the RD algorithm does not work well for the time-varying characteristic of Doppler frequency. To deal with this problem, many articles have developed different solution for different degree of maneuverability, and they can fall into two categories: linear frequency modulation (LFM) signal model6–9 for targets with smooth motion and cubic phase signal (CPS) model2,4,5,10–14 for targets with complex motion. For the slow-maneuvering target, the assumption that the azimuth echo is characterized as LFM signal is appropriate and practical, and has been demonstrated by Ref. 6 with experimental ship data. Many representative algorithms based on LFM signal model have been proposed, such as Radon Wigner transform,3 the modified Wigner-Ville distribution (WVD),6 the stretch keystone-Wigner transform,7 the Lv’s distribution,8 and the fractional Fourier transform.9 For targets with complex motion, in which rotational motion may cause higher-order phase terms, the traditional LFM signal model is not appropriate. For a more reasonable and accurate approximation, many articles have characterized the azimuth echo in a range bin as multicomponent CPSs, which is identified to be much more appropriate than the LFM model.10–14
With regard to the parameter estimation of CPS, numerous algorithms, including the cubic phase function (CPF),15,16 the higher-order ambiguity function,17 the product generalized CPF (PGCPF),18 the product high-order matched-phase transform,19 and the modified version of CPF (MCPF),14 have been addressed, but all of them involve multilinear transformation and discrete Fourier transform for the nonuniformly spaced signal sample, which result in extensive cross-terms under multi-CPSs and undesirable computational complexity.10–14,18,19–23 In Ref. 20, the parameter estimation algorithm based on the local polynomial Wigner distribution (LPWD) is proposed and has been successfully applied to ISAR imaging.21,22 However, after compensating the third-order term with the estimated parameter, the LPWD algorithm estimates the second-order coefficient using CPF method, which does lead to the heavy computational burden due to the Fourier transform with respect to the nonuniformly spaced data. In general, the performance and computational cost are conflicting requirements, e.g., the modified discrete chirp Fourier transform for CPS can obtain robust performance against noise and handle multi-CPSs, but it requires a multidimensional search.
Inspired by the proposed LPWD in Refs. 21 and 22 and considering the tradeoff between the computational cost and the performance of parameter estimation, an effective parameter estimation algorithm for multi-CPSs is introduced based on the local polynomial ambiguity function (LPAF) in this paper, which does not suffer from the considerable cross-term interference and heavy computational burden due to its moderate order of nonlinearity, coherent integration via fast Fourier transform (FFT). Then, the high-quality ISAR image can be achieved by the corresponding ISAR imaging algorithm based on LPAF.
This paper is organized as follows. ISAR imaging model with complex motion is given in Sec. 2. In Sec. 3, the parameter estimation algorithm for CPS is proposed based on LPAF. In Sec. 4, the performance of LPAF is discussed in detail, which mainly includes the computational cost and the robustness in the presence of noise. The LPAF-based ISAR imaging algorithm is presented in Sec. 5. Finally, the application of the LPAF-based ISAR imaging algorithm and the conclusion are given in Secs. 6 and 7, respectively.
ISAR imaging model on the imaging project plane is shown in Fig. 1, where the radar is supposed to locate at the same plane with the target, and to orient along the X axis. Pi(xi,yi) is the location on the i’th point scatterer in the Cartesian coordinates with its origin at the geometric center O of the target, and the i’th point scatterer can also be described as Pi(ri,θi), where xi=ri cos θi, yi=ri sin θi. Rt and θ(tm) are the translational distance and rotational displacement, respectively.
ISAR imaging model on the project plane.
The instant distance R(tm) between the radar and the i’th point scatterer at the instant slow time tm can be approximated as Display Formula
For targets with complex motion, e.g., the highly maneuvering airplane and ship fluctuation with oceanic waves, Rt and θ(tm) are time-varying and can be expressed in Taylor series as follows: Display Formula
Assume the transmitted LFM signal takes the following form: Display Formula
After the pulse compression by using the matched filter (MF) H(t^)=rect(t^/T)exp(jπμt^2), we have Display Formula
In Eq. (7), the complex motion will result in range migration and Doppler frequency shift. The translational range migration has the same feature for all scatterers and has no contribution for ISAR imaging. Hence, in this paper, the standard range alignment method24 and the phase gradient autofocus method25 are utilized to eliminate the effect caused by translation, and keystone transform to correct the rotational range migration. Here, we only concentrate on the Doppler frequency shift caused by rotation.
Consequently, after the motion compensation (the range migration compensation and the translational induced phase error correction), the azimuth echo of the i’th scatter in the range-azimuth domain can be rewritten as Display Formula
In realistic ISAR imaging applications, there will be multiple scatterers in a range bin simultaneously, which is more accurate and has been confirmed in many literatures. Thus, the azimuth echo signals in a range bin can be rewritten as Display Formula
LPWD,20 a special kind of time-frequency distribution proposed by Stankovic, is widely used for the parameters estimation of multi-CPSs. Then the ISAR imaging algorithm21,22 is also proposed based on LPWD. LPWD is defined in Ref. 20 as follows: Display Formula
Thereby, the instantaneous frequency estimation of polynomial phase signal (PPS) can be achieved with a significantly reduced distribution dimension.
Considering an analytic signal s(t), the ambiguity function (AF) is defined as Display Formula
Borrowing the idea from the definitions of LPWD and AF, the LPAF is defined as follows: Display Formula
Similar to the definition of LPWD, we have Display Formula
Consider mono-CPS with the following form: Display Formula
Substituting Eq. (15) in Eq. (16), we obtain Display Formula
From Eq. (19), it can be seen that the different lag-time instance τ0 indicates the different zoom factor, which can enlarge the estimation ranges of CR and QCR to avoid the spectrum aliasing and estimation error26 in some extreme ISAR imaging case, e.g., higher maneuverability. On the other hand, the symmetric instantaneous self-correlation operation will shorten the effective signal in Eq. (16) when τ0≠0. Based on the recommendation from Ref. 27, the optimal discrete lag is Nτ0/2opt=(N+1)/2 H, where H is the PPS order. For convenience to analyze, we choose τ0=1 so that the values of CR and QCR can be read directly from CRQCRD.
Thereby, we can estimate the parameters ϕ2 and ϕ3 simultaneously by searching the position of the maximum of |LPAF(ω1,ω2)|. Display Formula
With the estimated parameters (ϕ^2,ϕ^3) and dechirp technique, the other parameters can be estimated via FFT after compensating the high-order phase terms.
According to the aforementioned analyses on LPWD and LPAF for CPS, we know that LPAF is similar to LPWD except for the definition of θ(ω→,Δ) in Eqs. (11) and (14), where Δ=τ for LPWD and Δ=t for LPAF. However, instead of the special time instance as t=0,21 LPAF can choose the different lag-time instance τ0 to introduce a zoom factor adjusting the estimation ranges of CR and QCR in Eq. (19). It must be emphasized that such a difference circumvents the problems of LPWD including spectrum aliasing, error propagation, and computational burden, which will be discussed in Sec. 4.
In this section, the performance of LPAF for multi-CPSs is shown. As we know, correlation operation will result in the cross-term effects among multi-CPSs, which determines whether the proposed algorithm is suitable for multi-CPSs or not. For simplicity, we choose special case as τ0=1 and consider two-component CPSs in the following structure: Display Formula
Substituting Eq. (21) in Eq. (16), we obtain the LPAF for s(t) in Eq. (21) as follows: Display Formula
Note that we have two cross-terms LPAFpq,cross(ω1,ω2) and LPAFqp,cross(ω1,ω2), and we only take LPAFpq,cross(ω1,ω2) as an example to analyze. Obviously, LPAFpq,cross(ω1,ω2) is with the form Display Formula
As is obvious from Eq. (23), when ϕp,3≠ϕq,3, the auto-terms can be detected correctly by peak detection, and the phase coefficients can be estimated easily. Here, assume that the received signal contains the same QCR, i.e., ϕp,3=ϕq,3. Then, LPAFpq,cross(ω1,ω2) can be further simplified as Display Formula
From Eq. (24), we can find that the spurious peak will be located at Display Formula
It is worthwhile to mention that only when ϕp,3=ϕq,3, the LPAF for s(t) causes spurious peak, which affects the detection and estimation for multi-CPSs. Fortunately, the phase coefficients ϕj(j=0,1,2…) are not only related to target rotational parameter (e.g., CF corresponding to the initial rotational velocity, CR corresponding to the rotational acceleration, and QCR corresponding to the rotational acceleration rate), but also proximately proportional to the cross-range coordinates. Thus, for different scatterers with different coordinates in a range bin, we have ϕp3≠ϕq3(p≠q).14 Accordingly, the cross-terms of the LPAF induced by multi-CPSs can be avoided effectively.
Example 1. Consider a mono-component CPS and two-component CPSs, respectively. The sampling frequency and the effective length of signal are 256 Hz and 256, respectively. The signal parameters are shown in Table 1.
The WVD of CPS-A1 is shown in Fig. 2(a), in which the curves appear due to QCR of CPS. With the proposed parameter estimation algorithm, a sole peak appears in the CRQCRD in Figs. 2(b) and 2(c), where CR and QCR can be obtained by peak detection. Moreover, we obtain the CRQCRD of the two-component CPSs in Fig. 2(d). For easy visualization, the stereogram, the right view (XZ plane), and the left view (XZ plane) of CRQCRD are depicted in Figs. 2(d)–2(f), where two different peaks appear in CRQCRD. Therefore, the corresponding parameters of each component can be easily estimated via searching the peaks and CLEAN technique,7,8 sequentially.
Simulation results of example 1: (a) Wigner-Ville distribution of first CPS, (b) stereogram of CRQCRD in first CPS, (c) contour of (b), (d) stereogram of CRQCRD in two CPSs, (e) XZ plane of (d), and (f) YZ plane of (d).
In this section, we will discuss the performance of LPAF, mainly including the computational cost and the robustness in the presence of noise, which determine whether the proposed algorithm can be utilized for the parameter estimation of multi-CPSs in realistic ISAR imaging applications or not. In order to verify the advantages of the proposed estimation algorithm, other two representative algorithms, MCPF and LPWD, are chosen as reference to compare with the proposal.
For MCPF, the estimation procedures of the parameter ϕ3 require O(3Ntm2) complex multiply operations (CMs) for definition of multilinear function and O(Ntm3) operations for discrete Fourier transform (DFT) along the nonuniformly spaced lag-time axis. For estimation of the parameter ϕ2, the implementation procedures require O(Ntm2) operations for CPF and O(Ntm3) operations for DFT due to the two-dimensional (2-D) kernel. Thus, the computational complexity with MCPF is in the order of O(Ntm3). According to Refs. 11 and 13, the nonuniform fast Fourier transform (NUFFT)11,13 can be utilized to speed up the implementation of the Fourier transform along the nonuniformly spaced lag-time axis. Unfortunately, the interpolation necessary for NUFFT is not preferable in some realistic applications.
The LPAF algorithm, which can directly obtain CRQCRD, estimates CR and QCR with 2-D maximizations and can reduce error propagation effect to some extent. According to the aforementioned analyses and discussion on the modified discrete chirp-Fourier transform,28,29 which retains capability of efficient FFT implementation, the implementation procedures of parameters estimation based on LPAF include (1) the self-correlation operation in Eq. (16), after which the signal has the form of the linear frequency modulated signal, (2) FFT operation after matching the QCR ϕ3 in Eq. (18), and (3) the parameters estimation by peak detection technique. In this paper, the searching range of the parameter ϕ3 is [−ROUND(Ftm/2),ROUND(Ftm/2)], where ROUND(·) and Ftm denote the round-up operator and pulse repetition frequency (PRF), respectively. The resolution of the parameter ϕ3 is PRF/Ntm, where Ntm is the effective length of tm (the same as the number of effective echoes); thus, the searching point of the parameter ϕ3 is exactly Ntm in ISAR imaging. In conclusion, the computational complexity with this approach is about in the order of O(Ntm2 log2 Ntm).
For LPWD, the implementation procedures for the estimation of parameters ϕ1 and ϕ3 are analogous to the LPAF algorithm. However, the parameter ϕ2 is estimated by the integrated form of cubic phase function, which includes the defined bilinear function [O(Ntm2)] and Fourier transform [O(Ntm3)] along the nonuniformly spaced lag-time axis. Therefore, the computational complexity with the LPWD-based algorithm is also in the order of O(Ntm3), which is obviously greater than the LPAF algorithm.
It is obvious that, compared with MCPF and LPWD, the proposed algorithm is computationally more efficient than the other two algorithms since the latter involves nonuniformly spaced data computation. Furthermore, as mentioned in Ref. 21, MCPF estimates phase coefficients one by one via dechirp technique, which will reduce the estimation precision due to error propagation. On the contrary, the LPAF algorithm can obtain CR and QCR simultaneously. Table 2 gives the computational costs of the three algorithms.
In order to further illustrate the practicability of the proposed parameter estimation algorithm, the robustness of performance against noise, which is evaluated by input-output signal-to-noise ratio (SNR)9 and mean square error (MSE),30,31 will be analyzed in this section.
Example 2. Consider a mono-CPS as shown in Eq. (15); the sampling frequency and the effective length of signal are 256 Hz and 256, respectively. The signal parameters are set as follows: A=1, ϕ1=40 Hz, ϕ2=24 Hz/s, ϕ3=60 Hz/s2. Hereon, the CPS is embedded in complex white Gaussian noise, and the value of SNRin varies from −5 to 2 dB in steps of 1 dB, and 100 Monte Carlo simulations are performed for each SNRin value. Meanwhile, the MSEs for the parameters ϕ2 and ϕ3 with the corresponding Cramer-Rao bounds (CRBs)30 are also obtained and shown in Figs. 3(b) and 3(c), respectively.
Robustness of performance against noise: (a) input-output SNR, (b) MSE of CR estimation, and (c) MSE of QCR estimation.
In Fig. 3(a), MF, LPWD, and MCPF are chosen to compare with the proposed LPAF. For MCPF, its implementation is similar to PGCPF,18 which defines a multilinear function of fourth order and estimates unknown parameters one by one with two one-dimensional maximizations. Thus, the performance of MCPF closely matches that of MF until the input SNR falls below 0 dB, and the evidence for MCPF can be found in Ref. 10. For LPWD, its performance also closely matches that of MF until the input SNR falls below 0 dB, and we can find the corresponding evidence in Ref. 21. According to the analyses in Refs. 8 and 11, compared to the self-correlation function [s(τ/2)s*(−τ/2), for t=0] of LPWD, the self-correlation function [s(t+τ0/2)s*(t−τ0/2), for τ=τ0] of LPAF can reduce the influence of nonlinearity and benefit the performance in the presence of noise. Thereby, the proposed LPAF can obtain robust performance against noise and its threshold SNR is −3 dB. In addition, the MSEs for the estimations ϕ2 and ϕ3 are plotted in Figs. 3(b) and 3(c) as a function of SNR, respectively. Obviously, the MSEs of CR and QCR are close to the corresponding CRB when SNRin≥−3 dB, which conforms to the result in Fig. 3(a).
On the above analyses and confirmation, a new ISAR imaging algorithm for target with complex motion is proposed based on LPAF. Detail implementation procedures are given as follows:
In this section, we will utilize the synthetic ship target and aircraft target with the model shown in Figs. 4(a) and 6(a) to verify the proposed ISAR imaging algorithm successively. Motion parameters and corresponding radar parameters are listed in Tables 3 and 4, respectively.
(a) Ship target model, (b) result of the conventional RD algorithm, (c) result after the RD algorithm, (d) result after the Doppler spread compensation with LPAF, (e) result after the Doppler spread compensation with LPWD, and (f) result after the Doppler spread compensation with MCPF.
Similar to Refs. 10 and 13, a fluctuating ship target modeled as a set of ideal scatterers is provided to demonstrate the effectiveness of the proposed algorithm in this paper. After pulse compression and motion compensation, the received signal can be simplified as the 2-D reflectivity function shown in Eq. (8). Figure 4(b) gives the conventional RD imaging results with FFT operation along the azimuth slow time axis under the situation of SNRin=−3 dB, which is blurred severely by CR and QCR. In this section, we choose the special lag-time instance τ0=1 to achieve the LPAF so that CR and QCR can be read directly.
To demonstrate the accuracy of the parameters estimation by the proposed LPAF, the data of the 153rd range bin in Fig. 4(b) are extracted to analyze, and the corresponding CRs and QCRs of these scatterers are estimated as shown in Table 5. Then, with the estimated CRs and QCRs, we compensate the Doppler frequency shift and obtain Figs. 4(d)–4(f), where the ideal matched results are also shown to make a comparison with each mentioned algorithm. Compared with the results in Figs. 4(c), 4(e), and 4(f), the result obtained by LPAF in Fig. 4(d) is very close to the ideal results, which demonstrate that the accuracy of the parameters estimation by the proposed LPAF is high. Although some spurious peaks appear in Fig. 4(d), the same as in Table 5, the LPAF is also better than other algorithms in the estimation accuracy and suppression on cross-terms.
In addition, the comparison of the LPAF-based algorithm with the MCPF-based algorithm and the LPWD-based algorithm for ISAR imaging is discussed. Generally, to quantitatively evaluate the performance of the proposed algorithm, the MSEs of 2-D reconstructed coordinates are defined as MSE=‖Rest−R‖2, where R and Rest represent the original data and the reconstructed data obtained by the proposed algorithm, that is, the corresponding relation between R and Rest is determinate and known ahead of time. However, the number and the position of the scatterers on target are usually unknown and cannot be obtained in advance. Entropy can be used to measure the smoothness of a distribution function. The smoother a distribution function, the larger is its entropy.10,13 Owing to this property, we use the entropy of Eq. (30) as a criterion to measure the quality of the image in Table 6. Better focus results in a sharper image and thus smaller entropy. Meanwhile, simulation times of Fig. 5 (only including effective range bins) are also listed to certify the high efficiency of LPAF-based ISAR imaging algorithm. Display Formula
Simulation results: (a) result of LPAF-based ISAR imaging algorithm, (b) result of LPWD-based ISAR imaging algorithm, and (c) result of MCPF-based ISAR imaging algorithm.
From Fig. 5, we see that LPAF-based ISAR imaging algorithm has better suppression on cross-terms compared with the LPWD-based ISAR imaging algorithm and MCPF-based ISAR imaging algorithm. Evidently, most scatterers are relocated correctly in Fig. 5(a). Due to the four order of nonlinearity of MCPF and error propagation effect of LPWD (the estimation precision for CR is influenced by QCR), the spurious scatterers in Figs. 5(b) and 5(c) are more than those in Fig. 5(a). Furthermore, the entropy of Fig. 5(a) is also smaller than those of Figs. 5(b) and 5(c) in Table 6, which further shows that the LPAF-based ISAR imaging algorithm can obtain a well-focused image.
In this section, the synthetic aircraft target as shown in Fig. 6(a), which consists of 141 scatterers, is chosen to demonstrate the validity of the proposed ISAR imaging algorithm in this paper. We know that the CR and QCR of the azimuth echo in ISAR imaging are determined by the target’s size and the effective rotating vector. When the target’s size is a little larger or its maneuverability is a little higher, we have to adjust the parameter estimation range to meet the practical requirements. Here, we choose the special instance τ0=1/2 to enlarge the parameter estimation ranges and avoid spectrum aliasing and estimation error.
(a) Aircraft target model, (b) result of the conventional RD algorithm, (c) WVD of the extracted data, (d) simulation result with the proposed algorithm, (e) zoomed-in plot of ellipse area in (d), and (f) result after the Doppler spread compensation with LPAF.
Obviously, the image obtained by the traditional RD algorithm is blurred severely in Fig. 6(b). In order to confirm the validity of the proposed parameter estimation algorithm, when all scatterers have been corrected into right slant range bins after motion compensation, we extract the data of the 211st range bin to analyze in detail. Under the situation of SNRin=4 dB, the WVD as shown in Fig. 6(c) is computed for the radar echo signal of the 211st range bin. It is not difficult to find that the Doppler frequency is time-varying and nonlinear, which also demonstrates that the CPS model is very appropriate. Here we apply the LPAF-based parameter estimation algorithm to the extracted data. In Figs. 6(d) and 6(e), two peaks appear and can be easily picked out. Then the corresponding CR and QCR of two peaks can be correctly estimated via peak detection technique and CLEAN technique.
With the parameter estimation algorithm, we can estimate (ϕp,2,ϕp,3) and (ϕq,2,ϕq,3) as (4 Hz/s,8 Hz/s2) and (−4 Hz/s,−8 Hz/s2), respectively. After the Doppler spread compensation with the estimated parameters, we obtain the results in Fig. 6(f) by performing FFT, where the results of the conventional RD algorithm are also provided. Simulation shows that signal energy cannot be focused efficiently for Doppler frequency shift induced by CR and QCR. Consequently, the effective parameter estimation and compensation are quite essential.
Generally speaking, the better the focused quality of image, the smaller is the entropy.5,10,11,13,14 As before, the entropy is applied to measure the quality of the image. Compared with the LPWD-based imaging algorithm and the MCPF-based imaging algorithm, this proposed LPAF-based imaging algorithm can obtain robust performance against noise. Therefore, most scatterers have been reconstructed correctly in Fig. 7(a) and spurious scatterers in Fig. 7(a) are fewer than those in Figs. 7(b) and 7(c). Also, the entropy of image in Fig. 7(a) is smaller than those in other figures in Table 7, which indicates that the image quality for the LPAF-based ISAR imaging algorithm is better than those for other algorithms.
In order to further validate the effectiveness of the LPAF-based ISAR imaging algorithm for targets with complex motion, we provide the results of the LPAF-based ISAR imaging algorithm under a range of different SNRs in Fig. 8, where the scope of SNRin is from 2 to −8 dB with an increment of −2 dB. It is worthwhile to mention that the less spurious scatterers appear in Figs. 8(a)–8(c), and the entropies of Figs. 8(a)–8(c) are small than those of Figs. 8(d)–8(f) in Table 8. Especially, the images in Figs. 8(e) and 8(f) are blurred severely by spurious scatterers. Thus, we conclude that the proposed ISAR imaging algorithm fails when the input SNR falls below −4 dB.
Result of LPAF-based ISAR imaging algorithm under the situation of (a) SNRin=2 dB, (b) SNRin=0 dB, (c) SNRin=−2 dB, (d) SNRin=−4 dB, (e) SNRin=−6 dB, and (f) SNRin=−8 dB.
For ISAR imaging of targets with complex motion, the azimuth echo in a range bin can be modeled as multicomponent CPSs, which have been proved more appropriate than the LFM model. This paper presents a new ISAR imaging algorithm for targets with complex motion based on LPAF. LPAF can directly obtain CRQCRD only by FFT operation, where two parameters can be estimated simultaneously via peak detection. With the analyses and the simulations, we demonstrate that the LPAF outperforms MCPF and LPWD for the parameter estimation of CPS in the robustness of performance against noise, the computational cost, and the cross-term suppression. Finally, the effectiveness of the proposed ISAR imaging algorithm is validated by the experimental results of simulated data.
This work was supported in part by the National Natural Science Foundation of China under Grant Nos. 61271024 and 61201292 and in part by the Science and Technology Foundation of Shaanxi Province under Grant 2012JM8015. Qian Lv conceived the work in this paper that led to the submission, experiment design, and accomplishment of the writing of the manuscript. Jibin Zheng played an important role in revising the manuscript and providing English language support. Tao Su is mainly responsible for approving the final version.
Qian Lv received her BS degree in measuring and control techniques and instruments from Xi’an Shiyou University, Shaanxi, China, in 2013. She is currently working toward her PhD with the National Laboratory of Radar Signal Processing, Xidian University, Xi’an, China. Her research interests include synthetic aperture radar (SAR) and inverse SAR signal processing, time-frequency analysis, and interferometric inverse synthetic aperture radar (InISAR) image.
Tao Su received his BS degree in information theory, his MS degree in mobile communication, and his PhD in signal and information processing from Xidian University, Xi’an, China, in 1990, 1993, and 1999, respectively. Since 1993, he has been a professor with the National Laboratory of Radar Signal Processing, School of Electronic Engineering. His research interests include high-speed real-time signal processing on radar, sonar and telecommunications, digital signal processing, parallel processing system design, and FPGA IP design.
Jibin Zheng received his BS degree in electronic information science and technology from Shandong Normal University, Shandong, China, in 2009 and his PhD in signal and information processing from Xidian University, Xi’an, China, in 2015. From September 2012 to September 2014, he worked as a visiting PhD student at the Department of Electrical Engineering, Duke University, Durham, North Carolina. His research interests include SAR and inverse SAR signal processing, cognitive radar.
Qian Lv ; Tao Su and Jibin Zheng"Inverse synthetic aperture radar imaging of targets with complex motion based on the local polynomial ambiguity function", J. Appl. Remote Sens. 10(1), 015019 (Mar 18, 2016). ; http://dx.doi.org/10.1117/1.JRS.10.015019
Some tools below are only available to our subscribers or users with an online account.
Download citation file:
Customize your page view by dragging & repositioning the boxes below.
Shibboleth is an access management service that provides single sign-on to protected resources. It replaces the multiple user names and passwords necessary to access subscription-based content with a single user name and password that can be entered once per session. It operates independently of a user's location or IP address. If your institution uses Shibboleth authentication, please contact your site administrator to receive your user name and password.
Access to SPIE eBooks is limited to subscribing institutions and is not available
as part of a personal subscription. Print or electronic versions of individual SPIE
books may be purchased via SPIE.org.
Shibboleth is access management services that provide single sign-on to protected resources. It replace the multiple user names and passwords necessary to access subscription-based content with a single user name and password that can be entered once per session. It operates independently of a user's location or IP address. If your institution uses Shibboleth authentication, please contact your site administrator to receive your user name and password.