Efficient Parameter Estimation for Sparse SAR Imaging Based on Complex Image and Azimuth-Range Decouple

Sparse signal processing theory has been applied to synthetic aperture radar (SAR) imaging. In compressive sensing (CS), the sparsity is usually considered as a known parameter. However, it is unknown practically. For many functions of CS, we need to know this parameter. Therefore, the estimation of sparsity is crucial for sparse SAR imaging. The sparsity is determined by the size of regularization parameter. Several methods have been presented for automatically estimating the regularization parameter, and have been applied to sparse SAR imaging. However, these methods are deduced based on an observation matrix, which will entail huge computational and memory costs. In this paper, to enhance the computational efficiency, an efficient adaptive parameter estimation method for sparse SAR imaging is proposed. The complex image-based sparse SAR imaging method only considers the threshold operation of the complex image, which can reduce the computational costs significantly. By utilizing this feature, the parameter is pre-estimated based on a complex image. In order to estimate the sparsity accurately, adaptive parameter estimation is then processed in the raw data domain, combining with the pre-estimated parameter and azimuth-range decouple operators. The proposed method can reduce the computational complexity from a quadratic square order to a linear logarithm order, which can be used in the large-scale scene. Simulated and Gaofen-3 SAR data processing results demonstrate the validity of the proposed method.


Introduction
Synthetic aperture radar (SAR) is an important imaging technology that has been applied in environmental protection and marine observation [1,2]. In recent years, the sparse signal processing method based on CS [3] has been implemented in microwave imaging [4,5]. It can recover the scene by solving L q (0 < q ≤ 1) regularization.
In [6], Çetin et al. proposed a sparsity-driven SAR imaging model for achieving autofocusing and moving targets imaging. Zhang et al. [7] explored the principles and applications in sparse microwave imaging. Patel et al. [8] analyzed different azimuth sampling methods based on the CS model. Luo et al. [9] developed a multiple scatterers detection method for SAR tomography with CS approach. Hossein et al. [10] proposed a polarimetric SAR estimator under the frame of CS. In [11], Zhu reviewed the CS-based super-resolving algorithm. Zhang et al. [12] proposed a novel where y ∈ C N×1 is the SAR echo data vector, N = N a (azimuth) × N r (range), x ∈ C N×1 is the backscattered coefficient vector, Φ is the observation matrix, and n ∈ C N×1 is the noise vector.
For the data formation model Equation (1), if the considered scene x is sparse enough and the observation matrix Φ satisfies the restricted isometry property (RIP) [24], x can be reconstructed by solving the L 1 optimization problem: where λ is the regularization parameter. There are many algorithms to solve (e.g., Equation (2)), such as the convex optimization algorithm [25], Bayesian learning algorithm [26], nonconvex optimization algorithm [27,28], and greedy algorithm [29]. After reconstruction,x is transferred to a matrixX Φ .

Azimuth-Range Decouple-Based Sparse SAR Imaging
The azimuth-range decouple-based sparse SAR imaging method is proposed in [7,15], the echo simulation operator G(·) is used to replace the observation matrix Φ, which is the inverse of the imaging operator I(·), that is, G(·) = I −1 (·) ≈ Φ. Then the azimuth-range decouple-based sparse SAR data formation model is represented as follows: E a ∈ C N a ×N a is an azimuth downsampling matrix, E r ∈ C N r ×N r is a range downsampling matrix. Both E a and E r are the binary matrices to denote the downsampling strategy, which are no longer identity matrices, thus reducing the number of measurements. Y ∈ C N a ×N r is the SAR raw data matrix, X ∈ C N a ×N r is the backscattered coefficient matrix, and N ∈ C N a ×N r is the noise matrix.
For this data formation model, the considered scene can be reconstructed by solving the L 1 optimization problem:X where · 2 is the 2-norm of a matrix.

Complex Image-Based Sparse SAR Imaging
A complex image-based sparse SAR imaging method is proposed in [22]. This method first establishes the imaging model with the complex image after MF recovery as the input, then represents the reconstruction of sparse scene as an L 1 optimization problem, and finally utilizes the iterative recovery algorithm to get the focused high-resolution SAR imagery. The signal model is represented as follows: where X MF ∈ C N a ×N r is the MF-reconstructed SAR complex image and N 0 ∈ C N a ×N r is the noise matrix. For this model, the considered scene can also be reconstructed by solving the L 1 optimization problem:X

Automatic Parameter Estimation Method
Several methods have been presented for automatically estimating the regularization parameter. We choose the generalized cross-validation (GCV) method [14,30] as the parameter estimation method, which can estimate λ by minimizing the following cost function without knowing the noise variance: where Y ∈ C N a ×N r is the SAR raw data matrix, N is the scene size, and tr(·) is the trace operator of a matrix. H λ is given in Equation (8): In (8), and β is a small positive constant.

Efficient Adaptive Parameter Estimation for Sparse SAR Imaging
In this section, the parameter estimation method based on azimuth-range decouple and the parameter estimation method based on complex image are introduced. Next, we introduce the proposed method in detail. Finally, the computational complexity of these methods is analyzed.

The Adaptive Parameter Estimation Method Based on Azimuth-Range Decouple
Combining the azimuth-range decouple operators with GCV, we can get the adaptive parameter estimation method for sparse SAR imaging. Compared with the adaptive parameter estimation method based on observation matrix, this method can reduce the computational complexity. Considering that M(X λ , β) is a large diagonal matrix of N × N, the computational cost of the trace of it is also large, we replace the trace operator tr(·) with the sum operator. Equation (7) can be rewritten as follows: which is the cost function of the adaptive parameter estimation method based on azimuth-range decouple, where G(·) is the echo simulation operator andX λ1 is shown in Equation (4). There are several algorithms to achieve the sparse reconstruction, such as iterative soft thresholding (IST) [31] and complex approximated message passing (CAMP) [32,33]. In this paper, we choose CAMP as sparse reconstruction algorithm, which has been applied to constant false-alarm rate (CFAR) detection in sparse SAR imaging [34].
The optimal regularization parameter is estimated by minimizing Equation (9). However, considering that finding the optimal regularization parameter requires the iterative processing, the total computational cost of the adaptive parameter estimation method based on azimuth-range decouple is still large.

The Adaptive Parameter Estimation Method Based on Complex Image
Compared with the azimuth-range decouple-based sparse SAR imaging method, the complex image-based sparse SAR imaging method only considers the threshold operation, which can further reduce the computational and memory costs. Combining it with GCV, we can get the adaptive parameter estimation method for sparse SAR imaging based on complex image. Equation (7) can be rewritten as follows: which is the cost function of the adaptive parameter estimation method based on complex image, whereX λ2 is shown in Equation (6).

The Proposed Method
The proposed method is mainly for the case of the downsampled data. On the one hand, although the adaptive parameter estimation method based on azimuth-range decouple can estimate the sparsity accurately, as mentioned above, the total computational cost of this method is large. On the other hand, due to the energy dispersion and ambiguities, the estimated sparsity of the parameter estimation method based on complex image will be greater than the true value, and we cannot simply use the parameter estimation method based on complex image to replace the parameter estimation method based on azimuth-range decouple. Therefore, we need to find a method to adaptively estimate the sparsity accurately while having the lower computational complexity. A good solution is to combine these two adaptive methods together, utilizing the complex image to pre-estimate the parameter and reduce the iteration range, then estimating the accurate parameter with raw data.
The proposed method has three steps. First, set the iteration range of sparsity to [K min , K max ] and adaptively estimate the sparsity based on the complex SAR image which is reconstructed by the downsampled raw data. The pre-estimated sparsity is set to K mid , which is greater than the true value due to ambiguities and energy dispersion caused by downsampling. Second, update the iteration range from [K min , K max ] to [K min , K mid ]. Third, get the adaptive reconstructed image and the optimal adaptive result of sparsity K opt on the new range [K min , K mid ] based on raw data.
The flowchart is shown in Figure 1.

The Proposed Method
The proposed method is mainly for the case of the downsampled data. On the one hand, although the adaptive parameter estimation method based on azimuth-range decouple can estimate the sparsity accurately, as mentioned above, the total computational cost of this method is large. On the other hand, due to the energy dispersion and ambiguities, the estimated sparsity of the parameter estimation method based on complex image will be greater than the true value, and we cannot simply use the parameter estimation method based on complex image to replace the parameter estimation method based on azimuth-range decouple. Therefore, we need to find a method to adaptively estimate the sparsity accurately while having the lower computational complexity. A good solution is to combine these two adaptive methods together, utilizing the complex image to pre-estimate the parameter and reduce the iteration range, then estimating the accurate parameter with raw data.
The proposed method has three steps. First, set the iteration range of sparsity to The flowchart is shown in Figure 1.   The details of adaptive parameter estimation based on azimuth-range decouple are shown in Algorithm 1, where [K min , K mid ] is the range of the sparsity; η λ,µ,CAMP (·) is the threshold function of CAMP.

Analysis of Computational Complexity
The computational complexity of different adaptive parameter estimation methods is analyzed in this section. A common characteristic of the adaptive parameter estimation methods mentioned above is that regularization parameter iterations are required. The difference lies in the different sparse reconstruction algorithms.
The measure of the computational complexity is the floating point operation (FLOP). Each FLOP represents a real addition operation or a real multiplication operation. In the observation matrix-based sparse SAR imaging method and azimuth-range decouple-based sparse SAR imaging method, the main calculation includes the imaging process, the echo simulation process, and the threshold process. The computational complexity of the threshold process is (8 n + n log 2 n) FLOPs, where n = N a × N r is the scene size. In the observation matrix-based sparse SAR imaging method, the imaging process and echo simulation process are two matrix multiplications. The main computational complexity of a single-step iteration of the observation-matrix-based sparse SAR imaging method is (16 mn + 8n + n log 2 n) FLOPs, where m is the sampling number. This computational complexity is approximately proportional to the quadratic square of the scene size.
In this paper, the chirp scaling [35] operator is chosen as the imaging operator. Therefore, I(·) and G(·) can be expressed as follows: where F a and F −1 a are the azimuth Fourier transform (FFT) operators and azimuth inverse Fourier transform (IFFT) operators, F r and F −1 r are the range FFT operators and range IFFT operators, Θ sc , Θ rc and Θ ac are three complex phase matrix. Chirp scaling and inverse chirp scaling both contain two FFTs, two IFFTs, and three time complex phase multiplications. According to [2], the computational complexity of FFT and IFFT with length l 0 is (5l 0 log 2 l 0 ). FLOPs, and the computational complexity of a complex multiplication operation is six FLOPs. Assuming that the data are sampled in the manner of uniform/nonuniform downsampling, the main computational complexity of a single-step iteration of the azimuth-range decouple-based sparse SAR imaging method is (46 n + 2m + 21n log 2 n) FLOPs, which is approximately proportional to the product of the linear logarithm of the scene size.
The complex image-based sparse SAR imaging method includes only threshold process. The computational complexity of a single-step iteration of this method is (8 n + n log 2 n) FLOPs, which is much lower than the azimuth-range decouple-based sparse SAR imaging method.
Let I represent the iteration steps of the recovery for sparse reconstruction algorithms. Let J and J 2 denote the number of iteration steps required for regularization parameter convergence when the iteration ranges of sparsity are [K min , K max ] and [K min , K mid ], respectively. Assuming that I = 20, J = 16, J 2 = J/4, the scene size n = 4096 × 4096, and the downsampling rate m/n = 80%, the computational complexity of different adaptive parameter estimation methods is shown in Table 1. Since the proposed method utilizes the complex image as the prior information to pre-estimate the parameter, the iteration range of the sparsity is reduced when the adaptive parameter estimation is processed in the raw data domain. Therefore, the proposed method has the lower computational complexity compared with the parameter estimation method based on azimuth-range decouple. For example, if the scene size is 4096 × 4096 and the downsampling rate is 80%, the proposed method can increase the computational efficiency about 3-4-fold.

Experiments
In this section, both simulation and real data experiments have been carried out to validate the effectiveness of the proposed method. The 1D simulation experiments compare the performance and reconstruction precision of the parameter estimation method based on observation matrix, parameter estimation method based on complex image and the proposed method. The 2D simulation experiments compare the adaptive result and computational complexity of different adaptive parameter estimation methods. Airborne SAR data and Gaofen-3 SAR data experiments are done to validate the ability of the proposed method to suppress energy dispersion and ambiguities. At last, the computational complexity of different adaptive parameter estimation methods is compared for different scene size.

1D Simulation
To validate the effectiveness of the proposed method, 1D simulation experiments are carried out. We set five point targets. Figure 2a shows the reconstructed images obtained by MF, parameter estimation method based on observation matrix, parameter estimation method based on complex image and the proposed method. In Figure 2a, the signal-to-noise ratio (SNR) is 15 dB and the downsampling rate is 80%. The adaptive λ opt of the adaptive parameter estimation method based on the observation matrix is 0.10 and the adaptive λ opt of the proposed method is 0.09. Due to downsampling, the adaptive λ of the adaptive parameter estimation method based on a complex image is 0.06, which is smaller than other two methods. From Figure 2a, we can conclude that the proposed method can effectively suppress the sidelobes and energy dispersion, and can recover the positions of target accurately compared with the positions of the ground truth. L 1 regularization is known as a biased estimator [36,37], and the bias would underestimate the intensities of the targets. Therefore, in Figure 2a, the target amplitude of the proposed method is lower than the ground truth.
image is 0.06, which is smaller than other two methods. From Figure 2a, we can conclude that the proposed method can effectively suppress the sidelobes and energy dispersion, and can recover the positions of target accurately compared with the positions of the ground truth. L1 regularization is known as a biased estimator [36,37], and the bias would underestimate the intensities of the targets. Therefore, in Figure 2a, the target amplitude of the proposed method is lower than the ground truth. In order to explore the accuracy of different adaptive parameter estimation methods, Figure 2b shows the relative mean square error (RMSE) curves of three methods at different SNR and downsampling rate, where the downsampling rate are 50% and 80%, respectively. It can be seen from Figure 2b that the proposed method has the similar sparse recovery performance as the adaptive method based on an observation matrix. The reconstruction precision of the adaptive method based on a complex image is worse than other two methods when the raw data are downsampled.

2D Simulation
In order to further analyze the effectiveness of the proposed method, 2D simulation experiments are carried out. The major simulation parameters are given in Table 2. The imaging results of nine point targets are shown in Figure 3. In Figure 3, the signal-to-noise ratio (SNR) is 20 dB and the downsampling rate is 80%. Figure 3a shows the image reconstructed by MF. Figure 3b  To better compare the reconstruction results of different methods, Figure 3d shows the azimuth profile of the 2D simulation experiment. Due to the bias of L1 regularization, in Figure 3d, the target amplitude of the proposed method is lower than the ground truth. In order to explore the accuracy of different adaptive parameter estimation methods, Figure 2b shows the relative mean square error (RMSE) curves of three methods at different SNR and downsampling rate, where the downsampling rate are 50% and 80%, respectively. It can be seen from Figure 2b that the proposed method has the similar sparse recovery performance as the adaptive method based on an observation matrix. The reconstruction precision of the adaptive method based on a complex image is worse than other two methods when the raw data are downsampled.

2D Simulation
In order to further analyze the effectiveness of the proposed method, 2D simulation experiments are carried out. The major simulation parameters are given in Table 2. The imaging results of nine point targets are shown in Figure 3. In Figure 3, the signal-to-noise ratio (SNR) is 20 dB and the downsampling rate is 80%. Figure 3a shows the image reconstructed by MF. Figure 3b shows the image reconstructed by adaptive parameter estimation method based on complex image, with the adaptive result λ being 0.17. From Figure 3b, we can see that the adaptive result based on complex image is not accurate when the raw data are downsampled, with the sidelobes still existing. Figure 3c shows the image reconstructed by the proposed method, with the adaptive result λ opt being 0.32. To better compare the reconstruction results of different methods, Figure 3d shows the azimuth profile of the 2D simulation experiment. Due to the bias of L 1 regularization, in Figure 3d, the target amplitude of the proposed method is lower than the ground truth.    Table 3 shows the adaptive λ and RMSE of the parameter estimation method based on azimuthrange decouple, the parameter estimation method based on complex image, and the proposed method, respectively. According to Table 4, the adaptive λ of the adaptive parameter estimation method based on complex image is smaller than other two methods when the raw data are downsampled, and varies with the SNR and downsampling rate. We can also conclude that the adaptive parameter estimation method based on azimuth-range decouple and the proposed method have almost the same sparse recovery performance. With the decrease in the downsampling rate, the RMSE of different adaptive parameter estimation methods increases. Therefore, the downsampling rate is crucial for the reconstruction accuracy of the adaptive parameter estimation methods. In this  2D simulation experiments with different SNR and downsampling rates are also carried out. Table 3 shows the adaptive λ and RMSE of the parameter estimation method based on azimuth-range decouple, the parameter estimation method based on complex image, and the proposed method, respectively. According to Table 4, the adaptive λ of the adaptive parameter estimation method based on complex image is smaller than other two methods when the raw data are downsampled, and varies with the SNR and downsampling rate. We can also conclude that the adaptive parameter estimation method based on azimuth-range decouple and the proposed method have almost the same sparse recovery performance. With the decrease in the downsampling rate, the RMSE of different adaptive parameter estimation methods increases. Therefore, the downsampling rate is crucial for the reconstruction accuracy of the adaptive parameter estimation methods. In this experiment, when the downsampling rate is 80% and the SNR is 25 dB, the proposed method has the smallest RMSE, which is the best result. Next, we will analyze the computational complexity. To illustrate that the proposed method has lower computational complexity, the computational complexity of different adaptive parameter estimation methods is compared for different scene size is represented in Figure 4. experiment, when the downsampling rate is 80% and the SNR is 25 dB, the proposed method has the smallest RMSE, which is the best result. Next, we will analyze the computational complexity. To illustrate that the proposed method has lower computational complexity, the computational complexity of different adaptive parameter estimation methods is compared for different scene size is represented in Figure 4. Figure 4 illustrates the computational complexity of three different adaptive parameter estimation methods for different scene size clearly. If the size of scene is over 1024 1024 × , the computational complexity of the adaptive parameter estimation method based on azimuth-range decouple increases dramatically. Although the computational complexity of the adaptive parameter estimation method based on complex image is the lowest, the adaptive result of this method is not accurate when the raw data are downsampled, as shown in Table 4. The proposed method utilizes complex image as prior information, thus having the lower computational complexity compared with the adaptive parameter estimation method based on azimuth-range decouple.   Figure 4 illustrates the computational complexity of three different adaptive parameter estimation methods for different scene size clearly. If the size of scene is over 1024 × 1024, the computational complexity of the adaptive parameter estimation method based on azimuth-range decouple increases dramatically. Although the computational complexity of the adaptive parameter estimation method based on complex image is the lowest, the adaptive result of this method is not accurate when the raw data are downsampled, as shown in Table 4. The proposed method utilizes complex image as prior information, thus having the lower computational complexity compared with the adaptive parameter estimation method based on azimuth-range decouple.

Airborne Data
The airborne SAR data processing results are shown in Figure 5. The raw data are 80% randomly downsampled, received by the C-band airborne SAR system of Institute of Electronics, Chinese Academy of Sciences. The accurate sparsity of this scene is 0.02.

Airborne Data
The airborne SAR data processing results are shown in Figure 5. The raw data are 80% randomly downsampled, received by the C-band airborne SAR system of Institute of Electronics, Chinese Academy of Sciences. The accurate sparsity of this scene is 0.02. In order to better evaluate the performance of different adaptive methods, the integrated sidelobe ratio (ISLR) is chosen to quantitatively measure the ability to suppress the energy dispersion [2]: total main 10 main where main P is the main-lobe power, total P is the total power. Figure 5a shows the image reconstructed by MF, with the obvious energy dispersion, and Figure  5d is the azimuth profile of the imaging result of MF, with the ISLR being −6.59 dB. Figure 5b shows the imaging result of the adaptive parameter estimation method based on complex image, with the adaptive result of sparsity mid K = 0.21. Figure 5e is the corresponding azimuth profile, with ISLR being −9.14 dB. From Figure 5b,e, we can see that when the raw data are downsampled, the adaptive parameter estimation method based on complex image cannot obtain an accurate result, with energy dispersion still existing. Figure 5c is the imaging result of the proposed method, with the adaptive result of sparsity opt K = 0.02, which converges to the accurate sparsity of the scene. Figure 5f is the azimuth profile of the imaging result of the proposed method, with ISLR being −10.55 dB. The proposed method can accurately estimate the sparsity and effectively suppress the noise and energy dispersion. In order to better evaluate the performance of different adaptive methods, the integrated sidelobe ratio (ISLR) is chosen to quantitatively measure the ability to suppress the energy dispersion [2]: where P main is the main-lobe power, P total is the total power. Figure 5a shows the image reconstructed by MF, with the obvious energy dispersion, and Figure 5d is the azimuth profile of the imaging result of MF, with the ISLR being −6.59 dB. Figure 5b shows the imaging result of the adaptive parameter estimation method based on complex image, with the adaptive result of sparsity K mid = 0.21. Figure 5e is the corresponding azimuth profile, with ISLR being −9.14 dB. From Figure 5b,e, we can see that when the raw data are downsampled, the adaptive parameter estimation method based on complex image cannot obtain an accurate result, with energy dispersion still existing. Figure 5c is the imaging result of the proposed method, with the adaptive result of sparsity K opt = 0.02, which converges to the accurate sparsity of the scene. Figure 5f is the azimuth profile of the imaging result of the proposed method, with ISLR being −10.55 dB. The proposed method can accurately estimate the sparsity and effectively suppress the noise and energy dispersion.

Gaofen-3 Data
The proposed method is also applicable to the spaceborne data. The Gaofen-3 satellite is a remote sensing satellite of China's high-resolution special project, which was launched in August 2016. It is the first C-band multipolarized SAR imaging satellite with a resolution of 1 m. Gaofen-3 data are processed to verify the background clutter and noise suppressing ability and ambiguity suppressing ability of the proposed method. In this experiment, we perform 80% random downsampling for the fully sampled data. The Gaofen-3 data processing results are shown in Figure 6.

Gaofen-3 Data
The proposed method is also applicable to the spaceborne data. The Gaofen-3 satellite is a remote sensing satellite of China's high-resolution special project, which was launched in August 2016. It is the first C-band multipolarized SAR imaging satellite with a resolution of 1 m. Gaofen-3 data are processed to verify the background clutter and noise suppressing ability and ambiguity suppressing ability of the proposed method. In this experiment, we perform 80% random downsampling for the fully sampled data. The Gaofen-3 data processing results are shown in Figure 6.  Figure 6a gives the MF imaging results of the downsampled raw data, with the obvious energy dispersion and azimuth ambiguities. Figure 6b shows the imaging result of the adaptive parameter estimation method based on azimuth-range decouple, with the adaptive result of sparsity opt K = 0.3514. It can be seen that this method can reconstruct the scene successfully and suppress the noise, energy dispersion and ambiguities efficiently. Figure 6c shows the imaging result of the adaptive parameter estimation method based on complex image, with the adaptive result of sparsity mid K =  Figure 6a gives the MF imaging results of the downsampled raw data, with the obvious energy dispersion and azimuth ambiguities. Figure 6b shows the imaging result of the adaptive parameter estimation method based on azimuth-range decouple, with the adaptive result of sparsity K opt = 0.3514. It can be seen that this method can reconstruct the scene successfully and suppress the noise, energy dispersion and ambiguities efficiently. Figure 6c shows the imaging result of the adaptive parameter estimation method based on complex image, with the adaptive result of sparsity K mid = 0.46. From Figure 6c, we can see that the adaptive result based on complex image is not accurate, with energy dispersion still existing. These two experimental results prove that the parameter estimation method based on azimuth-range decouple and the parameter estimation method based on complex image are not equivalent when the raw data are downsampled. However, we can use this pre-estimated sparsity as prior information to reduce the iteration ranges. Figure 6d is the imaging result of the proposed method, with the adaptive result of sparsity K opt = 0.3522, which is basically the same with the adaptive parameter estimation method based on azimuth-range decouple.
To further evaluate the ability of different adaptive methods to suppress the noise and ambiguity, target-to-background ratio (TBR) [38] and azimuth ambiguity-to-signal ratio (AASR) [23] are selected as two evaluation indicators. Their discrete expressions are defined as follows: TBR(X) = 20 log 10 where B is the background area, N B is the pixel number in B, and T is the target region. AASR = 10 log 10 where A is target region, N a is the pixels number in A, M a is the ambiguity area, and N m is the pixel number in M a . In this experiment, we chose five ships as performance test regions, as shown in the corresponding red frames. These five ships are represented as Ship 1-5, from left to right. Their corresponding azimuth ambiguity areas are shown in the blue frames.
The TBR of these five ships reconstructed by different methods are shown in Table 4. It can be seen from Table 5 that the proposed method can suppress the noise and energy dispersion effectively when the raw data are downsampled. Table 5. AASR of target regions based on different methods with downsampled data (80% downsampling).

Methods
Azimuth Ambiguity-to-Signal Ratio (dB) The AASR of these five ships reconstructed by different methods are shown in Table 5. From Table 5, we can see that the adaptive parameter estimation method based on complex image cannot suppress the azimuth ambiguity effectively. As a contrast, the adaptive parameter estimation method based on azimuth-range decouple and the proposed method both have the ability to decrease the azimuth ambiguity-to-signal ratio.
It can be seen from Tables 4 and 5 that the adaptive parameter estimation method based on azimuth-range decouple and the proposed method have almost the same sparse recovery performance. According to the previous analysis and the simulation experiments, the proposed method has the lower computational complexity, which can be used in the large-scale scene.

Conclusions
In this paper, an efficient adaptive parameter estimation method for sparse SAR imaging based on complex image and azimuth-range decouple is proposed. The proposed method combines the advantages of the azimuth-range decouple-based sparse SAR imaging and the complex image-based sparse SAR imaging method. In the proposed method, the parameter is pre-estimated based on the complex image. Adaptive parameter estimation is then processed in the raw data domain combining with the pre-estimated parameter and azimuth-range decouple operators. Compared with the adaptive parameter estimation method based on complex image, the proposed method can estimate the sparsity accurately when the raw data are downsampled. Compared with the adaptive parameter estimation method based on azimuth-range decouple, the proposed method has the lower computational complexity, which can be used in the large-scale scene. The simulation, airborne SAR data and Gaofen-3 SAR data experiment results demonstrate its validity.