InTech uses cookies to offer you the best online experience. By continuing to use our site, you agree to our Privacy Policy.

Earth and Planetary Sciences » Geology and Geophysics » "Geoscience and Remote Sensing", book edited by Pei-Gee Peter Ho, ISBN 978-953-307-003-2, Published: October 1, 2009 under CC BY-NC-SA 3.0 license. © The Author(s).

Chapter 12

Moving Target Detection and Velocity Estimation in Multi-Channel AT-InSAR Systems from Amplitude and Phase Data

By Alessandra Budillon
DOI: 10.5772/8304

Article top

Overview

Along-Track Interferometry system single baseline geometry.
Figure 1. Along-Track Interferometry system single baseline geometry.
Inteferogram joint amplitude phase pdf in the hypothesis H0 for CNR=10 dB.
Figure 2. Inteferogram joint amplitude phase pdf in the hypothesis H0 for CNR=10 dB.
Inteferogram joint amplitude phase pdf in the hypothesis H1 for CNR= 10 dB, SCR= 10 dB, ur=ur,amb/2, corresponding to phase 1.5 rad.
Figure 3. Inteferogram joint amplitude phase pdf in the hypothesis H1 for CNR= 10 dB, SCR= 10 dB, ur=ur,amb/2, corresponding to phase 1.5 rad.
CRLB1/2 relative to the estimation of the radial velocity vr,Vs. the radial velocity and relevant to the model based on phase-only data (a) and to the model based on amplitude and phase data (b), for CNR=10 dB, and varying SCR (0,5,10,15,20 dB).
Figure 4. CRLB1/2 relative to the estimation of the radial velocity vr,Vs. the radial velocity and relevant to the model based on phase-only data (a) and to the model based on amplitude and phase data (b), for CNR=10 dB, and varying SCR (0,5,10,15,20 dB).
CRLB1/2 relative to the estimation of SCR,Vs. SCR and relevant to the model based on phase-only data (a) and to the model based on amplitude and phase data (b), for CNR=10 dB, and varying vr in the range (0, vr,amb).
Figure 5. CRLB1/2 relative to the estimation of SCR,Vs. SCR and relevant to the model based on phase-only data (a) and to the model based on amplitude and phase data (b), for CNR=10 dB, and varying vr in the range (0, vr,amb).
CRLB1/2 and RMSE relative to the estimation of the radial velocity vr,Vs. the radial velocity and relevant to the model based on phase-only data (a) and to the model based on amplitude and phase data (b), for CNR=10 dB, and SCR=10 dB.
Figure 6. CRLB1/2 and RMSE relative to the estimation of the radial velocity vr,Vs. the radial velocity and relevant to the model based on phase-only data (a) and to the model based on amplitude and phase data (b), for CNR=10 dB, and SCR=10 dB.
Mean value and RMSE relative to the estimation of the radial velocity vr,Vs. the radial velocity and relevant to the model based on phase-only data (a) and to the model based on amplitude and phase data (b), for CNR=10 dB, and SCR=10 dB.
Figure 7. Mean value and RMSE relative to the estimation of the radial velocity vr,Vs. the radial velocity and relevant to the model based on phase-only data (a) and to the model based on amplitude and phase data (b), for CNR=10 dB, and SCR=10 dB.
CRLB1/2 and RMSE relative to the estimation of SCR,Vs. the radial velocity and relevant to the model based on phase-only data (a) and to the model based on amplitude and phase data (b), for CNR=10 dB, and SCR=10 dB.
Figure 8. CRLB1/2 and RMSE relative to the estimation of SCR,Vs. the radial velocity and relevant to the model based on phase-only data (a) and to the model based on amplitude and phase data (b), for CNR=10 dB, and SCR=10 dB.
Inteferogram joint amplitude phase pdf in the hypothesis H0 and H1 for ur=ur,amb/2, CNR=10 dB, SCR=10 dB, and the separation curve.
Figure 9. Inteferogram joint amplitude phase pdf in the hypothesis H0 and H1 for ur=ur,amb/2, CNR=10 dB, SCR=10 dB, and the separation curve.
Inteferogram joint amplitude phase pdf in the hypothesis H0 and H1 for ur=ur,amb/2, CNR=10 dB, SCR=0 dB, and the separation curve.
Figure 10. Inteferogram joint amplitude phase pdf in the hypothesis H0 and H1 for ur=ur,amb/2, CNR=10 dB, SCR=0 dB, and the separation curve.
ROC in the amplitude phase approach (solid line) and in the phase-only approach (dashed line) for ur=ur,amb/2, CNR=10 dB, SCR=10dB (a) and SCR=0 dB (b).
Figure 11. ROC in the amplitude phase approach (solid line) and in the phase-only approach (dashed line) for ur=ur,amb/2, CNR=10 dB, SCR=10dB (a) and SCR=0 dB (b).
ROC in the amplitude phase approach for CNR=10 dB, SCR=10 dB, varying ur in the range (0, ur,amb) (a) and for ur=ur,amb/2, varying SCR (-5 dB, 0 dB, 5 dB, 10 dB) (b).
Figure 12. ROC in the amplitude phase approach for CNR=10 dB, SCR=10 dB, varying ur in the range (0, ur,amb) (a) and for ur=ur,amb/2, varying SCR (-5 dB, 0 dB, 5 dB, 10 dB) (b).
ROC in the multi-channel amplitude phase approach (solid line) for CNR=10 dB, SCR=10 dB, N=8 interferograms compared with the single channel amplitude phase approach (dashed line baseline b1 dotted line baseline b2).
Figure 13. ROC in the multi-channel amplitude phase approach (solid line) for CNR=10 dB, SCR=10 dB, N=8 interferograms compared with the single channel amplitude phase approach (dashed line baseline b1 dotted line baseline b2).

Moving Target Detection and Velocity Estimation in Multi-Channel AT-InSAR Systems from Amplitude and Phase Data

Alessandra Budillon1

1. Introduction

Recently, Along Track Interferometric Synthetic Aperture Radar systems (AT-InSAR) have been applied for traffic monitoring of ground vehicles (Meyer et al. 2006, Chapin & Chen, 2006, Hinz et al. 2007).

AT-InSAR systems are composed by more than one SAR antennas (typically two), mounted on the same platform and displaced along the platform moving direction. The separation distance between the antennas is denoted as baseline.

From the acquisitions of two or more image signals these systems are able to recover additional information about the observed scene: they allow the detection of moving targets on the ground and the estimation of their radial velocity (Raney, 1971). This is possible because the interferometric phase, i. e. the (-π,π] wrapped phase of the signal obtained from the point to point correlation between the complex images acquired from the two interferometric antennas, is related to the radial velocity through a known mapping. Then, after the so-called Phase Unwrapping (PhU) operation, a map of the target range velocity can be retrieved.

Detection and radial velocity estimation of a ground moving target are challenging problems, due to the difficulty of separating the moving target signal from the stationary background (clutter) (Chiu, 2003). Several methods, based on very different approaches, have been proposed in literature, such as Displaced Phase Centre Antennas (DPCA systems) techniques ( Gierull & Livingstone, 2004 , Chiu & Livingstone 2005), and Space-Time Adaptive Processing (STAP) (Ender, 1999, Klemm, 2002, Gierull & Livingstone, 2004 ). Interest in investigating AT-InSAR processor is motivated since such alternative techniques attempt to reject or cancel the stationary clutter but have the drawback that can attenuate slowly moving targets (Chiu & Livingstone 2005).

AT-InSAR systems usually use only interferometric phase information in order to estimate radial velocity (Chen, 2004, Budillon et al. 2005, Budillon et al. 2008a) while complex data in place of phase-only data have already been used in AT-InSAR systems detection applications ( Gierull, 2004 , Zhang et al. 2005, Budillon et al., 2008b) showing that detection performance improve when complex data are used.

In this paper it is proposed to consider both amplitude and phase of the interferometric SAR image since in the case of target velocity estimation of not extended targets, the exploitation of the image amplitude together with the image phase can add more information, even if the amplitude is influenced in a less sensitive way with respect to the phase. AT-InSAR approach can be considered clutter-limited since when the signal to clutter ratio decreases, the velocity estimation becomes gradually more and more critical, till it fails completely. Moreover there are solution ambiguities than can keep the velocity estimation from working correctly due to the wrapped phase measurements.

In this paper both above mentioned problems are solved by using statistical estimation methods, and exploiting multi-channel interferograms. The statistical estimation methods allow taking into account the correct statistics of the involved noise (likelihood model). The use of a multi-channel interferogram, that in this case can be obtained exploiting frequency diversity and/or baseline diversity, has a twofold effect: multi-channel interferograms can help to reduce the variance of the estimation, and, if properly chosen, can allow avoiding solution ambiguities (Budillon et al. 2005, Budillon et al. 2008c).

It is shown that combining the real and imaginary part of more than two acquired images (multi-channel approach) produce significative improvements in the velocity estimation accuracy and a sensitive reduction in the false alarm rate compared with AT-InSAR conventional systems using phase-only data.

In section 2 the AT-InSAR statistical model has been presented and the joint interferogram amplitude and phase distribution has been derived. Based on this distribution, in section 3 a radial velocity estimation maximum likelihood approach using more interferogram channels has been reported. Cramer Lower Bounds and Root Mean Squared Error show the method performance on simulated data using Terra SAR-X parameters and are evaluated in the case of phase-only data and amplitude and phase data. In section 4 a likelihood ratio test is adopted to detect the moving target and performance detection in terms of Probability of detection and false alarm have been examined comparing results obtained on phase-only data and on amplitude and phase data. Moreover a multi-channel detection strategy is proposed and compared with the one based on a single interferogram. Finally follow conclusions in section 5.

2. AT-InSAR statistical model

In this section is presented the statistical model of the AT-InSAR signal. Consider an AT-InSAR system constituted by two antennas moving along the direction x (azimuth) (see Figure 1), and suppose that the two antennas are separated by a baseline b along the azimuth direction x, such that b«H, where H is the platform quota. Assume a target on the ground moving with a constant velocity v T = v Tx x + v Tr r, where v Tx and v Tr are the velocity components along the azimuth and the line of sight direction (range) r, respectively. The azimuth velocity component v Tx produce a Doppler slope change causing a defocusing in the moving target image. The radial velocity component v Tr produce a Doppler history different from that of the stationary background, and an azimuth displacement of the target.

Suppose that |v Tx|, |v Tr|<<|v P|, where v P=v P x is the velocity of the flying platform and H»X and H»W, where X and W are the antenna footprint dimensions.

media/image1.png

Figure 1.

Along-Track Interferometry system single baseline geometry.

The SAR image signal formed by each antenna can be modeled as the superposition of the contributions of the moving target, of the stationary clutter, and of the additive noise. Then in a fixed image pixel we have:

Z1={Zc1+N1+ZT1 inpresenceofmovingtargetH1Zc1+N1  inabsenceofmovingtarget    H0     Z2={Zc2+N2+ZT2in presence of moving target   H1Zc2+N2in absence of moving target     H0
(1)

Z 1 and Z 2 are the computed image signals in the considered pixel, Z c1 and Z c2 are the clutter signals acquired by the two antennas, N 1 and N 2 represent the receiver thermal noise, and Z T1 and Z T2 denote the SAR images of the moving target produced by the two interferometric antennas, which will exhibit a phase factor related to the radial velocity:

ZT1=A1,           ZT2=A2ejφv 
(2)

where A1 and A2 are the target images obtained for zero velocity, and φv is given by (Raney, 1971):

φv=4πbλvr|vP|2π=4πbλur2π
(3)

where λ is the wavelength corresponding to the working frequency f=c/ λ of the SAR system, and < >represents the “modulo-2π” operation. In (3) the normalized radial velocity u r=v r/|v P| has been also introduced. Where the moving target is present, v r≠0 and consequently ϕ v≠0, otherwise the along-track inteferometric phase (3) is null. From (3) it is easy to derive that the ambiguity velocity value, such that the interferometric phase is equal to ± π, is given by u r,amb λ /(4b). For |u r|> λ /(4b) the interferometric phase wraps, as evidenced also by the “modulo-2 π” operation. Also disturbing effects have to be taken into account, they are related to different parameters such as the signal to clutter ratio (SCR), the clutter to thermal noise ratio (CNR), and the clutter coherence γ c. Since the time elapsing between the two interferometric acquisitions is very small (typically of the order of a millisecond) the clutter coherence can be considered equal to one. Then only the effect of SCR and CNR has to be considered.

To analyze the effect that the clutter and noise signals have on the velocity estimation accuracy, a statistical model for the involved signals has to be introduced. It is well known that the clutter signals Z c1 and Z c2 can be assumed random processes, whose real and imaginary parts are mutually uncorrelated Gaussian signals, with zero mean and same variance σc2 , since they are resulting from the superposition of the signals backscattered from many scattering centres lying in the resolution cell. N 1 and N 2 can be modelled as two additive (to the clutter) zero mean Gaussian complex processes independent of each other, independent on the clutter and with same variance 2 σN2

When the moving target is present, a deterministic model is applicable to the case of a target whose Radar Cross Section (RCS) can be expressed by a deterministic function of the incidence angle (Budillon et al., 2008a). This model applies to canonical scattering objects (such as corner reflectors, spheres, etc.), and to complex or extended targets whose RCS does not rapidly change between the interferometric acquisitions. An accurate knowledge of the average RCS values can be available only for accurately characterized targets (Palubinskas et al. 2004).

A Gaussian model for the target allows to take into account the lack of knowledge of the target RCS values (that can be described in terms of variance σT 2), and then of the SCR, and applies to complex or extended targets which can be considered to consist of a large number of isotropic scattering elements, randomly distributed in a region whose dimensions are large compared to the wavelength of the illuminating radiation, and all contributing to the overall signal with the same weight. In the following the target signals ZT1 and ZT2, have been modelled as zero mean (complex) Gaussian processes.

Then, when the moving target is absent (ZT1 = ZT2 = 0), the two processes Z1 and Z2 are Gaussian with zero-mean and correlation coefficient γH0 given by:

γH0=E[(Zc1+N1)(Zc2+N2)*]E[|Zc1+N1|2]E[|Zc2+N2|2]=γc(1+1CNR)
(4)

where E[] denotes the expectation operation, * denotes the conjugate, γc is the clutter coherence (real valued, in the ATI application can be considered equal to one), representing the correlation between images Zc1 and Zc2, and CNR = σc2 / σN2 , where 2 σc2 and 2 σN2 are the clutter and thermal noise powers respectively (the factor two is due to the sum of the powers of the real and imaginary parts).

Instead, when we are in presence of the moving target (ZT1 0, ZT2 0), the expression of correlation coefficient change with respect to (4) and γH1 is given by:

γH1=E[(Zc1+N1+ZT1)(Zc2+N2+ZT2)*]E[|Zc1+N1+ZT1|2] E[|Zc2+N2+ZT2|2]=     =γcσc2+σT2γTσc2+σn2+σT2=γc+γTSCR1+1CNR+SCR
(5)

where SCR = σT2 / σc2 , where 2 σT2 is the target power, and γT is the target (complex) coherence that depends on the target velocity through the nominal phase (3):

γT=E[ZT1ZT2*]E[|ZT1|2]E[|ZT2|2]=E[A1A2*]E[|A1|2]E[|A2|2]ejφv=γAejφv
(6)

where γA is the target coherence for zero radial velocity that is usually assumed equal to one.

The two processes Z1=Z1r+jZ1i and Z2=Z2r+jZ2i are Gaussian, then the joint probability density function of Z=[Z1r Z2r Z1i Z2i]T, is Gaussian with zero mean and covariance matrix

C={Cc+CN                    H0Cc+CN+CT            H1
(7)

where the matrices Cc, CN and CT are respectively the clutter, noise and target covariance matrix. It can be easily shown (Davenport & Root, 1958) that:

C=σ2[1Re(γ)0Im(γ)Re(γ)1Im(γ)00Im(γ)1Re(γ)Im(γ)0Re(γ)1]
(8)

where, in the hypothesis H0, has to be taken σ2=σc2+σN2 and γ=γH0, in the hypothesis H1, σ2=σc2+σN2+σT2 and γ=γH1.

Then the joint probability density function of Z=[Z1r Z2r Z1i Z2i]T, is Gaussian, i.e.

fZ(z)=fZ1r,Z2r,Z1i,Z2i(z1r,z2r,z1i,z2i)=1(2π)2|C|1/2exp{12zTC1z}
(9)

The SAR interferometric amplitude and phase distribution can be derived from (9) introducing the interferometric signal I:

Ι=Z1Z2*=Wexp(jΦ)
(10)

The joint W and Φ pdfs, in the hypothesis H1 and H0, derived from (9) via variables transformations (Davenport & Root, 1958), are respectively:

fWΦ(w,φ;H1)=w2πσc2(1+1CNR+SCR)2(1|γH1|2)  K0(wσc2(1+1CNR+SCR)(1|γH1|2))exp{|γH1|wcos(φφo)σc2(1+1CNR+SCR)(1|γH1|2)}fWΦ(wφH0)=w2πσc2(1+1CNR)2(1γH02)  K0(wσc2(1+1CNR)(1γH02))exp{γH0wcos(φ)σc2(1+1CNR)(1γH02)}               
(11)

where K0 denote the modified Bessel function of order zero, w0,0φ2π and ϕ0 is γH1 phase. The analytical expressions of the pdfs ( (9) and (11)) and of the coherences ( (4) and (5)), reveal their dependence on CNR, SCR, radial velocity (trough the phase ϕ v), clutter coherence γc and target coherence γA. As far as the coherence values are concerned, they are assumed equal to 1 in ATI applications.

In Figure 1 and Figure 2 the joint pdf of W and Φ, respectively in the hypothesis H0 and H1, with CNR= 10 dB, SCR= 10 dB), ur=ur,amb/2 corresponding to phase 1.5 rad, are shown.

media/image30.png

Figure 2.

Inteferogram joint amplitude phase pdf in the hypothesis H0 for CNR=10 dB.

media/image31.png

Figure 3.

Inteferogram joint amplitude phase pdf in the hypothesis H1 for CNR= 10 dB, SCR= 10 dB, ur=ur,amb/2, corresponding to phase 1.5 rad.

3. Multi-Channel AT-InSAR moving target velocity estimation

3.1. Joint estimation of velocity and SCR via Maximum likelihood approach

Since interferometric phase is measured in the interval (-π,π], then a Phase Unwrapping (PhU) operation is required to retrieve the target radial velocity. The PhU operation presents solution ambiguities when only one phase interferogram (single-channel) is used. It has already been shown in (Budillon et al. 2005, Budillon et al. 2008c) that the joint use of multi-channel configurations (deriving from the use of more than two interferometric images acquired with different baselines or at different working frequencies) and of classical statistical estimation techniques allows to obtain very accurate solutions and to overcome the limitations due to the presence of ambiguous solutions, intrinsic in the single-channel configurations.

Different baseline data sets (at least two) can be generated when the AT-InSAR system is constituted by more than two antennas (at least three). Different frequency data sets can be generated in two ways. In the first, we can suppose that the SAR sensors can operate at different working frequencies, for instance in X and C bands simultaneously. In the second, the multi-frequency interferograms can be obtained by sub-band filtering the interferometric images splitting the overall bandwidth (Budillon et al. 2008c). Note that while the use of a different working frequency or baseline does not affect the SCR values, the generation of adding frequencies by partitioning the available band reduces the SCR value. This SCR reduction is in inverse relation to the number of looks, as the spatial resolution worsens increasing the number of looks.

Likelihood function is easily derived from either pdf (9) or (11) in the H1 hypothesis. As discussed in the previous section, it depends on CNR, SCR and radial velocity since clutter and target coherence are assumed equal to 1 in ATI applications. CNR value can be easily computed from the data isolating an area where the target is absent. Then, the final estimation can be casted as a joint maximum likelihood estimation (Kay, 1993) of velocity and SCR:

[u^r,SC^R]=argmaxur,SCRL(ur,SCR)L(ur,SCR)=k=1,,NfΖ(Z(k)|ur,SCR,H1)single channel likelihood function
(12)

where the samples {Z(k)}k=1,N  represent the SAR signals acquired in N channels. The factorization in (12) comes from the assumed statistical independence of the multi-channel interferograms.

3.2. Performance assessment

Estimation performance evaluation has been carried out using Terra SAR-X parameters in Table 1, but in order to consider a multi-channel system a second baseline b2=1.8b1 [m] (b1=1.2 [m]), has been adjoined. By sub-band filtering the interferometric images also 4 azimuth looks have been considered, obtaining in total N=8 channels. The maximum normalized radial velocity value that can be unambiguously detected results |ur,max|=λ/(4b)=6.5×10-3, corresponding to a not normalized velocity |vr,amb| of about 49.4 m/sec (178 Km/h).

TerraSAR-X
Quota514.8 Km
Platform velocity7.6 Km/s
Along track antenna dimension4.8 m
Across track antenna dimension0.8 m
Along track baseline1.2 m
Working frequencyX band? f X =9.65 GHz
Wavelength3.12 cm
Range bandwidth150 MHz

Table 1.

Main parameters of Terra SAR X system.

To evaluate the performance of the ML estimator (12), the Cramer Rao Lower Bounds (CRLBs) (Kay, 1993) and the Root Mean Square Errors (RMSEs) for the unknown parameters ( vr,SCR ) have been computed. CRLBs depend on the data model and represent the maximum accuracy attainable with given data. In order to point out the advantages of taking into account amplitude and phase information they have been compared with the ones obtained using a phase-only approach, i.e. a maximum likelihood estimation based on the phase-only distribution (Budillon et al. 2008a).

The CRLB1/2 relative to the not normalized radial velocity vr and SCR are reported respectively in Figures 4 and 5. In Figure 4(a) it is shown the CRLB1/2 relative to the estimation of the not normalized radial velocity vr, Vs. the radial velocity and relevant to the model based on phase-only data. It is evident the significative improvements in the maximum accuracy attainable using the amplitude and phase model reported in Figure 4(b). The CRLB1/2 have been evaluated numerically and for different values of vr in the range (0, vr,amb), for CNR=10 dB, and varying SCR (0,5,10,15,20 dB). It can be appreciated that as expected the CRLB1/2 are lower for higher SCRs.

media/image35.png

Figure 4.

CRLB1/2 relative to the estimation of the radial velocity vr,Vs. the radial velocity and relevant to the model based on phase-only data (a) and to the model based on amplitude and phase data (b), for CNR=10 dB, and varying SCR (0,5,10,15,20 dB).

In Figure 5(a) the CRLB1/2 relative to the estimation of the SCR, Vs. the SCR and relevant to the model based on phase-only data is shown. Also for this parameter it is evident the significative improvements in the maximum accuracy attainable using the amplitude and phase model reported in Figure 5(b). Moreover it can be seen that using the amplitude and phase model the accuracy is slightly dependent on the SCR values, in both case the CRLB1/2 is smaller for higher values of radial velocities, i.e. as expected it is easier to estimate the SCR of a faster target. In Figure 6(a) are shown the CRLB1/2 and the RMSE relative to the estimation of the not normalized radial velocity vr, Vs. the radial velocity and relevant to the model based on phase-only data, for CNR=10 dB, and SCR=10 dB. They can be compared with the correspondent CRLB1/2 and the RMSE relevant to the model based on amplitude and phase data in Figure 6(b). It is noticeable that the performances are improved and also the RMSE is closer to the CRLB1/2 when the estimation is based on the amplitude and phase model. In Figures. 7(a) and 7(b) the statistical mean values and the RMSEs for different estimated velocities in the range (0, vr,amb), for CNR=10 dB and SCR=10 dB, respectively for the model based on phase-only data and to the model based on amplitude and phase data, are reported. Finally in Figures. 8(a) and 8(b) CRLB1/2 and the RMSE relative to the estimation of SCR,Vs. the radial velocity are shown respectively for the two models. It is evident again the improvement attainable in case it is assumed the amplitude and phase model.

media/image36.png

Figure 5.

CRLB1/2 relative to the estimation of SCR,Vs. SCR and relevant to the model based on phase-only data (a) and to the model based on amplitude and phase data (b), for CNR=10 dB, and varying vr in the range (0, vr,amb).

media/image37.png

Figure 6.

CRLB1/2 and RMSE relative to the estimation of the radial velocity vr,Vs. the radial velocity and relevant to the model based on phase-only data (a) and to the model based on amplitude and phase data (b), for CNR=10 dB, and SCR=10 dB.

media/image38.png

Figure 7.

Mean value and RMSE relative to the estimation of the radial velocity vr,Vs. the radial velocity and relevant to the model based on phase-only data (a) and to the model based on amplitude and phase data (b), for CNR=10 dB, and SCR=10 dB.

media/image39.png

Figure 8.

CRLB1/2 and RMSE relative to the estimation of SCR,Vs. the radial velocity and relevant to the model based on phase-only data (a) and to the model based on amplitude and phase data (b), for CNR=10 dB, and SCR=10 dB.

4. Multi-Channel AT-InSAR moving target detection

4.1. Likelihood ratio test

A moving target can be detected in the conventional way by comparing the interferometric phase ϕ with a threshold ηT in the interval (-π,π].

The performance of the detection process can be evaluated using the interferometric phase statistics (Budillon et al. 2088a). They are, as expected, better for high values of SCR, i.e. when the moving targets power is significantly larger than the clutter power. For moving targets mingling with the background clutter, the detection capability worsen, so that if one wants low values of PFA, the PD can decrease to very low values, not consistent with the applications.

As in the case of the velocity estimation (see Section 3) both amplitude and phase of the interferogram are considered instead of taking into account only the interferometric phase. Based on the pdfs (11) a constant false alarm rate (CFAR) detector can be designed.

In order to detect a moving target a likelihood ratio test is proposed, likelihood is derived from (11):

Λ(z)=fWΦ(w,φ;u^r,SC ^R,H1)fWΦ(w,φ;H0)H0H1η
(13)

Probability of false alarm PFA and Probability of detection PD are derived from (13)

PFA=Pr{Λη;H0}PD=Pr{Λη;H1}
(14)

The threshold η depends on a fixed PFA.

4.2. Performance assessment

The detection performance evaluation has been carried out using the same multi-channel system presented in section 3.2.

The proposed approach provides curves of separation between the two classes (see Figure 9-10). In Figure 9 the separation curve for the two hypothesis, presence and absence of a moving target, has been evaluated for CNR=10 dB, SCR=10 dB, ur=ur,amb/2=3.25×10-3 (corresponding to the nominal noise-free value ϕv=π/2), a threshold has been chosen such that PD=0.91 and PFA=0.001. Figure 10 shows the separation curve for the two hypothesis for CNR=10 dB, SCR= 0 dB, ur=ur,amb/2=3.25×10-3, a threshold has been chosen such that PD=0.7 and PFA=0.05.

For comparison with the conventional interferometric approach, the two Receiver Operating Characteristics (ROC) have been derived in both case SCR= 10 dB and SCR= 0 dB, for the amplitude and phase approach (solid line) and for the phase-only case (dashed line) (see Figure 11). It is clear the advantage in considering both amplitude and phase, for a fixed PFA a higher PD can be obtained.

In Figure 12 it can be appreciated the ROC dependence, in the amplitude and phase approach, on the radial velocity (Figure 12(a)), for CNR=10 dB, SCR=10 dB and on SCR (Figure 12 (b)), for CNR=10 dB and ur=ur,amb/2. As expected it is easier to detect a faster and stronger (in terms of reflectivity) target.

In order to exploit the multi-channel interferograms in the detection process, suppose that the detection probability of one of the channel corresponding to the first baseline is equal to PD1, and that the detection probability of one of the channel corresponding to the second baseline is equal to PD2.

The probability that the target is detected from (N/2+j) channels (j=1,…,N/2) on a total of N channels results:

Pj=k=jN/2[(N/2k)PD1k(1PD1)N/2k][(N/2N/2k+j)PD2N/2k+j(1PD2)kj]
(15)

The proposed multi-channel detection strategy consists in considering the moving target present when the majority of the interferogram values are above prefixed thresholds, so that the detection probability results:

PDN/2=j=1N/2Pj
(16)

For the estimation of the false alarm probability can be used the same reasoning.

In Figure 13 the ROC in the multi-channel amplitude phase approach (solid line) for CNR=10 dB, SCR=10 dB, N=8 interferograms compared with the single channel amplitude phase approach is reported (dashed lines). It is evident the advantages in considering a multi-channel approach that allows to keep low PFA and at the same time high PD.

media/image44.png

Figure 9.

Inteferogram joint amplitude phase pdf in the hypothesis H0 and H1 for ur=ur,amb/2, CNR=10 dB, SCR=10 dB, and the separation curve.

media/image45.png

Figure 10.

Inteferogram joint amplitude phase pdf in the hypothesis H0 and H1 for ur=ur,amb/2, CNR=10 dB, SCR=0 dB, and the separation curve.

media/image46.png

Figure 11.

ROC in the amplitude phase approach (solid line) and in the phase-only approach (dashed line) for ur=ur,amb/2, CNR=10 dB, SCR=10dB (a) and SCR=0 dB (b).

media/image47.png

Figure 12.

ROC in the amplitude phase approach for CNR=10 dB, SCR=10 dB, varying ur in the range (0, ur,amb) (a) and for ur=ur,amb/2, varying SCR (-5 dB, 0 dB, 5 dB, 10 dB) (b).

media/image48.png

Figure 13.

ROC in the multi-channel amplitude phase approach (solid line) for CNR=10 dB, SCR=10 dB, N=8 interferograms compared with the single channel amplitude phase approach (dashed line baseline b1 dotted line baseline b2).

5. Conclusion

In this paper it has been presented the performance evaluation of multi-channel AT-InSAR systems, exploiting both amplitude and phase interferogram information, in terms of target radial velocity estimation accuracy and moving target detection ability.

A Gaussian target response model has been considered and the amplitude and phase joint pdf of the inteferogram has been derived. Based on this model a maximum likelihood approach has been used to estimate the target radial velocity. A comparison of the proposed approach with a phase-only system has been reported considering different system parameters, such as radial velocity, SCR, CNR. It reveals the benefits in exploiting both amplitude and phase interferogram information in terms of CRLB and RMSE. The analysis has been performed considering a multi-channel system based on the Terra SAR-X parameters but with a second baseline.

A constant false alarm rate (CFAR) detector has been considered. This approach provides curves of separation between the two classes, hypothesis H0, and H1, in order to detect the pixels in the SAR images where a moving target is present. The proposed approach outperforms the conventional phase-only approach and in particular is able to detect target with low SCR. ROCs have been presented varying SCR and the normalized radial velocity.

Regarding the detection process, the use of multi-channel interferograms, after the application of a threshold stage to each channel, allows to adopt a binary integration to combine single-channel decisions. Such a strategy, compared with the one based on a single interferogram, provides better results in terms of simultaneous low values of PFA and high values of PD.

References

1 - A. Budillon, G. Ferraiuolo, V. Pascazio, G. Schirinzi, 2005 “Multi-Channel SAR Interferometry via Classical and Bayesian Estimation Techniques”, J. of Applied Signal Processing, 20 3180 3193 .
2 - A. Budillon, V. Pascazio, G. Schirinzi, 2008a “Estimation of Radial Velocity of Moving Targets by Along-Track Interferometric SAR Systems”, IEEE Geosci. Remote Sensing Letters, 5 349 353 .
3 - A. Budillon, V. Pascazio, G. Schirinzi, 2008b “Moving Target Detection in Along Track SAR Interferometry from In-Phase and Quadrature Components Data”, in Proc. of IEEE International Geoscience and Remote Sensing Symposium (IGARSS’08), III-1178 III-1181.
4 - A. Budillon, V. Pascazio, G. Schirinzi, 2008c “Multichannel Along-Track Interferometric SAR Systems: Moving targets Detection and Velocity Estimation“, International Journal of Navigation and Observation, 2008 16
5 - E. Chapin, C. W. Chen, 2006 “GMTI Along-Track Interferometry Experiment”, IEEE Aerosp. Electronic Syst. Mag., 21 15 20 .
6 - C. W. Chen, 2004 “Performance Assessment of Along-Track Interferometry for Detecting Ground Moving Targets”, in Proc. of 2004 IEEE Radar Conference, 99 104 .
7 - S. Chiu, 2003 “Clutter effects on ground moving target velocity estimation with SAR along-track interferometry”, in Proc. of 2003 IEEE International Geoscience and Remote Sensing Symposium (IGARSS’03), Toulouse (France), 2 1314 1319 .
8 - S. Chiu, C. Livingstone, 2005 “A comparison of displaced phase centre antenna and along-track interferometry techniques for RADARSAT-2 ground moving target indication”, Can. J. Remote Sensing, 31 1 37 51 .
9 - W. B. Davenport, W. L. Root, 1958 An Introduction to the Theory of Random Signals and Noise, McGraw-Hill, Kogakusha,
10 - J. H. G. Ender, 1999 “Space-Time Processing for Multichannel Synthetic Aperture Radar”, IEE Electronics & Communication Engineering Journal, 11 29 38 .
11 - C. H. Gierull, 2004 “Statistical Analysis of Multilook SAR Interferograms for CFAR Detection of Ground Moving Targets,” IEEE Trans. Geosci. Remote Sensing, 42 691 701 .
12 - C. H. Gierull, C. Livingstone, 2004 “SAR-GMTI concept for RADARSAT-2” in Applications of Space-Time Adaptive Processing, R. Klemm (Ed.), IEE Publishers, London, UK.
13 - S. Hinz, F. Meyer, M. Eineider, R. Bamler, 2007 “Traffic Monitoring with Spaceborne SAR-Theory, Simulations, and Experiments”, Computer Vision and Image Understanding, 106 231 244 .
14 - S. M. Kay, 1993 Fundamentals of Statistical Signal Processing: Vol. I, Estimation Theory, Prentice-Hall.
15 - R. K. Klemm, 2002 Principles of space-time adaptive processing, 2nd edn, London: IEE.
16 - F. Meyer, S. Hinz, A. Laika, D. Weihing, R. Bamler, 2006 "Performance Analysis of the TerraSAR-X Traffic Monitoring Concept", ISPRS J. of Photogr. and Remote Sensing, 61 3/4), 225 242 .
17 - G. Palubinskas, H. Runge, P. Reinartz, 2004 “Radar signatures of road vehicles”, In Proc. of 2004 IEEE International Geoscience and Remote Sensing Symposium (IGARSS’04), Anchorage (USA), 2 1498 1501 .
18 - R. Raney, 1971 “Synthetic Aperture Imaging Radar and Moving Targets”, IEEE Trans Aerosp. Electron. Syst, 7 499 505 .
19 - Y. Zhang, A. Hajjari, K. Kim, B. Himed, 2005 “A Dual-Threshold ATI-SAR Approach for Detecting Slow Moving Targets”, in Proc. of 2005 IEEE Radar Conference, 295 299 .