Open AccessArticle

A Method for Estimating the Polarimetric Scattering Matrix of Moving Target for Simultaneous Fully Polarimetric Radar

Fulai Wang

Chao Li

Chen Pang

Yongzhen Li

and

Xuesong Wang

College of Electronic Science, National University of Defense Technology (NUDT), Changsha 410073, China

Author to whom correspondence should be addressed.

Sensors 2018, 18(5), 1418; https://doi.org/10.3390/s18051418

Submission received: 22 March 2018 / Revised: 28 April 2018 / Accepted: 30 April 2018 / Published: 3 May 2018

(This article belongs to the Section Remote Sensors)

Download

Browse Figures

Versions Notes

Abstract

To precisely obtain the polarimetric scattering matrix (PSM) of moving target, a measurement model for the simultaneous fully polarimetric radar is formulated. The calibration errors and isolation of the transmitted waveforms are considered. To address the decline in performance of the traditional PSM estimation methods when the target moves, a novel method with measurement selection is proposed. Numerical experiments are conducted to demonstrate and validate the superiority of the proposed method, especially for the PSM estimation of the target with non-uniform motion.

Keywords:

polarization scattering matrix; moving target; simultaneous polarimetric radar; non-uniform motion

1. Introduction

Advancements in radar technology and theory have provided a better understanding of the polarimetric information contained in radar targets [1,2,3,4]. The polarimetric features, which can be described by a second order polarimetric scattering matrix (PSM) have been widely used in various fields, such as terrain observation, disaster surveillance and atmospheric remote sensing. To accurately obtain the PSM, two fully polarimetric measurement schemes, called the alternately transmitting and simultaneously receiving (ATSR) scheme and the simultaneously transmitting and simultaneously receiving (STSR) scheme, have been widely investigated since the 1980s [5,6,7]. The ATSR radar alternately transmits waveforms through horizontal (H) and vertical (V) polarizations while both polarizations are received simultaneously on reception. At least two pulses are required in this mode to obtain the four elements for the PSM. The ATSR in essence is a time-sharing polarimetric radar, hence, the target decorrelation may influence the measurements results. In contrast, for STSR radar, the two orthogonal polarization states are transmitted and received simultaneously. Thus, the PSM of the targets can be retrieved within one pulse recurrent time (PRT). In this case, the limitation caused due to the change of the transmitted polarization states in ATSR scheme can be overcome [8].

Estimation of the target scattering matrix is based on the measurements from the fully polarimetric radar. For the static target, its echoes are coherent, and the pulse-integration (PI) method can be used to estimate the target’s PSM. For moving targets, the Doppler velocity which refers to the radial velocity, and the PSM of the target are coupled. It is unclear whether the change of the radar-echo phase is caused by the target’s displacement or the PSM. Therefore, the true PSM of the moving target is usually difficult to obtain [9,10,11]. Fortunately, in real radar applications, the PSM of a slow-moving target, such as that for an unmanned aerial vehicle (UAV), is assumed to be deterministic over the observation duration [12]. Under this assumption, provided that the target Doppler velocity is known, the phase changes in the radar echoes due to the target motion can be compensated by the estimated velocity. After the compensation, the incoherent echoes become coherent, and the PSM of moving target can be estimated by the pulse-compensation (PC) method. Obviously, the estimated accuracy of the PC method is related to the accuracy of the estimated Doppler velocity. When the velocity cannot be determined exactly, the estimated PSM becomes inaccurate [13].

To obtain the precise PSM of a moving target, a method of measurement selection (MS) is proposed in this paper. Using the criterion based on the signal to noise ratio (SNR) of the integration echoes, partial measurements are selected to estimate the PSM of the target. After the selection, the influences of the target motion on the four polarization channels can be considered uniform, and the PSM with relative amplitude and phase can be estimated. The advantage of this method is that the MS does not require any prior information about the target velocity, and it can still be used even if the target exhibits non-uniform motion. The rest of this paper is organized as follows: Section 2 presents the problem formulation; in Section 3, three PSM estimation methods are introduced; the numerical simulations are provided in Section 4 to verify the performance of the proposed method, followed by the conclusions in Section 5.

N o t a t i o n

: In this paper, it is assumed that a lower-case letter (e.g., a) denotes a scalar; a boldface lowercase letter (e.g.,

a

) denotes a vector; and a boldface uppercase letter (e.g.,

A

) indicates a matrix. Additionally,

A^{T}

and

A^{H}

denote the transpose and the conjugate transpose of the matrix

A

, and the symbol

|\cdot|

denotes the modulus of a complex number.

2. Signal Model for Moving Target in STSR Radar

The simplified signal processing flow chart of the STSR radar is depicted in Figure 1 [8]. Suppose the STSR radar transmits a pair of opposite (up-going and down-going) slope of linear frequency-modulation (LFM) waveforms, which can be expressed as

\{\begin{matrix} s_{H} (t) = rect (\frac{t}{T_{r}}) \exp (j 2 π f_{0} t + j π γ t^{2}) \\ s_{V} (t) = rect (\frac{t}{T_{r}}) \exp (j 2 π f_{0} t - j π γ t^{2}) \end{matrix}; t \in [- \frac{T_{r}}{2}, \frac{T_{r}}{2}]

(1)

where

rect (\frac{t}{T_{r}}) = \{\begin{matrix} 1, |t| \leq \frac{T_{r}}{2} \\ 0, else \end{matrix}

(2)

and

T_{r}

is the pulse duration.

f_{0}

is the carrier frequency,

γ

is the modulation slope, and the radar bandwidth is

B = γ T_{r}

. To facilitate the discussion, the transmitted waveforms can be given in vector form as

s (t) = {[s_{H} (t), s_{V} (t)]}^{T}

(3)

For a point target, the received signals are the time-delayed version of the transmitted signals. Thus, for the

k th

pulse, the received signal is [14]

r_{k} (t) = [\begin{matrix} r_{H, k} (t) \\ r_{V, k} (t) \end{matrix}] = R^{T} \cdot S_{k} \cdot T \cdot s (t - \frac{2 d_{k}}{c}) + w_{k} (t); k = 1, 2, \dots, N

(4)

where N is the number of pulses in a coherent process interval (CPI), c is the speed of the light,

w_{k} (t) = {[w_{H, k} (t), w_{V, k} (t)]}^{T}

is the thermal noise of the receivers with variance

σ_{w}^{2}

d_{k}

is the radial distance from the point target to the radar, and

S_{k}

is the PSM of the target that can be described as

S_{k} = [\begin{matrix} S_{H H, k} S_{H V, k} \\ S_{V H, k} S_{V V, k} \end{matrix}]

(5)

where the corner marks

H H

H V

V H

and

V V

denote four polarization channels.

R

and

T

represent the effect of channels and antennas on the PSM during the reception and transmission, respectively. Here, the calibration errors of the STSR radar system, including the cross-polarization isolation of the antennas, the amplitude and phase difference of the channels, are considered. Suppose the same antennas and channels are used during the transmission and reception, then the

R

and

T

can be set as follows

R = T = [\begin{matrix} 1 α_{H V} \exp (j ϕ_{H V}) \\ α_{V H} \exp (j ϕ_{V H}) α_{V V} \exp (j ϕ_{V V}) \end{matrix}]

(6)

Processed by the mixer and matched filters shown in Figure 1, the high-resolution range profile (HRRP) can be obtained

G_{k} (t) = R^{T} S T [\begin{matrix} g_{H H, k} (t) g_{H V, k} (t) \\ g_{V H, k} (t) g_{V V, k} (t) \end{matrix}] + [\begin{matrix} w_{H H, k} (t) w_{H V, k} (t) \\ w_{V H, k} (t) w_{V V, k} (t) \end{matrix}]

(7)

and

g_{H H, k} (t) = A \cdot \sin c (π B (t - \frac{2 d_{k}}{c})) \cdot \exp (- j \frac{4 π d_{k}}{λ})

(8)

g_{V V, k} (t) = A \cdot \sin c (π B (t - \frac{2 d_{k}}{c})) \cdot \exp (- j \frac{4 π d_{k}}{λ})

(9)

where

A = B T_{r}

is the gain of the matched filtering (MF),

λ

is the wavelength, and

\sin c (\cdot)

is the Sinc function. One thing should be pointed out is that, generally, the scattering matrix is related to the shape, geometrical structure, reflectivity and orientation of the target. Meanwhile, the PSM may fluctuate whether or not the target is in motion. Fortunately, this effect can be controlled by limiting the CPI of the radar system. If the CPI is short enough, the PSM of the target can be assumed to be same for different pulses. Thus,

S

is used to replace

S_{k}

in Equation (7). With

t = 2 d_{k} / c

, Equations (8) and (9) can be written into

g_{H H, k} = A \cdot \exp (- j \frac{4 π d_{k}}{λ})

(10)

g_{V V, k} = A \cdot \exp (- j \frac{4 π d_{k}}{λ})

(11)

Additionally, in [15], the author has pointed out that the isolation of opposite slope of LFM waveforms, which is defined as

I = \max_{\forall t} \{20 \log_{10} \frac{|g_{H V, k} (t)|}{\max_{\forall t^{'}} |g_{H H, k} (t^{'})|}\}

(12)

is related to the time-bandwidth product of the waveforms and an approximate equation is given as follows

I \approx - 10 \log_{10} (B T_{r})

(13)

Obviously, when the time-bandwidth product is large enough, such as

B T_{r} = 10^{4}

, the I is equal to

- 40 dB

approximately, which means for arbitrary

t \in [- T_{r}, T_{r}]

, the modulus of

g_{H V, k} (t)

and

g_{V H, k} (t)

are much less than that of

g_{H H, k}

and

g_{V V, k}

. Therefore, when

t = 2 d_{k} / c

, Equation (7) can be further rewritten as

G_{k} = [\begin{matrix} G_{H H, k} G_{H V, k} \\ G_{V H, k} G_{V V, k} \end{matrix}] = R^{T} S T [\begin{matrix} g_{H H, k} I^{'} \\ I^{'} g_{V V, k} \end{matrix}] + [\begin{matrix} w_{H H, k} w_{H V, k} \\ w_{V H, k} w_{V V, k} \end{matrix}]

(14)

where

I^{'} = 10^{I / 20}

. If the second order small quantities are ignored, it comes

\{\begin{matrix} G_{H H, k} = g_{H H, k} (S_{H H} + α_{V H} \exp (j ϕ_{V H}) (S_{V H} + S_{H V})) + I^{'} \cdot α_{V V} S_{H V} \exp (j ϕ_{V V}) + w_{H H, k} \\ G_{V H, k} = g_{H H, k} (S_{H H} α_{H V} \exp (j ϕ_{H V}) + α_{V V} (S_{V H} \exp (j ϕ_{V V}) + S_{V V} α_{V H} \exp (j (ϕ_{V V} + ϕ_{V H})))) \\ + I^{'} \cdot S_{V V} α_{V V}^{2} \exp (2 j ϕ_{V V}) + w_{V H, k} \\ G_{H V, k} = g_{V V, k} (S_{H H} α_{H V} \exp (j ϕ_{H V}) + α_{V V} \exp (j ϕ_{V V}) (S_{H V} + α_{V H} S_{V V} \exp (j ϕ_{V H}))) \\ + I^{'} \cdot S_{H H} + w_{H V, k} \\ G_{V V, k} = g_{V V, k} α_{V V} (α_{H V} \exp (j (ϕ_{V V} + ϕ_{H V})) (S_{V H} + S_{H V}) + S_{V V} α_{V V} \exp (2 j ϕ_{V V})) \\ + I^{'} \cdot S_{V H} α_{V V} \exp (j ϕ_{V V}) + w_{V V, k} \end{matrix}

(15)

Our goal is to estimate the PSM of moving target from the

G_{k}

that includes the target echo and thermal noise of the receivers. As mentioned before, the variance of the thermal noise is assumed to be

σ_{w}^{2}

. Furthermore, the thermal noise is supposed to follow the Gaussian distribution. In [16], it has been proved that the noise components in the MF output, which are

w_{H H, k}

w_{V H, k}

w_{H V, k}

and

w_{V V, k}

, follow the Gaussian distribution similarly and the variance is

T_{r} σ_{w}^{2}

. To obtain the real PSM of the target, three estimation methods are analyzed in the next Section.

3. PSM Estimation for the Moving Target

Instead of estimating the real PSM of the target, the PSM with relative amplitude and phase is estimated. With

S_{H H}

as a reference, the relative PSM can be expressed as

\bar{S} = [\begin{matrix} 1 \frac{|S_{H V}|}{|S_{H H}|} \exp (j ϕ_{H V}) \\ \frac{|S_{V H}|}{|S_{H H}|} \exp (j ϕ_{V H}) \frac{|S_{V V}|}{|S_{H H}|} \exp (j ϕ_{V V}) \end{matrix}] = [\begin{matrix} 1 S_{H V / H H} \\ S_{V H / H H} S_{V V / H H} \end{matrix}]

(16)

where the

ϕ_{V H}

ϕ_{H V}

and

ϕ_{V V}

are the phases of different polarized channels, and the

S_{V H / H H}

S_{H V / H H}

and

S_{V V / H H}

are the normalized scattering matrix parameters. To facilitate the following discussion, some notations are defined as follows

\begin{matrix} A_{H H} = S_{H H} + α_{V H} \exp (j ϕ_{V H}) (S_{V H} + S_{H V}) \\ A_{V H} = S_{H H} α_{H V} \exp (j ϕ_{H V}) + α_{V V} (S_{V H} \exp (j ϕ_{V V}) + S_{V V} α_{V H} \exp (j (ϕ_{V V} + ϕ_{V H}))) \\ A_{H V} = S_{H H} α_{H V} \exp (j ϕ_{H V}) + α_{V V} \exp (j ϕ_{V V}) (S_{H V} + α_{V H} S_{V V} \exp (j ϕ_{V H})) \\ A_{V V} = α_{V V} (α_{H V} \exp (j (ϕ_{V V} + ϕ_{H V})) (S_{V H} + S_{H V}) + S_{V V} α_{V V} \exp (2 j ϕ_{V V})) \end{matrix}

(17)

and

\begin{matrix} B_{H H} = I^{'} \cdot α_{V V} S_{H V} \exp (j ϕ_{V V}) \\ B_{V H} = I^{'} \cdot S_{V V} α_{V V}^{2} \exp (2 j ϕ_{V V}) \\ B_{H V} = I^{'} \cdot S_{H H} \\ B_{V V} = I^{'} \cdot S_{V H} α_{V V} \exp (j ϕ_{V V}) \end{matrix}

(18)

Moreover, introducing the notation,

x_{p q, k} = A_{c} (A_{p q} \exp (j \frac{4 π v_{0} k T_{PRT}}{λ}) + \frac{B_{p q}}{A} \exp (j \frac{4 π d_{0}}{λ})); p, q = H, V

(19)

where

A_{c} = A \exp (- j 4 π d_{0} / λ)

d_{0}

is the target initial distance,

v_{0}

is the radial velocity and

T_{PRT}

is the pulse repetition time. In the rest of the article,

p, q

are used to denoted

H, V

. Then Equation (15) can be rewritten as

G_{p q, k} = x_{p q, k} + w_{p q, k}

(20)

For a static target,

x_{p q, 1} = x_{p q, 2} = \dots = x_{p q, N}

. Using the PI method,

S_{p q}

can be estimated by

{\hat{S}}_{p q} = \sum_{k = 1}^{N} (G_{p q, k} / N)

, and the estimation of the relative PSM becomes

{\hat{S}}_{p q / H H_PI} = \frac{\sum_{k = 1}^{N} G_{p q, k}}{\sum_{k = 1}^{N} G_{H H, k}}

(21)

However, when the target moves, its echoes are incoherent. That means the motion of the target has impact on the measurements of the PSM. Generally, the velocity of the target can be estimated by the relative phase change of the echoes. In [17], the standard deviation of the estimation error is given as follows

σ_{v} = \frac{λ}{2} \sqrt{6 / [{(2 π)}^{2} χ {(N T_{PRT})}^{2}]}

(22)

where

χ

is the SNR of the input signal. With the estimated velocity

{\hat{v}}_{0}

, the PC method can be used to estimate the parameters of the scattering matrix by

{\hat{S}}_{p q} = \sum_{k = 1}^{N} (\frac{G_{p q, k}}{A_{c} N} \cdot \exp (- j \frac{4 π {\hat{v}}_{0} k T_{PRT}}{λ}))

(23)

Therefore, the estimation of the relative PSM is

{\hat{S}}_{p q / H H_PC} = \frac{\sum_{k = 1}^{N} G_{p q, k} \exp (- j \frac{4 π {\hat{v}}_{0} k T_{PRT}}{λ})}{\sum_{k = 1}^{N} G_{H H, k} \exp (- j \frac{4 π {\hat{v}}_{0} k T_{PRT}}{λ})}

(24)

It can be observed from Equations (21) and (24) that the PI method is a special case of the PC method. When the velocity of the target reduces to

0 m / s

, PC has the same expression as PI. Additionally, for the PC method, accurate estimation of the PSM requires that the velocity of the target is estimated precisely. If the estimation error is large, the phase caused by the target’s motion cannot be compensated, leading to the inaccurate estimation of the PSM. Compared with the PC method, the estimation of the velocity is avoided in the MS method. The term

m_{p q, k}

is added to indicate the selection state of the measurements

m_{p q, k} = \{\begin{matrix} 1, selected \\ 0, unselected \end{matrix}, k = 1, 2, \dots, N

(25)

The integration results for different channels in MS method are

G_{p q} = X_{p q} + W_{q}

(26)

where

\begin{matrix} G_{p q} = \sum_{k = 1}^{N} m_{p q, k} G_{p q, k} \\ X_{p q} = \sum_{k = 1}^{N} m_{p q, k} x_{p q, k} \\ W_{q} = \sum_{k = 1}^{N} m_{p q, k} w_{p q, k} \end{matrix}

(27)

For arbitrary channel, the SNR of the integration terms is defined as:

{SNR}_{p q}^{integration} = \frac{E (X_{p q} X_{p q}^{*})}{T_{r} σ_{w}^{2} \sum_{k = 1}^{N} m_{p q, k}} = \frac{E (\sum_{k = 1}^{N} m_{p q, k} x_{p q, k} \sum_{k = 1}^{N} m_{p q, k} {(G_{p q, k} - w_{p q, k})}^{*})}{T_{r} σ_{w}^{2} \sum_{k = 1}^{N} m_{p q, k}}

(28)

where

E (\cdot)

represents the mathematical expectation. Since the thermal noise is assumed to follow the Gaussian distribution, and the mean and the variance are zero and

T_{r} σ_{w}^{2}

, respectively, Equation (28) can be simplified to

\begin{matrix} {SNR}_{p q}^{integration} & = \frac{E (\sum_{k = 1}^{N} m_{p q, k} x_{p q, k} \sum_{k = 1}^{N} m_{p q, k} G_{p q, k}^{*})}{T_{r} σ_{w}^{2} \sum_{k = 1}^{N} m_{p q, k}} \\ = \frac{E (\sum_{k = 1}^{N} m_{p q, k} (G_{p q, k} - w_{p q, k}) \sum_{k = 1}^{N} m_{p q, k} G_{p q, k}^{*})}{T_{r} σ_{w}^{2} \sum_{k = 1}^{N} m_{p q, k}} \\ = \frac{E (\sum_{k = 1}^{N} m_{p q, k} G_{p q, k} \sum_{k = 1}^{N} m_{p q, k} G_{p q, k}^{*} - \sum_{k = 1}^{N} m_{p q, k} w_{p q, k} {\sum_{k = 1}^{N} m_{p q, k} (x_{p q, k} + w_{p q, k})}^{*})}{T_{r} σ_{w}^{2} \sum_{k = 1}^{N} m_{p q, k}} \end{matrix}

(29)

Since for N observations, the measurement results are certain. Equation (29) can be further simplified to

\begin{matrix} {SNR}_{p q}^{integration} & = \frac{{|G_{p q}|}^{2} - E (\sum_{k = 1}^{N} m_{p q, k} w_{p q, k} \sum_{k = 1}^{N} m_{p q, k} w_{p q, k}^{*}) - E (\sum_{k = 1}^{N} m_{p q, k} w_{p q, k} \sum_{k = 1}^{N} m_{p q, k} x_{p q, k}^{*})}{T_{r} σ_{w}^{2} \sum_{k = 1}^{N} m_{p q, k}} \\ = \frac{{|G_{p q}|}^{2}}{T_{r} σ_{w}^{2} \sum_{k = 1}^{N} m_{p q, k}} - 1 \end{matrix}

(30)

It is clear that to improve the estimation performance, the SNR of the integration terms should be as high as possible. Besides, it should be pointed out that the term

m_{p q, k}

cannot be all zero. When

m_{p q, 1} = m_{p q, 2} = \dots = m_{p q, N} = 0

, the term

G_{p q} = 0

, and the SNR in (30) becomes meaningless. Another thing should be noticed is that to ensure the phase consistency of the selected measurements from different channels, the terms

m_{H H, k}

m_{V H, k}

m_{H V, k}

and

m_{V V, k}

should be equal, which means

m_{H H, 1} = m_{V H, 1} = m_{H V, 1} = m_{V V, 1}, \dots \dots, m_{H H, N} = m_{V H, N} = m_{H V, N} = m_{V V, N}

. Then the target term

X_{p q}

can be estimated by the observation term

G_{p q}

through minimizing the reciprocal of the first term of the SNR in Equation (30), and the criterion of the measurements selection can be expressed as

\{\begin{matrix} {\hat{X}}_{H H} = G_{H H}; {\hat{X}}_{V H} = G_{V H}; {\hat{X}}_{H V} = G_{H V}; {\hat{X}}_{V V} = G_{V V} \\ \min_{m_{H H, k}} \{\frac{T_{r} σ_{w}^{2} \sum_{k = 1}^{N} m_{H H, k}}{{|G_{H H}|}^{2}}, \frac{T_{r} σ_{w}^{2} \sum_{k = 1}^{N} m_{H H, k}}{{|G_{V H}|}^{2}}, \frac{T_{r} σ_{w}^{2} \sum_{k = 1}^{N} m_{H H, k}}{{|G_{H V}|}^{2}}, \frac{T_{r} σ_{w}^{2} \sum_{k = 1}^{N} m_{H H, k}}{{|G_{V V}|}^{2}}\} \end{matrix}

(31)

The sequence of

m_{p q, 1}, m_{p q, 2}, \dots, m_{p q, N}

has

2^{N}

combinations. The combination, which makes the SNR in Equation (30) maximum, can be obtained by enumerating. For instance, for

N = 3

, the value space of the sequence

m_{p q, 1}, m_{p q, 2}, \dots, m_{p q, N}

[\begin{matrix} m_{p q, 1} \\ m_{p q, 2} \\ m_{p q, 3} \end{matrix}] \in \{[\begin{matrix} 1 \\ 0 \\ 0 \end{matrix}], [\begin{matrix} 0 \\ 1 \\ 0 \end{matrix}], [\begin{matrix} 0 \\ 0 \\ 1 \end{matrix}], [\begin{matrix} 1 \\ 1 \\ 0 \end{matrix}], [\begin{matrix} 0 \\ 1 \\ 1 \end{matrix}], [\begin{matrix} 1 \\ 0 \\ 1 \end{matrix}], [\begin{matrix} 1 \\ 1 \\ 1 \end{matrix}]\}

(32)

By substituting each element of the value space into Equation (30) and calculating the

{SNR}_{p q}^{integration}

of each channel, then the combination which satisfies the criterion in Equation (31) can be obtained. With

x_{p q, k} = S_{p q / H H} x_{H H, k}

, the target integration results for the

p q

and

H H

channels satisfy

X_{p q} = \sum_{k = 1}^{N} m_{p q, k} x_{p q, k} = \sum_{k = 1}^{N} m_{p q, k} S_{p q / H H} x_{H H, k} = S_{p q / H H} X_{H H}

(33)

Thus, the estimation of

S_{p q / H H}

{\hat{S}}_{p q / H H_MS} = \frac{{\hat{X}}_{p q}}{{\hat{X}}_{H H}} = \frac{G_{p q}}{G_{H H}}

(34)

Based on Equation (34), the estimated target’s PSM is

\hat{\bar{S}} = [\begin{matrix} 1 \frac{G_{H V}}{G_{H H}} \\ \frac{G_{V H}}{G_{H H}} \frac{G_{V V}}{G_{H H}} \end{matrix}]

(35)

As is shown in Equation (35), no prior information are required in the MS method. Therefore, the PSM estimation of the target in non-uniform motion can also be solved. Performance of these three methods are analyzed in the next section.

4. Simulation Results and Discussion

In this section, the proposed method is evaluated using simulations. Subsequently, some metrics are defined herein. The polarized correlation coefficient (PCC) is selected to illustrate the performance of different methods. It can be defined as [18,19]

P C C = \frac{|{S^{'}}^{H} \cdot \hat{S^{'}}|}{\sqrt{{S^{'}}^{H} S^{'} {\hat{S^{'}}}^{H} \hat{S^{'}}}}

(36)

where

S^{'}

is the vector form of the relative PSM which is

S^{'} = {[1 S_{V H / H H} S_{H V / H H} S_{V V / H H}]}^{T}

(37)

and

{\hat{S}}^{'}

is the vector form of the estimated target’s PSM

\hat{\bar{S}}

P C C

indicates the similarity between two PSMs. In addition, based on Equations (19) and (20), the SNR of

p q

polarized channels can be defined as

{SNR}_{p q} = \frac{{|A_{c} A_{p q}|}^{2}}{T_{r} σ_{w}^{2}}

(38)

and the average SNR of the four polarized channels is

SN R_{a v} = \frac{SN R_{H H} + SN R_{V H} + SN R_{H V} + SN R_{V V}}{4}

(39)

4.1. PSM Estimation without System Errors

In this subsection, the system errors, including the cross-polarization isolation of the antennas, and the amplitude and phase difference of the channels, are supposed to be calibrated. Thus, the matrices

R

and

T

can be set to

R = T = [\begin{matrix} 1 0 \\ 0 1 \end{matrix}]

(40)

Radar and target parameters are shown in Table 1. To analyze the performance of different methods, a point target was simulated with a

{SNR}_{a v}

ranging from

- 10 dB

30 dB

. For each polarized channel, the peak of the target’s HRRP is used to compose the measured PSM. The velocity of the target can be obtained by adding a sinusoidal shift to the discrete time. The trajectories of the target in this paper do not represent the real, while they can be used to investigate the robustness of the method. The motion of the target in a particular Monte Carlo (MC) trial is

v_{u} (k) = v_{u, 0} + v_{u, 0}^{2} \sin (\frac{π k T_{P R T}}{P})

(41)

where k is the pulse number, the radial velocity

v_{u, 0}

is a uniform distribution on

[100 m / s, 150 m / s]

, which means the Doppler frequency belongs to

[6.67 kHz, 10 kHz]

, u is the index of an arbitrary MC trial, and the control factor

P = 10

. It can be observed from Equation (41) that when

P ≫ π k T_{P R T}

(e.g.,

P = 500

), the target’s motion becomes uniform. Besides, as this paper concentrates on the PSM estimation of moving targets, target detection is not discussed. With the simulated trajectories, measurements from a moving target can be obtained by the equations shown in Table 2, and the PSM estimation methods are carried out.

Utilizing the parameters shown in Table 1, 500 MC trials have been done at each SNR. The means and standard deviations (STD) of the

P C C

are plotted as curves and error bar, respectively. Figure 2 shows the performance of these three PSM estimation methods when the target experiences uniform motion. Since the target is not static, the echoes are not incoherent. Thus, direct sum of the echoes cannot obtain effective accumulation, leading to the worst performance of the PI method. For the PC and MS method, the

P C C

means increase to 1 with the increase of the

{SNR}_{a v}

, and the STD decreases owing to the same factor. As is mentioned before, the PC method uses the estimated velocity to compensate the phases of the measurements. However, it can be observed from Equation (22) that the accuracy of the velocity estimation is influenced by the SNR, and random errors in the estimation are inevitable. With a potentially inaccurate estimation, the phase generated by the target motion cannot be compensated completely, leading to the degradation in the PC results. By using partial measurements selected, the MS method gets the best performance. The essential reason is that by measurements selection, the echoes with similar phase are chosen, which means the selection coefficient

m_{p q, k} = 1

. The phases of these selected measurements are approximatively coherent, leading to the best performance of the MS method.

Furthermore, PSM estimation for a non-uniformly moving target is also analyzed. Figure 3 clearly demonstrates the performance of the PC method is highly affected by the target’s motion. For the non-uniform case, it is difficult to estimate the target velocity accurately, causing the incomplete compensation of the measurements. The performance of PC method significantly deteriorates compared with the results shown in Figure 2. Compared with PC, the MS method does not require any prior information about the target motion and thus it can be valid for the non-uniform motion.

Besides, as mentioned in Section 2, in order to assure the constant target PSM, the CPI should be limited. Thus, the analysis of the proposed method for different integrated pulses N is given as follows. Here, the control factor is set to be

P = 10

and other parameters are same as those in Table 1. It can be observed from Figure 4 that the number of N has almost no influence on the properties of these three methods for estimating the moving target PSM. The proposed method still performs better than the other two methods when the CPI is short, and the performance of the MS method is not sensitive to the number of N.

4.2. PSM Estimation with System Errors

In this subsection, the influence of the system errors on the PSM estimation is simulated. With the system errors of the well-known PARSAX radar system as a reference, the matrices are set as follows [20]

R = T = [\begin{matrix} 1 0.05 \exp (j \frac{π}{180}) \\ 0.05 \exp (j \frac{π}{180}) 1.05 \exp (j \frac{π}{180}) \end{matrix}]

(42)

Other simulation parameters are same as those in Table 1. Similarly, the PSMs of targets with different motion states are estimated and the results are shown in Figure 5 and Figure 6. The performance of the MS method is the best among these three PSM estimation methods under two kinds of motion states, and the reason has been given in the last subsection. Here, the difference is that, compared with the results shown in Figure 2 and Figure 3, even if the SNR is high, the maximum value of the

P C C

mean is lower than 1. The reason is that when the isolation I of the transmitted is ignored, the limitation of Equation (34) is

\lim_{SNR_a v \to \infty} {\hat{S}}_{p q / H H_MS} = \lim_{SNR_a v \to \infty} \frac{G_{p q}}{G_{H H}} = \frac{A_{p q}}{A_{H H}}

(43)

If the system errors do not exist,

A_{p q}

is equal to

S_{p q}

and the limitation of Equation (34) is

S_{p q / H H}

. However, in this subsection, the system errors are set as Equation (42), leading to the limitation is not equal to

S_{p q / H H}

. Therefore, the

P C C

mean is lower than 1.

5. Conclusions

For moving target PSM estimation, the measurement selection (MS) method is proposed in this paper. This method selects the measurements by the criterion based on the SNR of the integration echoes. Several numerical simulations and comparative analysis are conducted to demonstrate and validate the superior performance of the proposed method, especially when the target exhibits non-uniform motion, compared with the traditional methods. Furthermore, the authors plan to conduct research on the PSM estimation for the extended target.

Author Contributions

F.W. and C.L. proposed the idea of this paper, X.W. and C.P. conceived and designed the experiments; F.W. performed the experiments; F.W. and Y.L. analyzed the data; F.W. wrote the paper.

Funding

This work was funded by the National Natural Science Foundation of China, grant numbers 61490690, 61490694 and 61701512.

Conflicts of Interest

The authors declare no conflict of interest.

References

Sandeep, G.; Arye, N. Polarimetric MIMO Radar With Distributed Antennas for Target Detection. IEEE Trans. Signal Process. 2010, 58, 1689–1697. [Google Scholar]
Martorella, M.; Giusti, E.; Demi, L.; Zhou, Z.; Cacciamano, A.; Berizzi, F.; Bates, B. Target recognition by means of polarimetric ISAR images. IEEE Trans. Aerosp. Electron. Syst. 2011, 47, 225–239. [Google Scholar] [CrossRef]
Sharma, J.J.; Hajnsek, I.; Papathanassiou, K.P.; Moreira, A. Polarimetric decomposition over glacier ice using long-wavelength airborne PolSAR. IEEE Trans. Geosci. Remote Sens. 2011, 49, 519–535. [Google Scholar] [CrossRef] [Green Version]
Li, C.; Yang, Y.; Li, Y.; Wang, X. Frequency Calibration for the Scattering Matrix of the Simultaneous Polarimetric Radar. IEEE Geosci. Remote Sens. Lett. 2018, 15, 429–433. [Google Scholar] [CrossRef]
Nashashibi, A.Y.; Sarabandi, K.; Frantzis, P.; de Roo, R.D.; Ulaby, F.T. An ultrafast wide-band millimeter-wave (MMW) polarimetric radar for remote sensing applications. IEEE Trans. Geosci. Remote Sens. 2002, 40, 1777–1786. [Google Scholar] [CrossRef]
Yanovsky, F.J.; Russchenberg, H.W.J.; Unal, C.M.H. Retrieval of information about turbulence in rain by using Doppler-polarimetric radar. IEEE Trans. Microw. Theory Tech. 2005, 53, 444–450. [Google Scholar] [CrossRef]
Giuli, D.; Fossi, M.; Facheris, L. Radar target scattering matrix measurement through orthogonal signals. IEE Proc. F Radar Signal Process. 1993, 140, 233–242. [Google Scholar] [CrossRef]
Wang, F.L.; Pang, C.; Li, Y.Z.; Wang, X.S. Algorithms for Designing Unimodular Sequences with High Doppler Tolerance for Simultaneous Fully Polarimetric Radar. Sensors 2018, 18, 905. [Google Scholar] [CrossRef] [PubMed]
Xu, J.; Yu, J.; Peng, Y.N.; Xia, X.G. Radon-Fourier Transform for Radar Target Detection, I: Generalized Doppler Filter Bank. IEEE Trans. Aerosp. Electron. Syst. 2011, 47, 1186–1202. [Google Scholar] [CrossRef]
Li, M.; He, Z.S.; Li, W.C. Transient Interference Mitigation via Supervised Matrix Completion. IEEE Geosci. Remote Sens. Lett. 2016, 13, 907–911. [Google Scholar] [CrossRef]
Perry, R.P.; DiPietro, R.C.; Fante, R.L. SAR imaging of moving targets. IEEE Trans. Aerosp. Electron. Syst. 1999, 35, 188–200. [Google Scholar] [CrossRef]
Santalla, V.; Antar, Y.M.M. A comparison between different polarimetric measurement schemes. IEEE Trans. Geosci. Remote Sens. 2002, 40, 1007–1017. [Google Scholar] [CrossRef]
Niemeijer, R.J. Doppler-polarimetric radar signal processing. Electr. Eng. Math. Comput. Sci. 1996. [Google Scholar]
Xiao, J.J.; Nehorai, A. Joint Transmitter and Receiver Polarization Optimization for Scattering Estimation in Clutter. IEEE Trans. Signal Process. 2009, 57, 4142–4147. [Google Scholar] [CrossRef]
Pang, C. Study on Theory and Technology for Accurate Polarimetry Using Phased Array Radar. Ph.D. Thesis, Information and Communication Engineering, College of Electronic Science, National University of Defense Technology, Changsha, China, 2016. [Google Scholar]
Qian, L.C.; Xu, J.; Xia, X.G.; Sun, W.F.; Long, T.; Peng, Y.N. Wideband-scaled Radon-Fourier transform for high-speed radar target detection. IET Radar Sonar Navig. 2014, 8, 501–512. [Google Scholar] [CrossRef]
Richards, M.A. Fundamentals of Radar Signal Processing, Second Edition; McGraw-Hill: New York, NY, USA, 2014; pp. 71–73. [Google Scholar]
John, R.R.; Frederick, J.M. Foundations of Electromagnetic Theory; Iaddison-Wesley: Reading, MA, USA, 1960. [Google Scholar]
Harold, M. Remote Sensing with Polarimetric Radar; John Wiely and Sons, Inc.: Hoboken, NJ, USA, 2007. [Google Scholar]
Li, Z. Design, simulation and validation of dual-channel polaremetric agile radar technology. Ph.D. Thesis, Faculty of Electrical Engineering, Mathematics and Computer Science, Delft University of Technology, Delft, The Netherlands, 2011. [Google Scholar]

Figure 1. The signal processing flow chart of the STSR radar.

Figure 2. Results of three PSM estimation methods where

P = 500

Figure 2. Results of three PSM estimation methods where

P = 500

Figure 3. Results of three PSM estimation methods where

P = 10

Figure 3. Results of three PSM estimation methods where

P = 10

Figure 4. Results of three PSM estimation methods for different integrated pulses. (a) Results of three PSM estimation methods where N = 4; (b) Results of three PSM estimation methods where N = 6; (c) Results of three PSM estimation methods where N = 8; (d) Results of three PSM estimation methods where N = 10.

Figure 5. Results of three PSM estimation methods where

P = 500

Figure 5. Results of three PSM estimation methods where

P = 500

Figure 6. Results of three PSM estimation methods where

P = 10

Figure 6. Results of three PSM estimation methods where

P = 10

Table 1. Parameters for numerical experiments.

Component	Value	Component	Value
Carrier frequency	10 GHz	Target PSM	$[\begin{matrix} 1 0.5 j \\ 0.5 j 0.9 \end{matrix}]$
Mode	STSR	Control factor	$P = 500, 10$
Bandwidth	1 GHz	Initial position	5000 m
Pulse duration	100 $μ s$	${SNR}_{a v}$	$- 10 \sim 30 dB$
PRF	1000 Hz	Noise	Gaussian
Measurements	12	Trials of each SNR	500

Table 2. Methods for PSM Estimation.

	PI	PC	MS
HH	1	1	1
VH	$\frac{\sum_{k = 1}^{N} G_{V H, k}}{\sum_{k = 1}^{N} G_{H H, k}}$	$\frac{\sum_{k = 1}^{N} G_{V H, k} \exp (- j \frac{4 π {\hat{v}}_{0} k T_{PRT}}{λ})}{\sum_{k = 1}^{N} G_{H H, k} \exp (- j \frac{4 π {\hat{v}}_{0} k T_{PRT}}{λ})}$	$\frac{G_{V H}}{G_{H H}}$
HV	$\frac{\sum_{k = 1}^{N} G_{H V, k}}{\sum_{k = 1}^{N} G_{H H, k}}$	$\frac{\sum_{k = 1}^{N} G_{H V, k} \exp (- j \frac{4 π {\hat{v}}_{0} k T_{PRT}}{λ})}{\sum_{k = 1}^{N} G_{H H, k} \exp (- j \frac{4 π {\hat{v}}_{0} k T_{PRT}}{λ})}$	$\frac{G_{H V}}{G_{H H}}$
VV	$\frac{\sum_{k = 1}^{N} G_{V V, k}}{\sum_{k = 1}^{N} G_{H H, k}}$	$\frac{\sum_{k = 1}^{N} G_{V V, k} \exp (- j \frac{4 π {\hat{v}}_{0} k T_{PRT}}{λ})}{\sum_{k = 1}^{N} G_{H H, k} \exp (- j \frac{4 π {\hat{v}}_{0} k T_{PRT}}{λ})}$	$\frac{G_{V V}}{G_{H H}}$

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

Share and Cite

MDPI and ACS Style

Wang, F.; Li, C.; Pang, C.; Li, Y.; Wang, X. A Method for Estimating the Polarimetric Scattering Matrix of Moving Target for Simultaneous Fully Polarimetric Radar. Sensors 2018, 18, 1418. https://doi.org/10.3390/s18051418

AMA Style

Wang F, Li C, Pang C, Li Y, Wang X. A Method for Estimating the Polarimetric Scattering Matrix of Moving Target for Simultaneous Fully Polarimetric Radar. Sensors. 2018; 18(5):1418. https://doi.org/10.3390/s18051418

Chicago/Turabian Style

Wang, Fulai, Chao Li, Chen Pang, Yongzhen Li, and Xuesong Wang. 2018. "A Method for Estimating the Polarimetric Scattering Matrix of Moving Target for Simultaneous Fully Polarimetric Radar" Sensors 18, no. 5: 1418. https://doi.org/10.3390/s18051418

APA Style

Wang, F., Li, C., Pang, C., Li, Y., & Wang, X. (2018). A Method for Estimating the Polarimetric Scattering Matrix of Moving Target for Simultaneous Fully Polarimetric Radar. Sensors, 18(5), 1418. https://doi.org/10.3390/s18051418

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

Article Menu

A Method for Estimating the Polarimetric Scattering Matrix of Moving Target for Simultaneous Fully Polarimetric Radar

Abstract

1. Introduction

2. Signal Model for Moving Target in STSR Radar

3. PSM Estimation for the Moving Target

4. Simulation Results and Discussion

4.1. PSM Estimation without System Errors

4.2. PSM Estimation with System Errors

5. Conclusions

Author Contributions

Funding

Conflicts of Interest

References

Share and Cite

Article Metrics

Article Access Statistics

Further Information

Guidelines

MDPI Initiatives

Follow MDPI