Fringe 2017 > Session details
Paper 11 - Session title: Thematic mapping, vegetation and DEMs
09:40 The Global TanDEM-X Digital Elevation Model: Final Performance Assessment
Rizzoli, Paola; Martone, Michele; Gonzales, Carolina; Wecklich, Christopher; Borla Tridon, Daniela; Bachmann, Markus; Fritz, Thomas; Wessel, Birgit; Krieger, Gerhard; Zink, Manfred German Aerospace Center, Germany
Digital elevation models are of fundamental importance for a large variety of scientific and commercial applications. For example, precise and up-to date information about the Earth’s topography is required in many geoscience areas, such as geology, forestry, glaciology, oceanography, and hydrology. For reliable navigation using global positioning systems, such as GPS and Galileo, digital elevation maps are of essential importance as well. Up to now, the primary source of elevation data on an almost global scale has been provided by the Shuttle Radar Topography Mission (SRTM), which is characterized by a spatial sampling of 30 m between [-60°,+60°] latitudes. For higher latitudes and over Antarctica, only low resolution DEMs are available, such as GTOPO, RAMP, and GLAS/ICESat.
With the main goal of acquiring a global and consistent DEM with unprecedented accuracy, the TanDEM-X mission (TerraSARX add-on for Digital Elevation Measurements) opened a new era in spaceborne synthetic aperture radar (SAR). Developed in a public-private partnership between the German Aerospace Center (DLR) and Airbus Defence and Space, it is comprised of two almost identical satellites, TerraSAR-X and TanDEM-X, mounting a synthetic aperture radar antenna operating at X-band. Since October 2010, both satellites have been flying in a close orbit configuration at an altitude of around 500 kilometers, acting as a single-pass SAR interferometer and allowing for a flexible selection of baselines and acquisition geometries. Images are nominally
acquired in bistatic configuration, where one satellite transmits and both simultaneously receive the backscattered signal from the Earth’s surface, enabling the acquisition of highly accurate interferograms, which do not suffer from temporal and atmospheric decorrelation.
A dedicated acquisition strategy has been developed and optimized throughout the years, in order to achieve a homogenous performance globally. The satellite formation has been adjusted accordingly during different acquisition phases, depending on the target height performance to be achieved.
This paper reports the final performance of the TanDEM-X global DEM, presenting the developed acquisition and processing strategy, and a detailed analysis of all characteristic parameters. The global performance is then assessed in terms of vertical and horizontal errors, and coverage statistics. The obtained results confirm the outstanding quality of the delivered product, which can be now globally exploited for both scientific and commercial applications.
Paper 339 - Session title: Thematic mapping, vegetation and DEMs
09:00 High Precision DSM Generation in Densely Vegetated Mountainous Areas with Dual-Baseline InSAR Assisted by StereoSAR
Dong, Yuting (1); Jiang, Houjun (2); Zhang, Lu (1,3); Liao, Mingsheng (1,3) 1: State Key Laboratory of Information Engineering in Surveying, Mapping and Remote Sensing, Wuhan University; 2: Nanjing University of Posts and Telecommunications; 3: Collaborative Innovation Center for Geospatial Technology
Synthetic Aperture Radar Interferometry (InSAR) is a powerful tool for large-area topographic mapping due to its capability of all-time all-weather imaging and high sensitivity to terrain relief. However, there is an inherent contradiction between geometric decorrelation and sensitivity of height measurement for topographic mapping with a single InSAR pair. A normal baseline of proper length is required to keep a balance between the two issues. A promising solution to this problem is the so-called multi-baseline InSAR analysis. The basic principle of multi-baseline InSAR is to derive an optimal height estimate by joint analysis of multiple phase measurements from a few interferograms with different normal baselines. Compared with single-baseline InSAR, the major benefit of using multi-baseline observations is the possibility of exploiting redundant topographic phase observations with different height of ambiguities to improve the accuracy of phase unwrapping, or even avoid phase unwrapping.
Stereo radargrammetry (StereoSAR) is also widely used to reconstruct digital surface models (DSMs) Compared with InSAR, it is less affected by temporal decorrelation. Therefore, StereoSAR is often used for DSM generation in heavily vegetated areas but with lower accuracy. As two major approaches to use SAR remote sensing data for topographic mapping, both StereoSAR and InSAR techniques have their own advantages and disadvantages. Consequently, it is not easy to obtain accurate and reliable height information in terrain relief areas with vegetation cover by adopting only either technique. Aiming at such a problem, we are going to carry out a study on joint utilization of StereoSAR and dual/multi-baseline InSAR to make full use of advantages of the two techniques to substantially improve our capability of height measurement.
In this study, we developed a maximum a posteriori (MAP) estimation method for multi-baseline InSAR assisted by StereoSAR. According to Bayesian theory, the combination of StereoSAR and InSAR for topographic mapping can be viewed as update of the StereoSAR DSM with InSAR phase observations. At the same time, StereoSAR DSM is also a constraint to InSAR phase observations, which can solve the problem of elevation ambiguity and avoid phase unwrapping process. Under the condition of distributed scattering, the likelihood function of height could be derived from InSAR phase. The prior probability of height is acquired form StereoSAR DSM.
In order to verify the proposed method, we choose Mount Song as the test area, which is one of the five sacred mountains of China. Although the mountain peaks reach only around 1500 m, the slopes are very steep and they are densely vegetated, making Mount Song an difficult area for InSAR DSM generation. In our study two interferometric data pairs acquired by TerraSAR-X with different normal baselines and one stereo pair of TerraSAR-X data in Stripmap mode are used together to reconstruct the DSM of 10 m spatial resolution. In order to evaluate the accuracy of generated DSM, we use a 1 m resolution DSM created by airborne photogrammetry as reference data. The experimental result shows that there is neither systematic error nor large data voids in MAP estimated DSM and the standard deviation of height error σh of MAP estimated DSM is less than 10 m with respect to the photogrammetric DSM for the whole area, while in plain areas σh is about 5 m.
Paper 448 - Session title: Thematic mapping, vegetation and DEMs
10:00 On the Use of TanDEM-X Bistatic InSAR Images for Scene Recognition
Cagatay, Nazli Deniz; Datcu, Mihai German Aerospace Center, Germany
Research on SAR interferometry is mainly focused on the two main application areas, namely the generation of accurate digital elevation models and change detection, and also on the crucial processing steps like interferogram filtering, phase unwrapping, co-registration etc. This study rather aims to make use of interferometric SAR (InSAR) images, once they are generated, for scene recognition purposes.
In literature, limited research is available on the use of InSAR images for object recognition and scene classification. For those studies, the main trend is in the direction of using the interferometric coherence for mostly binary classification such as forest/non-forest, urban/non-urban or change/no-change classification. On the other hand, available research on multi-class classification is mostly based on the temporal variation of interferometric coherence and/or backscatter intensity.
However, in our studies, we emphasize the use of whole complex-valued InSAR image for multi-class classification, i.e., recognition of more complex scenes such as forest, agricultural fields, water body, different kinds of residential and industrial areas, etc., and also their combinations. Furthermore, a new complex-valued phase-gradient InSAR (PGInSAR) image is defined whose phase represents the magnitude of the phase gradient of InSAR in range and azimuth directions. Interferometric phase being related to the terrain height, phase of PGInSAR image can be considered as a measure of the terrain slope, i.e., how fast the interferometric phase changes over the image.
This review study serves as a comparative assessment of various feature descriptors such as Gabor-based, FrFT-based, BoVW-based and partial derivatives based features extracted from the complex-valued InSAR and PGInSAR images, and used for patch-based classification. For this purpose, an image patch database is generated from bistatic interferometric pairs acquired by the TanDEM-X mission over the test site Toulouse, France. In order to investigate the impact of effective baseline on the classification, 3 datasets are constructed from the bistatic acquisitions with 3 different effective baselines over the same area. Supervised KNN classification is performed on the image patches of 200 x 200 pixels from 8 different scene classes representing different natural and man-made structures on a 400m x 340m terrain (Figure 1).
The classification results can be summarized as follows:
The features extracted from the complex-valued InSAR and PGInSAR images improve the mean accuracy by 15% and 27% compared to detected or complex-valued SAR images, respectively .
The use of InSAR and PGInSAR images improves the individual class accuracies for scenes with natural structures, such as agricultural area, forest and mixed vegetation .
PGInSAR images are found to be more robust to effective baseline changes than InSAR images, as the phase discontinuities are reduced compared to InSAR images .
Although the global Gabor features yield better accuracies than the global FrFT features for InSAR images, by the implementation of BoVW model, FrFT outperforms Gabor features .
The spatial discriminability of BoVW model together with the strong interferometric signature of urban structures for larger baseline results in better classification performance .
The use of partial derivatives based scale-space representation improves the classification performance for scenes like industrial and urban areas. The presence of many local variations in such scenes is represented quite successfully by these features .
Our previous work:
 N. D. Cagatay and M. Datcu, “Scene Recognition Based on Phase Gradient InSAR Images,” IEEE International Conference on Image Processing (ICIP), Paris, France, 27-30 October 2014.
 N. D. Cagatay and M. Datcu, “FrFT Based Scene Classification of Phase Gradient InSAR Images and Effective Baseline Dependence,” IEEE Geoscience and Remote Sensing Letters, vol. 12, no. 5, pp. 1131-1135, May 2015.
 N. D. Cagatay and M. Datcu, “Bag-of-Visual-Words Model for Classification of Interferometric SAR Images,” European Conference on Synthetic Aperture Radar (EUSAR), Hamburg, Germany, 6-9 June 2016.
 N. D. Cagatay and M. Datcu, “Multi-Scale Feature Extraction Approaches for Classification of InSAR and Phase Gradient InSAR Images,” IEEE International Conference on Image Processing (ICIP), Phoenix, Arizona, USA, 25-28 September 2016.
Paper 505 - Session title: Thematic mapping, vegetation and DEMs
09:20 TanDEM-X and Cartosat-1 DEM: Height accuracy evaluation and weighted–average fusion for DEM quality improvement
Bagheri, Hossein (1); Schmitt, Michael (1); Zhu, Xiao Xiang (1,2) 1: Signal Processing in Earth Observation, Technical University of Munich (TUM); 2: Remote Sensing Technology Institute (IMF), German Aerospace Center (DLR)
Digital Elevation Models (DEMs) as a representation of terrain surface with different resolution and coverage are produced for diverse applications. Special interest lies in the generation of global DEMs, which provide homogeneous information about the world’s topography for (almost) the complete globe. One of the most recent and most prominent examples of global terrain datasets is the TanDEM-X DEM, which was produced by the German TanDEM-X mission, consisting of two SAR satellites flying in close orbit configuration. The mission goals aimed at a DEM of HRTI‑3 standard, i.e. with a relative height accuracy better than 2 m for slopes lower than 20 percent. However, initial visual inspection of TDX DEM data demonstrates undesirable horizontal resolution and height accuracy in morphologically complex terrains like urban areas and in consequence the need for quality enhancement in these areas. One possibility for improvement of TDX elevation data in problematic terrain is the fusion with other available elevation data which do not suffer from SAR-inherent imaging effects such as layover and shadowing. Examples for these alternative data are DEMs derived from optical stereo imagery. Optical DEMs generally reach HRTI-3 standard (with a relative height accuracy of 2-3 m) when they are produced using high-resolution optical imagery like acquired by Cartosat-1.
The objective of this paper is two-fold: First, the height accuracies of TDX and Cartosat-1 elevation data over different land types is empirically evaluated in order to analyze the potential of TDX-Cartosat-1 DEM data fusion. After the quality assessment, a preliminary DEM fusion is carried out over urban areas aiming at TDX DEM enhancement.
For height accuracy evaluation, six subsets of different land types are selected:
A) Industrial area. This land type is characterized by flat, simple and rectangular buildings which often are not very high.
B) Inner city area. This land type includes many dense buildings and possibly also some higher buildings.
C) High buildings are skyscrapers and other high rise buildings.
D) Residential area. This type contains many small one-family homes.
E) Forested and agricultural area.
F) Agricultural areas and lake.
The Cartosat-1 DEM used in this study has been produced by the XDibias image processing system of (DLR) and has a nominal resolution of 5 m. The TDX raw DEM used in this study was produced by DLR’s Integrated TanDEM-X Processor (ITP) from data with an incidence angle of about 38° and a height of ambiguity of about 46 m and has a pixel size of 0.2 arc seconds. All these test areas are located in the area of Munich, Bavaria so that for each of them a highly-accurate LiDAR point cloud (with a density of 1 point per m²) provided by the Bavarian surveying administration is available as a reference.
The uncertainty assessment confirmed that the Cartosat-1 DEM is more precise than the TDX DEM over urban areas, especially when high buildings are present while in non-urban subsets both DEMs have almost identical accuracy. On the other hand, Cartosat-1 DEM suffers from misalignment and TDX DEM with appropriate absolute accuracy can be employed as a great reference for absolute orientation, bias compensation, and RPC refinement.
As explained earlier, a feasible solution for improvement of TDX elevation data over urban areas is a fusion with Cartosat-1 DEM data, which provides much better accuracy in urban areas. The simplest method for DEM fusion is weighted averaging (WA). The main issue of WA is to employ weights which are proportional to the expected height errors. In the following, the weight maps are generated based on the effective factors and features such as height error map of TDX DEM and matching standard deviation derived from optical imagery in the process of Cartosat-1 DEM generation. The consequences of fusion confirm the improvement of TDX in all subsets.
Paper 530 - Session title: Thematic mapping, vegetation and DEMs
10:20 Exploitation of Sentinel-1 Interferometric Coherence for Land Cover and Vegetation Mapping (SInCohMap project)
Vicente-Guijalba, F. (1); Duro, J. (1); Lopez-Martinez, C. (2); Lopez-Sanchez, J. M. (3); Notarnicola, C. (4); Jacob, A. (4); Sonnenschein, R. (4); Dabrowska-Zielinska, K. (5); Hoscilo, A. (5); Pottier, E. (6); Lavalle, M. (7); Engdahl, M. (8) 1: Dares Technology, Spain; 2: UPC, Spain; 3: UA, Spain; 4: EURAC, Italy; 5: IGIK, Poland; 6: University of Rennes 1, France; 7: JPL, USA; 8: ESA-ESRIN, Italy
The Sentinel-1 mission represents a challenging opportunity to develop operational applications for classification and vegetation mapping using periodical interferometric SAR acquisitions. The SinCohMap project under the ESA SEOM framework is intended to develop novel methodologies exploiting the evolution of the interferometric products. To compare with state-of-the-art classification strategies, within the project context a round-robin consultation is defined to involve both internal and external expert groups to test their algorithms using a common set of well-known and in detail monitored scenarios.
The main objective of this research is to develop, analyze and validate novel methodologies for land cover and vegetation mapping using time series of Sentinel-1 data and by exploiting the temporal evolution of the interferometric coherence. Further the project aims on quantifying the impact and possible benefit of using Sentinel-1 InSAR (Interferometric Synthetic Aperture Radar) data relative to traditional land cover and vegetation mapping using optical data (especially Sentinel-2) and traditional intensity-based SAR (Synthetic Aperture Radar) approaches. A careful revision of state-of-the-art classifiers along with the definition of novel methodologies specifically designed for this scenario is being addressed within the project framework.
The most innovative concept is the creation of a robust library based on a forward physical model, which can help identify or separate temporal decorrelation due to physical changes of the scatterers. In such a way, the models can be derived from expected physical changes towards time (e.g. vegetation growth, crops evolution) to identify the classes and to separate from spike measurements or outliers. Machine learning algorithms are also being analyzed due to their high performance and their ability to deal with complex data such as high-dimensional imagery and multi-source data sets.
The main classes sought after are Forests, Agricultural areas (e.g. Crops), Artificial surfaces (e.g. Urban), Water Bodies, Scrub and Herbaceous Vegetation, Open or bare land with little to no vegetation and Wetlands. Three different reference test site areas are defined with very accurate ground truth data are defined for performing quantitative assessment and validation in Spain, Italy and Poland. Additional evaluations to detect the ability to extrapolate the extracted ideas will be performed over boreal forest.