Sparse SAR Imaging and Quantitative Evaluation Based on Nonconvex and TV Regularization

: Sparse signal processing has been used in synthetic aperture radar (SAR) imaging due to the maturity of compressed sensing theory. As a typical sparse reconstruction method, L 1 regularization generally causes bias effects as well as ignoring region-based features. Our team has proposed to linearly combine the nonconvex penalty and the total variation (TV)-norm penalty as a compound regularizer in the imaging model, called nonconvex and TV regularization, which can not only reduce the bias caused by L 1 regularization but also enhance point-based and region-based features. In this paper, we use the variable splitting scheme and modify the alternating direction method of multipliers (ADMM), generating a novel algorithm to solve the above optimization problem. Moreover, we analyze the radiometric properties of sparse-signal-processing-based SAR imaging results and introduce three indexes suitable for sparse SAR imaging for quantitative evaluation. In experiments, we process the Gaofen-3 (GF-3) data utilizing the proposed method, and quantitatively evaluate the reconstructed SAR image quality. Experimental results and image quality analysis verify the effectiveness of the proposed method in improving the reconstruction accuracy and the radiometric resolution without sacriﬁcing the spatial resolution.


Introduction
SAR is the major modern microwave imaging technology in remote sensing and is capable of producing high-resolution images of the Earth's surface. Compared with optical sensing technology, SAR is an active sensing technology and has all-time and all-weather observing ability. In recent years, with the development of compressed sensing, sparse signal processing has been widely used in SAR imaging, forming a novel strategy of SAR imaging named sparse SAR imaging. The sparse SAR imaging system can reconstruct a sparse signal with far fewer samples than that required by Nyquist-Shannon theory; in addition, it can also improve the reconstruction performance under full sampling conditions [1][2][3]. Researchers have applied sparse SAR imaging to a variety of SAR imaging modes, including stripmap SAR, ScanSAR, spotlight SAR and TOPS SAR modes [2,4,5].
L 1 regularization is a common reconstruction model in compressed sensing theory. The L 1 regularization-based SAR imaging model was first proposed in [6]. Then the 1dimension (1-D) and 2-dimension (2-D) L 1 regularization-based SAR imaging models were proposed in [7,8]. When implemented in sparse SAR imaging, L 1 regularization can suppress the sidelobes as well as improving the target background ratio. However, there exist two problems in L 1 regularization. Firstly, the L 1 -norm penalty is convex, and the analytic solution of the sparse signal reconstruction algorithm based on L 1 regularization is a biased estimation according to the related theory in statistical signal processing [9,10], which leads to the underestimation of targets' reflectivity, affecting the reconstruction accuracy. Secondly, L 1 regularization can only enhance point-based features in SAR images while ignoring region-based features, which means its reconstruction results will suffer from speckles due to the SAR imaging mechanism just the same as matched filtering [11]. In this paper, the region-based features refer to the characteristics of distributed targets, such as uniformity, continuity, edge preservation, etc.
The graphics and properties of nonconvex penalties are more similar to those of the L 0 -norm penalty compared with the L 1 -norm penalty, therefore, several works have used nonconvex penalties in the model, including the minimax concave (MC) penalty [12,13], the smoothly clipped absolute deviation (SCAD), the log sum penalty (LSP), and the L q -norm (0 < q < 1) penalty [14], to reduce the bias caused by L 1 regularization while ensuring the advantages of sparse signal processing. In addition to point-based features, region-based features are also important for many applications such as target classification and image segmentation, for which the TV-norm of the image amplitude is introduced into the SAR imaging model as a constraint. In [15], the authors propose a novel way of integrating both TV penalty and phase error into the cost function of the sparse SAR image reconstruction problem to improve the overall reconstruction quality. In [16], a TV regularization model for SAR image despeckling is introduced. Then [11] presents a solution to the complex SAR imaging problem within a constrained optimization formulation where the objective function includes a combination of the L 1 -norm and the TV of the magnitude of the complex-valued reflectivity field. In [17], the authors analyze the performance of the compressive sensing-TV method in SAR imaging and autofocus. The above research shows that TV regularization can result in reconstructions with emphasized piecewise-constant features. Therefore, in a recently submitted contribution, we have proposed to linearly combine the nonconvex penalty and the TV-norm penalty as a compound regularizer in the imaging model, called nonconvex and TV regularization, which can not only reduce the bias caused by L 1 regularization but also enhance point-based and region-based features [18,19].
In this paper, based on the existing research work, we further optimize the algorithm and analyze the reconstruction performance for the nonconvex and TV regularization model. Firstly, we use the variable splitting scheme and modify the alternating direction method of multipliers (ADMM) in the algorithm, generating a novel algorithm to solve the nonconvex and TV regularization problem [11,20]. The algorithm can avoid the intermediate minimization problems becoming increasingly ill-conditioned, which appears in the quadratic penalty approach used in the previous paper [21]. Then, to quantitatively evaluate the reconstructed SAR image quality based on nonconvex and TV regularization, this paper analyzes the radiometric properties of sparse-signal-processing-based SAR imaging results, and introduces three indexes suitable for sparse SAR imaging [22,23]. Finally, as the first C-band multi-polarization spaceborne SAR in China, GF-3 plays an important role in providing high-quality data for scientific experiments, therefore, it has aroused our great interest to process data from the GF-3 satellite utilizing the nonconvex and TV regularization method. Experimental results and image quality analysis verify the effectiveness and advantages of the method: compared with L 1 regularization, the method can not only improve the reconstruction accuracy but also enhance region-based features represented by the radiometric resolution; compared with the matched filtering method, the method can suppress speckles as well as sidelobes and additive noise; compared with the method of speckle removal in the image domain including multilook processing [24,25], the method can process the raw echo data and generate SAR images after speckle reduction, without sacrificing the spatial resolution.
The rest of the paper is organized as follows. Section 2 constructs the imaging model, introduces the nonconvex and TV regularization method, and proposes the variable splitting-modified ADMM and sparse SAR imaging evaluation indexes. Section 3 shows the experimental results and image quality analysis based on GF-3 data processing. The discussions and conclusions are presented in Sections 4 and 5.

Sparse SAR Imaging Formation
In SAR imaging, X ∈ C N P ×N Q denotes the complex-valued reflectivity matrix of the surveillance region and Y ∈ C N η ×N τ represents the sampled 2-D echo data. Then where the vectorization operation vec(·) stacks the columns one after the other. In this way, the relationship between the SAR measurements and the reflectivity of the scene to be reconstructed can be expressed as [2]: where n ∈ C M×1 is the additive noise vector and Φ ∈ C M×N is the corresponding system measurement matrix between y and α, as shown in Figure 1. The SAR imaging process can be considered as a linear inverse problem (LIP), in which we try to solve α with y and Φ. The most commonly used method for solving a LIP is to minimize the regularized linear least-squares cost function: where p(·) is the penalty function or regularizer and λ is the regularization parameter. If raw data y are the down-sampled data, the model in (1) turns out to be an under-determined linear system, which has either no solution or infinitely many solutions. Downsampling can be divided into azimuth downsampling and range downsampling. If Y ∈ C N η ×N τ represents the full sampled 2-D echo data, the azimuth down-sampled data can be written as D a Y and the range down-sampled data are Y D r , where the azimuth downsampling matrix D a extracts some rows of the identity matrix I ∈ R N η ×N τ as its own row elements and the range downsampling matrix D r extracts some columns of the identity matrix I ∈ R N η ×N τ as its own column elements. For the ill-conditioned inverse problem, when x is sparse enough and Φ satisfies the condition of restricted isometry properties, we can obtain the unique solution of (1) via the L 1 regularization-based sparse SAR imaging method: where α i is the ith element of α. However, there are two main problems in the above imaging model: firstly, the L 1 -norm penalty is convex and the analytic solution of the sparse signal reconstruction algorithm based on L 1 regularization is the biased estimation; secondly, L 1 regularization can only enhance point-based features in SAR images while ignoring region-based features, which means its reconstruction results will suffer from speckles due to the SAR imaging mechanism just the same as matched filtering. Therefore, it is necessary to find a suitable penalty function substituting for the L 1 -norm penalty.

Nonconvex and TV Regularization
Given the above two problems, we linearly combine the nonconvex penalty and the TV-norm penalty as a compound regularizer, generating the following nonconvex and TV regularization model [18]: In the above formula, λ 1 and λ 2 are regularization parameters and p NC (α) represents the nonconvex penalty of image vector α. In this paper, the minimax concave (MC) penalty is used as the representative of nonconvex penalties, as shown below [12,13]: And p TV (|α|) means the TV-norm penalty of image vector α. For two-dimensional scenes, the discrete isotropic TV-norm of α is defined as follows [11]: where A is the two-dimensional complex-valued matrix corresponding to image vector α, the operator |·| represents element-wise magnitude calculation and ∇(|A| i,j ) is the gradient vector of a pixel in the ith row and jth column, which is defined as: In the nonconvex and TV regularization-based sparse SAR imaging model, the nonconvex penalty, playing the role of sparsity inducing regularizer, can enhance point-based features as well as reducing the bias caused by L 1 regularization, and the TV-norm penalty will enhance region-based features, maintaining the continuity of the backscattering coefficient of distributed targets within a certain area. Therefore, the proposed method can improve the reconstruction accuracy of targets as well as protecting the reconstruction result from speckles, improving the radiometric resolution of SAR images.

Variable Splitting and Modified ADMM
For solving the nonconvex and TV regularization problem described in (4), we use the variable splitting scheme and modify the ADMM, generating a novel algorithm as shown in Figure 2. Firstly, under the guidance of the variable splitting scheme, we convert the unconstrained optimization problem (4) to the following equivalent constrained optimization problem [11]: where z 1 , z 2 are two auxiliary variables and I is the identity matrix. Then, according to the principle of ADMM, the above optimization problem can be solved via alternating minimization as below: In the above optimization problems, d is a variant of the Lagrange multiplier vector and γ ≥ 0 is called the augmented Lagrangian penalty parameter.
For ease of understanding, the objective function of sub-problem 1 in (11) can be written as a strictly convex quadratic function: which leads to the following optimal solution according to the principle of gradient descent, and 1 L is the step size of the iteration: Sub-problem 2 in (11) is a nonconvex optimization problem, which cannot be solved via the traditional ADMM. Therefore, we modify the ADMM by borrowing the idea of the general iterative shrinkage-thresholding algorithm, in which the author made some constraints on the nonconvex penalty so that the algorithm could solve nonconvex optimization problems [14]. According to the reference, when the nonconvex penalty p NC (α) satisfies the following three assumptions: • p NC (α) is a continuous function that is possibly nonsmooth and nonconvex, and can be rewritten as the difference of two convex functions; there exists a simple closed-form solution to sub-problem 2 in (11). The MC penalty used in this paper obviously meets the above conditions, and the closed-form solution corresponding to the MC penalty can be expressed as(θ > 1): Sub-problem 3 in (11) is the TV regularization problem and we can obtain its closedform solution using the Chambolle algorithm: for more details of the Chambolle algorithm, please refer to the related literature [26].
There still exists a problem in the above algorithm: it has to store and process the large-scale measurement matrix, which would bring huge memory and computational cost to SAR imaging. Thus, we adopt the azimuth-range decouple scheme in the real data processing [27], in which the measurement matrix Φ and its Hermitian transpose Φ H can be approximately replaced by the raw data generation operator Ge(·) and the imaging operator Im(·). Applying the azimuth-range decouple scheme to iterations, we can obtain the optimal solution of the problem (11) through the following iterative algorithm, as shown in Algorithm 1. In practical applications, the value of λ 1 needs to be determined according to the sparsity of the scene, and the value of λ 2 needs to be determined according to the uniformity of distributed targets in the scene.

Algorithm 1 Variable splitting and modified ADMM for nonconvex and total variation regularization
Input: The echo data y, initial sub-solution z the azimuthrange decouple operators Ge(α) ∼ = Φα and Im(y) ∼ = Φ H y, error parameter , maximum number of iterations T max . 1: while res > and t < T max do 2:

Sparse SAR Imaging Evaluation Index
• Physical Significance of Sparse SAR Imaging Results As shown in Equation (1) and Figure 1, the observation model of the SAR can be regarded as a linear system. Likewise, the imaging results of the matched filtering method can be written as [28]: In the above expression, the weighted complex reflectivity α not only represents the radar cross-section (RCS) of the target (σ), but also includes the energy weighting of the radar system (C 1 ), as shown in Figure 3: where P t is the peak transmitted power, G is the one-way power gain of the antenna pattern, R is the distance between radar and target, and L allows for system losses. In addition, C 2 = G PR × G PA is the processing gain due to the range compression (G PR ) and azimuth compression (G PA ), corresponding to Φ H Φ in Equation (16). n is additive noise, and it also contains ambiguity under the condition of downsampling.
The theoretical result of measurement-matrix-based sparse SAR imaging can be expressed as:α while the azimuth-range decouple-operators-based sparse SAR imaging results have the following expression:α Therefore, we can think of the azimuth-range decouple-operators-based sparse SAR imaging resultsα S O as the noiseless version of the matched-filtering-method-based SAR imaging resultsα MF . It should be noted that the above analysis is carried out for the point target. For distributed targets, we need to replace σ with σ = σ 0 ρ g ρ a in Equations (16)-(19), where ρ g denotes the ground resolution and ρ a denotes the azimuth resolution.

• Radiometric Accuracy
Through the above analysis, we can see that in the ideal case, for targets in the scene, the reconstruction results of azimuth-range decouple-operators-based sparse SAR imaging (α S O ,Tar ) should be almost the same as that of the matched filtering method (α MF,Tar ). However, some sparse SAR imaging methods cannot reconstruct targets accurately due to the bias effects caused by L 1 regularization. Here, we take the reconstruction result of the matched filtering method as the standard to measure the reconstruction accuracy of azimuth-range decouple-operatorsbased sparse SAR imaging. The more similarα S O ,Tar is toα MF,Tar , the higher the reconstruction accuracy of azimuth-range decouple-operatorsbased sparse SAR imaging is, and the more accurate the estimation of RCS is, which means higher radiometric accuracy. To quantitatively analyze the reconstruction accuracy, we define relative bias (RB) as the evaluation criteria: where RB p means the relative bias of point targets and RB d is the relative bias of distributed targets. Correspondingly, α is the imaging result of the point target based on the matched filtering method, which is regarded as the ground truth here, andα is the estimation of α using the azimuth-range decouple-operators-based sparse SAR imaging method. Because the noise of coherent speckle exists in the SAR image, we use the mean value of the reconstruction results (α µ -matched filtering,α µ -sparse signal processing) in the definition of relative bias for distributed targets here.

• Radiometric Resolution
In addition to impulse response properties, image quality assessment must consider the response of the system to distributed targets. Generally, the equivalent number of looks (ENL) and radiometric resolution (RR) are used to measure the reconstruction quality of distributed targets. The definition of ENL is derived from multilook processing, which can be operated based on an intensity or amplitude SAR image. According to the relevant reference [29], amplitude-averaging-based ENL has the following expression: where µ(A) is the mean of amplitude SAR images and σ(A) is the standard deviation. Furthermore, we can get the definition of radiometric resolution as below [22,23]: ).
As nonconvex and TV-regularization-based sparse SAR imaging generates the amplitude SAR image, which is the same as operating multilook processing based on the amplitude SAR image, it is reasonable for us to use the above definition to evaluate the reconstruction quality of distributed targets.

• Spatial Resolution
Main lobe width (MLW), also called the spatial resolution, is measured by the distance between the points with intensities 3 dB below the maximum intensity of the main lobe peak in the azimuth and range direction, standing for the distance between two objects on the ground to the extent that they can be distinct or separate. In the radar image obtained by matched filtering, the system impulse response of the ideal point target is the sinc(·) function, so we can use MLW to measure the spatial resolution. However, the ideal reconstruction result of the point target in measurement-matrix-based sparse microwave imaging is the impulse function δ(·), so it is no longer suitable to calculate MLW to measure spatial resolution [2]. Fortunately, as shown in Equation (19), the azimuth-range decouple-operators-based sparse SAR imaging resultsα S O can be regarded as the noiseless version of the matched-filtering-method-based SAR imaging resultsα MF , which means the system impulse response of the ideal point target is approximately a sinc(·) function, so we can calculate MLW to measure the spatial resolution ofα S O .

GF-3 SAR Data Description
GF-3 is the first C-band remote sensing satellite with quad-polarization and multiangle capability in the civilian field of China, and operates in 12 imaging modes with spatial resolution up to 1 m. In this paper, we select two pieces of stripmap mode GF-3 SAR raw echo data for experiments. Data 1: To verify the ability of nonconvex and TV regularization in improving reconstruction accuracy, we choose an area that includes many farmland regions with different backscattering coefficients as shown in Figure 4a. Considering that the backscattering coefficient of the farmland region is relatively uniform, we can think of them as distributed targets. In SAR images, the gray scale of farmland regions is different, which indicates the difference in backscattering coefficients between them. This scene can be used to analyze the reconstruction accuracy of distributed targets with different backscattering coefficients, and to study the maintenance of the relative relationship between them.

Data 2:
The second piece of data was collected in the Amazon rainforest area, which has a quite homogeneous backscattering coefficient around a wide area, as shown in Figure 4b. The rainforest area in the scene is divided into several pieces by the river, but they still have almost the same backscattering coefficient; the scene also shows a certain degree of sparsity. Therefore, we will use it to explore the performance of nonconvex and TV regularization in improving radiometric resolution. In addition, we have added a piece of simulated echo data of a point target to the raw SAR data, for verifying the ability of imaging methods to maintain point-based features when enhancing region-based features.

• Experiment 1
Firstly, we utilize the chirp scaling algorithm (CSA) to process the data 1 [30]. As shown in Figure 5a, there are obvious speckles in farmland regions. The simplest approach to despeckling is to average the amplitude over several pixels within a window centered on a specific pixel, known as multilook (ML) processing [24]. The result of applying this method to the original image in Figure 5a, using a 5 by 5 window, is shown in Figure 5b. Secondly, L 1 regularization is used to process the data, with results shown in Figure 5c. Compared with Figure 5a,b, the figure has less noise and clutter. However, as L 1 regularization-based sparse SAR imaging methods just aim at enhancing point-based features, their reconstruction results also suffer from speckles due to the SAR imaging mechanism just the same as matched filtering. Finally, we reconstruct the experimental scene via L 1 and TV regularization and nonconvex and TV regularization. As shown in Figure 5d,e, the two methods can not only suppress noise and clutter but also maintain the uniformity and continuity of distributed targets. All the images in Figure 5 are plotted with the same color map. For quantitative analysis of reconstruction accuracy, we select four areas as experimental distributed targets, which are surrounded by the red rectangle and labeled with A1/A2/A3/A4 in Figure 5a-e, and calculate the relative bias (RB) of the four targets according to Equation (14). Table 1 shows that the average reconstruction RB of nonconvex and TV regularization is 2.42% when we regard the result of the matched filtering method as the ground truth, while the average reconstruction RBs of L 1 regularization and L 1 and TV regularization are 23.4% and 32.9% respectively. The statistical results show that nonconvex and TV regularization has higher reconstruction accuracy than L 1 regularization and L 1 and TV regularization. Furthermore, we also calculate the ratio of the mean amplitude of target A2, A3, A4 to target A1; the results in Table 1 can demonstrate that nonconvex and TV regularization can better maintain the relative relationship between different targets compared with L 1 regularization and L 1 and TV regularization. In this part, we process data 2 using the same five methods as above: CSA, CSA + multilook processing using a 9 by 9 window, L 1 regularization, L 1 and TV regularization, and the proposed nonconvex and TV regularization. As shown in Figure 6, the latter three sparse SAR imaging methods can obviously suppress clutter and noise, meanwhile, nonconvex and TV regularization and L 1 and TV regularization can also maintain the uniformity and continuity of the backscattering coefficient in rainforest areas. All the images in Figure 6 are plotted with the same color map. To evaluate the radiometric resolution of the reconstructed SAR image, as in experiment 1, we also select four areas as experimental distributed targets, which are surrounded by the red rectangle and labeled with F1/F2/F3/F4 in Figure 6. In addition, there is a simulated point target in the scene, surrounded by the red dotted circle and labeled with P1, for assessment of the spatial resolution. Then we calculate the mean and variance of the four distributed targets for getting ENL and radiometric resolution. As shown in Table 2, the average radiometric resolution (γ(dB)) of reconstruction results based on the matched filtering method, L 1 regularization and nonconvex and TV regularization are 3.13 dB, 4.25 dB and 0.46 dB, respectively. The lower the value of γ, the higher the radiometric resolution, therefore, it can be further calculated that nonconvex and TV regularization can improve the radiometric resolution by an average of 85.28%, 89.19% compared with the matched filtering method and L 1 regularization, which verifies the effectiveness of nonconvex and TV regularization in improving the radiometric resolution on the one hand, and that the bias effect caused by L 1 regularization has a certain degree of adverse impact on the radiometric resolution on the other hand. Finally, we slice the edge of the rainforest marked with a blue line in Figure 6 and the simulated point along the azimuth direction to compare the spatial resolution between the multilook processing and nonconvex and TV regularization. Figure 7a shows that the edge in the imaging results of nonconvex and TV regularization is sharper than that of the multilook processing, which means that nonconvex and TV regularization has better edge-preserving ability than the multilook processing, and further proves the ability of nonconvex and TV regularization to maintain the spatial resolution. In Figure 7b, it is obvious that nonconvex and TV regularization has nearly the same spatial resolution as chirp scaling while the multilook processing widens the main lobe, which also proves that nonconvex and TV regularization can maintain the same spatial resolution as the matched filtering method.

Discussion
In this section, we will explore the effectiveness of nonconvex and TV regularization under the condition of downsampling. We first perform 80% random downsampling in the azimuth direction for the fully sampled GF-3 echo data 1, and the scene corresponding to data 1 is shown in Figure 4a. Figure 8 shows the reconstructed images from downsampled echo data by CSA, CSA + multilook processing using a 5 by 5 window, L 1 regularization, L 1 and TV regularization and the proposed nonconvex and TV regularization, respectively. All the images in Figure 8 are plotted with the same color map. It can be seen that there are obvious ambiguities and energy dispersion along the azimuth direction due to the data downsampling, especially for targets with strong scattering intensity such as the target surrounded by a blue dotted rectangle in Figure 8. However, the latter three sparse SAR imaging methods reconstruct the considered scene successfully and simultaneously obtain an image with lower noise and clutter. Besides, it is obvious that nonconvex and TV regularization and L 1 and TV regularization can also maintain the uniformity and continuity of the backscattering coefficient in farmland regions. To quantitatively analyze the reconstruction accuracy and radiometric resolution, we select four regions as experimental distributed targets, which are the same as those in Experiment 1, and calculate the mean value, relative bias, variance, ENL, and radiometric resolution of the four targets. Statistical results in Table 3 show that under the condition of downsampling, nonconvex and TV regularization can also reconstruct targets at almost the same amplitude as the matched filtering method, as well as improve the radiometric resolution. Finally, we design a set of simulation experiments to verify the reconstruction performance of the proposed method under different downsampling ratios. The size of the simulated scene is 1024 × 1024 pixels. We place a distributed target that occupies 101 × 101 pixels at the center of the scene. According to Oliver and Quegan [24], we set an equivalent phase center in each pixel cell of the target, whose amplitudes are i.i.d. Rayleigh distributed (the mean value: µ = √ πσ 0 /2, the variance: σ 2 = (1 − π/4)σ 0 , σ 0 : the backscattering coefficient) and phases are i.i.d. uniformly distributed U[−π, +π). In this experiment, the SNR is set to 20 dB, the range of the azimuth downsampling ratio is from 80% to 30%, and the interval is 10%. Under each downsampling ratio condition, we use the chirp scaling algorithm and the proposed method to reconstruct the target scene, and the imaging results are shown in Figure 9. Based on the above experimental results, we calculate the relative bias of the reconstruction results of the proposed method, as shown in Table 4, to quantitatively analyze the influence of the downsampling ratio on imaging performance. According to the calculation results in Table 4 and the imaging results in Figure 9, we can find the following: first, the amount of missing data will cause degradation of image quality, as the amount of missing data increases (the downsampling ratio decreases), the reconstruction performance becomes worse and the reconstruction bias of the proposed method gradually increases; second, as the amount of missing data increases, the probability of successful reconstruction gradually decreases, and under the conditions of this experiment, the proposed method will not be able to reconstruct the target scene successfully when the downsampling ratio is less than 40%.

Conclusions
In this paper, we linearly combined the nonconvex penalty and the TV-norm penalty as a compound regularizer in the imaging model, called nonconvex and TV regularization, and proposed a novel algorithm combining the variable splitting scheme and the modified ADMM to solve it. To quantitatively evaluate the reconstructed SAR image quality based on nonconvex and TV regularization, we analyzed the radiometric property of sparsesignal-processing-based SAR imaging results and introduced three indexes suitable for sparse SAR imaging. Experimental results and image quality analysis based on GF-3 data verified the effectiveness of the proposed method in improving the reconstruction accuracy and the radiometric resolution without sacrificing the spatial resolution, which is beneficial to the quantitative use of SAR data, such as RCS estimation, target classification, and image segmentation. Based on the content of this paper, our next research work includes two aspects. The first is to optimize the TV-norm constraints (such as weighting processing) to improve the reconstruction performance of the method proposed in this paper. The second is to extend the method to the data processing of 3-D SAR imaging such as downward-looking array 3-D SAR and tomoSAR to improve the performance of 3-D reconstruction.