performance analysis of the tandem-x dem generated · pdf fileperformance analysis of the...

21
Performance analysis of the TanDEM-X DEM generated with DIAPASON XTI_LAND1554 Oberpfaffenhofen 12/06/2013 Javier Duro 1 , Fifamè. N. Koudogbo 1 , Palma Blonda 2 , Richard Lucas 3 , Sander Mucher 4 , Patrick Ordoqui 1 , Javier García Robles 1 TanDEM-X science team meeting 10-14 June 2013 DLR Oberpfaffenhofen, Germany 1. Altamira Information 2. CNR Italy 3. Aberystwyth University 4. Alterra Netherlands

Upload: doanduong

Post on 07-Feb-2018

224 views

Category:

Documents


1 download

TRANSCRIPT

Page 1: Performance analysis of the TanDEM-X DEM generated  · PDF filePerformance analysis of the TanDEM-X DEM generated with DIAPASON ... Blonda2, Richard Lucas3, Sander Mucher4,

Performance analysis of the TanDEM-X DEM generated with DIAPASON

XTI_LAND1554Oberpfaffenhofen 12/06/2013

Javier Duro1, Fifamè. N. Koudogbo1, Palma Blonda2, Richard Lucas3, Sander Mucher4,

Patrick Ordoqui1, Javier García Robles1

TanDEM-X science team meeting10-14 June 2013

DLR Oberpfaffenhofen, Germany

1. Altamira Information2. CNR Italy3. Aberystwyth University4. Alterra Netherlands

Page 2: Performance analysis of the TanDEM-X DEM generated  · PDF filePerformance analysis of the TanDEM-X DEM generated with DIAPASON ... Blonda2, Richard Lucas3, Sander Mucher4,

BIO_SOS - FP7-SPACE-2010-1 project

DIAPASON adaptation to TanDEM-X data

TanDEM–X for HR DEM generation

Veluwe (The Netherlands) HR DSM generation

Fors Cohno (Wales) HR DSM generation

Conclusion

Agenda

Page 3: Performance analysis of the TanDEM-X DEM generated  · PDF filePerformance analysis of the TanDEM-X DEM generated with DIAPASON ... Blonda2, Richard Lucas3, Sander Mucher4,

BIO_SOS aims to monitor the induced human effects on Biodiversity within and around Natura 2000 sites by detecting VHR Land Cover/Use (LCLU) changes.

The main objective is the development of a pre-operational multimodular ecological modeling system (EODHaM) based on the integration of EO and in-situ data.Outputs: VHR LCLU and habitat maps. Biodiversity Indicators

Vegetation height is used for habitat classification and habitat management monitoring.

[email protected]

BIO_SOS GA 263435, FP7-SPACE-2010-1 Project BIOdiversity Multi-SOurce Monitoring System: From Space To Species

www.biosos.eu

Page 4: Performance analysis of the TanDEM-X DEM generated  · PDF filePerformance analysis of the TanDEM-X DEM generated with DIAPASON ... Blonda2, Richard Lucas3, Sander Mucher4,

Evaluation of capacity of TanDEM-X to measure vegetation height: Digital Surface Model (DSM) over Canopy.

Evaluation of the TanDEM-X coherence for supporting Landcoverclassification procedures.

Adapt DIAPASON InSAR processor to the particularities of TanDEM-X:• Bi-static acquisition• CoSSCs products

On going project by now.

Veluwe (The Netherlands) DEM Generation

Cors Fochno (Wales) DEM Generation

BIO_SOS FP7-SPACE-2010-1 ProjectContribution of TanDEM-X data to the Project

Page 5: Performance analysis of the TanDEM-X DEM generated  · PDF filePerformance analysis of the TanDEM-X DEM generated with DIAPASON ... Blonda2, Richard Lucas3, Sander Mucher4,

Data extraction: Annotation QCs check, orbit state vectors, image timing annotation, Id Master scene

RAW interferogram generation: Hermitian product master-slave

Reference interferogram generation (with SRTM at 3” if available, if not WGS84)

Phase Multilooking, definition of GCPs & definition of ground projected geometry

Phase unwrapping

Phase to Height conversion with geocoding

Phase UnwrappingPhase Unwrapping

Height ConversionHeight Conversion

Output DEM

CoSSC images

MLML

RAW Interferogram RAW Interferogram

Data ExtractionData Extraction

Reference Interferogram

Reference Interferogram

DIAPASON adaptation to TanDEM-X dataBi-static InSAR processing adaptation

GCPGCP Proj.Proj.

ML interf.

Page 6: Performance analysis of the TanDEM-X DEM generated  · PDF filePerformance analysis of the TanDEM-X DEM generated with DIAPASON ... Blonda2, Richard Lucas3, Sander Mucher4,

Simulation of the SAR synthetic amplitude image due to the ground topography (SRTM at 3”) and comparison with the product amplitude (ml55)

DIAPASON adaptation to TanDEM-X dataImage timing annotation QCs

Page 7: Performance analysis of the TanDEM-X DEM generated  · PDF filePerformance analysis of the TanDEM-X DEM generated with DIAPASON ... Blonda2, Richard Lucas3, Sander Mucher4,

DIAPASON adaptation to TanDEM-X dataImage timing annotation QCs

Simulation of the SAR synthetic amplitude image due to the ground topography (SRTM at 3”) and comparison with the product amplitude (ml55)

Offsets of 7 samples and 15 samples in range and azimuth

Page 8: Performance analysis of the TanDEM-X DEM generated  · PDF filePerformance analysis of the TanDEM-X DEM generated with DIAPASON ... Blonda2, Richard Lucas3, Sander Mucher4,

After the timing correction the differential interferogram present less residual fringes

DIAPASON adaptation to TanDEM-X dataImage timing annotation Correction

Before Correction After Correction

Before Correction

Page 9: Performance analysis of the TanDEM-X DEM generated  · PDF filePerformance analysis of the TanDEM-X DEM generated with DIAPASON ... Blonda2, Richard Lucas3, Sander Mucher4,

After the timing correction the differential interferogram present less residual fringes

DIAPASON adaptation to TanDEM-X dataImage timing annotation Correction

After Correcton

Page 10: Performance analysis of the TanDEM-X DEM generated  · PDF filePerformance analysis of the TanDEM-X DEM generated with DIAPASON ... Blonda2, Richard Lucas3, Sander Mucher4,

S2

S1

R0

R0 + Rref + RDifPhase

zref

b

H

Solving an equitation system with 3 unknowns

DIAPASON adaptation to TanDEM-X dataPhase to Height conversion with geocoding

Calibration of the slant-range InSAR phase based on mean offset retrieval (GCPs)

Direct transformation of the unwrapped interferometric phase to geocoded heights for each SAR pixel

Initial guess P0 provided based on SRTM to speed up the solution

Output geometry is set according to the product needs (& ML of the processing)

Zref+measured

| 1 |2 12

12

12 02 

| 2 |2 22 2

22

2 02

2

 

2∙ 1

| 1| 

Page 11: Performance analysis of the TanDEM-X DEM generated  · PDF filePerformance analysis of the TanDEM-X DEM generated with DIAPASON ... Blonda2, Richard Lucas3, Sander Mucher4,

DIAPASON adaptation to TanDEM-X dataGlobal comparison TDX with LiDAR DSM

LiDAR – TDX DSM Global Image:Mean = -0.35 mStddev = 7.5 m

LiDAR DSM 15m. res.

TDX DSM 15m. res.

Page 12: Performance analysis of the TanDEM-X DEM generated  · PDF filePerformance analysis of the TanDEM-X DEM generated with DIAPASON ... Blonda2, Richard Lucas3, Sander Mucher4,

DIAPASON adaptation to TanDEM-X dataLiDAR and TDX DSMs over fat areas

LiDAR DSM 15m. res.

TDX DSM 15m. res.

Page 13: Performance analysis of the TanDEM-X DEM generated  · PDF filePerformance analysis of the TanDEM-X DEM generated with DIAPASON ... Blonda2, Richard Lucas3, Sander Mucher4,

DIAPASON adaptation to TanDEM-X dataLiDAR and TDX DSMs over fat areas

LiDAR – TDX DSM Flat areas:

Mean = -0.66 mStddev = 0.96 m

TDX provides more details of the groundtopography

LiDAR DSM 15m. res.

TDX DSM 15m. res.

Page 14: Performance analysis of the TanDEM-X DEM generated  · PDF filePerformance analysis of the TanDEM-X DEM generated with DIAPASON ... Blonda2, Richard Lucas3, Sander Mucher4,

Dates (November/December 2012) with normally low to moderate vegetation (post-flush).

On-going processing of those datasets ─ Preliminary DEMs have been generated for

Veluwe.─ Comparison with ground truth data (GCPs,

orthophoto, LiDAR DEM for Cors Fochno site) + comparison with SRTM and ASTER DEMs

New acquisitions have been obtained in beginning of June 2013 (peak flush).

TanDEM-X for HR DEM generationStatus of the work

Page 15: Performance analysis of the TanDEM-X DEM generated  · PDF filePerformance analysis of the TanDEM-X DEM generated with DIAPASON ... Blonda2, Richard Lucas3, Sander Mucher4,

• Clear forested zones in mountain area

• River and water reservoir

• LiDAR DSM at 15m

• Stack Worldview-2 for classification

TanDEM-X for HR DEM generationFors Cohno (Wales) area & data

Page 16: Performance analysis of the TanDEM-X DEM generated  · PDF filePerformance analysis of the TanDEM-X DEM generated with DIAPASON ... Blonda2, Richard Lucas3, Sander Mucher4,

Preliminary Analysis of the bi-static InSAR with landcover Acq date: 2012/12/25 hamb=37.2m

Low coherence areas are tied with forested areas and water surfaces

TanDEM-X for HR DEM generationFors Cohno (Wales) TDX data

Page 17: Performance analysis of the TanDEM-X DEM generated  · PDF filePerformance analysis of the TanDEM-X DEM generated with DIAPASON ... Blonda2, Richard Lucas3, Sander Mucher4,

Exploitation of the phase for trees height measurement

• Phase profiles captured crossing the trees induced discontinuity

• Relative phase measurements in the discontinuity used to retrieve possible trees heights from (7 to 14 m this case)

0

1

2

3

4

5

6

7 Wrapped phase section

discontinuity > π

TanDEM-X for HR DEM generationFors Cohno (Wales) InSAR phase over trees

Page 18: Performance analysis of the TanDEM-X DEM generated  · PDF filePerformance analysis of the TanDEM-X DEM generated with DIAPASON ... Blonda2, Richard Lucas3, Sander Mucher4,

TanDEM-X for HR DEM generationFors Cohno (Wales) HR DSM generation

Exploitation of the TDX DSM for trees height measurement

Page 19: Performance analysis of the TanDEM-X DEM generated  · PDF filePerformance analysis of the TanDEM-X DEM generated with DIAPASON ... Blonda2, Richard Lucas3, Sander Mucher4,

Exploitation of the TDX DSM for trees height measurement

TanDEM-X for HR DEM generationFors Cohno (Wales) HR DSM generation

Phase unwrapping cannot resolve the phase ambiguity in forest areas since there is no transition

LiDAR –TDX DSM differences

-18m

+18m

Page 20: Performance analysis of the TanDEM-X DEM generated  · PDF filePerformance analysis of the TanDEM-X DEM generated with DIAPASON ... Blonda2, Richard Lucas3, Sander Mucher4,

• Produced TSX DSM have a very high level of details of the ground topographic surface

• Trees volumetric backscattering seems to reduce scene interferometric coherence

• Coherence might be used combined with height information as land cover classifier given its sensitivity to forested and water areas

• Phase unwrapping present difficulties to retrieve trees height Trees height sensitivity is constrained by height of ambiguity

• Interferometric phase might be used to retrieve trees height with DInSAR interferograms at different seasons of the year

• Future and sustainability of BIO_SOS project: Evaluate the costs and feasibility of the elaboration of a TD-X HR DSM vs. LIDAR DSM

Preliminary Analysis and Key Facts

Page 21: Performance analysis of the TanDEM-X DEM generated  · PDF filePerformance analysis of the TanDEM-X DEM generated with DIAPASON ... Blonda2, Richard Lucas3, Sander Mucher4,

ALTAMIRA INFORMATIONC/ Còrsega, 381-387E-08037 Barcelona

Tel.: +34 93 183 57 50Fax: +34 93 183 57 59

{javier.duro, fifame.koudogbo}@altamira-information.com

www.biosos.wur.nl