Extraction of Anisotropic Characteristics of Scattering Centers and Feature Enhancement in Wide-Angle SAR Imagery Based on the Iterative Re-Weighted Tikhonov Regularization

The anisotropic characteristic reflects discriminating features of the geometry of a scattering center. In this study, we propose a novel method based on the iterative re-weighted Tikhonov regularization (IRWTR) to achieve the extraction of anisotropic characteristics of scattering centers in the wide-angle SAR synthetic aperture radar (SAR) imaging. Moreover, based on the extracted anisotropic scattering behaviors, the incomplete imaging results of the distributed scattering centers are restored. In this paper, we first discussed the scattering property in SAR imagery from the perspective of attributed scattering center model (ASCM). The reason for the incomplete imaging results of the distributed scattering centers was also discussed based on the ASCM. Subsequently, we modeled the aspect-dependent amplitude responses of a scattering center as a linear combination of a set of orthogonal basis. Consequently, the extraction of anisotropic characteristics can be transformed into an inverse problem, which was solved by the proposed IRWTR with high efficiency and accuracy. After the extraction, we attempted to restore the complete image of a distributed scattering center, which consisted of only two points. The enhanced SAR image can reveal the actual shape of a target. Processing results of electromagnetic computation data validated that the proposed method is effective and efficient.


Introduction
Wide-angle SAR synthetic aperture radar (SAR) imagery can provide SAR images with higher azimuth resolution in comparison with the traditional small-angle SAR. In addition to the finer cross-range resolution, wide-angle SAR also benefits the detection and extraction of anisotropic characteristics of scatterers. Anisotropic characteristics, which can be identified in wide-angle SAR imaging, provide us with more interesting and important information about targets. Anisotropic characteristic implies the amplitude response of a scattering center can vary with the look angle [1][2][3]. We can infer geometric features of a scatterer from its anisotropic characteristic. Therefore, it is well worth studying the method to extract anisotropy characteristics from radar measurements to broaden the understanding of targets in wide-angle SAR image.
imaging scene. The radar transmitted signals, which illuminated the imaging scene. The signals reflected from the imaging scene were collected by the radar, which were determined by the scattering characteristics and the positions of the targets in the imaging scene. The radar transmitted and collected the signals at different positions on the flight path. The radar antenna always illuminated the imaging scene during the procedure. With the assumption of far-field and plane-wave propagation,  denotes the look angle of the radar that is identical for all the scattering centers in the imaging scene. In wide-angle SAR imaging, the dependence of the amplitude response on the look angle  should be taken into consideration. In wide-angle SAR imaging, the amplitude response of the p th scattering center can be expressed as ( ) , , ,   (2) for a distributed scattering center. Trihedral, top hat, and sphere are canonical local scattering centers. Distributed scattering centers include dihedral, flat plate, cylinder, and broadside edge diffraction. Generally, the local scattering center shows isotropic characteristic. In contrast, distributed scattering center presents anisotropic characteristic. In (1) and (2)  In wide-angle SAR imaging, the dependence of the amplitude response on the look angle ϕ should be taken into consideration. In wide-angle SAR imaging, the amplitude response of the pth scattering center can be expressed as E p f , ϕ, x p , y p , where f is the frequency of the transmitted signal. x p , y p denotes the position of the scattering center. From the perspective of ASCM, E p f , ϕ, x p , y p can be expressed as: for a local scattering center or: · exp −j 4π f c · x p cos ϕ + y p sin ϕ (2) for a distributed scattering center. Trihedral, top hat, and sphere are canonical local scattering centers. Distributed scattering centers include dihedral, flat plate, cylinder, and broadside edge diffraction. Generally, the local scattering center shows isotropic characteristic. In contrast, distributed scattering center presents anisotropic characteristic. In (1) and (2), sin c(x) = sin x x and j = √ −1. A p is the relative amplitude of the pth scatterer; α p reflects the frequency dependence; γ p , L p , ϕ 0p reflects the aspect dependence. γ p denotes the tiny aspect dependence of a local scattering center, whose influence can be omitted in practical. L p , ϕ 0p are the length and the orientation angle of a distributed scattering center. The orientation angle is the angle between the flank of a distributed scattering center and X-axis of the imaging scene. The SAR image of a local scattering center is a point. The images of a distributed scattering center with different orientations are quite different. For a distributed scattering center, if ϕ = ϕ 0p at a look angle, the image will be a line whose length equals the length of the distributed scattering center. However, if ϕ = ϕ 0p in the acquiring procedure, the imaging result of the distributed Remote Sens. 2018, 10, 2066 4 of 19 scattering center will be composed of two isolated points. If ϕ = ϕ 0p in the acquiring procedure, sinc 2π f L p c · sin ϕ − ϕ 0p in (1) can be written as: From (3), a distributed scattering center can be regarded as two isolated local scattering centers and their responses can be expressed as: and: According to (4) and (5), a distributed scattering center becomes two local scattering centers, and they are at the ends of the distributed scattering center. Consequently, the imaging result of the distributed scattering center only comprises of two points. To illustrate this phenomenon, the backscattered data of the distributed scattering centers with different orientation angles are simulated by FEKO [11] and the corresponding images are shown in Figures 2 and 3. The relative information of this simulation is given in Section 4.2.
( ) From (3), a distributed scattering center can be regarded as two isolated local scattering centers and their responses can be expressed as: and: According to (4) and (5), a distributed scattering center becomes two local scattering centers, and they are at the ends of the distributed scattering center. Consequently, the imaging result of the distributed scattering center only comprises of two points. To illustrate this phenomenon, the backscattered data of the distributed scattering centers with different orientation angles are simulated by FEKO [11] and the corresponding images are shown in

Signal Model of Return Signals
The , where c f is the center frequency. The approximations are based on the fact that the bandwidth of a normal radar system is quite small in comparison with the center frequency [1,4]. Then, the simplified ASCM can be expressed as: , exp cos sin The extraction of anisotropic characteristics can be achieved by estimating ( ) , pc sf  .
Assuming there are P scattering centers in the scene and the radar measurements are collected at K discrete frequencies and N discrete look angles. The received signal ( ) , kn yf , which is collected by the radar at the frequency k f and the look angle n  , can be expressed as: where ( ) , kn nf is the zero mean additive complex white Gaussian noise and 12 , , , P A A A are the relative amplitudes of the P scattering centers, respectively. We used a vector set   12 , , , , ,  (8) The matrix B plays an important role in designing the processing algorithm of anisotropic characteristic extraction. We will discuss B in Section 3. Subsequently, we can use (8) to express (7) in matrix form: =+ y Φan (9)

Signal Model of Return Signals
The frequency-dependent factor g( f ) = (j · f / f c ) α p and aspect-dependent factor (for the distributed scattering centers) or s p ( f , ϕ) = exp −2πγ p f sin ϕ (for the local scattering centers) can be simplified as g( f ) ≈ 1 and where f c is the center frequency. The approximations are based on the fact that the bandwidth of a normal radar system is quite small in comparison with the center frequency [1,4]. Then, the simplified ASCM can be expressed as: The extraction of anisotropic characteristics can be achieved by estimating s p ( f c , ϕ).
Assuming there are P scattering centers in the scene and the radar measurements are collected at K discrete frequencies and N discrete look angles. The received signal y( f k , ϕ n ), which is collected by the radar at the frequency f k and the look angle ϕ n , can be expressed as: where n( f k , ϕ n ) is the zero mean additive complex white Gaussian noise and A 1 , A 2 , · · · , A P are the relative amplitudes of the P scattering centers, respectively.
We used a vector set T . b m is a column vector containing N elements. Then, s p can be represented as: The matrix B plays an important role in designing the processing algorithm of anisotropic characteristic extraction. We will discuss B in Section 3. Subsequently, we can use (8) to express (7) in matrix form: Remote Sens. 2018, 10, 2066 The variables in (9) are listed below: where represents Hadamard product and 1 M is a M × 1 vector of ones. According to (9), we can get a by solving an inverse problem. After obtaining a, the azimuth-dependent amplitudes of scattering centers can be retrieved by B and (8).

The IRWTR Algorithm and Feature Enhancement of Wide-Angle SAR Images
In this section, we firstly present the IRWTR algorithm for the extraction of anisotropic characteristics. Secondly, based on the anisotropic characteristics extracted by the IRWTR and the minimization of the least squares error criterion, we propose a method to restore the incomplete imaging results of certain distributed scattering centers.
A diagram to illustrate the proposed algorithm is shown in Figure 4. Figure 4a shows the imaging scene, which includes a metal rod and metal ball. The corresponding imaging result is shown in Figure 4b. There are three points in the image, which cannot reflect the actual imaging scene. We used the IRWTR algorithm to extract the aspect-dependent amplitudes of the scatterers that are shown in Figure 4c. From the aspect-dependent amplitudes, we restored the incomplete imaging result of the metal rod, which is shown in Figure 4d. The variables in (9) are listed below:

The IRWTR Algorithm and Feature Enhancement of Wide-Angle SAR Images
In this section, we firstly present the IRWTR algorithm for the extraction of anisotropic characteristics. Secondly, based on the anisotropic characteristics extracted by the IRWTR and the minimization of the least squares error criterion, we propose a method to restore the incomplete imaging results of certain distributed scattering centers.
A diagram to illustrate the proposed algorithm is shown in Figure 4. Figure 4a shows the imaging scene, which includes a metal rod and metal ball. The corresponding imaging result is shown in Figure 4b. There are three points in the image, which cannot reflect the actual imaging scene. We used the IRWTR algorithm to extract the aspect-dependent amplitudes of the scatterers that are shown in Figure 4c. From the aspect-dependent amplitudes, we restored the incomplete imaging result of the metal rod, which is shown in Figure 4d.
In the following subsections, the main steps of the proposed method are introduced.  In the following subsections, the main steps of the proposed method are introduced.

Selection of Scattering Centers and Estimation of the Relative Amplitudes and Positions
In the first part of the proposed scheme, the positions and relative amplitudes of scatterers were estimated. Firstly, we used the Backprojection (BP) algorithm [12] to achieve the SAR image. Based on the imaging result, we could perform the estimation of the positions of scatterers. After the position estimation, the RELAX method [4] was utilized to accurately estimate the relative amplitudes of the scatterers. The positions and relative amplitudes were required in forming Φ.

Anisotropic Characteristic Estimation
Apart from the positions and the relative amplitudes of the scattering centers, we still needed an appropriate matrix B to form matrix Φ. In this paper, we proposed a parsimony B which is an identity matrix with the dimensions N × N and can be expressed as: The physical meaning of B is clear and sample. Its nth column is the basis that only represents the amplitude of a scattering center at the nth look angle. Then, the estimation of the scattering centers' aspect-dependent amplitudes can be achieved by solving a from y, which can be expressed as the following problem:â = argmin a y − Φa 2 2 (11) Problem (11) is a least squares problem, and its analytic solution is: Equation (12) only provides a simple and reasonable solution. Under the condition of noise and other perturbations, the performance of (12) may be poor.
Because of the poor performance of LS method, it was necessary to impose a regularization term on the problem in (11), which could make the solution more accurate and suppress noise's influence. We used the Tikhonov regularization term to make the solution more accurate. The rationality of using Tikhonov regularization was based on the fact that the aspect-dependent amplitude function of look angle of a scatterer was normally continuous and smooth. The objective function with the Tikhonov regularization can be expressed as:â where F is a gradient operator and can be expressed as: D is the matrix of first-order derivatives with the dimensions (M − 1) × M and can be expressed as: λ is the regularization parameter to control the tradeoff between y − Φa 2 2 and Fa 2 2 . With a relative smaller regularization parameter, the solution will vary rapidly. Consequently, a smaller regularization parameter was suitable to retrieve the anisotropic scattering. In contrast, a larger regularization parameter was appropriate for the inversion of the isotropic scattering. To use a uniform regularization parameter was not appropriate. It was necessary to choose specific λ for each scattering center based on its own scattering property. By taking the individual scattering attributes into consideration, we proposed the MRP Tikhonov regularization. The problem with the MRP Tikhonov regularization can be expressed as: where F w is expressed as: In (17), λ p is the regularization parameter for the pth scattering center. Since λ = λ 1 , λ 2 , · · · , λ P T has a significant influence on the solution, the selection of it is crucial. In order to accurately select a suitable λ, we discussed the exact physical meaning of λ, based on which selection method was introduced. Similar to [13], we also revealed the regularization parameter's meaning from the perspective of the maximum a posteriori (MAP) estimation. If we use the MAP theory to understand (16), the procedure of estimating a is:â where P(y − Φa) can be expressed as: where P(n) is the probability density function (PDF) of the noise which is a complex Gaussian PDF and σ 2 n is the noise variance. σ 2 n can be estimated from the SAR imaging result [13]. To estimate σ 2 n from the data, we selected an area of the imaging result of the data which contained no targets. Then, the estimation of σ 2 n can achieved by calculating the mean power of this area. If Ith pixels of the image are used, the mean square error of the estimator is σ 2 n I . P(a) is the PDF of the amplitudes of scatterers and can be expressed as: since the distributions of the amplitudes of different scattering centers were independent. In a p , a p,n is the amplitude response of the pth scattering center at the nth look angle, which indicates that adjacent elements in a p show correlation. Consequently, it is reasonable to model a p as the first-order Gaussian-Markov process [14,15]: a p,n = a p,n−1 + w p,n where w p,n is a random variable and w p,n ∼ CN 0, σ 2 p . For a different n, w p,n is independent and identically distributed (i.i.d). The real part and imaginary part of a p,1 are two random variables with the independent uniform distributions whose PDF are a constant C p . With a larger σ 2 p , the variation between two adjacent elements of a p would be more bigger, which indicates that the anisotropic feature of pth scattering center is more significant. In Figure 5, we show an example to illustrate the influence of σ 2 p . In Figure 5, it is clear that the elements in a p varies more significantly as σ 2 p increases. w is independent and identically distributed (i.i.d). The real part and imaginary part of ,1 p a are two random variables with the independent uniform distributions whose PDF are a constant p C . With a larger 2 p  , the variation between two adjacent elements of p a would be more bigger, which indicates that the anisotropic feature of p th scattering center is more significant. In Figure 5, we show an example to illustrate the influence of 2 p  . In Figure 5, it is clear that the elements in From (16), (17) and (23), λ can be expressed as: , , , , , From (24), λ is related to the anisotropic characteristics of scattering centers, which we cannot know before the problem (16) is solved. Thus, we needed to jointly estimate a and λ . Based on [16], Gamma distribution was used as the prior of p  and its PDF can be expressed as: Based on (21) and the PDF of a p,1 , P a p could be written as: P a p = P a p,1 , a p,2 , · · · a p,N = P a p,1 P w p,2 P w p,3 · · · P w p,N According to (18) From (16), (17) and (23), λ can be expressed as: From (24), λ is related to the anisotropic characteristics of scattering centers, which we cannot know before the problem (16) is solved. Thus, we needed to jointly estimate a and λ. Based on [16], Gamma distribution was used as the prior of λ p and its PDF can be expressed as: where Γ(·) is the Gamma function and α and β are the hyper-parameters. The MAP estimation of a and λ can be achieved by solving the following problem: As σ 2 n can be estimated from the SAR image, we mainly focused on the estimation of σ = σ 2 1 , · · · , σ 2 p , · · · σ 2 P T . Based on (25) and (26), the MAP estimation of a and σ can be expressed as: It is difficult to solve the problem (27) directly. In this study, the EM algorithm was employed to solve (27), which yielded an iterative algorithm. We use σ t−1 (t is the iteration number) to calculateâ t by the following equation:â where F w,t−1 is expressed as Then, σ t = σ 2 1,t , · · · , σ 2 p,t , · · · , σ 2 P,t T is calculated by the following equation: where var Da p,t denotes the variance of Da p,t . After σ t is obtained, the matrix F w,t−1 is updated for the next iteration. Problem (28) has an analytic solution. In the following, we derived the solution with the closed form. The conjugate gradient of J(a t ) = y − Φa t 2 2 + F w,t−1 a t 2 2 can be expressed as By setting ∂J(a t ) ∂a * t = 0, the minimizer of (28) can be expressed as Equation (32) is the closed form of the solution to (28), so (28) can be solved directly. Consequently, the proposed method can be implemented with high efficiency.
To sum up, the IRWTR algorithm consists of two steps. In the first step, we calculate a t by (28) and σ t−1 . In the second step, a new σ t is calculated via (30) to update the regularization parameters. This procedure will be repeated until the difference between σ t−1 and σ t in two successive iterations is smaller than a small threshold.

The Restoration of the Incomplete Images of Certain Distributed Scattering Centers
In this subsection, we discuss the method to restore the incomplete imaging results of the distributed scattering centers. We have mentioned that some points in the SAR image are actually the ends of distributed scattering centers. With the aspect-dependent amplitudes extracted by the IRWTR, we can distinguish these ends from the local scattering centers. Furthermore, we need to know which two points belong to the same distributed scattering center and restore the incomplete imaging result.
Firstly, we distinguished the ends of the distributed scattering centers from the local scattering centers based on the variance of the aspect-dependent amplitudes. A threshold was set and the scattering centers' aspect-dependent-amplitude variances below the threshold will be identified to be local scattering centers, while the rest were identified to be the ends of distributed scattering centers.
Secondly, these ends were arranged in pairs, of which the procedure is referred to as the matching procedure. The matching procedure exploits the ASCM and the minimization of least squares error criterion. With the aspect-dependent amplitudes extracted by the IRWTR, the position and the relative amplitude, the estimated backscatter of the pth scatterer at f c can be expressed as: where E est p ( f c , ϕ n ) is calculated by (6). If the i and jth scattering center are the two ends of the same distributed scattering center, the least squares error e LS i,j , expressed as the following equation, will be minimized: where E est i,j ( f c ) represents the backscatter of the distributed scattering center whose ends are at the ith and jth scatterers' positions. The length L est i,j , orientation angle ϕ est 0,i,j , position x est i,j , y est i,j and relative amplitude A est i,j of this distributed scattering center are calculated by the following equations: are the estimated positions of the ith and jth scatterer. A est i and A est j are the estimated relative amplitides of the ith and jth scatterer. The backscatter E est i,j ( f c ) can be calculated by (1)

and the parameters in (35)-(39).
At the beginning of the matching procedure, we chose one end and found the other end that was paired with it by minimizing (34). Then, this procedure was repeated among the residual ends. The matching procedure was stopped until all the ends were arranged in pairs.
Thirdly, we restored the incomplete imaging results of the distributed scattering centers based on the matching result. We added a signal to the original data in the frequency-aspect domain. The imaging result of this signal was the flank of a distributed scattering center. For the distributed scattering center whose ends are the ith and jth scatterer, we added the signal expressed as the following equation to the original data: where the parameters are calculated by (35)-(39). After we fixed all the incomplete imaging results of the distributed scattering centers by using (40), a BP algorithm was performed to yield the improved SAR image. Based on the aforementioned contents, we present the flowchart of the whole algorithm in Figure 6. In the first part, the extraction of the scattering centers' positions and relative amplitudes is performed. The relative amplitudes and positions are used to form the matrix Φ in (9). In the second part, the IRWTR is utilized to estimate the aspect-dependent amplitudes of the scattering centers. In the third part, the incomplete imaging results of the distributed scattering centers are restored based on the extracted anisotropic characteristics.
At the beginning of the matching procedure, we chose one end and found the other end that was paired with it by minimizing (34). Then, this procedure was repeated among the residual ends. The matching procedure was stopped until all the ends were arranged in pairs. Thirdly, we restored the incomplete imaging results of the distributed scattering centers based on the matching result. We added a signal to the original data in the frequency-aspect domain. The imaging result of this signal was the flank of a distributed scattering center. For the distributed scattering center whose ends are the i th and j th scatterer, we added the signal expressed as the following equation to the original data: ,  est  est  est  est  est  ,  ,  ,  ,  0, ,  0, ,   2  , , , sinc sin sin where the parameters are calculated by (35)-(39). After we fixed all the incomplete imaging results of the distributed scattering centers by using (40), a BP algorithm was performed to yield the improved SAR image.
Based on the aforementioned contents, we present the flowchart of the whole algorithm in Figure 6. In the first part, the extraction of the scattering centers' positions and relative amplitudes is performed. The relative amplitudes and positions are used to form the matrix Φ in (9). In the second part, the IRWTR is utilized to estimate the aspect-dependent amplitudes of the scattering centers. In the third part, the incomplete imaging results of the distributed scattering centers are restored based on the extracted anisotropic characteristics.
In the next section, the processing results by the proposed method will be shown.  Figure 6. Flowchart of the proposed method.

Experimental Results
In this section, we apply the proposed algorithm to the simulated data and the corresponding results are shown. In Section 4.1, we show the IRWTR processing result of the data simulated by MATLAB. In the next section, the processing results by the proposed method will be shown.

Experimental Results
In this section, we apply the proposed algorithm to the simulated data and the corresponding results are shown. In Section 4.1, we show the IRWTR processing result of the data simulated by MATLAB.

Results of Matlab Simulated Data and Analysis
The processing result of Matlab simulated data will be shown in this section. There are three scattering centers whose positions are (0 m, 1 m) (0 m, 0 m) and (1.5 m, 0 m), respectively, in the imaging scene. Their amplitude functions of azimuth were a constant, Sinc function and hamming window function. The radar transmitted step sweep frequency signal with a start frequency of 9.5 GHz and bandwidth of 1 GHz. The number of frequency signals was 32. The look angles ranged from −22.5 • to 22.5 • . The number of the look angles was 64.
Two methods were used to process the simulated data: solving (13) and the IRWTR algorithm. The processing was undertaken under different noise conditions (10 dB and 0 dB) and the results are shown in Figures 7 and 8, where the blue dashed line shows the real value and the red solid line shows the estimation.
In Figure 7, the performance of the estimation of the Sinc-like amplitude function was the best. In contrast, the other estimations were inaccurate, for the regularization parameter was not large enough. This phenomenon is more significant in Figure 8. According to estimation of the Sinc-like amplitude function in Figure 8b, the regularization parameter was larger than the requirement, for the estimation was too smooth and deviated from the real values. In contrast, the regularization parameters for the rest amplitude functions were smaller than the requirement, which made the estimations oscillate. Compared with the regularization with a single regularization parameter, the MRP regularization provided better estimations. Both in Figures 7 and 8, the estimations obtained via IRWTR were more close to the real values.
To evaluate the performance of different methods under different conditions of signal-to-noise ratio (SNR) SNR, a Monte Carlo simulation was used, in which 50 trials were performed for each SNR. In each trial, the relative mean squared error (RMSE) [17] was calculated. The statistical analysis is shown in Figure 9. In Figure 9, the performance of the IRWTR was more superior as SNR decreased.  In Figure 7, the performance of the estimation of the Sinc-like amplitude function was the best. In contrast, the other estimations were inaccurate, for the regularization parameter was not large enough. This phenomenon is more significant in Figure 8. According to estimation of the Sinc-like amplitude function in Figure 8b, the regularization parameter was larger than the requirement, for the estimation was too smooth and deviated from the real values. In contrast, the regularization parameters for the rest amplitude functions were smaller than the requirement, which made the estimations oscillate. Compared with the regularization with a single regularization parameter, the MRP regularization provided better estimations. Both in Figures 7 and 8, the estimations obtained via IRWTR were more close to the real values.

Azimuth
To evaluate the performance of different methods under different conditions of signal-to-noise ratio (SNR) SNR, a Monte Carlo simulation was used, in which 50 trials were performed for each SNR. In each trial, the relative mean squared error (RMSE) [17] was calculated. The statistical analysis is shown in Figure 9. In Figure 9, the performance of the IRWTR was more superior as SNR decreased.

The Results of Electromagnetic Simulation Data and Analysis
In this section, we apply the proposed processing procedure to the data simulated by FEKO software and the processing results will be shown.
The imaging scene is shown in Figure 10, which included two metal rods and three metal balls. The radar transmitted a step sweep frequency signal with a start frequency of 8 GHz, bandwidth 4 GHz. The azimuth start and end angle were −22.5° and 22.5°. We applied the BP imaging algorithm to the simulation data and the result is shown in Figure 11.
In Figure 11, there are only seven points. In fact, we knew that there were two distributed scattering centers in the scene from Figure 10. However, because their specific orientations, we could not recognize these two distributed scattering centers from Figure. 11. We applied the proposed method in Section 3.2 to the data and extract the aspect-dependent amplitudes of all points. The estimated aspect-dependent amplitudes are presented in Figure 12. According to Figure 12, we can find that the amplitudes of certain scatterers were almost independent of the look angle. These scatterers with slowly varying amplitudes were the imaging results of the metal balls, which were the actual local scattering centers. The variations in the amplitude responses of the rest were more

The Results of Electromagnetic Simulation Data and Analysis
In this section, we apply the proposed processing procedure to the data simulated by FEKO software and the processing results will be shown.
The imaging scene is shown in Figure 10, which included two metal rods and three metal balls. The radar transmitted a step sweep frequency signal with a start frequency of 8 GHz, bandwidth 4 GHz. The azimuth start and end angle were −22.5 • and 22.5 • . We applied the BP imaging algorithm to the simulation data and the result is shown in Figure 11. remarkable. The degrees of variation in amplitude responses could help us distinguish the ends of distributed scattering centers from the local scattering centers. Based on the azimuth-dependent scattering behaviors shown in Figure 12, the ends of distributed scattering centers are indicated by red circles in Figure 11.  The estimated aspect-dependent amplitudes of the scatterer 1, 2, 4, and 7 were compared with the ones calculated by FEKO. The comparisons are shown in Figure 13. In Figure 11, there are only seven points. In fact, we knew that there were two distributed scattering centers in the scene from Figure 10. However, because their specific orientations, we could not recognize these two distributed scattering centers from Figure 11. We applied the proposed method in Section 3.2 to the data and extract the aspect-dependent amplitudes of all points. The estimated aspect-dependent amplitudes are presented in Figure 12. According to Figure 12, we can find that the amplitudes of certain scatterers were almost independent of the look angle. These scatterers with slowly varying amplitudes were the imaging results of the metal balls, which were the actual local scattering centers. The variations in the amplitude responses of the rest were more remarkable. The degrees of variation in amplitude responses could help us distinguish the ends of distributed scattering centers from the local scattering centers. Based on the azimuth-dependent scattering behaviors shown in Figure 12, the ends of distributed scattering centers are indicated by red circles in Figure 11. distributed scattering centers from the local scattering centers. Based on the azimuth-dependent scattering behaviors shown in Figure 12, the ends of distributed scattering centers are indicated by red circles in Figure 11. Observing angle(°) Amplitude Figure 12. The aspect-dependent amplitudes of the scattering centers.
The estimated aspect-dependent amplitudes of the scatterer 1, 2, 4, and 7 were compared with the ones calculated by FEKO. The comparisons are shown in Figure 13. The estimated aspect-dependent amplitudes of the scatterer 1, 2, 4, and 7 were compared with the ones calculated by FEKO. The comparisons are shown in Figure 13.
From Figure 13, it can be found that the estimated amplitudes corresponded closely with the ones calculated by FEKO. Moreover, we arranged the ends in pairs based on the extracted amplitudes. The signals to form the images of the flanks were constructed based on the arrangement and added to the original return signals. The imaging result of the reproduced data is shown in Figure 14. In Figure 14, the distributed scattering centers were easier to be distinguished from the local scattering centers. From Figure 13, it can be found that the estimated amplitudes corresponded closely with the ones calculated by FEKO. Moreover, we arranged the ends in pairs based on the extracted amplitudes. The signals to form the images of the flanks were constructed based on the arrangement and added to the original return signals. The imaging result of the reproduced data is shown in Figure 14. In Figure 14, the distributed scattering centers were easier to be distinguished from the local scattering centers.

Range(m)
The restored images of the distributed scattering centers  From Figure 13, it can be found that the estimated amplitudes corresponded closely with the ones calculated by FEKO. Moreover, we arranged the ends in pairs based on the extracted amplitudes. The signals to form the images of the flanks were constructed based on the arrangement and added to the original return signals. The imaging result of the reproduced data is shown in Figure 14. In Figure 14, the distributed scattering centers were easier to be distinguished from the local scattering centers.

The Results of the Backhoe Data and Discussion
In this subsection, we apply the proposed method to the backhoe data. The backhoe model is shown in Figure 15.
The corresponding processing results by the proposed method are presented in Figure 16b,e. The original imaging result of the backhoe with the look angles ranging from 25 • to 70 • is presented in Figure 16a. The portion within the white box in Figure 16a is shown in Figure 16d. The imaging result with the look angles ranging from −10 • to 35 • is shown in Figure 16c and the corresponding enlargement is shown in Figure 16f. In Figure 16a,d, some distributed scattering centers cannot be seen, which can be found in Figure 16c,f. The proposed algorithm was applied to the data with look angles ranging from 25 • to 70 • to restore the incomplete images of distributed scattering centers. The restored image is shown in Figure 16b,e. According to Figure 16b,e, some incomplete images of distributed scattering centers were restored.

The Results of the Backhoe Data and Discussion
In this subsection, we apply the proposed method to the backhoe data. The backhoe model is shown in Figure 15. The corresponding processing results by the proposed method are presented in Figure 16b,e. The original imaging result of the backhoe with the look angles ranging from 25°to 70° is presented in Figure 16a. The portion within the white box in Figure 16a is shown in Figure 16d. The imaging result with the look angles ranging from −10°to 35° is shown in Figure 16c and the corresponding enlargement is shown in Figure 16f. In Figure 16a,d, some distributed scattering centers cannot be seen, which can be found in Figure 16c,f. The proposed algorithm was applied to the data with look angles ranging from 25° to 70° to restore the incomplete images of distributed scattering centers. The restored image is shown in Figure 16b,e. According to Figure 16b,e, some incomplete images of distributed scattering centers were restored.
However, in the restored image, we also found that some images of the distributed scattering centers were not restored correctly. According to the inaccurate restorations, we discussed some problems and limitations of the proposed method. In the discussion part, we mainly focus on the processing result of the bucket of the backhoe, whose images are within the white boxes in Figure 16.
The black and red lines in Figure 15 mark the distributed scatterers whose images are not complete in Figure 16d. The red lines in Figure 15 denote the distributed scatterers with the correctly restored images in Figure 16e, corresponding to the areas within the red dashed boxes. The distributed scatterers that were not restored correctly were marked by the black lines on the backhoe model in Figure 15. These scatterers correspond to the areas within the black dashed boxes. The two distributed scatterers marked by the black lines on the model were inside the bucket under the condition that the look angles ranged from 25° to 70°, since we can only find the images of their ends on one side in Figure 16e. Consequently, although the algorithm detected the two ends, it cannot perform the correct pairing, since the ends that can be correctly paired with them cannot be illuminated by the signals. However, the proposed algorithm has to find the ends to be paired with them. As a result, the wrong pairing occured. Another phenomenon we did not consider in this algorithm is the interactions [10] among the scatterers. Interactions among the scatterers are very common in a complex scene, especially for the bucket of the backhoe which can be regarded as a cavity. The interactions may give rise to the images of scatterers with wrong positions. This may influence the anisotropic characteristics extraction and the pairing.

5.Conclusions
A novel method for extracting anisotropic characteristics of scattering centers and the feature enhancement in the wide-angle SAR imagery was proposed in this paper. The main target of this study was to restore the incomplete images of the distributed scattering centers in a wide-angle SAR image based on the anisotropic characteristics. In the proposed scheme, firstly, the anisotropic However, in the restored image, we also found that some images of the distributed scattering centers were not restored correctly. According to the inaccurate restorations, we discussed some problems and limitations of the proposed method. In the discussion part, we mainly focus on the processing result of the bucket of the backhoe, whose images are within the white boxes in Figure 16.
The black and red lines in Figure 15 mark the distributed scatterers whose images are not complete in Figure 16d. The red lines in Figure 15 denote the distributed scatterers with the correctly restored images in Figure 16e, corresponding to the areas within the red dashed boxes. The distributed scatterers that were not restored correctly were marked by the black lines on the backhoe model in Figure 15. These scatterers correspond to the areas within the black dashed boxes. The two distributed scatterers marked by the black lines on the model were inside the bucket under the condition that the look angles ranged from 25 • to 70 • , since we can only find the images of their ends on one side in Figure 16e. Consequently, although the algorithm detected the two ends, it cannot perform the correct pairing, since the ends that can be correctly paired with them cannot be illuminated by the signals. However, the proposed algorithm has to find the ends to be paired with them. As a result, the wrong pairing occured. Another phenomenon we did not consider in this algorithm is the interactions [10] among the scatterers. Interactions among the scatterers are very common in a complex scene, especially for the bucket of the backhoe which can be regarded as a cavity. The interactions may give rise to the images of scatterers with wrong positions. This may influence the anisotropic characteristics extraction and the pairing.

Conclusions
A novel method for extracting anisotropic characteristics of scattering centers and the feature enhancement in the wide-angle SAR imagery was proposed in this paper. The main target of this study was to restore the incomplete images of the distributed scattering centers in a wide-angle SAR image based on the anisotropic characteristics. In the proposed scheme, firstly, the anisotropic characteristics were extracted by the proposed IRWTR algorithm. Secondly, the anisotropic characteristics were used to determine if a point in the SAR image was actually the image of an end of a distributed scattering center. Thirdly, these point-like images of the ends of distributed scattering centers will be lined up in pairs based on their anisotropic characteristics extracted by IRWTR.
Author Contributions: Y.G. conceived the main idea; M.X., L.G., and Z.Z. conceived and designed the experiments; Y.G. and M.X. analyzed the data and wrote the paper.