Journal of Sensors

Volume 2015, Article ID 293141, 11 pages

http://dx.doi.org/10.1155/2015/293141

## Superresolution Polarimetric ISAR Imaging Based on 2D CP-GTD Model

College of Electronic Science and Engineering, National University of Defense Technology, Changsha, Hunan 410073, China

Received 24 December 2014; Accepted 20 April 2015

Academic Editor: Jong-Sen Lee

Copyright © 2015 Da-hai Dai et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

This paper presented a new approach to superresolution ISAR imaging based on a scattering model called coherent polarized geometrical theory of diffraction (CP-GTD) which is better matched to the physical scattering mechanism. The algorithm is a joint processing between polarization and superresolution essentially. It can also estimate the number, position, frequency dependence, span, and normalized scattering matrix of scattering centers instantaneously for each channel rather than the one which extracts parameters from each channel separately, and its performance is better than the latter because the fully polarized information is used. The superiority of the CP-GTD is verified by experiment results based on simulated and real data.

#### 1. Introduction

Since the 1990s there has been a boost in radar target recognition (RTR) based on fully polarized information. Specifically, RTR based on fully polarized high-resolution radars has been taken as a very promising approach [1, 2]. That is because that the association of high resolution and full polarization can afford more information about targets’ structure [3–6]. Backscattering from manmade objects is strongly dependent on polarization, so using only a single channel some important scatterers might be totally missed. The appropriate joint processing of the polarization channels provides higher SNR level which is required to achieve high superresolution gains and makes the superresolution algorithm more robust. Furthermore, the scattering matrix will present very useful information for RTR [7–11].

In the high-frequency region, the electric field backscattered from a complex radar target can be modeled approximately by a sum of fields scattered from some dominant scattering centers on the target. How to realize superresolution is always a key problem for ISAR [12]. There are two main approaches: one is the conventional FFT-based processing; the other is model-based processing. The common scattering center models include point-scattering model, damped exponential (DE) [5] model, and the model based on geometry theory of diffraction (GTD) [6]. The GTD model can accurately match the high-frequency electromagnetic scattering mechanism and approximates the physical mechanism of high-frequency electromagnetic scattering more closely. The existing GTD-based scattering model is already mature for singly polarized modeling. Although it was mentioned that polarization effects can be included for 1D scattering center extraction in [6], it had not explained how to do it in detail and it cannot fit for full-polarization ISAR.

The remainder of the paper is organized as follows. In Section 2, for full-polarization ISAR system, a new model, called 2D coherent polarization GTD (CP-GTD) model, which is capable of accurately describing the high-frequency polarization scattering mechanism of the radar target, is presented. In Section 3, a novel method based on CP-GTD model for full-polarization scattering center extraction and parameter estimation is proposed in detail. In Section 4, the validity is proved by the experimental results based on simulated and real data. Some conclusions are given in Section 5.

#### 2. 2D CP-GTD Signal Model for Fully Polarized ISAR

A fully polarized radar measures the backscattered electric field with two linearly independent transmit and receive polarizations also called dual-channel orthogonal polarization measurements. Without regard to the translational motion, the raw data of ISAR can be represented by a 2D signal of frequency-angular domain. Suppose horizontal and vertical linear polarization as the two linearly independent transmit and receive polarizations. We adopt a full-polarization parametric scattering model based on geometrical theory of diffraction (GTD) which is called 2D coherent polarization GTD (CP-GTD) model in this paper.

Taking the coherent stepped-frequency full-polarization radar as an example, fully polarized ISAR can measure a set of scattering matrices at a set of stepped frequencies at an angular domain. The backscattered signals from scattering centers are then given byThat is,where is the incident signal; is the backscattered signal; is the scattering matrix of the th scattering center of the radar target; is the wavenumber; is the wavelength; is the frequency; is the aspect angle. is a frequency dependence factor of the scattering center which characterizes the geometry and hence is named type factor as well. According to the geometrical theory of diffraction, the type factor is an integer multiple of 1/2. For canonical scattering geometries, the type parameters are as follows: is for flat plate at broadside or dihedral; is for singly curved surface reflection; is for point scatterer or doubly curved surface reflection or straight edge specular; is for edge diffraction; is for corner diffraction.

Thus, the raw data of every polarization channel can be represented aswhere, or , is the received noise.

Generally, the aspect angle is small; varies very little. Here suppose it is invariable in the imaging processing. Then,

Taking Stolt interpolation the data in frequency-angular domain [6], let

After the interpolation, (4) can be written in the Cartesian spectrum domain as

Rewrite (6) in the matrix form aswhere is a matrix corresponding to the full-polarization measurement information; is a matrix corresponding to the scattering matrix information of the scattering centers of radar target; is a matrix; is a measurement noise matrix. Their concrete expressing forms are given aswhere denotes the single-polarization measurement information; and denote the vectors of the single-polarization scattering information and the measurement noise, respectively; is the steering vector corresponding to the th scattering center. Their detailed expressions are

What is worthy to point out is that the CP-GTD model is not just a simple extension of the traditional GTD model. It is distinguished from the existing models by considering that each scattering center has an exclusive coherent polarization scattering matrix which can model full-polarization scattering as a whole. Its advantage will be verified in the following two sections.

#### 3. Superresolution and Feature Extraction for Fully Polarized ISAR

In this section, the method of superresolution and feature extraction for fully polarized ISAR based on 2D CP-GTD model will be proposed using a modified 2D Multiple Signal Classification (MUSIC) algorithm [13].

In this modified MUSIC algorithm, using the raw data of full-polarization ISAR, a fully polarized correlation matrix can be obtained which is defined as , where superscript “” denotes conjugate transpose.

Make eigenvector-based decomposition of :

The eigenvectors of the correlation matrix are divided into two subspaces: the signal subspace and the noise subspace .

Define two-dimensional full-polarization MUSIC spectrum for a radar target as

Theoretically speaking, through 3D search, the position () and the type factor of each scattering center can be estimated. Since the CP-GTD model is adopted here, the estimations for the position and the type are consistent in every polarization channel, which is very useful because it can guarantee the coherence of different polarization channels.

In order to decrease the computation burden in 3D search, let and estimate the position of the scattering centers firstly:

When searching , a roughly searching method can be used firstly and thereby the accurate estimation by a local searching method based on the rough estimation is suggested.

Once has been derived, searching is only needed for the discrete values of the frequency-dependent factors. The signal steering vector’s projection on the signal subspace will reach the maximum when the estimation of is consistent with the true value at the right range; namely,

After the position and type estimation of all scattering centers, the estimation of can be derived by (9). Then the estimation of matrix can be obtained using least square (LS) criterion as follows:

Rearrange th row of as , with as the fully polarized scattering matrix of the th scattering center. It can be figured out that the span (total power) and normalized scattering matrix are

The other fully polarized characteristic just like polarization invariables can also be calculated for RTR [10], where denotes the phase angle.

In fact, although the algorithm is studied aiming at the fully polarized stepping frequency signal, it is also applicable for linear frequency modulated (LFM) signal, and it is compatible with the case of single or double polarization processing (they are just some simple forms of the CP-GTD method). It is a joint processing of polarization and superresolution. It can estimate the number, type factor, position, span, and normalized scattering matrix of scattering centers instantaneously for each channel rather than extracting parameters from every channel separately. Its performance is better than the latter because the fully polarized information is used. The number and the location of scattering centers in each channel are consistent. This characteristic saves complexity in the scattering centers association in different polarizations.

Before applying the modified 2D MUSIC algorithm presented above, the number of the scattering centers and the fully polarized correlation matrix must be estimated. If the estimation is less than the true value, some real scattering center will be lost. Contrarily, if the estimation is less than the true value, some false scattering centers will appear because of noise. Generally, when the SNR (signal-to-noise ratio) is high enough, the number of scattering centers can be estimated by observing the eigenvalues of . However, this method is subjective and sensitive to the SNR. The Gerschgorin Disk method is adopted because it is not sensitive to the SNR which is fit for the situation in this paper [14].

To estimate the fully polarized correlation matrix , a preprocessing scheme known as SSP (spatial smoothing preprocessing) can be used. For idea point scattering condition, Odendaal et al. [15] applied the MSSP (Modified SSP) to estimate the singly polarized correlation matrix and Kim et al. [5] extended and generalized it in order to estimate the fully polarized correlation matrix. However, considering that the MSSP will disturb the estimation of type, SSP is still adopted here.

A set of overlapping subarrays with dimensions is defined as

Applying the full-polarization SSP, the fully polarized correlation matrix is obtained aswhere and . In order to obtain an exclusive LS solution of in (16), we should make sure that is not singular; that is to say, .

To sum up, the flow chart on the CP-GTD algorithm is as follows.

*Step 1. *Estimate the full-polarization covariance matrix. Choose the forward smoothing algorithm (see (19)).

*Step 2. *With the eigendecomposition of the full-polarization covariance matrix, use the Gerschgorin Disks to estimate the number of the scattering centers (see (12)).

*Step 3. *Estimate the ranges of the scattering centers using the full-polarization MUSIC spectrum search (see (14)).

*Step 4. *Estimate the frequency-dependent factor of each scattering center (see (15)).

*Step 5. *The estimation of matrix can be obtained using least square (LS) criterion (see (16)).

*Step 6. *Calculate the span and normalized scattering matrix (see (17)).

*Step 7. *Classify the attribute of each scattering center and apply the result to target recognition.

#### 4. Performance Analysis

To evaluate the performance of the CP-GTD method, we choose the conventional FFT method and the singly polarized GTD method as comparisons. The last two kinds of method present the progress of extracting the scattering centers in each polarization channel and, then, associating the scattering centers in multiple polarization channels and gathering the associated scattering coefficient in each channel together as the estimation of scattering matrix of the scattering center.

##### 4.1. Evaluation Criteria

The evaluation indexes are classified into two types: one is used to evaluate accuracy of the noncoherent parameters, including the number, type, range, and span of scattering centers. Here assuming that both the traditional GTD method and CP-GTD method are capable of estimating the number of scatterers correctly, we choose the mean square roots of the last three estimated errors as the evaluation indicators corresponding to the estimating accuracy of the noncoherent parameters. The criterion is defined by

Similarly, and , where denotes the assembly average.

The other type is used to evaluate accuracy of the coherent parameters. As we all know, the scattering matrix is a complex matrix, and the phase of each element in the scattering matrix, especially the relative phase difference, is the important information which exactly reflects the polarimetric scattering mechanism of the radar target. Under the condition of symmetry, the coherent polarization scattering matrix includes three relative phases, namely, , and (where is the phase difference of the elements and in the scattering matrix). Choose the mean square roots of the above three difference phases’ estimated errors as the evaluation criteria corresponding to the estimating accuracy of the coherent parameters; namely,

and are defined similarly.

##### 4.2. Simulation Results

An aircraft type target whose size was about 6 m × 8 m defined by a set of eight point scatterers as shown in Figure 3 was used to generate raw data of ISAR over a frequency-angular domain. The initial frequency is 8.75 GHz; the frequency step is 10 MHz, the number of frequency points is 32; the total angle is ; , ; the Monte-Carlo number is 1000; and the signal-to-noise Ratio (SNR) is 25 dB. The SNR is defined as the ratio of the target span and noise intensity; namely, . The results of full-polarization ISAR are given in Figure 1.