Maximum Likelihood Multi-Baseline SAR Interferometry
|
|
- Roy Elliott
- 6 years ago
- Views:
Transcription
1 Maximum Likelihood Multi-Baseline SAR Interferometry G. Fornaro(*), A. Monti Guarnieri(**), A. Pauciullo(*), F. Rocca(**) (*) Istituto per il Rilevamento Elettromagnetico dell Ambiente (IREA), Consiglio Nazionale delle Ricerche (CNR), Via Diocleziano 38, 824 Napoli Italy, (**) Dipartimento di Elettronica e Informazione Politecnico di Milano Piazza Leonardo da Vinci 32, 233 Milano Italy, monti@elet.polimi.it. May 7, 25 Abstract We propose a technique to provide interferometry by combining multiple images of the same area. This technique exploits all the images jointly and performs an optimal spectral shift pre-processing to remove most of the decorrelation for distributed targets. It s applications are mainly for DEM generation at centimetric accuracy, and for differential interferometry. The major requirement is that targets are coherent over all the images: this may be the case of current multi-pass over desert areas, or better the case of images coming from future short revisit time systems (constellations, cart-wheel, geosynchronous SAR etc.). Introduction Interferometric SAR (InSAR) surveys have been exploited for terrain mapping and DEM generation since more than two decades. The enhanced performances that would have been achieved by combining multiple interferometric acquisitions have been studied in the early papers in SAR interferometry [], even though a large amount of literature has been built up since the early sixties for holographic applications. In papers [2], [3] a combination was proposed of many interferograms by exploiting techniques similar to the Chinese remainder theorem to reduce phase unwrapping problems. A Maximum Likelihood (ML) approach was proposed in [4] for solving the same problem, applied to Tandem ERS acquisition. In that case, the problem was made simple as only pairs of interferogram, assumed each one independent upon the other, were combined. The same technique was extended to the case of 3 images by exploiting jointly all the information in [5], or in the case of multi-frequency acquisitions in [6]. In all these papers, the model assumed for SAR reflectivity resembles a collection of point-scatterers. Then, in that case, a complete and assessed methodology has recently developed under the name of Permanent Scatterer [7]. On the other side, when the SAR reflectivity is modeled by a distributed target (speckle), these approaches becomes far from optimal, as they do not properly handle the frequency domain incorrelation of the sources. Indeed, it was shown in [8] that, a sensible increment in quality could be obtained in 2-pass interferometry if a proper prefiltering, tuned to the baseline and the terrain slope is exploited. In [9] and [] this spectral shift (SS) filtering was extended to the case of non-constant slopes, and in [] this approach was shown to be optimal in the MMSE sense, and extended to the combination of two images coming very different SAR modes (ScanSAR, SPOT etc.). The purpose of this paper is to extend such approach to a multi-channel environment, where multiple acquisitions are made available either with different frequencies or with different baselines. The purpose is to establish This work has been committed under ASI contract I/R/69/2 Multiple pass SAR Interferometry for very accurate Earth modeling 23-24
2 a theoretical framework based on a Maximum Likelihood estimate, in which the best estimate of the underlying geometry, linked to the terrain topography, is provided. Although we focus mainly on DEM generation, we will show that the approach is equally suited to retrieve motion fringes in multi-temporal acquisition, thus for mapping landslides, earthquakes etc. The major scope of the paper is to jointly exploit the information provided by distributed targets. Hence, the technique would give the best result with a set of coherent images. This would be the case of long temporal baselines in long term correlated areas (they should be non-vegetated in C or X band surveys), or using short temporal baseline. For such reasons the main applications will be framed in the next L band missions (ALOS, TERRASAR-L), or in the forthcoming constellations (Cartwheel etc.). 2 Multi-Channel model Let us refer to the multibaseline geometry in Fig., an extension of the conventional SAR Interferometry [2, 3] to the multi channel case. In the figure, the target P corresponds to a distributed scatterer, that is imaged using three different positions. Let us define the normal baseline B n as the distance between the master and the slave s track, measured normal to the slant range, azimuth plane. In the case shown in the picture, we could compute up to three interferograms, as the Hermitian products between two images, all of them coregistered in the same master reference. The contribute of the target P to the interferogram s phase would be proportional to the travel path difference between the two acquisitions ϕ = 4π λ (R s(p ) R M (P )) () where R S and R M are respectively the slave-target and the master-target travel path. Eventually we approximate the interferogram phase as a slow varying linear term[4]: θ φ(t) = f tan(θ α(t)) t (2) B n f R tan(θ α(t)) t (3) α(t) being the local slope (range-varying), R the closest approach, f the carrier frequency, and θ the look angle difference. The approximation in (3) holds for range intervals narrow enough to assume θ constant. In (3) the interferogram phase scales linearly both with the B n and with f, therefore we may get different interferograms either by exploiting multi-baseline or multi-frequency systems. In this paper we will refer to both cases as multi-channel, however we will address mostly the multi baseline case, that is much more common in spaceborne SARs. Let us assume to have N channels, including one master and (N ) slave images, and M samples (range bins) from each acquisition. Following reference [], we express each single acquisition as a filtered version of the wide-band reflectivity: y i (P (t)) = γ(t) exp (jφ(t)) + w i (t) (4) We limit our study to the D case, as the slant range, or fast-time, t, is the direction mostly affected by phase aliasing and baseline decorrelation due to the geometric deformations in the ground-to-range mapping [5]. In (4), γ is the source reflectivity, y i the focused signal of the i-th acquisition, and w i a noise source that cumulates all the decorrelation terms [6]. In order to proper process the data, we need to up-sample the received signal, so that the final range sampling would accommodate for all the possible spectral shifts, e.g. by computing (3) for the largest baseline. After up sampling each channel in 4, we end up in a discrete-time model that can be expressed by the following vector formulation [], represented schematically in Fig. 2: y i = F i Φ i γ + w i = H i γ + w i (5) where we assume bold notation for matrixes and vectors. In particular, we assume that the impulse response on each channel is of length L, and we to have N channels, and M samples out of each channel. We require then M+L- samples of the source. The vectors and matrixes involved in (5) are the following: 2
3 S B n B n2 Master S 2 P Figure : Multi-baseline Interferometric SAR: a master and two slave sensors (may be the same sensor following the same orbit) are shown. 3
4 y i is the column vector [M,] that corresponds to the complex SAR image, the data, coregistered in the reference of the Master: y i = [ y y M ] T where the suffix T stands for Hermitian matrix transposition, and we assume row and column indexes starting from. γ is the column vector [D,] that represents the source reflectivity. This is white over a large band. Φ i is a diagonal modulation matrix, [D,D] that express the topographic-dependent (or motion dependent) contributions: φ i ()... Φ i = φ i () φ i (D ) its element on the diagonal being just the sampled version of the topographic-dependent phase, φ(t) in 3. Notice that for different channel the phases simply scales according to the baseline. Henceforth, we will gather all the elements on the diagonal of Φ i on a vector of unknown parameters (to be estimated): { } φ = f R tan(θ α(t i )) t i (6) t i being the sampling times. φ is computed for a normalized Baseline (B n = ). The model further simplifies when we assume first that the terrain has constant slope. In that case there is only parameter to retrieve: φ(α). F i is the filter matrix [M,D] where D=M+L-, that is Toeplitz and contains the impulse response of the equivalent SAR end-to-end channel (that can be well approximated by an ideal bandpass), f L f L... f F i = f L f f L f H i is a matrix [M,D] H i = F i Φ i (7) that represents the channel to be estimated. Notice that the channel is now linear, but space-variant due to the modulation matrix Φ i, hence, H i is not block Toeplitz as usually assumed in literature; w i is the noise in the specific acquisition, it has the same size as y i. The model (5) frames in the Single-Input-Multiple-Output (SIMO) channel blind estimation, where the unknown channels, H i are to be retrieved from the outputs, y i, given only a few information on the inputs. Such model is widely studied in literature, as it is shared by many fields: estimates of Direction Of Arrivals (DOA), wireless cellular networks, tomography, etc. The reader could refer to paper [7] for a summary on blind deconvolution techniques for SIMO problems. Following the trend in literature, we provide a compact matrix formulation for the model (5), where the N input vectors, the N outputs and the N channel matrixes are stacked one upon the other: where: y = Hγ + w (8) y = [ y T... yi T... yn T w = [ w T... wi T... wn T ] T ] T 4
5 y and w are then column vector of size [NM,]. The channel matrix H in (8) is also a block matrix of size [NM,M+L-]: h H =... h i h N = φ ()f L φ (D )f M+L... φ ()f... h i h N w Φ F. γ Φ i F i. Φ N- F N- x w i x i w N- x N- y y i y N- Figure 2: Multi-Baseline Interferometric SAR system, modeled as a SIMO system. 3 ML estimate The straight and simplest derivation of the Maximum Likelihood estimate of the topographic dependent contributes, the vector φ, starts from the observation that the data, y, belong to a mutlivariate Normal process, therefore the Log-Likelihood maximization result in the following expression: ( )) φ ML = arg max log f (y; φ i bφ i = arg min ( y C y y ) (9) The inverse covariance matrix in (9) accounts for both the superimposed noise, w, and the so-called baseline or geometric decorrelation, that may arise when different spectral components of the same source, γ are imaged by the different channels. The covariance matrix C y, in (9) can be expressed as the following block notation: E [yy ] E [y y]... E [ ] y N yn E [yy ] E [yy ] E [yi y i]... C y = E [ yn y ] E [ yn y ] N () 5
6 Each block can be estimated basing on the model (8) and (7) (see Fig. 2): E [y i y i ] = F i φ i C γ φ i F i + C wi () E [y i y j ] = i j F i φ i C γ φ j F j + E [w i w j ] (2) = F i φ i C γ φ j F j (3) where C wi is the covariance matrix of the noise term in the i-th channel, and where we have introduced the definition: C γ = E [γγ ] The need of the covariance matrix inversion makes the ML estimator unfeasible, but for simple cases requiring one or two parameters to be retrieved and involving few images. The most interesting case could be the estimate of slopes basing on a very local scale (a few pixels) in a conventional two channel interferometric system [8]. 4 DML estimate The Deterministic Maximum Likelihood (DML) method approaches the same problem by including in the vector of parameters the samples of the source, γ [7], [9]. This choice provides a simpler solution, at the price of increasing the number of parameters to be estimated and potentially worsening the accuracy of the estimates. The problem is linear in both the channel and the sources, therefore we can cascade the optimizations with respect to γ and φ: ( ( { φdml, γ} = arg min min (y Hγ) C (y Hγ) ) ) (4) φ b bγ where C is the noise covariance matrix. Eventually we define the SIMO channel covariance: We first minimize with respect to the sources: R H = H C H = Φ F C FΦ (5) that leads to the solution: min bγ ( y C γh + H γ C y H γ C γh ) (6) γ DML = R H H C y = R H Φ F C y (7) where R H is the pseudo-inverse of the channel covariance. We do not need the explicit computation of γ DML if we are not interested, we just need to substitute its expression back into (4), to retrieve the estimate of φ: where P H is a projector in the norm of the channel H : φ DML = max (y P H y) (8) P H = C HR H H C (9) = C FΦR H Φ F C (2). The DML estimate φ DML appears as complex as the ML approach in (8), as one matrix inversion is required to compute R H at each iteration. However, a sub-optimal implementation is known as the Two-Step ML (TSML), introduced by Hua in [9]. This approach makes a first iteration by assuming that the channel covariance matrix is diagonal: R H = I. In that case the DML estimate would be ( φ DMLI = max y C FΦΦ F C y ) (2) bφ 6
7 this leads to an estimate of φ DMLI that is shown to be optimal as for SNR [9]. An approximate value of R H is computed by combining (5) with the first-step estimate φ DMLI in (2). Thereafter, the final estimate φ DML is retrieved from (2). In a similar way, in papers [] and [] a preliminary information of the geometry and the topography was used to derive the proper whitening operator. 5 Single baseline interferometry Let us compare the two estimates of the topographic-dependent phase in the simplest case, the single baseline interferometry. The model of the Interferometric SAR channel is represented in Fig. 3.a, whereas Fig. 3.b draws the linear estimate of the parameters φ. The linear estimate is justified by the Normal statistics of both sources and noise. We will however derive it, and show that, by a proper selection of the coefficients, the same scheme in Fig. 3.b can represent either the ML or the DML approach. In the figure, we have conventionally attributed half of the interferometric phase, Φ /2, to the master, and the other half to the slave. Furthermore, we changed the conventional SIMO model (see Fig. 2) by moving the noise contribution in front of the blocks F (that represents the SAR acquisition and focusing overall transfer function). This choice is closer to the actual SAR acquisition. However, the two models are both correct, provided that the proper noise covariance matrix is used. w Master Slave synthesis γ Φ F y H w 2 Slave G 2 y 2 y 2 Φ 2 F 2 y 2 G 2 y 2 y 2 H 2 (a) Δφ Master synthesis (b) Figure 3: Single baseline interferometry: (a) the acquisition model, and (b) the linear estimator of the interferometric phase. This estimator applies for both the ML and the DML technique. 5. DML estimate Let us first find a suitable approximation of the inverse of the noise covariance matrix involved in the DML estimate (4). Let us assume the same kernel for both channels, F = F = F that is usually the case in multi-pass interferometry. The model in Fig. 3.a allows us to express the noise covariance as the following Kronecker product: [ ] σ 2 C = σ 2 R F where σ 2 and σ 2 are the noise variances in the two channels, and R F expresses the channel cross-correlation. Eventually we take advantage of the inversion expression for Kronecker products: (A B) = A B 7
8 that leads to the inverse of the channel covariance matrix: [ ] C σ = σ R F [ σ = R F σ R F ] = [ ] C C Having expressed the inverse noise covariance as a block matrix, with two non null blocks, we may write the DML estimator, (8), (2), as follows: where the channel covariance matrix is: M = y C H R H H C y + y C H R H H C y (23) (22) + y C H R H H C y + y C H R H H C y (24) R H = H C H = σ H R F H + σ H R F H (25) 5.. TSML: Fist Step Let us now follow the TSML approach by providing a first estimate of φ by assuming R H = I. The two addends in (23), involve auto-correlations of each of outputs separately, therefore are blind to the interferometric phases and can be dropped. In fact the following expression holds for both i = and i = : y i H i H i y i = y i F i Φ /2 Φ /2 F i y i = y i F i F i y i that no longer depends in φ. The first step of the DML is then the maximization of the summation in (24): M = y C H H C y + y C H H C y (26) = 2 Re(y C H H C y ) (27) = 2 Re(s s ) (28) According to the scheme in Fig. 3.b, we interpret the DML (st step) as the maximization of the cross-energy of two signals, s and s obtained by applying the kernels: G = H C = σ Φ/2 F R F (29) G = H C = σ Φ /2 F R F (3) to the outputs y and y. We may assume that the channel kernels have no phase bias, i.e. F and F are Toeplitz, hence the products P = F R F and P = F R F will still be two Toeplitz matrixes, P and P. We can now derive the DML estimate form (27): φ DMLI = arg max (Re (yp ΦP y )) ( = arg max Re ) exp(jφ n )yp P y n (3) that is maximized by imposing: φ DMLI (n) = arg(y(n)y (n)) = arg I (32) I being the complex interferogram. Not surprisingly, the first step of the DML corresponds to the usual approach, that retrieves the phase, as the argument of the complex interferogram. Notice that this approach is 8
9 optimal as SNR. In the case of constant sloped terrain (3), would be the Fourier Transform of the complex interferogram, evaluated at the frequency φ n /(2π), and the maximization of the likelihood corresponds to the conventional periodogram-based frequency estimation TS-DML: Second Step The second step is performed by maximizing the likelihood: M = 2 Re(y C H R H H C y ) (33) where R H is the pseudo inverse of the channel autocorrelation (25) computed basing on the parameter φ DMLI from (32). Apart from this term, we would get the same result as in (32). The linear estimator in Fig. 3.b still holds, but we now need to modify (29,3) by including R /2 H to both terms. In order to have a reasonable idea of the role played by the term R H, let us assume () a constant sloped terrain, and (2) a stationary case. The system can now be studied as in terms of linear, time variant convolutions. The forward model in Fig. 3.a becomes: y (t) = {γ(t) exp( jωt) + w (t)} f (t) (34) y (t) = {γ(t) exp(jωt) + w (t)} f (t) (35) that is then sampled at a frequency consistent with the bandwidth of the on-board filters, f (t) and f (t). The frequency domain representation of the spectra Y (f) and Y (f) are shown in Fig. 4 (the two upper plots). Figure 4: Single baseline model in the frequency domain, for a constant sloped terrain. Notice the spectral shift, f, that depends on the local slope and the baseline. We have assumed ideal boxcar bandpass spectra, centered on the carrier frequency. Under this assumption, the autocorrelation of the two channels will have the same bandpass shape, then we will approximate the inverse noise covariance (22) as follows: [ C σ = F σ ] (36) 9
10 Eventually, we get from (25): R H = σ ( ) ( Φ /2 F F F F Φ /2) + σ ( ) ( Φ /2 F F F F Φ /2) (37) In the constant slope case, the inverse kernel R H will perform a whitening, e.g. the equalization in the frequency domain, shown in Fig. 4, last plot. Even if in this specific case, we do not expect any difference with the result achieved at the first step, this equalization would be effective in the Multi-Baseline (SIMO) environment, when many channels could contribute to the same spectrum. 5.2 ML estimate Let us now find a suitable approximation for the ML estimate in (9). We first introduce a matrix in the complementary space of the signal covariance: G = I C y, the ML is now the following maximization The inverse of the matrix G can be expressed as the block matrix: [ ] G G G = G G The ML estimate is then: φ ML = arg max ( y G y ) (38) Φ = arg max (y G y + y G y + y G y + y G y ) = arg max (y G y + 2 Re (y G y ) + y G y ) (39) The likelihood to be maximized is similar in topology (but non identical in weights) to the one in the DML estimate, (23). Instead of inverting the covariance matrix, we can directly derive the ML by finding the weights and G ij that maximizes (39). We observe that the solution is already known in the literature, in fact, in [] the linear model in Fig. 3.b: s = G y and s = G y (4) was exploited to minimize the squared error: [ min y G y 2 + y G y 2] (4) G,G [ ] y = max G G y + yg y + yg y yg G y + yg y + yg y that is formally equivalent to the likelihood (39). The minimization of (4) gives the result (see []) : G + SNR F Φ /2 F (F F ) (42) + SNR F Φ /2 F (43) G = + SNR F Φ /2 F (44) where SNR, SNR are the ratio between the noise power (in the signal bandwidth) and the signal power in each channel, and (42) is an approximation that gives almost the same results if the channels F, F have flat frequency response in their pass-band. 5.3 Deterministic and statistic likelihood
11 We have just shown that both DML and ML estimates are performed by means of linear kernels, whose expressions (29,3) and (43,44) appears quite similar. The most outstanding difference is due to the inclusion, in the ML case, of the weighting factor: ρ = + SNR + SNR that corresponds to the definition of the coherence of the interferometric pair [6]. The channel decorrelation is not handled in the DML approach, as this one regards the sources as set of parameters to be estimated rather than the realization of a statistical process. A qualitative comparison of the two techniques is provided in Fig. 5. In the figure, the upper plots show the squared absolute value of the Fourier Transform of the interferogram, that we have shown to be equivalent to the DML approach. ML DML Coherence Frequency Figure 5: Estimate of the fringe frequency made by () the DML technique (upper plot), (2) the ML technique (middle plot). All frequencies are normalized wrt 2f s. The correct frequency (.35) is retrieved in both cases, however, the ML technique, that exploits the coherence (lower plot) in its Likelihood, appears to be more robust. The ML estimator was applied following the idea shown in Fig. 3.b. In theory, a loop would be required where frequencies f are evaluated on a grid as fine as the one provided by the FFT (the DML approach). For each frequency, both the master and the slave are filtered according to the spectral shift filtering (43,44), to keep their common band contributions (see Fig.4). The likelihood is then estimated as the cross-energy of the filtered signals multiplied times the coherence. We can estimate the coherence by the sampled estimator: ρ bφ = s Φs (s s ) (s s ) where Φ is the modulation matrix evaluated at the current frequency, f. In practice, this loop can be made much faster by using a coarser grid for the spectral shift filtering, and refining the estimate by performing a FFT on the interferogram: I(n) = s (n) s (n). As an example, Fig. 5 (middle plot) shows the Fourier Transform of the interferogram after filtering the two images in the case of the correct spectra shift, e.g. f = f. The lower plot shows the coherence estimated at the various spectra shifts. A comparison of the plots in Fig. 5, motivates the following observations: both the FT of the interferogram and the spectral shift approach give the same results if the search is limited to a narrow interval of the correct value: f = f ± ɛ
12 the spectral shift approach is more robust due to the use of the coherence in the likelihood. We then notice that the inclusion of coherence in the likelihood allows to better discriminate which of the spectral peaks in Fig. 3 is the true one, whereas the coherence peak is too flat to provide an improvement in the fine estimate f. We once again remind that the weakness of the DML is due to the higher number of parameters to be estimated. We than expect that, at the highest SNR, and for large estimation windows, the two techniques would provide the same results, whereas different results would be achieved at the low SNR. A practical idea of this limit value is provided in Fig. 5. The figure shows the standard deviation of the estimate of f as a function of the coherence and for different window sizes, L, achieved by means of a Monte- Carlo simulation, with the DML approach. For high coherence/snr, all the estimates gives the same accuracy, once that they are scaled times the square root of the window length, L. However, for the low SNR, the performances drops and the estimated frequency is biased towards, a fact well known in literature [2]. In that case, the ML estimate would get significant improvements. σ f L=5 L = L=2 L=5 L= Coherence Figure 6: Standard deviation of the frequency estimate made by the DML approach, for different size of the estimation window, as a function of the coherence. Notice the threshold behavior for coherences ranging form.35 (L=5) to.95 (L=). 6 Multi-Baseline Slope Estimate The multi-baseline case is more complicated than the single baseline case. In fact, although we are still able to extend the DML approach, we have no closed form solution for the ML in (9) that would avoid the inversion of the covariance matrix. We than propose a straight extension of (26), by simply adding coherence-based weights to get close to the ML results. We need then to maximize the following figure of merit: φ = argmax N j= N i= ρ ij (φ)y i G ij y j (45) where the weights are computed basing on (29), and assuming the noiseless case (as we already accounts for noise with the coherence): G i = Φ /2 ij F i The diagonal matrixes Φ ij keeps the elements of φ in (6) properly scaled by the baseline of the pair i, j. The maximum is then found by iterating for all the values of the normalized frequency shift, f. Furthermore the co-channel terms (for i = j), are useless for the estimate of the phase term φ, and they can be dropped from the summation (45). This can be shown by writing the Fourier Transform of the SAR channel (34): Y i (f) = {Γ(f) δ(f f(α(t)))} F i (f) = Γ(f f) F i (f) (46) 2
13 where [ f(α(t) is the spectral shift, that depends on the slope. The power spectrum of the reflectivity, E Γ(f) 2], is white over a large bandwith, therefore the power spectrum of the output, E [ Y i (f) 2] carries no information on f. Nor we can retrieve any information from the phase of Y i (f), as the sources are Gaussian. Having dropped the co-channel contribution, (45) reduces to the summation of averaged interferograms (i j): ρ ij (φ)y i G ij y j Attention should be taken in order to exclude, at each value of the parameter f, the contributions that would be aliased. The constraint we have to impose is indeed more stringent, as we want to avoid biasing in the coherence estimate. Given a spectral shift f, the common bandwidth would be B c = B f, where B i is the SAR channel bandwidth. If we exploit patches on the ground of L pixels, sampled at a rate f s, the number of effectively independent samples would be: N i = L B c f s = L ( ) B n B max f f s R tan(θ α) and we need N i to avoid biasing of the coherence estimate ρ ij. Finally, we remark that in the straight extension from the single baseline, (42), to the MB case (45), we have not accounted for a spectral equalization like the one provided by the operator R H (see Fig. 4), however this topic could be a subject for future analysis. 6. Experimental results The proposed MB InSAR estimation algorithm has been tested first with synthetic and then with real data. The estimation algorithm has been applied with respect to the pixel-to-pixel phase differences in range, in order to remove any potential problem coming from tropospheric phase screen or other unknown phase offsets. We selected a 5 samples window length as a compromise between accuracy and resolution. The algorithms tested were based on the combination (45) in the two following implementations: - DML unweighed, - ML, including the spectral shift filtering and weights derived by the a-posteriori estimate of the coherence. In presenting the results we have also added a single interferogram image for comparison. 6.2 Results from simulations A set of 6 images was simulated by assuming the parameters of ENVISAT-ERS systems and the topography of Mt. Vesuvius. The baselines were chosen to cover a range of m with uniform sampling but taking care to avoid periodicity. The baseline vector is [ ] m with respect to the master. We did not add any further decorrelation, apart from the unavoidable baseline decorrelation. Having the perfect knowledge of the topography, we could perform an accurate evaluation of the estimator performances for different slopes. These results are shown in Fig. 7. Fig. 7.a maps the local phase gradient in range expressed in radians and measured on the theoretical unwrapped phase that would come from the 47 m baseline interferogram.the other graphics in the same figure show the natural logartihm of the counts of measured slopes with respect to the true slope, a sort of scatter-plot. The estimates are based on a single interferogram, in Fig. 7.(b), on the DML estimator, Fig. 7.(c) and on the proposed ML-derived estimator, Fig. 7.(d). The single interferogram is clearly limited in the trade-offs between altimetric resolution and phase aliasing. In the case shown we selected a small baseline ( m), nonetheless there is a sensible loss of the contribution coming from the higher slopes. We also notice the large variance, due to the fact that only two images were exploited instead of the five available. The DML estimate, in Fig. 7.(c) is almost useless, due to the strong biasing of the estimate, as discussed in section 5.3. On the other hand, the estimate achieved by the proposed technique, shown in Fig. 7.d, provides the higher accuracy for the low values of the gradient, due to the cooperation of all the images, and the capability to solve the highest gradient with no aliasing, due to the combinations with the lowest baselines. 3
14 φ e (a) φ 2 (b) φ e 3 2 φ e φ (c) φ (d) 4 6 Figure 7: (a) Phase gradient φ, along slant range, of a 5 m interferogram on Mt. Vesuvius. (b-d) Scatter diagram of φ estimated versus true value for different estimator. (b) single-baseline interferogram, (c) MB based on the unweighted combination, (d) MB achieved by both exploiting coherence weights and spectral shift filtering. 4
15 6.3 Real data processing The same technique has been exploited in a real ENVISAT dataset acquired over the area of Las Vegas. In particular we have considered the data relative to track 356, frame Acquisition parameters, in terms of perpendicular baseline and temporal separation, are collected in Table. Name Data Baseline Normal [m] Baseline Parallel [m] Delta Days M 25-Oct-22 S 29-Nov S2 3-Jan S3 7-Feb S4 27-Jun S5 -Aug Table : ENVISAT datasets, Las Vegas test site. The estimation algorithm has been applied, with respect to the pixel-to-pixel phase differences in range, over an image patch of approximately pixels 2 km (range, azimuth). An external DEM of the SRTM mission have been used with different purposes: - to remove aliasing in the single baseline interferogram, - to avoid threshold effects, and slope biasing in the DML approach, - to reduce the search range, hence speed up computations in the ML estimate. A smoothing window of 3 5 pixels in range and azimuth was used. The reference baseline used for referring phase gradient measures is 237 m. The slope measured on the reference SRTM DEM are reported in Fig. 8.a, together with the single baseline estimate, 8.b and the proposed ML combination, 8.c. Figure 8: Phase gradient, referred to a 237 m baseline image and estimated from: (left) an SRTM DEM, (center) a single baseline interferogram, (right) the proposed ML combination. 5
16 Here again we find the same features as in the simulated results, although the altimetric resolution is less evident due to the smaller baseline span. Notice the ML combination gets result similar to the SRTM DEM for both robustness with respect to ambiguity and the accuracy. 6.4 Acknowledgments The author would like to thank the Italian Space Agency for sponsoring of the work (ASI contract I/R/69/2 Multiple pass SAR Interferometry for very accurate Earth modeling). Thanks also to Dr. Ing. F. De Zan for the proficuous discussions on the statistic of the ML approach. 7 Conclusions In the paper we have studied two different Maximum Likelihood techniques to provide estimates of the terrain topography in a Multi-Baseline SAR interferometric case, namely the Deterministic ML and the statistic ML. The first approach includes the sources in the vector of parameters, and thus it is less robust than the second (particularly when few samples are available). However, its quite efficient and suited for applications when a preliminary knowledge of the slope exists. The second estimator, the ML approach, was derived as an extension of the MB case the exact ML estimated for conventional interferometry. A preliminary validation of this technique has demonstrated the capabilities of solving aliasing for the higher slopes, due to the use of the lower baselines, and to attain at the same time the high altimetric resolution that would come from the higher baselines. The major limit of the technique, the computational complexity, leaves open opportunities for further researches. References [] F. K. Li and R. M. Goldstein, Studies of multibaseline spaceborne interferometric synthetic aperture radars, IEEE Transactions on Geoscience and Remote Sensing, vol. 28, no., pp , Jan. 99. [2] D. Massonnet, H. Vadon, and M. Rossi, Reduction of the need for phase unwrapping in radar interferometry, IEEE Transactions on Geoscience and Remote Sensing, vol. 34, no. 2, pp , Mar [3] R. Lanari, G. Fornaro, D. Riccio, M. Migliaccio, K. P. Papathanassiou, J. ao R Moreira, M. Schwäbisch, L. Dutra, G. Puglisi, G. Franceschetti, and M. Coltelli, Generation of digital elevation models by using SIR-C/X-SAR multifrequency two-pass interferometry: The Etna case study, IEEE Transactions on Geoscience and Remote Sensing, vol. 34, no. 5, pp. 97 4, Sept [4] A. Ferretti, C. Prati, F. Rocca, and A. Monti Guarnieri, Multibaseline SAR interferometry for automatic DEM reconstruction, in Third ERS Symposium Space at the Service of our Environment, Florence, Italy, 7 2 March 997, ser. ESA SP-44, 997, pp [Online]. Available: [5] F. Lombardini, Optimal absolute phase retrieval in three-element SAR interferometry, Electronics Letters, vol. 34, no. 5, pp , July 998. [6] V. Pascazio and G. Schirinzi, Estimation of terrain elevation by multifrequency interferometric wide band SAR data, IEEE Signal Processing Lett., vol. 8, pp. 7 9, 2. [7] A. Ferretti, C. Prati, and F. Rocca, Permanent scatterers in SAR interferometry, IEEE Transactions on Geoscience and Remote Sensing, vol. 39, no., pp. 8 2, Jan. 2. [8] F. Gatelli, A. Monti Guarnieri, F. Parizzi, P. Pasquali, C. Prati, and F. Rocca, The wavenumber shift in SAR interferometry, IEEE Transactions on Geoscience and Remote Sensing, vol. 32, no. 4, pp , July 994. [9] G. W. Davidson and R. Bamler, Multiresolution phase unwrapping for SAR interferometry, IEEE Transactions on Geoscience and Remote Sensing, vol. 37, no., pp , Jan
17 [] A. M. Guarnieri and F. Rocca, Combination of low- and high-resolution SAR images for differential interferometry, IEEE Transactions on Geoscience and Remote Sensing, vol. 37, no. 4, pp , July 999. [] G. Fornarno and A. M. Guarnieri, Minimum mean square error space-varying filtering of interferometric SAR data, IEEE Transactions on Geoscience and Remote Sensing, vol. 4, no., pp. 2, July 22. [2] R. Bamler and P. Hartl, Synthetic aperture radar interferometry, Inverse Problems, vol. 4, pp. R R54, 998. [3] P. Rosen, S. Hensley, I. R. Joughin, F. K. Li, S. Madsen, E. Rodríguez, and R. Goldstein, Synthetic aperture radar interferometry, Proceedings of the IEEE, vol. 88, no. 3, pp , Mar. 2. [4] A. M. Guarnieri, SAR interferometry and statistical topography, IEEE Trans. Geosc. Remote Sens., vol. 4, no. 2, pp , Dec. 22. [5] H. A. Zebker and J. Villasenor, Decorrelation in interferometric radar echoes, IEEE Transactions on Geoscience and Remote Sensing, vol. 3, no. 5, pp , sept 992. [6] D. Just and R. Bamler, Phase statistics of interferograms with applications to synthetic aperture radar, Applied Optics, vol. 33, no. 2, pp , 994. [7] L. Tong and S. Perreau, Multichannel blind identification: From subspace to maximum likelihood methods, Proc. IEEE, vol. 86, no., pp , Oct [8] F. De-Zan, Stima ottima della fase nell interferometria SAR, Ph.D. dissertation, Politecnico di Milano, July 24. [9] Y. Hua, Fast maximum likelihood for blind identification of multiple FIR channels, IEEE Transactions on Signal Processing, vol. 44, pp , Mar [2] R. Bamler, N. Adam, G. W. Davidson, and D. Just, Noise-induced slope distortion in 2-d phase unwrapping by linear estimators with application to SAR interferometry, IEEE Transactions on Geoscience and Remote Sensing, vol. 36, no. 3, pp , May
Multi-Baseline SAR interferometry
ulti-baseline SAR interferometry A. onti Guarnieri, S. Tebaldini Dipartimento di Elettronica e Informazione - Politecnico di ilano G. Fornaro, A. Pauciullo IREA-CNR ESA cat-1 project, 3173. Page 1 Classical
More informationMODELING INTERFEROGRAM STACKS FOR SENTINEL - 1
MODELING INTERFEROGRAM STACKS FOR SENTINEL - 1 Fabio Rocca (1) (1) Politecnico di Milano, Via Ponzio 34/5, 20133 Milano, Italy, Email: rocca@elet.polimi.it ABSTRACT The dispersion of the optimal estimate
More informationDEM GENERATION AND ANALYSIS ON RUGGED TERRAIN USING ENVISAT/ASAR ENVISAT/ASAR MULTI-ANGLE INSAR DATA
DEM GENERATION AND ANALYSIS ON RUGGED TERRAIN USING ENVISAT/ASAR ENVISAT/ASAR MULTI-ANGLE INSAR DATA Li xinwu Guo Huadong Li Zhen State Key Laboratory of Remote Sensing Science, Institute of Remote Sensing
More informationERS-ENVISAT Cross-interferometry for Coastal DEM Construction
ERS-ENVISAT Cross-interferometry for Coastal DEM Construction Sang-Hoon Hong and Joong-Sun Won Department of Earth System Sciences, Yonsei University, 134 Shinchon-dong, Seodaemun-gu, 120-749, Seoul, Korea
More informationPHASE UNWRAPPING. Sept. 3, 2007 Lecture D1Lb4 Interferometry: Phase unwrapping Rocca
PHASE UNWRAPPING 1 Phase unwrapping 2 1D Phase unwrapping Problem: given the wrapped phase ψ=w(φ) find the unwrapped one ψ. The wrapping operator: W(φ)=angle(exp(j φ)), gives always solution -π π and is
More informationSAR Data Analysis: An Useful Tool for Urban Areas Applications
SAR Data Analysis: An Useful Tool for Urban Areas Applications M. Ferri, A. Fanelli, A. Siciliano, A. Vitale Dipartimento di Scienza e Ingegneria dello Spazio Luigi G. Napolitano Università degli Studi
More informationA Radar Eye on the Moon: Potentials and Limitations for Earth Imaging
PIERS ONLINE, VOL. 6, NO. 4, 2010 330 A Radar Eye on the Moon: Potentials and Limitations for Earth Imaging M. Calamia 1, G. Fornaro 2, G. Franceschetti 3, F. Lombardini 4, and A. Mori 1 1 Dipartimento
More informationOn the Exploitation of Target Statistics for SAR Interferometry Applications
On the Exploitation of Target Statistics for SAR Interferometry Applications A. Monti Guarnieri, S. Tebaldini Politecnico di Milano Abstract This paper focuses on multi-image Synthetic Aperture Radar Interferometry
More informationJournal of Geodynamics
Journal of Geodynamics 49 (2010) 161 170 Contents lists available at ScienceDirect Journal of Geodynamics journal homepage: http://www.elsevier.com/locate/jog Recent advances on surface ground deformation
More informationERS-ENVISAT CROSS-INTERFEROMETRY SIGNATURES OVER DESERTS. Urs Wegmüller, Maurizio Santoro and Christian Mätzler
ERS-ENVISAT CROSS-INTERFEROMETRY SIGNATURES OVER DESERTS Urs Wegmüller, Maurizio Santoro and Christian Mätzler Gamma Remote Sensing AG, Worbstrasse 225, CH-3073 Gümligen, Switzerland, http://www.gamma-rs.ch,
More informationCHAPTER-7 INTERFEROMETRIC ANALYSIS OF SPACEBORNE ENVISAT-ASAR DATA FOR VEGETATION CLASSIFICATION
147 CHAPTER-7 INTERFEROMETRIC ANALYSIS OF SPACEBORNE ENVISAT-ASAR DATA FOR VEGETATION CLASSIFICATION 7.1 INTRODUCTION: Interferometric synthetic aperture radar (InSAR) is a rapidly evolving SAR remote
More informationSAR interferometry Status and future directions. Rüdiger Gens
SAR interferometry Status and future directions Rüdiger Gens Polarimetric InSAR Polarimetric InSAR InSAR - Status and future directions sensitivity to changes in surface scattering, even in the presence
More informationDeformation measurement using SAR interferometry: quantitative aspects
Deformation measurement using SAR interferometry: quantitative aspects Michele Crosetto (1), Erlinda Biescas (1), Ismael Fernández (1), Ivan Torrobella (1), Bruno Crippa (2) (1) (2) Institute of Geomatics,
More informationSlow Deformation of Mt. Baekdu Stratovolcano Observed by Satellite Radar Interferometry
Slow Deformation of Mt. Baekdu Stratovolcano Observed by Satellite Radar Interferometry Sang-Wan Kim and Joong-Sun Won Department of Earth System Sciences, Yonsei University 134 Shinchon-dong, Seodaemun-gu,
More informationModeling of Atmospheric Effects on InSAR Measurements With the Method of Stochastic Simulation
Modeling of Atmospheric Effects on InSAR Measurements With the Method of Stochastic Simulation Z. W. LI, X. L. DING Department of Land Surveying and Geo-Informatics, Hong Kong Polytechnic University, Hung
More informationP079 First Results from Spaceborne Radar Interferometry for the Study of Ground Displacements in Urban Areas SUMMARY
P079 First Results from Spaceborne Radar Interferometry for the Study of Ground Displacements in Urban Areas C.M. Crosetto (Instituto de Geomatica), C.A. Casas (University of Barcelona), R.G. Ranieri (University
More informationMethods and performances for multi-pass SAR Interferometry
Methods and performances for multi-pass SAR Interferometry Methods and performances for multi-pass SAR Interferometry Stefano Tebaldini and Andrea Monti Guarnieri Dipartimento di Elettronica e Informazione
More informationDIFFERENTIAL INSAR STUDIES IN THE BOREAL FOREST ZONE IN FINLAND
DIFFERENTIAL INSAR STUDIES IN THE BOREAL FOREST ZONE IN FINLAND Kirsi Karila (1,2), Mika Karjalainen (1), Juha Hyyppä (1) (1) Finnish Geodetic Institute, P.O. Box 15, FIN-02431 Masala, Finland, Email:
More informationALOS PI Symposium 2009, 9-13 Nov 2009 Hawaii MOTION MONITORING FOR ETNA USING ALOS PALSAR TIME SERIES
ALOS PI Symposium 2009, 9-13 Nov 2009 Hawaii ALOS Data Nodes: ALOS RA-094 and RA-175 (JAXA) MOTION MONITORING FOR ETNA USING ALOS PALSAR TIME SERIES Urs Wegmüller, Charles Werner and Maurizio Santoro Gamma
More informationAtmospheric Phase Screen (APS) estimation and modeling for radar interferometry
Atmospheric Phase Screen (APS) estimation and modeling for radar interferometry Ramon Hanssen¹, Alessandro Ferretti², Marco Bianchi², Rossen Grebenitcharsky¹, Frank Kleijer¹, Ayman Elawar¹ ¹ Delft University
More informationMULTILOOKING is an essential procedure for improving
IEEE GEOSCIENCE AND REMOTE SENSING LETTERS, VOL. 8, NO. 3, MAY 011 441 Adaptive InSAR Stack Multilooking Exploiting Amplitude Statistics: A Comparison Between Different Techniques and Practical Results
More informationBUILDING HEIGHT ESTIMATION USING MULTIBASELINE L-BAND SAR DATA AND POLARIMETRIC WEIGHTED SUBSPACE FITTING METHODS
BUILDING HEIGHT ESTIMATION USING MULTIBASELINE L-BAND SAR DATA AND POLARIMETRIC WEIGHTED SUBSPACE FITTING METHODS Yue Huang, Laurent Ferro-Famil University of Rennes 1, Institute of Telecommunication and
More informationVALIDATION OF THE PERMANENT SCATTERERS TECHNIQUE IN URBAN AREAS
VALIDATION OF THE PERMANENT SCATTERERS TECHNIQUE IN URBAN AREAS Alessandro Ferretti, Claudio Prati, Fabio Rocca, Carlo Colesanti Dipartimento di Elettonica e Informazione Politecnico di Milano Piazza L.
More informationMeasuring rock glacier surface deformation using SAR interferometry
Permafrost, Phillips, Springman & Arenson (eds) 2003 Swets & Zeitlinger, Lisse, ISBN 90 5809 582 7 Measuring rock glacier surface deformation using SAR interferometry L.W. Kenyi Institute for Digital Image
More informationModel based forest height estimation with ALOS/PalSAR: A first study.
Model based forest height estimation with ALOS/PalSAR: A first study. K.P. Papathanassiou*, I. Hajnsek*, T.Mette*, S.R. Cloude** and A. Moreira* * (DLR) (DLR-HR) Oberpfaffenhofen, Germany ** AEL Consultants
More informationAPPEARANCE OF PERSISTENT SCATTERERS FOR DIFFERENT TERRASAR-X ACQUISITION MODES
APPEARANCE OF PERSISTENT SCATTERERS FOR DIFFERENT TERRASAR-X ACQUISITION MODES Stefan Gernhardt a, Nico Adam b, Stefan Hinz c, Richard Bamler a,b a Remote Sensing Technology, TU München, Germany b Remote
More informationMONITORING THE CAMPI FLEGREI CALDERA BY EXPLOITING SAR AND GEODETICAL DATA: RECENT RESULTS AND FUTURE APPLICATIONS
MONITORING THE CAMPI FLEGREI CALDERA BY EXPLOITING SAR AND GEODETICAL DATA: RECENT RESULTS AND FUTURE APPLICATIONS P. Berardino (1), S. Borgström (2), I. Aquino (2), F. Casu (1), C. Del Gaudio (2), R.
More informationSNOW MASS RETRIEVAL BY MEANS OF SAR INTERFEROMETRY
SNOW MASS RETRIEVAL BY MEANS OF SAR INTERFEROMETRY Helmut Rott (1), Thomas Nagler (1), Rolf Scheiber (2) (1) ENVEO, Environmental Earth Observation OEG, Exlgasse 39, A-6020 Innsbruck, Austria E-mail: Helmut.Rott@enveo.at
More informationGeneration and Validation of Digital Elevation Model using ERS - SAR Interferometry Remote Sensing Data
Jour. Agric. Physics, Vol. 7, pp. 8-13 (2007) Generation and Validation of Digital Elevation Model using ERS - SAR Interferometry Remote Sensing Data SHELTON PADUA 1, VINAY K. SEHGAL 2 AND K.S. SUNDARA
More informationERS Track 98 SAR Data and InSAR Pairs Used in the Analysis
ERS Track 98 SAR Data and InSAR Pairs Used in the Analysis Date 1 Date 2 Date 1 Date 2 Date 1 Date 2 Date 1 Date 2 7/17/1992 6/19/2000 7/17/1992 7/2/1993 9/10/1993 10/28/1996 9/3/1995 10/18/1999 9/25/1992
More informationGround deformation monitoring in Pearl River Delta region with Stacking D-InSAR technique
Ground deformation monitoring in Pearl River Delta region with Stacking DInSAR technique Zhao Qing*, LIN Hui, JIANG Liming Institute of Space and Earth Information Science, Room 615, Esther Lee Building,
More informationGLOBAL FOREST CLASSIFICATION FROM TANDEM-X INTERFEROMETRIC DATA: POTENTIALS AND FIRST RESULTS
GLOBAL FOREST CLASSIFICATION FROM TANDEM-X INTERFEROMETRIC DATA: POTENTIALS AND FIRST RESULTS Michele Martone, Paola Rizzoli, Benjamin Bräutigam, Gerhard Krieger Microwaves and Radar Institute German Aerospace
More informationDEMONSTRATION OF TERRASAR-X SCANSAR PERSISTENT SCATTERER INTERFEROMETRY
DEMONSTRATION OF TERRASAR-X SCANSAR PERSISTENT SCATTERER INTERFEROMETRY Fernando Rodriguez Gonzalez (1) Ramon Brcic (1) Nestor Yague-Martinez (1) Robert Shau (1) Alessandro Parizzi (1) Nico Adam (1) (1)
More informationThe Potential of High Resolution Satellite Interferometry for Monitoring Enhanced Oil Recovery
The Potential of High Resolution Satellite Interferometry for Monitoring Enhanced Oil Recovery Urs Wegmüller a Lutz Petrat b Karsten Zimmermann c Issa al Quseimi d 1 Introduction Over the last years land
More informationTHREE DIMENSIONAL DETECTION OF VOLCANIC DEPOSIT ON MOUNT MAYON USING SAR INTERFEROMETRY
ABSTRACT THREE DIMENSIONAL DETECTION OF VOLCANIC DEPOSIT ON MOUNT MAYON USING SAR INTERFEROMETRY Francis X.J. Canisius, Kiyoshi Honda, Mitsuharu Tokunaga and Shunji Murai Space Technology Application and
More informationSAR Monitoring of Progressive and Seasonal Ground Deformation Using the Permanent Scatterers Technique
IEEE TRANSACTIONS ON GEOSCIENCE AND REMOTE SENSING, VOL. 41, NO. 7, JULY 2003 1685 SAR Monitoring of Progressive and Seasonal Ground Deformation Using the Permanent Scatterers Technique Carlo Colesanti,
More informationMAPPING DEFORMATION OF MAN-MADE LINEAR FEATURES USING DINSAR TECHNIQUE
MAPPING DEFORMATION OF MAN-MADE LINEAR FEATURES USING DINSAR TECHNIQUE H. Wu a, *, Y. Zhang a, J. Zhang a, X. Chen b a Key Laboratory of Mapping from Space of State Bureau of Surveying and Mapping, Chinese
More informationSnow Water Equivalent (SWE) of dry snow derived from InSAR -theory and results from ERS Tandem SAR data
Snow Water Equivalent (SWE) of dry snow derived from InSAR -theory and results from ERS Tandem SAR data Tore Guneriussen, Kjell Arild Høgda, Harald Johnsen and Inge Lauknes NORUT IT Ltd., Tromsø Science
More informationApplication of differential SAR interferometry for studying eruptive event of 22 July 1998 at Mt. Etna. Abstract
Application of differential SAR interferometry for studying eruptive event of 22 July 1998 at Mt. Etna Coltelli M. 1, Puglisi G. 1, Guglielmino F. 1, Palano M. 2 1 Istituto Nazionale di Geofisica e Vulcanologia,
More informationDETECTING ICE MOTION IN GROVE MOUNTAINS, EAST ANTARCTICA WITH ALOS/PALSAR AND ENVISAT/ASAR DATA
DETECTING ICE MOTION IN GROVE MOUNTAINS, EAST ANTARCTICA WITH ALOS/PALSAR AND ENVISAT/ASAR DATA TIAN Xin (1), LIAO Mingsheng (1), ZHOU Chunxia (2), ZHOU Yu (3) (1) State Key Laboratory of Information Engineering
More informationNoise covariance model for time-series InSAR analysis. Piyush Agram, Mark Simons
Noise covariance model for time-series InSAR analysis Piyush Agram, Mark Simons Motivation Noise covariance in individual interferograms has been well studied. (Hanssen, 2001) Extend the noise model to
More informationGround deformation monitoring at the Phlegrean Fields (Naples, Italy) from the exploitation of SAR data in the framework of CAT-1 and DUP activities
Ground deformation monitoring at the Phlegrean Fields (Naples, Italy) from the exploitation of SAR data in the framework of CAT-1 and DUP activities Borgström S., Aquino I., Del Gaudio C., Ricco C., Siniscalchi
More informationGEOMORPHOLOGICAL MAPPING WITH RESPECT TO AMPLITUDE, COHERENCEAND PHASE INFORMATION OF ERS SAR TANDEM PAIR
GEOMORPHOLOGICAL MAPPING WITH RESPECT TO AMPLITUDE, COHERENCEAND PHASE INFORMATION OF ERS SAR TANDEM PAIR AUNG LWIN Assistant Researcher Remote Sensing Department Mandalay Technological University, Myanmar
More informationRadar Remote Sensing: Monitoring Ground Deformations and Geohazards from Space
Radar Remote Sensing: Monitoring Ground Deformations and Geohazards from Space Xiaoli Ding Department of Land Surveying and Geo-Informatics The Hong Kong Polytechnic University A Question 100 km 100 km
More informationSubsidence-induced fault
Surveying Monitoring underground coal mining-induced subsidence by Yaobin Sheng, Linlin Ge, Chris Rizos, University of New South Wales, and Yunjia Wang, China University of Mining and Technology This paper
More informationTEMPORAL VARIABILITY OF ICE FLOW ON HOFSJÖKULL, ICELAND, OBSERVED BY ERS SAR INTERFEROMETRY
TEMPORAL VARIABILITY OF ICE FLOW ON HOFSJÖKULL, ICELAND, OBSERVED BY ERS SAR INTERFEROMETRY Florian Müller (1), Helmut Rott (2) (1) ENVEO IT, Environmental Earth Observation GmbH, Technikerstrasse 21a,
More informationAnalysis of the Temporal Behavior of Coherent Scatterers (CSs) in ALOS PalSAR Data
Analysis of the Temporal Behavior of Coherent Scatterers (CSs) in ALOS PalSAR Data L. Marotti, R. Zandona-Schneider & K.P. Papathanassiou German Aerospace Center (DLR) Microwaves and Radar Institute0 PO.BOX
More informationMitigation of Atmospheric Water-vapour Effects on Spaceborne Interferometric SAR Imaging through the MM5 Numerical Model
PIERS ONLINE, VOL. 6, NO. 3, 2010 262 Mitigation of Atmospheric Water-vapour Effects on Spaceborne Interferometric SAR Imaging through the MM5 Numerical Model D. Perissin 1, E. Pichelli 2, R. Ferretti
More informationUniversità di Roma La Sapienza, Facoltà d Ingegneria, Dipartimento di Idraulica,Trasporti e Strade (D.I.T.S.) 00185, Roma, Italy
DORIS FP7-EU PROJECT: EXPLOITATION OF 20 YEARS DINSAR DATA ARCHIVE FOR LANDSLIDE MONITORING Manunta M. (1), Calò F. (1), Paglia L., (1), Bonano M. (1,2), Lanari R. (1) (1) IREA-CNR, Via Diocleziano 328,
More informationCLOSE FORMATION FLIGHT OF PASSIVE RECEIVING MICRO-SATELLITES
CLOSE FORMATION FLIGHT OF PASSIVE RECEIVING MICRO-SATELLITES Hauke Fiedler (1) and Gerhard Krieger (2) (1) DLR, HR, Weßling, Germany, E-mail:hauke.fiedler@dlr.de (2) DLR, HR, Weßling, Germany, E-mail:gerhard.krieger@dlr.de
More informationA New Approach to Estimate Forest Parameters Using Dual-Baseline POL-InSAR Data
Jan 26-30, 2009 Frascati, Italy A New Approach to Estimate Forest Parameters Using Dual-Baseline POL-InSAR Data Lu Bai, Wen ong, Fang Cao, Yongsheng Zhou bailu8@gmail.com fcao@mail.ie.ac.cn National Key
More informationSupplementary information. Analytical Techniques and Measurement Uncertainty
Supplementary information Analytical Techniques and Measurement Uncertainty Interferomeric SAR (InSAR) allows measurement of millimetre-level surface displacements, including land subsidence, in the radar
More informationIn order to obtain a long term monitoring result for the Kilauea Volcano, ALOS PALSAR images taken on Track 287, Frame 38, ascending orbit with 21.5 d
ALOS PALSAR OBSERVATION OF KILAUEA VOLCANO ACTIVITIES FROM 2006 TO 2009 Zhe Hu, Linlin Ge, Xiaojing Li, Kui Zhang, Alex Hay-Man NG and Chris Rizos Cooperative Research Centre for Spatial Information &
More informationHigh-resolution temporal imaging of. Howard Zebker
High-resolution temporal imaging of crustal deformation using InSAR Howard Zebker Stanford University InSAR Prehistory SEASAT Topographic Fringes SEASAT Deformation ERS Earthquake Image Accurate imaging
More informationInvestigating landslides with space-borne Synthetic Aperture Radar (SAR) interferometry
Engineering Geology 88 (2006) 173 199 www.elsevier.com/locate/enggeo Investigating landslides with space-borne Synthetic Aperture Radar (SAR) interferometry Carlo Colesanti 1, Janusz Wasowski a, a CNR-IRPI,
More informationUrban land and infrastructure deformation monitoring by satellite radar interferometry
Urban land and infrastructure deformation monitoring by satellite radar interferometry Lei Zhang and Xiaoli Ding Department of Land Surveying and Geo-Informatics (LSGI) The Hong Kong Polytechnic University
More informationSpatiotemporal analysis of ground deformation at Campi Flegrei and Mt Vesuvius, Italy, observed by Envisat and Radarsat-2 InSAR during
Spatiotemporal analysis of ground deformation at Campi Flegrei and Mt Vesuvius, Italy, observed by Envisat and Radarsat InSAR during 233 Sergey V. Samsonov, Pablo J. González, Kristy F. Tiampo, Antonio
More informationTowards an Optimal Noise Versus Resolution Trade-off in Wind Scatterometry
Towards an Optimal Noise Versus Resolution Trade-off in Wind Scatterometry Brent Williams Jet Propulsion Lab, California Institute of Technology IOWVST Meeting Utrecht Netherlands June 12, 2012 Copyright
More informationAvailable online at GHGT-9. Detection of surface deformation related with CO 2 injection by DInSAR at In Salah, Algeria
Available online at www.sciencedirect.com Energy Procedia 100 (2009) (2008) 2177 2184 000 000 Energy Procedia www.elsevier.com/locate/procedia www.elsevier.com/locate/xxx GHGT-9 Detection of surface deformation
More informationAnalysis of ERS Tandem SAR Coherence From Glaciers, Valleys, and Fjord Ice on Svalbard
IEEE TRANSACTIONS ON GEOSCIENCE AND REMOTE SENSING, VOL. 39, NO. 9, SEPTEMBER 2001 2029 Analysis of ERS Tandem SAR Coherence From Glaciers, Valleys, and Fjord Ice on Svalbard Dan Johan Weydahl Abstract
More informationTHE USE OF DIFFERENT REMOTE SENSING TECHNIQUES FOR LANDSLIDE CHARACTERIZATION
THE USE OF DIFFERENT REMOTE SENSING TECHNIQUES FOR LANDSLIDE CHARACTERIZATION Ph.D. candidate: ANNA DE AGOSTINI Tutor: Dr. M. FLORIS Cycle: XXVI Abstract The aim of my research is to investigate the potential
More informationATMOSPHERIC ERROR, PHASE TREND AND DECORRELATION NOISE IN TERRASAR-X DIFFERENTIAL INTERFEROGRAMS
ATMOSPHERIC ERROR, PHASE TREND AND DECORRELATION NOISE IN TERRASAR-X DIFFERENTIAL INTERFEROGRAMS Steffen Knospe () () Institute of Geotechnical Engineering and Mine Surveying, Clausthal University of Technology,
More informationImproved PSI Performance for Landslide Monitoring Applications. J. Duro, R. Iglesias, P. Blanco-Sánchez, F. Sánchez and D. Albiol
Improved PSI Performance for Landslide Monitoring Applications J. Duro, R. Iglesias, P. Blanco-Sánchez, F. Sánchez and D. Albiol Outline Area of Study Previous PSI (and others) studies Main conclusions
More informationIEEE GEOSCIENCE AND REMOTE SENSING LETTERS, VOL. 1, NO. 2, APRIL
IEEE GEOSCIENCE AND REMOTE SENSING LETTERS, VOL. 1, NO. 2, APRIL 2004 57 Delta-K Interferometric SAR Technique for Snow Water Equivalent (SWE) Retrieval Geir Engen, Tore Guneriussen, and Øyvind Overrein
More information[1] Integrated Global Observing Strategy, Theme Report 2003 Monitoring of our Environment from Earth and Space.
Bibliography [1] Integrated Global Observing Strategy, Theme Report 2003 Monitoring of our Environment from Earth and Space. [2] Y. Chen, G. Zang, X. Ding, Z. Li, "Monitoring Earth Surface Deformations
More informationAPPLICABILITY OF PSINSAR FOR BUILDING HAZARD IDENTIFICATION
APPLICABILITY OF PSINSAR FOR BUILDING HAZARD IDENTIFICATION. STUDY OF THE 29 JANUARY 2006 KATOWICE EXHIBITION HALL COLLAPSE AND THE 24 FEBRUARY 2006 MOSCOW BASMANNY MARKET COLLAPSE Zbigniew Perski (1),
More informationAdvances in Geosciences
Advances in Geosciences (2003) 1: 137 142 c European Geosciences Union 2003 Advances in Geosciences From the time-wise to space-wise GOCE observables M. Reguzzoni DIIAR, Sez. Rilevamento, Politecnico di
More informationStatistical and Adaptive Signal Processing
r Statistical and Adaptive Signal Processing Spectral Estimation, Signal Modeling, Adaptive Filtering and Array Processing Dimitris G. Manolakis Massachusetts Institute of Technology Lincoln Laboratory
More informationNavneet Agrawal Deptt. Of Electronics & Communication Engineering CTAE,MPUAT,Udaipur,India
Navneet Agrawal et al / (IJCSE) International Journal on Computer Science and Engineering, Saturation Adaptive Quantizer Design for Synthetic Aperture Radar Data Compression Navneet Agrawal Deptt. Of Electronics
More informationInSAR techniques and applications for monitoring landslides and subsidence
Geoinformation for European-wide Integration, Benes (ed.) 2003 Millpress, Rotterdam, ISBN 90-77017-71-2 InSAR techniques and applications for monitoring landslides and subsidence H. Rott & T. Nagler Institut
More informationSignal Modeling Techniques in Speech Recognition. Hassan A. Kingravi
Signal Modeling Techniques in Speech Recognition Hassan A. Kingravi Outline Introduction Spectral Shaping Spectral Analysis Parameter Transforms Statistical Modeling Discussion Conclusions 1: Introduction
More informationAirborne Holographic SAR Tomography at L- and P-band
Airborne Holographic SAR Tomography at L- and P-band O. Ponce, A. Reigber and A. Moreira. Microwaves and Radar Institute (HR), German Aerospace Center (DLR). 1 Outline Introduction to 3-D SAR Holographic
More informationMonitoring long-term ground movements and Deep Seated Gravitational
Monitoring long-term ground movements and Deep Seated Gravitational Slope Deformations by InSAR time series: cases studies in Italy Salvatore Stramondo (1), M. Saroli (1, 2), M. Moro (1, 2), S. Atzori
More informationIvana Zinno, Francesco Casu, Claudio De Luca, Riccardo Lanari, Michele Manunta. CNR IREA, Napoli, Italy
An Unsupervised Implementation of the P-SBAS DiNSAR Algorithm for Processing Large Data Volumes through Distributed Computing Infrastructures within Operational Environments Ivana Zinno, Francesco Casu,
More informationLet us consider a typical Michelson interferometer, where a broadband source is used for illumination (Fig. 1a).
7.1. Low-Coherence Interferometry (LCI) Let us consider a typical Michelson interferometer, where a broadband source is used for illumination (Fig. 1a). The light is split by the beam splitter (BS) and
More informationRADAR Remote Sensing Application Examples
RADAR Remote Sensing Application Examples! All-weather capability: Microwave penetrates clouds! Construction of short-interval time series through cloud cover - crop-growth cycle! Roughness - Land cover,
More informationEasyChair Preprint. Height Accuracy and Data Coverage for the Final Global TanDEM-X DEM Data
EasyChair Preprint 123 Height Accuracy and Data Coverage for the Final Global TanDEM-X DEM Data Christopher Wecklich, Carolina Gonzalez and Paola Rizzoli EasyChair preprints are intended for rapid dissemination
More information3-Dimension Deformation Mapping from InSAR & Multiaperture. Hyung-Sup Jung The Univ. of Seoul, Korea Zhong Lu U.S. Geological Survey, U.S.A.
3-Dimension Deformation Mapping from InSAR & Multiaperture InSAR Hyung-Sup Jung The Univ. of Seoul, Korea Zhong Lu U.S. Geological Survey, U.S.A. Outline Introduction to multiple-aperture InSAR (MAI) 3-D
More informationThe financial and communal impact of a catastrophe instantiated by. volcanoes endlessly impact on lives and damage expensive infrastructure every
Chapter 1 Introduction The financial and communal impact of a catastrophe instantiated by geophysical activity is significant. Landslides, subsidence, earthquakes and volcanoes endlessly impact on lives
More informationSentinel-1A SAR Interferometry Verification
Sentinel-1A SAR Interferometry Verification Dirk Geudtner1, Pau Prats2, Nestor Yaguee-Martinez2, Andrea Monti Guarnieri3, Itziar Barat1, Björn Rommen1 and Ramón Torres1 1ESA ESTEC 2DLR, Microwave and Radar
More informationsensors ISSN
Sensors 008, 8, 546-5448; DOI: 10.3390/s809546 OPEN ACCESS sensors ISSN 144-80 www.mdpi.org/sensors Review Atmospheric Effects on InSAR Measurements and Their Mitigation Xiao-li Ding 1, *, Zhi-wei Li,
More informationAtmospheric Effects on InSAR Measurements in Southern China and Australia: A Comparative Study
Atmospheric Effects on InSAR Measurements in Southern China and Australia: A Comparative Study X. L. Ding a, *, L. Ge b, Z. W. Li a, C. Rizos b a Dept. of Land Surveying and Geo-Informatics, Hong Kong
More informationPSI Precision, accuracy and validation aspects
PSI Precision, accuracy and validation aspects Urs Wegmüller Gamma Remote Sensing AG, Gümligen, Switzerland, wegmuller@gamma-rs.ch Contents - Precision - Accuracy - Systematic errors - Atmospheric effects
More informationMERIS and OSCAR: Online Services for Correcting Atmosphere in Radar
National Aeronautics and Space Administration MERIS and OSCAR: Online Services for Correcting Atmosphere in Radar Eric Fielding and Evan Fishbein Jet Propulsion Laboratory, California Inst. of Tech. Zhenhong
More informationInterferometric Synthetic Aperture Radar (InSAR): Its Past, Present and Future
Interferometric Synthetic Aperture Radar (InSAR): Its Past, Present and Future by Zhong Lu, Ohig Kwoun, and Russell Rykhus Introduction Very simply, interferometric synthetic aperture radar (InSAR) involves
More informationBlind Identification of FIR Systems and Deconvolution of White Input Sequences
Blind Identification of FIR Systems and Deconvolution of White Input Sequences U. SOVERINI, P. CASTALDI, R. DIVERSI and R. GUIDORZI Dipartimento di Elettronica, Informatica e Sistemistica Università di
More informationACHIEVING THE ERS-2 ENVISAT INTER-SATELLITE INTERFEROMETRY TANDEM CONSTELLATION.
ACHIEVING THE ERS-2 ENVISAT INTER-SATELLITE INTERFEROMETRY TANDEM CONSTELLATION M. A. Martín Serrano (1), M. A. García Matatoros (2), M. E. Engdahl (3) (1) VCS-SciSys at ESA/ESOC, Robert-Bosch-Strasse
More informationTwo-pass DInSAR uses an interferometric image pair and an external digital elevation model (DEM). Among the two singlelook complex (SLC) images, one i
DIFFERENTIAL RADAR INTERFEROMETRY AND ITS APPLICATION IN MONITORING UNDERGROUND COAL MINING-INDUCED SUBSIDENCE Yaobin Sheng a, b, c, *, Yunjia Wang a, b, Linlin Ge c, Chris Rizos c a Jiangsu Key Laboratory
More informationAn Iterative Blind Source Separation Method for Convolutive Mixtures of Images
An Iterative Blind Source Separation Method for Convolutive Mixtures of Images Marc Castella and Jean-Christophe Pesquet Université de Marne-la-Vallée / UMR-CNRS 8049 5 bd Descartes, Champs-sur-Marne 77454
More informationStudy and Applications of POLSAR Data Time-Frequency Correlation Properties
Study and Applications of POLSAR Data Time-Frequency Correlation Properties L. Ferro-Famil 1, A. Reigber 2 and E. Pottier 1 1 University of Rennes 1, Institute of Electronics and Telecommunications of
More informationDUAL-POLARIZED COSMO SKYMED SAR DATA TO OBSERVE METALLIC TARGETS AT SEA
DUAL-POLARIZED COSMO SKYMED SAR DATA TO OBSERVE METALLIC TARGETS AT SEA F. Nunziata, M. Montuori and M. Migliaccio Università degli Studi di Napoli Parthenope Dipartimento per le Tecnologie Centro Direzionale,
More informationPERSISTENT SCATTERER INTERFEROMETRY: POTENTIAL AND LIMITS
PERSISTENT SCATTERER INTERFEROMETRY: POTENTIAL AND LIMITS M. Crosetto a, O. Monserrat a, A. Jungner, B. Crippa b a Institute of Geomatics, Av. del Canal Olímpic, s/n, Castelldefels, E-08860, Spain (michele.crosetto,
More informationSurface Deformation Measurements Scientific Requirements & Challenges
Surface Deformation Measurements Scientific Requirements & Challenges 1st Science and Application Workshop for Germany-Japan Next-Generation SAR M. Eineder, C. Minet, A. Parizzi Tokyo, 27.6.2013 Tandem-L
More informationCNES R&D and available software for Space Images based risk and disaster management
CNES R&D and available software for Space Images based risk and disaster management 1/21 Contributors: CNES (Centre National d Etudes Spatiales), Toulouse, France Hélène Vadon Jordi Inglada 2/21 Content
More informationCFAR TARGET DETECTION IN TREE SCATTERING INTERFERENCE
CFAR TARGET DETECTION IN TREE SCATTERING INTERFERENCE Anshul Sharma and Randolph L. Moses Department of Electrical Engineering, The Ohio State University, Columbus, OH 43210 ABSTRACT We have developed
More informationEfficient Semi-Blind Channel Estimation and Equalization Based on a Parametric Channel Representation
Efficient Semi-Blind Channel Estimation and Equalization Based on a Parametric Channel Representation Presenter: Kostas Berberidis University of Patras Computer Engineering & Informatics Department Signal
More information3D temporal evolution of displacements recorded on Mt. Etna from the 2007 to 2010 through the SISTEM method
3D temporal evolution of displacements recorded on Mt. Etna from the 2007 to 2010 through the SISTEM method Bonforte A., Guglielmino F.,, Puglisi G. INGV Istituto Nazionale di Gofisica e vulcanologia Osservatorio
More informationApplication of PSI technique to slope stability monitoring in the Daunia mountains, Italy
ESA ESRIN - Frascati 28 November- 2 December Application of PSI technique to slope stability monitoring in the Daunia mountains, Italy F. Bovenga (1) (fabio.bovenga@ba.infn.it) L. Guerriero (1) R. Nutricato
More informationAdaptive beamforming for uniform linear arrays with unknown mutual coupling. IEEE Antennas and Wireless Propagation Letters.
Title Adaptive beamforming for uniform linear arrays with unknown mutual coupling Author(s) Liao, B; Chan, SC Citation IEEE Antennas And Wireless Propagation Letters, 2012, v. 11, p. 464-467 Issued Date
More informationMOST synthetic aperture radar (SAR) satellites operate in
IEEE GEOSCIENCE AND REMOTE SENSING LETTERS, VOL. 9, NO. 2, MARCH 2012 257 MERIS Atmospheric Water Vapor Correction Model for Wide Swath Interferometric Synthetic Aperture Radar Zhenhong Li, Member, IEEE,
More information