Open AccessTechnical Note

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

Zhongqiu Xu

^1,2,*

Bingchen Zhang

¹,

Guoru Zhou

^1,2,

Lihua Zhong

¹ and

Yirong Wu

Key Laboratory of Spatial Information Processing and Application System Technology, Chinese Academy of Sciences, Beijing 100190, China

University of Chinese Academy of Sciences, Beijing 100190, China

Author to whom correspondence should be addressed.

Remote Sens. 2021, 13(9), 1643; https://doi.org/10.3390/rs13091643

Submission received: 22 March 2021 / Revised: 11 April 2021 / Accepted: 20 April 2021 / Published: 22 April 2021

(This article belongs to the Section Remote Sensing Communications)

Download

Browse Figures

Graphical abstract
"> Figure 1
Schematic diagram of sparse SAR imaging model. "> Figure 2
Flow chart of the proposed algorithm combining variable splitting and modified ADMM. "> Figure 3
Various factors affecting SAR imaging results. "> Figure 4
Optical remote sensing image of experimental areas (provided by Google Earth). (a) Scene 1 corresponding to data 1. (b) Scene 2 corresponding to data 2. The size of scene 1 is 8673 m (vertical direction) × 8867 m (horizontal direction). The size of scene 2 is 19,228 m (vertical direction) × 4935 m (horizontal direction). "> Figure 5
Experimental results for data 1 under full sampling condition. (a) Chirp scaling. (b) CSA + multilook processing using a 5 by 5 window. (c) <math display="inline"><semantics> <msub> <mi>L</mi> <mn>1</mn> </msub> </semantics></math> regularization. (d) <math display="inline"><semantics> <msub> <mi>L</mi> <mn>1</mn> </msub> </semantics></math> and TV regularization. (e) Nonconvex and TV regularization. The vertical direction of the image is the azimuth direction, and the horizontal direction is the range direction. SAR data 1 occupies 1404 (azimuth) × 3941 (range) sampling points, and the size of the corresponding scene is 8673 m (azimuth) × 8867 m (range). All the images are plotted with the same color map. "> Figure 6
Experimental results for data 2 under full sampling condition. (a) Chirp scaling. (b) CSA + multilook processing using a 9 by 9 window. (c) <math display="inline"><semantics> <msub> <mi>L</mi> <mn>1</mn> </msub> </semantics></math> regularization. (d) <math display="inline"><semantics> <msub> <mi>L</mi> <mn>1</mn> </msub> </semantics></math> and TV regularization. (e) Nonconvex and TV regularization. The vertical direction of the image is the azimuth direction, and the horizontal direction is the range direction. SAR data 2 occupies 3115 (azimuth) × 4387 (range) sampling points, and the size of the corresponding scene is 19,228 m (azimuth) × 4935 m (range). All the images are plotted with the same color map. "> Figure 7
The slice along the azimuth direction of data 2. (a) The edge of the rainforest. (b) The simulated point target. "> Figure 8
Experimental results for data 1 under downsampling condition. (a) Chirp scaling. (b) CSA + Multilook processing using a 5 by 5 window. (c) <math display="inline"><semantics> <msub> <mi>L</mi> <mn>1</mn> </msub> </semantics></math> regularization. (d) <math display="inline"><semantics> <msub> <mi>L</mi> <mn>1</mn> </msub> </semantics></math> and TV regularization. (e) Nonconvex and TV regularization. The vertical direction of the image is the azimuth direction, and the horizontal direction is the range direction. SAR data 1 occupies 1404 (azimuth) × 3941 (range) sampling points, and the size of the corresponding scene is 8673 m (azimuth) × 8867 m (range). All the images are plotted with the same color map. "> Figure 9
Simulation results under different downsampling conditions. (a–c,g–i): the reconstruction results of the chirp scaling algorithm. (d–f,j–l): the reconstruction results of the proposed method. The top two rows: from left to right the ratios of azimuth downsampling are 80%, 70% and 60%. The bottom two rows: from left to right the ratios of azimuth downsampling are 50%, 40% and 30%. ">

Review Reports Versions Notes

Abstract

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 sacrificing the spatial resolution.

Keywords:

sparse SAR imaging; nonconvex and TV regularization; image quality assessment; radiometric property; spatial resolution

Graphical Abstract

1. 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 1-dimension (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 Section 4 and Section 5.

2. Materials and Methods

2.1. Sparse SAR Imaging Formation

In SAR imaging,

X \in C^{N_{P} \times N_{Q}}

denotes the complex-valued reflectivity matrix of the surveillance region and

Y \in C^{N_{η} \times N_{τ}}

represents the sampled 2-D echo data. Then

x = v e c (X) \in C^{N \times 1}

with

N = N_{P} \times N_{Q}

, and

y = v e c (Y) \in C^{M \times 1}

with

M = N_{η} \times N_{τ}

, where the vectorization operation

v e c (\cdot)

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]:

y = Φ α + n

(1)

where

n \in C^{M \times 1}

is the additive noise vector and

Φ \in C^{M \times 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:

m i n_{α} {∥y - Φ α∥}_{2}^{2} + λ p (α)

(2)

where

p (\cdot)

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 \in C^{N_{η} \times 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 \in R^{N_{η} \times N_{τ}}

as its own row elements and the range downsampling matrix

D_{r}

extracts some columns of the identity matrix

I \in R^{N_{η} \times 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:

m i n_{α} {∥y - Φ α∥}_{2}^{2} + λ {∥α∥}_{1}, {∥α∥}_{1} = \sum_{i = 1}^{N} (|α_{i}|)

(3)

where

α_{i}

is the

i

th 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.

2.2. 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]:

m i n_{α} {∥y - Φ α∥}_{2}^{2} + λ_{1} p_{N C} (α) + λ_{2} p_{T V} (|α|) .

(4)

In the above formula,

λ_{1}

and

λ_{2}

are regularization parameters and

p_{N C} (α)

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]:

p_{M C} (α) = \sum_{i = 1}^{N} p_{M C} (α_{i}), w i t h p_{M C} (α_{i}) = \{\begin{matrix} (|α_{i}| - \frac{α_{i}^{2}}{2 θ}), & |α_{i}| \leq θ \\ θ / 2, & |α_{i}| > θ \end{matrix} .

(5)

And

p_{T V} (|α|)

means the TV-norm penalty of image vector

α

. For two-dimensional scenes, the discrete isotropic TV-norm of

α

is defined as follows [11]:

p_{T V} (|α|) = \sum_{i, j} |\nabla ({|A|}_{i, j})| = {∥\nabla (|A|)∥}_{1}

(6)

where

A

is the two-dimensional complex-valued matrix corresponding to image vector

α

, the operator

|\cdot|

represents element-wise magnitude calculation and

\nabla ({|A|}_{i, j})

is the gradient vector of a pixel in the

i

th row and

j

th column, which is defined as:

\nabla ({|A|}_{i, j}) = (D_{h} {|A|}_{i, j}, D_{v} {|A|}_{i, j}), |\nabla ({|A|}_{i, j})| = \sqrt{{(D_{h} {|A|}_{i, j})}^{2} + {(D_{v} {|A|}_{i, j})}^{2}}

(7)

D_{h} {|A|}_{i, j} = |A [i + 1, j]| - |A [i, j]|,

(8)

D_{v} {|A|}_{i, j} = |A [i, j + 1]| - |A [i, j]| .

(9)

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.

2.3. 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]:

\begin{matrix} m i n_{α, z_{1}, z_{2}} {∥y - Φ α∥}_{2}^{2} + λ_{1} p_{M C} (z_{1}) + λ_{2} p_{T V} (|z_{2}|) \\ s . t . G α = z, w i t h G = {[I^{T}, I^{T}]}^{T}, z = {[z_{1}^{T}, z_{2}^{T}]}^{T} \end{matrix}

(10)

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:

\begin{matrix} s u b - p 1 : & α^{(t + 1)} = a r g m i n_{α} \frac{γ}{2} {∥G α - z^{(t)} - d^{(t)}∥}_{2}^{2} + {∥y - Φ α∥}_{2}^{2}, \\ s u b - p 2 : & z_{1}^{(t + 1)} = a r g m i n_{z} \frac{γ}{2} ({∥α^{(t + 1)} - d_{1}^{(t)} - z_{1}∥}_{2}^{2}) + λ_{1} p_{M C} (z_{1}), \\ s u b - p 3 : & z_{2}^{(t + 1)} = a r g m i n_{z} \frac{γ}{2} ({∥α^{(t + 1)} - d_{2}^{(t)} - z_{2}∥}_{2}^{2}) + λ_{2} p_{T V} (|z_{2}|), \\ d^{(t + 1)} = d^{(t)} - G α^{(t + 1)} + z^{(t + 1)}, w i t h d = {[d_{1}^{T}, d_{2}^{T}]}^{T} . \end{matrix}

(11)

In the above optimization problems,

d

is a variant of the Lagrange multiplier vector and

γ \geq 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:

α^{(t + 1)} = a r g m i n_{α} {∥[\begin{matrix} \sqrt{\frac{γ}{2}} G \\ Φ \end{matrix}] α - [\begin{matrix} \sqrt{\frac{γ}{2}} (z^{(t)} + d^{(t)}) \\ y \end{matrix}]∥}_{2}^{2} = a r g m i n_{α} {∥H α - β∥}_{2}^{2},

(12)

which leads to the following optimal solution according to the principle of gradient descent, and

\frac{1}{L}

is the step size of the iteration:

α^{(t + 1)} = α^{(t)} + \frac{1}{L} H^{H} (β - H α^{(t)}) = α^{(t)} + \frac{1}{L} [\frac{γ}{2} \sum_{i = 1}^{2} (d_{i}^{(t)} + z_{i}^{(t)}) - γ α^{(t)} + Φ^{H} (y - Φ α^{(t)})] .

(13)

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_{N C} (α)

satisfies the following three assumptions:

$p_{N C} (α)$ is a continuous function that is possibly nonsmooth and nonconvex, and can be rewritten as the difference of two convex functions;
$p_{N C} (α)$ can be written as $p_{N C} (α) = \sum_{i = 1}^{N} p_{N C} (α_{i})$ ;
the objective function is lower-bounded;

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

\begin{matrix} z_{1}^{(t + 1)} = Ψ_{M C} (α^{(t + 1)} - d_{1}^{(t)}, \frac{λ_{1}}{γ}), w i t h \\ Ψ_{M C} (α_{i}, λ) = \{\begin{matrix} 0, & |α_{i}| < λ \\ θ λ \frac{|α_{i}| - λ}{θ λ - λ} s i g n (α_{i}), & λ \leq |α_{i}| \leq θ λ \\ α_{i}, & |α_{i}| > θ λ . \end{matrix} \end{matrix}

(14)

Sub-problem 3 in (11) is the TV regularization problem and we can obtain its closed-form solution using the Chambolle algorithm:

z_{2}^{(t + 1)} = Ψ_{T V} (α^{(t + 1)} - d_{2}^{(t)}, \frac{λ_{2}}{γ}),

(15)

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

G e (\cdot)

and the imaging operator

I m (\cdot)

. 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_{1}^{(0)} = z_{2}^{(0)} = 0$ , $d_{1}^{(0)} = d_{2}^{(0)} = 0$ , the azimuth-range decouple operators $G e (α) ≅ Φ α$ and $I m (y) ≅ Φ^{H} y$ , error parameter $ϵ$ , maximum number of iterations $T_{m a x}$ .
1:: while $r e s > ϵ$ and $t < T_{m a x}$ do
2:: $α^{(t + 1)} = α^{(t)} + \frac{1}{L} [\frac{γ}{2} \sum_{i = 1}^{2} (d_{i}^{(t)} + z_{i}^{(t)}) - γ α^{(t)} + I m (y - G e (α^{(t)}))]$
3:: $z_{1}^{(t + 1)} = Ψ_{M C} (α^{(t + 1)} - d_{1}^{(t)}, λ_{1} / γ)$
4:: $z_{2}^{(t + 1)} = Ψ_{T V} (α^{(t + 1)} - d_{2}^{(t)}, λ_{2} / γ)$
5:: $d^{(t + 1)} = d^{(t)} - G α^{(t + 1)} + z^{(t + 1)}$
6:: $r e s = {∥α^{(t + 1)} - α^{(t)}∥}_{2} / {∥α^{(t)}∥}_{2}$
7:: $t = t + 1$
8:: end while
Output:: $\hat{α} = α^{(t)}$

2.4. 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]:

{\hat{α}}_{M F} = Φ^{H} y = Φ^{H} Φ \cdot α + n^{^{'}} = C_{2} \cdot (C_{1} \sqrt{σ} e^{j ϕ (α)}) + n^{^{'}} .

(16)

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:

{|C_{1}|}^{2} = \frac{P_{t} G^{2} λ^{2}}{{(4 π)}^{3} R^{4} L} .

(17)

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_{P R} \times G_{P A}

is the processing gain due to the range compression (

G_{P R}

) and azimuth compression (

G_{P A}

), 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:

{\hat{α}}_{S_{M}} = α = C_{1} \sqrt{σ} e^{j ϕ (α)},

(18)

while the azimuth-range decouple-operators-based sparse SAR imaging results have the following expression:

{\hat{α}}_{S_{O}} = Φ^{H} Φ \cdot α = C_{2} \cdot (C_{1} \sqrt{σ} e^{j ϕ (α)}) .

(19)

Therefore, we can think of the azimuth-range decouple-operators-based sparse SAR imaging results

{\hat{α}}_{S_{O}}

as the noiseless version of the matched-filtering-method-based SAR imaging results

{\hat{α}}_{M F}

. 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 (

{\hat{α}}_{S_{O}, T a r}

) should be almost the same as that of the matched filtering method (

{\hat{α}}_{M F, T a r}

). 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

{\hat{α}}_{S_{O}, T a r}

is to

{\hat{α}}_{M F, T a r}

, the higher the reconstruction accuracy of azimuth-range decouple-operators-based 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:

\begin{matrix} R B_{p} = |\hat{α} - α| / |α|, \\ R B_{d} = |\hat{α_{μ}} - α_{μ}| / |α_{μ}| \end{matrix}

(20)

where

R B_{p}

means the relative bias of point targets and

R B_{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

\hat{α}

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,

\hat{α_{μ}}

—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:

E N L = 0 . 5227^{2} \times \frac{μ {(A)}^{2}}{σ {(A)}^{2}}

(21)

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]:

\begin{matrix} γ (d B) = 10 \cdot l g (1 + \frac{1}{\sqrt{E N L}}) = 10 \cdot l g (1 + \frac{1}{0.5227} \times \frac{σ (A)}{μ (A)}) . \end{matrix}

(22)

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

s i n c (\cdot)

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

δ (\cdot)

, 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

{\hat{α}}_{S_{O}}

can be regarded as the noiseless version of the matched-filtering-method-based SAR imaging results

{\hat{α}}_{M F}

, which means the system impulse response of the ideal point target is approximately a

s i n c (\cdot)

function, so we can calculate MLW to measure the spatial resolution of

{\hat{α}}_{S_{O}}

3. Results

3.1. 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.

3.2. Data Processing and Image Quality Assessment

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

R B

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

R B

s 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.

Experiment 2

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.

4. 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:

μ = \sqrt{π σ_{0}} / 2

, the variance:

σ^{2} = \sqrt{(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%.

5. 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 sparse-signal-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.

Author Contributions

Conceptualization, Z.X., B.Z. and Y.W.; methodology, Z.X.; validation, Z.X., G.Z. and L.Z.; formal analysis, Z.X.; investigation, Z.X.; resources, L.Z.; data curation, Z.X.; writing—original draft preparation, Z.X.; writing—review and editing, Z.X.; visualization, Z.X.; supervision, Z.X. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China Grant No. 61991421.

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Acknowledgments

The authors would like to thank the China Center for Resources Satellite Data and Application for providing the GF-3 data.

Conflicts of Interest

The authors declare no conflict of interest.

References

Ender, J.H. On compressive sensing applied to radar. Signal Process. 2010, 90, 1402–1414. [Google Scholar] [CrossRef]
Zhang, B.; Hong, W.; Wu, Y. Sparse microwave imaging: Principles and applications. Sci. China Inf. Sci. 2012, 55, 1722–1754. [Google Scholar] [CrossRef] [Green Version]
Cetin, M.; Stojanovic, I.; Onhon, O.; Varshney, K.; Samadi, S.; Karl, W.C.; Willsky, A.S. Sparsity-driven synthetic aperture radar imaging: Reconstruction, autofocusing, moving targets, and compressed sensing. IEEE Signal Process. Mag. 2014, 31, 27–40. [Google Scholar] [CrossRef]
Ao, D.; Wang, R.; Hu, C.; Li, Y. A Sparse SAR Imaging Method Based on Multiple Measurement Vectors Model. Remote Sens. 2017, 9, 297. [Google Scholar] [CrossRef] [Green Version]
Bi, H.; Zhang, B.; Zhu, X.X.; Jiang, C.; Hong, W. Extended Chirp Scaling-Baseband Azimuth Scaling-Based Azimuth-Range Decouple L₁ Regularization for TOPS SAR Imaging via CAMP. IEEE Trans. Geosci. Remote Sens. 2017, 55, 3748–3763. [Google Scholar] [CrossRef]
Baraniuk, R.; Steeghs, P. Compressive radar imaging. In Proceedings of the 2007 IEEE Radar Conference, Waltham, MA, USA, 17–20 April 2007; pp. 128–133. [Google Scholar]
Alonso, M.T.; López-Dekker, P.; Mallorquí, J.J. A novel strategy for radar imaging based on compressive sensing. IEEE Trans. Geosci. Remote Sens. 2010, 48, 4285–4295. [Google Scholar] [CrossRef] [Green Version]
Patel, V.M.; Easley, G.R.; Healy, D.M., Jr.; Chellappa, R. Compressed synthetic aperture radar. IEEE J. Sel. Top. Signal Process. 2010, 4, 244–254. [Google Scholar] [CrossRef]
Osher, S.; Ruan, F.; Xiong, J.; Yao, Y.; Yin, W. Sparse recovery via differential inclusions. Appl. Comput. Harmon. Anal. 2016, 41, 436–469. [Google Scholar] [CrossRef] [Green Version]
Candes, E.; Tao, T. The Dantzig selector: Statistical estimation when p is much larger than n. Ann. Stat. 2007, 35, 2313–2351. [Google Scholar]
Güven, H.E.; Güngör, A.; Cetin, M. An augmented lagrangian method for complex-valued compressed sar imaging. IEEE Trans. Comput. Imaging 2016, 2, 235–250. [Google Scholar] [CrossRef]
Selesnick, I. Sparse regularization via convex analysis. IEEE Trans. Signal Process. 2017, 65, 4481–4494. [Google Scholar] [CrossRef]
Wei, Z.; Zhang, B.; Wu, Y. A SAR imaging method based on generalized minimax-concave penalty. Sci. China Inf. Sci. 2019, 62, 1–3. [Google Scholar] [CrossRef] [Green Version]
Gong, P.; Zhang, C.; Lu, Z.; Huang, J.; Ye, J. A general iterative shrinkage and thresholding algorithm for non-convex regularized optimization problems. In Proceedings of the International Conference on Machine Learning, Atlanta, GA, USA, 16–21 June 2013; pp. 37–45. [Google Scholar]
Uḡur, S.; Arıkan, O. SAR image reconstruction and autofocus by compressed sensing. Digit. Signal Process. 2012, 22, 923–932. [Google Scholar] [CrossRef]
Zhao, Y.; Liu, J.G.; Zhang, B.; Hong, W.; Wu, Y.R. Adaptive total variation regularization based SAR image despeckling and despeckling evaluation index. IEEE Trans. Geosci. Remote Sens. 2014, 53, 2765–2774. [Google Scholar] [CrossRef] [Green Version]
Kang, M.S.; Kim, K.T. Compressive sensing based SAR imaging and autofocus using improved Tikhonov regularization. IEEE Sens. J. 2019, 19, 5529–5540. [Google Scholar] [CrossRef]
Xu, Z.; Liu, M.; Zhou, G.; Wei, Z.; Zhang, B.; Wu, Y. An Accurate Sparse SAR Imaging Method for Enhancing Region-Based Features Via Nonconvex and TV Regularization. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2021, 14, 350–363. [Google Scholar] [CrossRef]
Zhang, Y.; Miao, W.; Lin, Z.; Gao, H.; Shi, S. Millimeter-Wave InSAR Image Reconstruction Approach by Total Variation Regularized Matrix Completion. Remote Sens. 2018, 10, 1053. [Google Scholar] [CrossRef] [Green Version]
Afonso, M.V.; Bioucas-Dias, J.M.; Figueiredo, M.A.T. Fast Image Recovery Using Variable Splitting and Constrained Optimization. IEEE Trans. Image Process. 2010, 19, 2345–2356. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Afonso, M.V.; Bioucas-Dias, J.M.; Figueiredo, M.A.T. An Augmented Lagrangian Approach to the Constrained Optimization Formulation of Imaging Inverse Problems. IEEE Trans. Image Process. 2011, 20, 681–695. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Chen, Q.; Li, Z.; Zhang, P.; Tao, H.; Zeng, J. A preliminary evaluation of the GaoFen-3 SAR radiation characteristics in land surface and compared with Radarsat-2 and Sentinel-1A. IEEE Geosci. Remote Sens. Lett. 2018, 15, 1040–1044. [Google Scholar] [CrossRef]
Mittermayer, J.; Younis, M.; Metzig, R.; Wollstadt, S.; Martínez, J.M.; Meta, A. TerraSAR-X system performance characterization and verification. IEEE Trans. Geosci. Remote Sens. 2009, 48, 660–676. [Google Scholar] [CrossRef]
Oliver, C.; Quegan, S. Understanding Synthetic Aperture Radar Images; SciTech Publishing: Raleigh, NC, USA, 2004. [Google Scholar]
Cozzolino, D.; Verdoliva, L.; Scarpa, G.; Poggi, G. Nonlocal CNN SAR Image Despeckling. Remote Sens. 2020, 12, 1006. [Google Scholar] [CrossRef] [Green Version]
Chambolle, A. An algorithm for total variation minimization and applications. J. Math. Imaging Vis. 2004, 20, 89–97. [Google Scholar]
Fang, J.; Xu, Z.; Zhang, B.; Hong, W.; Wu, Y. Fast compressed sensing SAR imaging based on approximated observation. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2013, 7, 352–363. [Google Scholar] [CrossRef] [Green Version]
Potter, L.C.; Ertin, E.; Parker, J.T.; Cetin, M. Sparsity and compressed sensing in radar imaging. Proc. IEEE 2010, 98, 1006–1020. [Google Scholar] [CrossRef]
Lee, J.S.; Jurkevich, L.; Dewaele, P.; Wambacq, P.; Oosterlinck, A. Speckle filtering of synthetic aperture radar images: A review. Remote Sens. Rev. 1994, 8, 313–340. [Google Scholar] [CrossRef]
Raney, R.K.; Runge, H.; Bamler, R.; Cumming, I.G.; Wong, F.H. Precision SAR processing using chirp scaling. IEEE Trans. Geosci. Remote Sens. 1994, 32, 786–799. [Google Scholar] [CrossRef]

Figure 1. Schematic diagram of sparse SAR imaging model.

Figure 2. Flow chart of the proposed algorithm combining variable splitting and modified ADMM.

Figure 3. Various factors affecting SAR imaging results.

Figure 4. Optical remote sensing image of experimental areas (provided by Google Earth). (a) Scene 1 corresponding to data 1. (b) Scene 2 corresponding to data 2. The size of scene 1 is 8673 m (vertical direction) × 8867 m (horizontal direction). The size of scene 2 is 19,228 m (vertical direction) × 4935 m (horizontal direction).

Figure 5. Experimental results for data 1 under full sampling condition. (a) Chirp scaling. (b) CSA + multilook processing using a 5 by 5 window. (c)

L_{1}

regularization. (d)

L_{1}

and TV regularization. (e) Nonconvex and TV regularization. The vertical direction of the image is the azimuth direction, and the horizontal direction is the range direction. SAR data 1 occupies 1404 (azimuth) × 3941 (range) sampling points, and the size of the corresponding scene is 8673 m (azimuth) × 8867 m (range). All the images are plotted with the same color map.

Figure 5. Experimental results for data 1 under full sampling condition. (a) Chirp scaling. (b) CSA + multilook processing using a 5 by 5 window. (c)

L_{1}

regularization. (d)

L_{1}

Figure 6. Experimental results for data 2 under full sampling condition. (a) Chirp scaling. (b) CSA + multilook processing using a 9 by 9 window. (c)

L_{1}

regularization. (d)

L_{1}

and TV regularization. (e) Nonconvex and TV regularization. The vertical direction of the image is the azimuth direction, and the horizontal direction is the range direction. SAR data 2 occupies 3115 (azimuth) × 4387 (range) sampling points, and the size of the corresponding scene is 19,228 m (azimuth) × 4935 m (range). All the images are plotted with the same color map.

Figure 6. Experimental results for data 2 under full sampling condition. (a) Chirp scaling. (b) CSA + multilook processing using a 9 by 9 window. (c)

L_{1}

regularization. (d)

L_{1}

and TV regularization. (e) Nonconvex and TV regularization. The vertical direction of the image is the azimuth direction, and the horizontal direction is the range direction. SAR data 2 occupies 3115 (azimuth) × 4387 (range) sampling points, and the size of the corresponding scene is 19,228 m (azimuth) × 4935 m (range). All the images are plotted with the same color map.

Figure 7. The slice along the azimuth direction of data 2. (a) The edge of the rainforest. (b) The simulated point target.

Figure 8. Experimental results for data 1 under downsampling condition. (a) Chirp scaling. (b) CSA + Multilook processing using a 5 by 5 window. (c)

L_{1}

regularization. (d)

L_{1}

Figure 8. Experimental results for data 1 under downsampling condition. (a) Chirp scaling. (b) CSA + Multilook processing using a 5 by 5 window. (c)

L_{1}

regularization. (d)

L_{1}

Figure 9. Simulation results under different downsampling conditions. (a–c,g–i): the reconstruction results of the chirp scaling algorithm. (d–f,j–l): the reconstruction results of the proposed method. The top two rows: from left to right the ratios of azimuth downsampling are 80%, 70% and 60%. The bottom two rows: from left to right the ratios of azimuth downsampling are 50%, 40% and 30%.

Table 1. The relative bias and relative relationship of four areas after relative normalization in experiment 1 for the five methods.

Method	A1		A2			A3			A4
Method	$μ_{A 1}$	${RB}_{A 1}$	$μ_{A 2}$	${RB}_{A 2}$	$μ_{A 2} / μ_{A 1}$	$μ_{A 3}$	${RB}_{A 3}$	$μ_{A 3} / μ_{A 1}$	$μ_{A 4}$	${RB}_{A 4}$	$μ_{A 4} / μ_{A 1}$
CSA	10.912	/	11.383	/	1.0432	11.664	/	1.0689	12.374	/	1.1340
CSA + ML	10.906	/	11.375	/	1.0430	11.671	/	1.0701	12.367	/	1.1340
$L_{1}$	8.2309	0.2457	8.6820	0.2373	1.0548	8.9710	0.2309	1.0899	9.6036	0.2239	1.1668
$L_{1}$ and TV	7.1261	0.3469	7.5889	0.3333	1.0649	7.8861	0.3239	1.1067	8.5145	0.3119	1.1948
MC and TV	10.557	0.0325	11.077	0.0269	1.0493	11.447	0.0186	1.0843	12.141	0.0188	1.1500

Table 2. The radiometric resolution of four distributed targets in experiment 2 for the five methods.

Method	$μ$	$σ^{2}$	ENL	$γ$ (dB)	Target
CSA CSA + ML $L_{1}$ $L_{1}$ and TV MC and TV	2.3668 2.3661 1.4625 1.4029 2.3665	1.7217 0.1104 1.5837 0.0255 0.0193	0.8889 13.8518 0.3690 21.080 79.340	3.1401 1.0335 4.2263 0.8558 0.4621	F1
CSA CSA + ML $L_{1}$ $L_{1}$ and TV MC and TV	2.2768 2.2763 1.3752 1.3124 2.2759	1.5575 0.0871 1.4180 0.0217 0.0170	0.9093 16.2510 0.3644 21.667 83.403	3.1148 0.9624 4.2434 0.8452 0.4513	F2
CSA CSA + ML $L_{1}$ $L_{1}$ and TV MC and TV	2.2863 2.2864 1.3844 1.3218 2.2853	1.5989 0.1027 1.4597 0.0244 0.0183	0.8932 13.9013 0.3587 19.548 78.162	3.1347 1.0319 4.2646 0.8855 0.4654	F3
CSA CSA + ML $L_{1}$ $L_{1}$ and TV MC and TV	2.2835 2.2834 1.3823 1.3192 2.2827	1.6114 0.1056 1.4701 0.0240 0.0183	0.8941 13.4943 0.3551 19.774 77.959	3.1461 1.0456 4.2782 0.8809 0.4660	F4

Table 3. The reconstruction accuracy and radiometric resolution of four distributed targets in downsampling experiments based on data 1 for the five methods.

Method	$μ$	${RB}_{d}$	$σ^{2}$	ENL	$γ$ (dB)	Target
CSA CSA + ML $L_{1}$ $L_{1}$ & TV MC & TV	9.8115 9.8021 6.9433 6.0151 9.1418	/ / 0.2923 0.3869 0.0683	27.237 1.8626 37.382 1.1548 1.4164	0.9657 14.094 0.3524 8.5602 16.120	3.0484 1.0256 4.2889 1.2768 0.9659	A1
CSA CSA + ML $L_{1}$ $L_{1}$ & TV MC & TV	10.073 10.068 7.3991 6.5702 9.7939	/ / 0.2655 0.3477 0.0277	29.682 2.5807 42.103 1.8127 1.7385	0.9339 10.730 0.3553 6.5064 15.075	3.0852 1.1570 4.2777 1.4365 0.9953	A2
CSA CSA + ML $L_{1}$ $L_{1}$ & TV MC & TV	10.340 10.334 7.8030 6.7917 9.9845	/ / 0.2454 0.3432 0.0344	29.943 2.1585 43.340 1.5232 1.3465	0.9755 13.516 0.3838 8.2740 20.227	3.0373 1.0449 4.1732 1.2958 0.8719	A3
CSA CSA + ML $L_{1}$ $L_{1}$ & TV MC & TV	10.726 10.725 8.1305 7.2285 10.557	/ / 0.2420 0.3261 0.0158	32.826 2.6600 48.634 1.8624 1.7436	0.9576 11.815 0.3714 7.6652 17.463	3.0576 1.1090 4.2176 1.3392 0.9318	A4

Table 4. The relative bias of the simulated distributed target under different downsampling ratios for the proposed method.

MC & TV	Downsampling Ratio
MC & TV	80%	70%	60%	50%	40%	30%
$R B_{d}$	0.0241	0.0276	0.0357	0.0704	0.1195	0.2402

Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

© 2021 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).

Share and Cite

MDPI and ACS Style

Xu, Z.; Zhang, B.; Zhou, G.; Zhong, L.; Wu, Y. Sparse SAR Imaging and Quantitative Evaluation Based on Nonconvex and TV Regularization. Remote Sens. 2021, 13, 1643. https://doi.org/10.3390/rs13091643

AMA Style

Xu Z, Zhang B, Zhou G, Zhong L, Wu Y. Sparse SAR Imaging and Quantitative Evaluation Based on Nonconvex and TV Regularization. Remote Sensing. 2021; 13(9):1643. https://doi.org/10.3390/rs13091643

Chicago/Turabian Style

Xu, Zhongqiu, Bingchen Zhang, Guoru Zhou, Lihua Zhong, and Yirong Wu. 2021. "Sparse SAR Imaging and Quantitative Evaluation Based on Nonconvex and TV Regularization" Remote Sensing 13, no. 9: 1643. https://doi.org/10.3390/rs13091643

APA Style

Xu, Z., Zhang, B., Zhou, G., Zhong, L., & Wu, Y. (2021). Sparse SAR Imaging and Quantitative Evaluation Based on Nonconvex and TV Regularization. Remote Sensing, 13(9), 1643. https://doi.org/10.3390/rs13091643

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Menu

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

Abstract

1. Introduction

2. Materials and Methods

2.1. Sparse SAR Imaging Formation

2.2. Nonconvex and TV Regularization

2.3. Variable Splitting and Modified ADMM

2.4. Sparse SAR Imaging Evaluation Index

3. Results

3.1. GF-3 SAR Data Description

3.2. Data Processing and Image Quality Assessment

4. Discussion

5. Conclusions

Author Contributions

Funding

Informed Consent Statement

Acknowledgments

Conflicts of Interest

References

Share and Cite

Article Metrics

Article Access Statistics

Further Information

Guidelines

MDPI Initiatives

Follow MDPI