## 1.

## Introduction

Due to advantages over optical sensing tools, such as robust performance under all-time and all-weather circumstances, radar plays an important role in such realms as earth observation, oceanic monitoring, and military reconnaissance.^{1}^{,}^{2} To realize these applications, two-dimensional high resolution is usually required. The high range resolution can be obtained by transmitting the wideband signal and using the pulse compression technique. The high resolution in azimuth/angular can be realized using synthetical aperture technique or Doppler beam sharpening approaches. However, these techniques are unable to realize radar imaging in the forward-looking area of the radar platform. The reason is that the Doppler resolution and the range resolution in the forward-looking area are parallel.^{3}

Scanning radar works as a noncoherent sensor, which is suitable for any geometrical situation and is able to realize two-dimensional radar imaging in the forward-looking area of the radar platform. In the same way, the pulse compression provides the high range resolution, while the angular resolution is usually poor when compared with the range resolution. For scanning radar, the upper limit of angular resolution is determined by the effective wavelength and the size of the antenna. Moreover, the angular resolution suffers and degrades as a function of target distance. One method to improve the angular resolution in scanning radar image can be accomplished by increasing the physical size of the antenna aperture. However, there may be insufficient space to accommodate such a large size antenna in the application. Therefore, an improvement in the angular resolution of the scanning radar with fixed parameters using signal processing methods is considered.

According to Refs. 4 and 5, the echo of scanning radar can be modeled as the convolution of the transmitted signal with the scattering coefficients of the targets. Therefore, in theory, a deconvolution method can enhance the angular resolution of the scanning radar.^{6} Deconvolution aims at removing the transfer function, system noise, and other nonidealities of the imaging process, thereby improving the imaging resolution.^{7} The deconvolution approach has been widely used for video restoration,^{8} image processing,^{9} sparse signal restoration,^{10} and through-wall radar imaging,^{11} while less work on the deconvolution approach for angular super-resolution imaging in scanning radar has been reported.

This paper presents a deconvolution method for improving the angular resolution of scanning radar. We first convert the angular super-resolution imaging task in scanning radar into an equivalent deconvolution problem. It is well-known that the deconvolution problem is ill-posed, which means the solution to the deconvolution problem is unstable and sensitive to the noise. To stabilize the solution of the deconvolution problem, we must utilize some prior information of the solution to the deconvolution problem. This paper focuses on some applications such as the detection of detecting and city area imaging, and the case where a small number of dominant scatterers are distributed in the scene. In these applications, large coefficients of the scatterers can capture most of the information of the scene.^{12} Inspired by the Bayesian compressive sensing,^{13} the statistic characteristics of the dominant scatterers in the illuminated scene can be represented by ${l}_{1}$-norm. Then the deconvolution technique for dominant scatterers in angular super-resolution imaging can be converted into the corresponding ${l}_{1}$ optimization problem. Therefore, the ${l}_{1}$ regularization is adopted and the deconvolution problem can be written with a constrained optimization form. Recently, applications of ${l}_{1}$-norm minimization can be found in Ref. 14 for ISAR imaging,^{15} for feature enhanced, and for through-wall radar image.^{16} In this paper, we are focused on solving the deconvolution problem, which corresponds to the task of angular super-resolution imaging in the forward-looking area of radar platform. Finally, the constrained optimization problem is solved using an augmented Lagrangian method. The rationale behind the augmented Lagrangian method is to find a saddle point of the associated augmented Lagrangian function, which combines the function to be optimized with a Lagrange multiplier term and a penalty term for the constraints.^{17} This leads to an alternating direction method for finding the saddle point of the augmented Lagrangian function. Compared with the traditional Lagrangian multiplier method, the augmented Lagrangian method overcomes numerical difficulties due to the penalty going to infinity.^{18} Based on simulation and experimental results, we demonstrate that the proposed deconvolution algorithm has better angular super-resolution performances.

The rest of this paper is organized as follows: Sec. 2 presents the signal model for scanning radar and converts the angular super-resolution problem into an equivalent deconvolution problem. Section 3 presents an augmented Lagrangian method to solve the corresponding deconvolution problem. In Sec. 4, simulation and experimental results are given to demonstrate the validity of the proposed method. Finally, the conclusions are given in Sec. 5.

## 2.

## Signal Model for Scanning Radar in Azimuth

In this section, we introduce the signal model for scanning radar and convert the angular super-resolution problem into an equivalent deconvolution problem. The details of the deconvolution algorithm for the angular super-resolution problem will be subsequently provided.

Suppose that the radar works in the scanning mode. The geometrical model of scanning radar is illustrated in Fig. 1. The radar platform at the altitude $H$ is moving along the $Y$-axis corresponding to the range direction, while the antenna scans the scene along the $X$-axis corresponding to the azimuth direction with constant angular velocity $\omega $ and then receives the echo data from the observed scene. Assume that the target in the azimuth can be considered as an ideal point. For an isolated target scatter, the received signal at the receiver output is proportional to the antenna pattern. The left part of Fig. 1 shows that the scanning radar collects echo data, while it scans in the azimuth across the target field. It is noted that the adjacent echoes mostly cover the same target features on the observed scene but with different contributions to the overall signal, which is shown in the bottom part of Fig. 1.

At first, the conventional range compression and range migration correction are applied to the echo data with the current approaches.^{2}^{,}^{19} The range compression and range migration to correct data are denoted by $g(\theta ,\phi ,r)$, where $r$ is the slant range from the radar to the target and $\theta $, $\phi $ denote the angle between the direction of the antenna to the target and flight direction and the incident angle of the beam, respectively, which are smaller than 10 deg. On the other hand, it is assumed that the antenna pattern is of relative isotropy to the imaging scene. Then the received signal of the scanning radar along the range profile can be recognized as a convolution kernel comprising the antenna power pattern in the angle coordinates and the pulse modulation function in the range dimension,^{20} which could be expressed as

Without loss of generality, we consider the variation in scattering with the azimuth angle $\theta $ for a fixed range ${r}_{0}$ and a fixed angle $\phi $. Note that with the limited range resolution, the range averaging scattering would be exactly equal to the effective scattering evaluated at the range of the interested ${r}_{0}$. Thus, we only consider the azimuth variation in Eq. (1) that can be written as

## (2)

$$g({\phi}_{0},\theta ,{r}_{0})=f({\phi}_{0},\theta ,{r}_{0})*[a({\phi}_{0},\theta )x(\frac{2{r}_{0}}{c}\left)\right].$$For mathematical simplicity, we only consider the azimuth variation. Then we can rewrite Eq. (2) as

In addition, the presence of the noise contaminates the data $g(\theta )$. The formation of the data $g(\theta )$ can be described as

where $n(\theta )$ represents the noise.In the case of observing a certain scene $\mathrm{\Omega}$, the echo returned to the radar is the whole scene. The reflectivity of $\mathrm{\Omega}$ can be denoted as a two-dimensional matrix

## (5)

$$\mathrm{\Omega}={\left[\begin{array}{cccc}f({\phi}_{1},{\theta}_{1}),& f({\phi}_{1},{\theta}_{2}),& \cdots & f({\phi}_{1},{\theta}_{A})\\ f({\phi}_{2},{\theta}_{1}),& f({\phi}_{2},{\theta}_{2}),& \cdots & f({\phi}_{2},{\theta}_{A})\\ \vdots & \hspace{1em}& \ddots & \vdots \\ f({\phi}_{R},{\theta}_{1}),& f({\phi}_{R},{\theta}_{2}),& \cdots & f({\phi}_{R},{\theta}_{A})\end{array}\right]}_{R\times A},$$To attain the corresponding discrete convolution equation, we represent the scene $\mathrm{\Omega}$ as a vector by stacking the columns of the matrix:

## (6)

$$f=\mathbf{vec}(\mathrm{\Omega})\phantom{\rule{0ex}{0ex}}={[\begin{array}{cccccccccc}f({\phi}_{1},{\theta}_{1})& \dots & f({\phi}_{1},{\theta}_{A}),& f({\phi}_{2},{\theta}_{1})& \dots & f({\phi}_{2},{\theta}_{A}),& \cdots ,& f({\phi}_{R},{\theta}_{1})& \dots & f({\phi}_{R},{\theta}_{A})\end{array}]}^{T},$$**vec**maps ${\mathrm{\Omega}}^{R\times A}$ to ${\mathrm{\Omega}}^{RA\times 1}$ by stacking the rows of a matrix in a vector, and ${[\xb7]}^{T}$ is used for the matrix-vector transposition. Likewise,

## (7)

$$g={[\begin{array}{cccccccccc}g({\phi}_{1},{\theta}_{1})& \dots & g({\phi}_{1},{\theta}_{A}),& g({\phi}_{2},{\theta}_{1})& \dots & g({\phi}_{2},{\theta}_{A}),& \cdots ,& g({\phi}_{R},{\theta}_{1})& \dots & g({\phi}_{R},{\theta}_{A})\end{array}]}^{T}$$## (8)

$$n={[\begin{array}{cccccccccc}n({\phi}_{1},{\theta}_{1})& \dots & n({\phi}_{1},{\theta}_{A}),& n({\phi}_{2},{\theta}_{1})& \dots & n({\phi}_{2},{\theta}_{A}),& \cdots ,& n({\phi}_{R},{\theta}_{1})& \dots & n({\phi}_{R},{\theta}_{A})\end{array}]}^{T}.$$According to Ref. 21, a discrete convolution equation corresponding to Eq. (4) can be written as

where $g,f$, and $n$ are of dimension $AR\times 1$, and $H$ is of dimension $AR\times AR$ and is ordered as## (10)

$$H={\left[\begin{array}{cccc}{H}_{1}& {H}_{R}& \cdots & {H}_{2}\\ {H}_{2}& {H}_{1}& \cdots & {H}_{3}\\ \vdots & \vdots & \ddots & \vdots \\ {H}_{R}& {H}_{R-1}& \cdots & {H}_{1}\end{array}\right]}_{RA\times RA}.$$The matrix $H$ is a block circulant consisting of circulant blocks with dimensions $RA\times RA$. Each partition of ${H}_{i}(i=\mathrm{1,2},\cdots R)$ has the size $A\times A$, which is constructed from the $i$’th scanning antenna pattern. Each partition ${H}_{i}$ can be written as

## (11)

$${H}_{i}={\left[\begin{array}{cccc}a({\theta}_{1},{\phi}_{i})& a({\theta}_{A},{\phi}_{i})& \cdots & a({\theta}_{2},{\phi}_{i})\\ a({\theta}_{2},{\phi}_{i})& a({\theta}_{1},{\phi}_{i})& \cdots & a({\theta}_{3},{\phi}_{i})\\ a({\theta}_{3},{\phi}_{i})& a({\theta}_{2},{\phi}_{i})& \cdots & \hspace{1em}\\ \vdots & \vdots & \ddots & \vdots \\ a({\theta}_{A},{\phi}_{i})& a({\theta}_{A-1},{\phi}_{i})& \cdots & a({\theta}_{1},{\phi}_{i})\end{array}\right]}_{A\times A}.$$The goal of angular super-resolution in scanning radar is to infer $f$ from the $g$, as accurately as possible. This task is called the deconvolution problem in this paper. The challenge is that the data $g$ recorded at the output of the radar system are a low-pass-filtered version of the original data $f$ due to the “low-pass” filtering effect of the antenna. So the high-frequency parts of the $g$ are lost by the antenna, which leads to the decrease of the azimuth resolution of the scanning radar system.

## 3.

## Augmented Lagrangian Solver for Angular Super-Resolution Imaging

In this section, we first analyze the phenomenon of noise amplification in the conventional inverse filter methods for solving a deconvolution problem. The deconvolution problem is converted into a constrained optimization task. Then the augmented Lagrangian method is presented as the foundation of the following deconvolution method. Finally, the angular super-resolution is realized by solving the constrained optimization via augmented Lagrangian method.

## 3.1.

### Deconvolution Problem and Constrained Optimization Problem

The conventional linear deconvolution method would be a simple division of data $g$ by the transferring function in the Fourier domain. Using the Fourier transform, Eq. (4) can be written as

where $G(w)$, $H(w)$, $F(w)$, and $N(w)$ are the Fourier transforms of $g$, $h$, $f$, and $n$, respectively. The linear deconvolution approach can be stated as the task of finding a linear operator $T$ such asThe challenge consists in that the convolution in the Fourier corresponds to multiplication while the deconvolution is Fourier division. For a radar system, the multipliers are often small for high frequencies, and the inverse filter $1/T(w)$ is large when $T(w)$ is very small. This results in large noise amplification, thus the angular super-resolution performance might be degraded. In order to solve the deconvolution problem, we employ the framework of regularization theory which provides the most likely solution given the recorded data and a reasonable assumption of the true scene. The details of the proposed methods are presented next.

To perform the deconvolution, we employ the framework of regularization theory, which not only suppresses the noise amplification but also allows for solutions that reflect reality. Regularization theory allows us to incorporate the prior information about the targets. In order to apply the regularization theory to the deconvolution problem, we first transform the deconvolution problem into an equivalent constrained optimization problem by adding the prior information of the targets.

In some applications of radar imaging, such as detection of ocean ships and aircraft monitoring, the number of dominant scattering targets is much smaller than the number of the overall samples. Therefore, the few large coefficients of the scatterers can capture most of the information of the scene while the weak scattering center is regarded as the noise in the radar imaging.^{12}

Considering the points that lie on the constant ${l}_{1}$-norm curve, the ones with two coordinates having equal magnitudes are favored by minimization of the ${l}_{k}$-norm.^{22} The effect of the ${l}_{1}$-norm favors a field with smaller number of dominant scatterers and better preserves the scatterers and their magnitudes. The deconvolution problem can be solved by constrained optimization based on ${l}_{1}$-norm. Finally, the regularization theory for angular super-resolution by solving the corresponding deconvolution problem leads to the following constrained optimization problem:

## (14)

$$\underset{f}{\mathrm{arg}\text{\hspace{0.17em}}\mathrm{min}}{\Vert Lf\Vert}_{1}\phantom{\rule{0ex}{0ex}}\text{subject to}:{\Vert g-Hf\Vert}_{2}^{2}\le \epsilon ,$$Equation (14) can be written as an equivalent unconstrained optimization problem, which yields the following equation:

## (15)

$$\underset{f}{\mathrm{arg}\text{\hspace{0.17em}}\mathrm{min}}\lambda {\Vert Lf\Vert}_{1}+\frac{1}{2}{\Vert g-Hf\Vert}_{2}^{2}.$$It should be noted that the problem in Eq. (15) can be seen as the Lagrangian problem [Eq. (14)]. From the optimization, the solution of Eq. (15) well approximates that of Eq. (14) only when $\lambda $ becomes large, which results in difficulty in numerically solving the equation. To avoid $\lambda $ going to infinity, the augmented Lagrangian method is used to solve Eq. (15), which avoids numerical instabilities that would occur as $\lambda \to \infty $. The details of this approach are presented next.

## 3.2.

### General Framework of Augmented Lagrangian Method

In this section, the general framework of the augmented Lagrangian method is presented as the foundation of the following discussion. The basic idea of the augmented Lagrangian method goes back to the work of Gabay and Mercier^{23} and Glowinski and Marocco.^{24}

Considering an unconstrained optimization problem of the form

## (16)

$$\underset{f}{\mathrm{arg}\text{\hspace{0.17em}}\mathrm{min}}\{{\theta}_{1}(Lf)+{\theta}_{2}(Hf)\},$$## (17)

$$\underset{f,z}{\mathrm{arg}\text{\hspace{0.17em}}\mathrm{min}}\{{\theta}_{1}(v)+{\theta}_{2}(u)\}\phantom{\rule{0ex}{0ex}}\text{subject to}:\left[\begin{array}{c}u\\ v\end{array}\right]=\left[\begin{array}{c}H\\ L\end{array}\right]f.$$Therefore, the augmented Lagrangian function associated with Eq. (17) is defined as

## (18)

$${L}_{A}(f,u,v,{\gamma}_{u},{\gamma}_{v})={\theta}_{1}(v)+{\theta}_{2}(u)+{\gamma}_{v}^{T}(v-Lf)+{\gamma}_{u}^{T}(u-Hf)\phantom{\rule{0ex}{0ex}}+\frac{{\rho}_{v}}{2}{\Vert v-Lf\Vert}_{2}^{2}+\frac{{\rho}_{u}}{2}{\Vert u-Hf\Vert}_{2}^{2},$$^{25}

The idea of the augmented Lagrangian method is to find a saddle point of ${L}_{A}(f,u,v,{\gamma}_{u},{\gamma}_{v})$, which is also the solution to the original problem in Eq. (16). Then the saddle point of the augmented Lagrangian function [Eq. (18)] can be found using the alternating direction method, which results in the iterative scheme in Table 1.

## Table 1

Augmented Lagrangian method.

Augmented Lagrangian method for Eq. (18) |
---|

Input: Set k=0, started at ρu, ρv, uk, vk, (γu)k and (γv)k |

Repeat |

Step 1: fk+1=argminf{[(γv)k(γu)k]T[vk−Lfuk−Hf]+‖[ρv2(vk−Lf)ρu2(uk−Hf)]‖22} |

Step 2: vk+1=arg minv{θ1(v)+[(γv)k]T(vk−Lfk+1)+ρv2‖v−Lfk+1‖22}. |

Step 3: uk+1=arg minu{θ2(u)+[(γu)k]T(uk−Hfk+1)+ρu2‖u−Hfk+1‖22}. |

Step 4: (γv)k+1=(γv)k−ρv(vk+1−Lfk+1). |

Step 5: (γu)k+1=(γu)k−ρu(uk+1−Hfk+1). |

Step 6: k=k+1. |

Until stopping criterion is satisfied. |

The convergence of the augmented Lagrangian method is guaranteed by the theorem in Ref. 26. The augmented Lagrangian method has been recently used for an imaging inverse problem.^{25}^{,}^{27} In Ref. 8, the augmented Lagrangian method is used to restore video. It is also studied in computed tomography.^{28} In the following, we apply the idea of the augmented Lagrangian method to solve Eq. (14).

## 3.3.

### Applying Augmented Lagrangian Method to Deconvolution Problem

Starting from Eq. (14), we introduce auxiliary constraint variables $u$ and $v$. The unconstrained optimization problem [Eq. (14)] can be written in the form of Eq. (17) as

## (19)

$${\theta}_{1}(f)=\lambda {\Vert Lf\Vert}_{1}\phantom{\rule{0ex}{0ex}}{\theta}_{2}(Hf)=\frac{1}{2}{\Vert g-Hf\Vert}_{2}^{2}\phantom{\rule{0ex}{0ex}}u=Hf,v=Lf.$$The unconstrained optimization problem is

## (20)

$$\underset{f,u,v}{\mathrm{arg}\text{\hspace{0.17em}}\mathrm{min}}\{E(f,u,v)=\lambda {\Vert v\Vert}_{1}+\frac{1}{2}{\Vert g-u\Vert}_{2}^{2}\}\phantom{\rule{0ex}{0ex}}\text{subject to}:u=Hf,\phantom{\rule[-0.0ex]{1em}{0.0ex}}v=Lf.$$The authors of Refs. 27 and 29 have also utilized this method in the context of imaging restoration. However, our emphasis here is on angular super-resolution in scanning radar. Before processing, we rewrite Eq. (17) concisely as

## (21)

$$\underset{z}{\mathrm{arg}\text{\hspace{0.17em}}\mathrm{min}}\text{\hspace{0.17em}}E(z)\phantom{\rule{0ex}{0ex}}\text{subject to}:z=Cf,$$The augmented Lagrangian equation associated with Eq. (21) is

## (22)

$${L}_{A}(f,z,\gamma )=E(z)+{\gamma}^{T}(z-Cf)+\frac{\mu}{2}{\parallel z-Cf\parallel}_{2}^{2},$$^{28}For easier manipulation, the multiplier term can be absorbed into the penalty term in Eq. (22), which leads to the following equation:

## (23)

$$\begin{array}{c}{L}_{A}(f,z,\gamma )=E(z)+\frac{\mu}{2}{\Vert z-Cf-\eta \Vert}_{2}^{2}+\text{const}\end{array},$$Next, we find the saddle point of Eq. (23) using the alternating direction method shown in Table 1, which results in the following equations:

## (24)

$${f}_{k+1}=\underset{f}{\mathrm{arg}\text{\hspace{0.17em}}\mathrm{min}}{\parallel {z}_{k}-Cf-{\eta}_{k}\parallel}_{2}^{2},$$## (25)

$${z}_{k+1}=\underset{z}{\mathrm{arg}\text{\hspace{0.17em}}\mathrm{min}}\text{\hspace{0.17em}}E(z)+\frac{\mu}{2}{\parallel z-C{f}_{k+1}-{\eta}_{k}\parallel}_{2}^{2},$$We now solve these equations one by one.

## 3.3.1.

#### f-subproblem

By dropping the index $k$, the solution of Eq. (24) is found by considering the normal equation as follows:

or, equivalently,The convolution matrix $H$ in Eq. (28) is a Block Toeplitz matrice, with each element of it being the Toeplitz block structure. This leads to an efficient solution for Eq. (28) using the preconditioned conjugate gradient method. Therefore, Eq. (28) has the following solution:

## (29)

$$f={\mathfrak{F}}^{-1}\left\{\frac{\mathfrak{F}{(H)}^{T}\mathfrak{F}[{u}_{k}+{({\eta}_{u})}_{k}]+\mathfrak{F}[{v}_{k}+{({\eta}_{v})}_{k}]}{\mathfrak{F}{(H)}^{T}\mathfrak{F}(H)+\mathfrak{F}{(L)}^{T}\mathfrak{F}(L)}\right\},$$## 3.3.2.

#### z-subproblem

We see that due to the structure of $E(z)$ and $C$, Eq. (25) can be separated into the following parts:

## (30)

$${u}_{k+1}=\underset{u}{\mathrm{arg}\text{\hspace{0.17em}}\mathrm{min}}\{\frac{1}{2}{\Vert g-u\Vert}_{2}^{2}+{({\gamma}_{u}^{T})}_{k}(u-H{f}_{k+1})+\frac{\mu}{2}{\Vert u-H{f}_{k+1}\Vert}_{2}^{2}\},$$## (31)

$${v}_{k+1}=\underset{v}{\mathrm{arg}\text{\hspace{0.17em}}\mathrm{min}}\{\lambda \Vert v{\Vert}_{1}+{({\gamma}_{v}^{T})}_{k}(v-{f}_{k+1})+\frac{\mu}{2}{\parallel v-{f}_{k+1}\parallel}_{2}^{2}\}.$$Equations (30) and (31) are independent of each other and can be solved simultaneously. Equation (30) is quadratic and has a closed solution.

According to Ref. 30, Eq. (31) can be solved using the shrinkage equation. Then, Eq. (31) becomes

## (33)

$${v}_{k+1}={\text{shrink}}_{\frac{\lambda}{\mu}}[{v}_{k}-{f}_{k+1}-{({\eta}_{v})}_{k}]\phantom{\rule{0ex}{0ex}}=\frac{{v}_{k}-{f}_{k+1}-{({\eta}_{v})}_{k}}{{\Vert {v}_{k}-{f}_{k+1}-{({\eta}_{v})}_{k}\Vert}_{2}}\times \mathrm{max}\{{\Vert {v}_{k}-{f}_{k+1}-{({\eta}_{v})}_{k}\Vert}_{2}-\frac{\lambda}{\mu},0\},$$Finally, the parameter ${\eta}_{u}$ and ${\eta}_{v}$ are updated as

To make this clear, the main flowchart of the proposed decovolution algorithm for angular super-resolution in scanning radar is shown in Fig. 2.

## 4.

## Experimental Results

In order to show the validity of the proposed deconvolution method for angular super-resolution in scanning radar, simulation and experimental results are presented in this section.

## 4.1.

### Simulations

In this section, the simulation results are presented to demonstrate the performance of the proposed deconvolution method for angular super-resolution in scanning radar. The simulation scene is shown in Fig. 3. The results are compared with those obtained by the Wiener filter method and the operator $L$ is chosen as the $I$, where $I$ represents the identity matrix. In the tests, the iteration was manually adjusted to yield the best angular super-resolution result for each experimental condition. To get the final result of the solution, the parameter $\mu $ is essential. The issue of how to optimally select $\mu $ is important, but is outside of the scope of this work. In this paper, the parameter is set as 10 for all simulations and experiments according to the method from the authors of Ref. 31.

Some related scanning radar system parameters are set as follows: the radar carrier frequency is 10 GHz and the pulse repetition frequency is 4000 Hz with an antenna scanning speed $20\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{deg}/\mathrm{s}$. The bandwidth of the transmitted signal is 2 MHz and the 3-dB width of the real beam is about 3 deg.

In order to quantify and compare the angular super-resolution performance of the proposed deconvolution method on the simulation data, relative error (ReErr),^{32} structure similarity (SSIM),^{33} and the peak-to-valley point difference in decibels are used in the simulations. They are defined as follows:

## (36)

$$\mathrm{ReErr}=\frac{{\parallel \widehat{f}-f\parallel}_{2}}{{\parallel f\parallel}_{2}};\phantom{\rule[-0.0ex]{2em}{0.0ex}}\mathrm{SSIM}=\frac{2{\rho}_{(\widehat{f},f)}\xb7(2{\mu}_{\widehat{f}}\xb7{\mu}_{f})}{({\mu}_{\widehat{f}}^{2}+{\mu}_{f}^{2})({\sigma}_{\widehat{f}}^{2}+{\sigma}_{f}^{2})},$$In order to show the angular super-resolution performance of the proposed method, we simulate the case of different signal-to-noise ratios (SNRs). Figure 5 shows the angular super-resolution results by using the proposed method and Wiener filter approach under different noise levels. Figure 6 shows the evolution of the ReErr and the SSIM with different SNRs. Some observations can be made based on the results in Figs. 5 and 6. For the isolated target, the echo is proportional to the antenna pattern. The echo of two adjacent targets is added and form a single peak as illustrated in Fig. 5. Hence, the individual response cannot be distinguished. The second row of Fig. 5 is the angular super-resolution results obtained by the Wiener filter method, and the third row of Fig. 5 represents the angular super-resolution results given by the proposed deconvolution method. In Fig. 6, the ReErr, whose ideal value is 0, is shown. When the value of the ReErr gets smaller, the precision of the angular super-resolution method increases. The metric SSIM is a quality measurement between the angular super-resolution result and the original scene. The value of SSIM is between $-1$ and 1, where 1 means the angular super-resolution result that fully identifies with the original scene. This leads to the larger value of the SSIM and a better performance of the proposed angular super-resolution method in terms of quality. It is clear that the proposed deconvolution method gives better results than the Wiener filter method. For example, in the case of $\mathrm{SNR}=10\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{dB}$, the angular super-resolution result of our method not only seems more appealing, but also has a better performance in terms of ReErr and SSIMs.

## 4.2.

### Experimental Results

In order to show the effectiveness of the proposed deconvolution algorithm for angular super-resolution in scanning radar, the real data results are presented. Figure 7 shows the three buildings and their distributions. The range between the scene center and the radar is about 594 m and the distance of each building is 45 m. The transmitted antenna pattern is obtained through the measurement result in the microwave anechoic chamber in the University of Electronic Science and Technology of China. However, due to the noise and energy degradation in actual situations, the antenna pattern used in deconvolution is not consistent with the transmitted pattern. To achieve nice results, we use a pretreatment for the antenna pattern which changes the length of the sidelobe while retaining the width of the main beam of the original pattern. The width of the sidelobe is changed step by step and is chosen by the best experiment result for random data in one scan. In our experiment for scene data processing, this method for choosing the appropriate antenna pattern is effective. The echo from the three buildings is much stronger than that from other areas so that the scene can be considered as the combination of three point targets. Employing the scanning radar system in Fig. 8, the data are acquired. The scanning radar system parameters for the experiment are listed in Table 2.

## Table 2

Experimental parameters for real data.

Parameters | Value | Units |
---|---|---|

Carrier frequency | 10 | GHz |

Bandwidth | 30 | MHz |

Pulse duration | 2 | μs |

Pulse repetition frequency | 200 | Hz |

Antenna scanning velocity | 70 | deg/s |

Antenna scanning area | −45 to +45 | deg |

Main-lobe beam width | 5 | deg |

Compared with the Wiener filter result, the proposed method yields a greater radar image quality with increased resolution, as shown in Fig. 9. Different rows in Fig. 9 correspond to the imaging result with the azimuth profile of the target scene. The first row of Fig. 9 shows the scanning radar image. The image cannot distinguish the positions of the buildings. Compared with the second row of Fig. 9, the third row shows the improved angular resolution and higher visual quality. The reason is that the prior information of the target is used in the proposed method, while the Wiener filter method uses no prior information of targets. Therefore, the proposed method is more effective compared with the Wiener filter method.

## 5.

## Conclusion

This paper presents a deconvolution algorithm for angular super-resolution in scanning radar. After presenting and analyzing the signal model of the scanning radar, we first convert the angular super-resolution problem into an equivalent deconvolution problem. In order to overcome the ill-posed nature of the deconvolution problem, the deconvolution problem is converted into a constrained optimization problem by incorporating the prior information about the targets. Finally, the constrained optimization problem is solved in the framework of the augmented Lagrangian method. Compared with the conventional Wiener filter method, the presented deconvolution algorithm improves the angular super-resolution performance of the scanning radar in terms of precision. Simulations and experimental results show a better performance of the presented method compared with the traditional Wiener filter method. However, there are still some issues that need to be studied in angular super-resolution based on the deconvolution method. Future work will study the robustness of the deconvolution algorithm for angular super-resolution under low SNR levels.

## References

D. L. Mensa, High Resolution Radar Cross-Section Imaging, Vol. 1, p. 280, Artech House, Boston, MA (1991).Google Scholar

D. R. Wehner, High Resolution Radar, Vol. 1, p. 484, Artech House, Norwood, MA (1987).Google Scholar

X. Qiu, D. Hu and C. Ding, “Some reflections on bistatic SAR of forward-looking configuration,” IEEE Geosci. Remote Sens. Lett. 5(4), 735–739 (2008).IGRSBY1545-598Xhttp://dx.doi.org/10.1109/LGRS.2008.2004506Google Scholar

M. A. Richards, “Iterative noncoherent angular superresolution [radar],” in Proc. 1988 IEEE National Radar Conf., 1988, pp. 100–105, IEEE (1988).Google Scholar

M. T. Alonso, P. López-Dekker and J. J. Mallorqui, “A novel strategy for radar imaging based on compressive sensing,” IEEE Trans. Geosci. Remote Sens. 48(12), 4285–4295 (2010).IGRSD20196-2892http://dx.doi.org/10.1109/TGRS.2010.2051231Google Scholar

M. Migliaccio and A. Gambardella, “Microwave radiometer spatial resolution enhancement,” IEEE Trans. Geosci. Remote Sens. 43(5), 1159–1169 (2005).IGRSD20196-2892http://dx.doi.org/10.1109/TGRS.2005.844099Google Scholar

S. Shah, “Deconvolution algorithms for fluorescence and electron microscopy,” PhD Thesis, The University of Michigan (2006).Google Scholar

S. H. Chan et al., “An augmented Lagrangian method for total variation video restoration,” IEEE Trans. Image Process. 20(11), 3097–3111 (2011).IIPRE41057-7149http://dx.doi.org/10.1109/TIP.2011.2158229Google Scholar

M. S. Almeida and M. A. Figueiredo, “Deconvolving images with unknown boundaries using the alternating direction method of multipliers,” IEEE Trans. Image Process. 22(8), 3074–3086 (2013).IIPRE41057-7149http://dx.doi.org/10.1109/TIP.2013.2258354Google Scholar

C. Soussen et al., “From Bernoulli–Gaussian deconvolution to sparse signal restoration,” IEEE Trans. Signal Process. 59(10), 4572–4584 (2011).ITPRED1053-587Xhttp://dx.doi.org/10.1109/TSP.2011.2160633Google Scholar

K. Browne, R. Burkholder and J. Volakis, “Fast optimization of through-wall radar images via the method of Lagrange multipliers,” IEEE Trans. Antenna Propag. 61(1), 320–328 (2013).IETPAK0018-926Xhttp://dx.doi.org/10.1109/TAP.2012.2220321Google Scholar

S.-J. Wei et al., “Sparse reconstruction for SAR imaging based on compressed sensing,” Prog. Electromagn. Res. 109, 63–81 (2010).PELREX1043-626Xhttp://dx.doi.org/10.2528/PIER10080805Google Scholar

S. Ji, Y. Xue and L. Carin, “Bayesian compressive sensing,” IEEE Trans. Signal Process. 56(6), 2346–2356 (2008).ITPRED1053-587Xhttp://dx.doi.org/10.1109/TSP.2007.914345Google Scholar

L. Zhang et al., “High-resolution ISAR imaging by exploiting sparse apertures,” IEEE Trans. Antenna Propag. 60(2), 997–1008 (2012).IETPAK0018-926Xhttp://dx.doi.org/10.1109/TAP.2011.2173130Google Scholar

M. Çetin and W. C. Karl, “Feature-enhanced synthetic aperture radar image formation based on nonquadratic regularization,” IEEE Trans. Image Process. 10(4), 623–631 (2001).IIPRE41057-7149http://dx.doi.org/10.1109/83.913596Google Scholar

K. E. Browne and R. J. Burkholder, “Nonlinear optimization of radar images from a through-wall sensing system via the Lagrange multiplier method,” IEEE Geosci. Remote Sens. Lett. 9(5), 803–807 (2012).IGRSBY1545-598Xhttp://dx.doi.org/10.1109/LGRS.2012.2188093Google Scholar

D. P. Bertsekas, “Multiplier methods: a survey,” Automatica 12(2), 133–145 (1976).ATCAA90005-1098http://dx.doi.org/10.1016/0005-1098(76)90077-7Google Scholar

M. Tao, J. Yang and B. He, “Alternating direction algorithms for total variation deconvolution in image reconstruction,” TR0918, Department of Mathematics, Nanjing University (2009).Google Scholar

W. Li, J. Yang and Y. Huang, “Keystone transform-based space-variant range migration correction for airborne forward-looking scanning radar,” Electron. Lett. 48(2), 121–122 (2012).ELLEAK0013-5194http://dx.doi.org/10.1049/el.2011.2774Google Scholar

M. A. Richards, Fundamentals of Radar Signal Processing, Tata McGraw-Hill Education, New York (2005).Google Scholar

M. Piles et al., “Spatial-resolution enhancement of SMOS data: a deconvolution-based approach,” IEEE Trans. Geosci. Remote Sens. 47(7), 2182–2192 (2009).IGRSD20196-2892http://dx.doi.org/10.1109/TGRS.2009.2013635Google Scholar

M. Cetin, “Feature-enhanced synthetic aperture radar imaging,” PhD Thesis, University of Salford, Manchester, UK (2001).Google Scholar

D. Gabay and B. Mercier, “A dual algorithm for the solution of nonlinear variational problems via finite element approximation,” Comput. Math. Appl. 2(1), 17–40 (1976).CMAPDK0898-1221http://dx.doi.org/10.1016/0898-1221(76)90003-1Google Scholar

R. Glowinski and A. Marroco, “Sur l’approximation, par éléments finis d’ordre un, et la résolution, par pénalisation-dualité d’une classe de problèmes de dirichlet non linéaires,” ESAIM: Math. Modell. Numer. Anal. 9(R2), 41–76 (1975).Google Scholar

M. V. Afonso, J. M. Bioucas-Dias and M. A. Figueiredo, “Fast image recovery using variable splitting and constrained optimization,” IEEE Trans. Image Process. 19(9), 2345–2356 (2010).IIPRE41057-7149http://dx.doi.org/10.1109/TIP.2010.2047910Google Scholar

J. Eckstein and D. P. Bertsekas, “On the Douglas-Rachford splitting method and the proximal point algorithm for maximal monotone operators,” Math. Program. 55(1–3), 293–318 (1992).http://dx.doi.org/10.1007/BF01581204Google Scholar

M. V. Afonso, J. M. Bioucas-Dias and M. A. Figueiredo, “An augmented Lagrangian approach to the constrained optimization formulation of imaging inverse problems,” IEEE Trans. Image Process. 20(3), 681–695 (2011).IIPRE41057-7149http://dx.doi.org/10.1109/TIP.2010.2076294Google Scholar

S. Ramani and J. A. Fessler, “A splitting-based iterative algorithm for accelerated statistical x-ray CT reconstruction,” IEEE Trans. Med. Imaging 31(3), 677–688 (2012).ITMID40278-0062http://dx.doi.org/10.1109/TMI.2011.2175233Google Scholar

M. A. Figueiredo and J. M. Bioucas-Dias, “Restoration of Poissonian images using alternating direction optimization,” IEEE Trans. Image Process. 19(12), 3133–3145 (2010).IIPRE41057-7149http://dx.doi.org/10.1109/TIP.2010.2053941Google Scholar

P. L. Combettes and V. R. Wajs, “Signal recovery by proximal forward-backward splitting,” Multiscale Model. Simul. 4(4), 1168–1200 (2005).MMSUBT1540-3459http://dx.doi.org/10.1137/050626090Google Scholar

W. Yin et al., “Practical compressive sensing with Toeplitz and circulant matrices,” in Visual Communications and Image Processing 2010, pp. 77440K, SPIE Press, Bellingham, Washington (2010).Google Scholar

Y.-M. Huang, M. K. Ng and Y.-W. Wen, “A new total variation method for multiplicative noise removal,” SIAM J. Imaging Sci. 2(1), 20–40 (2009).1936-4954http://dx.doi.org/10.1137/080712593Google Scholar

Z. Wang and A. C. Bovik, “Mean squared error: love it or leave it? A new look at signal fidelity measures,” IEEE Signal Process. Mag. 26(1), 98–117 (2009).ISPRE61053-5888http://dx.doi.org/10.1109/MSP.2008.930649Google Scholar

## Biography

**Yuebo Zha** received his BS degree in applied mathematics and his MS degree in operation research from the University of Electronic Science and Technology of China in 2011, where he is currently working toward his PhD in signal and information processing.

**Yulin Huang** received his BS and PhD degrees in electronic engineering from the University of Electronic Science and Technology of China (UESTC), Chengdu, China, in 2002 and 2008, respectively. He is currently an associate professor with UESTC. His research interests include signal processing and radar imaging.

**Jianyu Yang** received his BS degree in electronic engineering from the National University of Defense Technology, Changsha, China, in 1984 and his MS and PhD degrees in electronic engineering from the University of Electronic Science and Technology of China (UESTC), Chengdu, China, in 1987 and 1991, respectively. He is currently a professor with the School of Electronic Engineering, UESTC. His research interests are mainly in synthetic aperture radar and signal processing.