Digital Canopy Model Estimation from TanDEM-X Interferometry using High-Resolution Lidar DEM

Interferometric TanDEM-X data are used together with high-resolution, airborne lidar-derived digital elevation models (DEMs) to produce digital canopy models (DCMs) for the boreal forests of Remningstorp and Krycklan, situated in southern and northern Sweden, respectively. An overview of interferometric data processing is given. First results showing the potential of TanDEM-X-based forest canopy mapping are presented. It is concluded that baselines giving height-of-ambiguity values in the order of 50-80 meters are preferable, although factors such as angle of incidence and along-track baseline are also of importance. Clear-cuts can easily be detected thanks to the high resolution of TanDEM-X imagery. Seasonal variations of scattering height are most visible for deciduous trees, where the scattering height is significantly lower in the winter, probably due to the lack of leaves. \ua9 2013 IEEE

Utilization of bistatic TanDEM-X data to derive land cover information

Forests have significance as carbon sink in climate change. Therefore, it is of high importance to track land use changes as well as to estimate the state as carbon sink. This is useful for sustainable forest management, land use planning, carbon modelling, and support to implement international initiatives like REDD+ (Reducing Emissions from Deforestation and Degradation). A combination of field measurements and remote sensing seems most suitable to monitor forests. Radar sensors are considered as high potential due to the weather and daytime independence. TanDEM-X is a interferometric SAR (synthetic aperture radar) mission in space and can be used for land use monitoring as well as estimation of biophysical parameters. TanDEM-X is a X-band system resulting in low penetration depth into the forest canopy. Interferometric information can be useful, whereas the low penetration can be considered as an advantage. The interferometric height is assumable as canopy height, which is correlated with forest biomass. Furthermore, the interferometric coherence is mainly governed by volume decorrelation, whereas temporal decorrelation is minimized. This information can be valuable for quantitative estimations and land use monitoring. The interferometric coherence improved results in comparison to land use classifications without coherence of about 10% (75% vs. 85%). Especially the differentiation between forest classes profited from coherence. The coherence correlated with aboveground biomass in a R² of about 0.5 and resulted in a root mean square error (RSME) of 14%. The interferometric height achieved an even higher correlation with the biomass (R²=0.68) resulting in cross-validated RMSE of 7.5%. These results indicated that TanDEM-X can be considered as valuable and consistent data source for forest monitoring. Especially interferometric information seemed suitable for biomass estimation

Estimation of Forest Biomass From Two-Level Model Inversion of Single-Pass InSAR Data

A model for aboveground biomass estimation from single-pass interferometric synthetic aperture radar (InSAR) data is presented. Forest height and canopy density estimates Delta h and eta(0), respectively, obtained from two-level model (TLM) inversion, are used as biomass predictors. Eighteen bistatic VV-polarized TanDEM-X (TDM) acquisitions are used, made over two Swedish test sites in the summers of 2011, 2012, and 2013 (nominal incidence angle: 41 degrees; height-of-ambiguity: 32-63 m). Remningstorp features a hemiboreal forest in southern Sweden, with flat topography and where 32 circular plots have been sampled between 2010 and 2011 (area: 0.5 ha; biomass: 42-242 t/ha; height: 14-32 m). Krycklan features a boreal forest in northern Sweden, 720-km north-northeast from Remningstorp, with significant topography and where 31 stands have been sampled in 2008 (area: 2.4-26.3 ha; biomass: 23-183 t/ha; height: 7-21 m). A high-resolution digital terrain model has been used as ground reference during InSAR processing. For the aforementioned plots and stands and if the same acquisition is used for model training and validation, the new model explains 65%-89% of the observed variance, with root-mean-square error (RMSE) of 12%-19% (median: 15%). By fixing two of the three model parameters, accurate biomass estimation can also be done when different acquisitions or different test sites are used for model training and validation, with RMSE of 12%-56% (median: 17%). Compared with a simple scaling model computing biomass from the phase center elevation above ground, the proposed model shows significantly better performance in Remningstorp, as it accounts for the large canopy density variations caused by active management. In Krycklan, the two models show similar performance

Rice Plant Height Monitoring from Space with Bistatic Interferometry

This chapter provides an overview of the possibility to derive paddy rice plant heights with spaceborne bistatic SAR interferometry (InSAR). By using the only available interferometer in space, TanDEM-X, an investigation of rice crops located in Turkey is performed. Before analyzing the main outcomes, an introduction to the generation of elevation models with InSAR is provided, with a special focus on the agricultural land cover. The processing chain and the modifications foreseen to properly produce plant elevations and a roadmap for the quality assessment are described. The results obtained, with a very high interferometric coherence supporting an accurate estimation due to a limited electromagnetic wave penetration into the canopy, support a temporal change analysis on a field-by-field basis. For the purpose, an automatic approach to segment the fields without external auxiliary data is also provided. The study is concluded with an analysis of the impact of the wave polarization in the results

Vertical Accuracy of Freely Available Global Digital Elevation Models (ASTER, AW3D30, MERIT, TanDEM-X, SRTM, and NASADEM)

https://www.mdpi.com/2072-4292/12/21/348

Biomass Representation in Synthetic Aperture Radar Interferometry Data Sets

This work makes an attempt to explain the origin, features and potential applications of the elevation bias of the synthetic aperture radar interferometry (InSAR) datasets over areas covered by vegetation. The rapid development of radar-based remote sensing methods, such as synthetic aperture radar (SAR) and InSAR, has provided an alternative to the photogrammetry and LiDAR for determining the third dimension of topographic surfaces. The InSAR method has proved to be so effective and productive that it allowed, within eleven days of the space shuttle mission, for acquisition of data to develop a three-dimensional model of almost the entire land surface of our planet. This mission is known as the Shuttle Radar Topography Mission (SRTM). Scientists across the geosciences were able to access the great benefits of uniformity, high resolution and the most precise digital elevation model (DEM) of the Earth like never before for their a wide variety of scientific and practical inquiries. Unfortunately, InSAR elevations misrepresent the surface of the Earth in places where there is substantial vegetation cover. This is a systematic error of unknown, yet limited (by the vertical extension of vegetation) magnitude. Up to now, only a limited number of attempts to model this error source have been made. However, none offer a robust remedy, but rather partial or case-based solutions. More work in this area of research is needed as the number of airborne and space-based InSAR elevation models has been steadily increasing over the last few years, despite strong competition from LiDAR and optical methods. From another perspective, however, this elevation bias, termed here as the “biomass impenetrability”, creates a great opportunity to learn about the biomass. This may be achieved due to the fact that the impenetrability can be considered a collective response to a few factors originating in 3D space that encompass the outermost boundaries of vegetation. The biomass, presence in InSAR datasets or simply the biomass impenetrability, is the focus of this research. The report, presented in a sequence of sections, gradually introduces terminology, physical and mathematical fundamentals commonly used in describing the propagation of electromagnetic waves, including the Maxwell equations. The synthetic aperture radar (SAR) and InSAR as active remote sensing methods are summarised. In subsequent steps, the major InSAR data sources and data acquisition systems, past and present, are outlined. Various examples of the InSAR datasets, including the SRTM C- and X-band elevation products and INTERMAP Inc. IFSAR digital terrain/surface models (DTM/DSM), representing diverse test sites in the world are used to demonstrate the presence and/or magnitude of the biomass impenetrability in the context of different types of vegetation – usually forest. Also, results of investigations carried out by selected researchers on the elevation bias in InSAR datasets and their attempts at mathematical modelling are reviewed. In recent years, a few researchers have suggested that the magnitude of the biomass impenetrability is linked to gaps in the vegetation cover. Based on these hints, a mathematical model of the tree and the forest has been developed. Three types of gaps were identified; gaps in the landscape-scale forest areas (Type 1), e.g. forest fire scares and logging areas; a gap between three trees forming a triangle (Type 2), e.g. depending on the shape of tree crowns; and gaps within a tree itself (Type 3). Experiments have demonstrated that Type 1 gaps follow the power-law density distribution function. One of the most useful features of the power-law distributed phenomena is their scale-independent property. This property was also used to model Type 3 gaps (within the tree crown) by assuming that these gaps follow the same distribution as the Type 1 gaps. A hypothesis was formulated regarding the penetration depth of the radar waves within the canopy. It claims that the depth of penetration is simply related to the quantisation level of the radar backscattered signal. A higher level of bits per pixels allows for capturing weaker signals arriving from the lower levels of the tree crown. Assuming certain generic and simplified shapes of tree crowns including cone, paraboloid, sphere and spherical cap, it was possible to model analytically Type 2 gaps. The Monte Carlo simulation method was used to investigate relationships between the impenetrability and various configurations of a modelled forest. One of the most important findings is that impenetrability is largely explainable by the gaps between trees. A much less important role is played by the penetrability into the crown cover. Another important finding is that the impenetrability strongly correlates with the vegetation density. Using this feature, a method for vegetation density mapping called the mean maximum impenetrability (MMI) method is proposed. Unlike the traditional methods of forest inventories, the MMI method allows for a much more realistic inventory of vegetation cover, because it is able to capture an in situ or current situation on the ground, but not for areas that are nominally classified as a “forest-to-be”. The MMI method also allows for the mapping of landscape variation in the forest or vegetation density, which is a novel and exciting feature of the new 3D remote sensing (3DRS) technique. Besides the inventory-type applications, the MMI method can be used as a forest change detection method. For maximum effectiveness of the MMI method, an object-based change detection approach is preferred. A minimum requirement for the MMI method is a time-lapsed reference dataset in the form, for example, of an existing forest map of the area of interest, or a vegetation density map prepared using InSAR datasets. Preliminary tests aimed at finding a degree of correlation between the impenetrability and other types of passive and active remote sensing data sources, including TerraSAR-X, NDVI and PALSAR, proved that the method most sensitive to vegetation density was the Japanese PALSAR - L-band SAR system. Unfortunately, PALSAR backscattered signals become very noisy for impenetrability below 15 m. This means that PALSAR has severe limitations for low loadings of the biomass per unit area. The proposed applications of the InSAR data will remain indispensable wherever cloud cover obscures the sky in a persistent manner, which makes suitable optical data acquisition extremely time-consuming or nearly impossible. A limitation of the MMI method is due to the fact that the impenetrability is calculated using a reference DTM, which must be available beforehand. In many countries around the world, appropriate quality DTMs are still unavailable. A possible solution to this obstacle is to use a DEM that was derived using P-band InSAR elevations or LiDAR. It must be noted, however, that in many cases, two InSAR datasets separated by time of the same area are sufficient for forest change detection or similar applications

Estimating mangrove canopy height and above-ground biomass in the Everglades National Park with airborne LiDAR and TanDEM-X data

Mangrove forests are important natural ecosystems due to their ability to capture and store large amounts of carbon. Forest structural parameters, such as canopy height and above-ground biomass (AGB), provide a good measure for monitoring temporal changes in carbon content. The protected coastal mangrove forest of the Everglades National Park (ENP) provides an ideal location for studying these processes, as harmful human activities are minimal. We estimated mangrove canopy height and AGB in the ENP using Airborne LiDAR/Laser (ALS) and TanDEM-X (TDX) datasets acquired between 2011 and 2013. Analysis of both datasets revealed that mangrove canopy height can reach up to ~25 m and AGB can reach up to ~250 Mg·ha-1. In general, mangroves ranging from 9 m to 12 m in stature dominate the forest canopy. The comparison of ALS and TDX canopy height observations yielded an R2 = 0.85 and Root Mean Square Error (RMSE) = 1.96 m. Compared to a previous study based on data acquired during 2000-2004, our analysis shows an increase in mangrove stature and AGB, suggesting that ENP mangrove forests are continuing to accumulate biomass. Our results suggest that ENP mangrove forests have managed to recover from natural disturbances, such as HurricaneWilma

Elevation and Deformation Extraction from TomoSAR

3D SAR tomography (TomoSAR) and 4D SAR differential tomography (Diff-TomoSAR) exploit multi-baseline SAR data stacks to provide an essential innovation of SAR Interferometry for many applications, sensing complex scenes with multiple scatterers mapped into the same SAR pixel cell. However, these are still influenced by DEM uncertainty, temporal decorrelation, orbital, tropospheric and ionospheric phase distortion and height blurring. In this thesis, these techniques are explored. As part of this exploration, the systematic procedures for DEM generation, DEM quality assessment, DEM quality improvement and DEM applications are first studied. Besides, this thesis focuses on the whole cycle of systematic methods for 3D & 4D TomoSAR imaging for height and deformation retrieval, from the problem formation phase, through the development of methods to testing on real SAR data. After DEM generation introduction from spaceborne bistatic InSAR (TanDEM-X) and airborne photogrammetry (Bluesky), a new DEM co-registration method with line feature validation (river network line, ridgeline, valley line, crater boundary feature and so on) is developed and demonstrated to assist the study of a wide area DEM data quality. This DEM co-registration method aligns two DEMs irrespective of the linear distortion model, which improves the quality of DEM vertical comparison accuracy significantly and is suitable and helpful for DEM quality assessment. A systematic TomoSAR algorithm and method have been established, tested, analysed and demonstrated for various applications (urban buildings, bridges, dams) to achieve better 3D & 4D tomographic SAR imaging results. These include applying Cosmo-Skymed X band single-polarisation data over the Zipingpu dam, Dujiangyan, Sichuan, China, to map topography; and using ALOS L band data in the San Francisco Bay region to map urban building and bridge. A new ionospheric correction method based on the tile method employing IGS TEC data, a split-spectrum and an ionospheric model via least squares are developed to correct ionospheric distortion to improve the accuracy of 3D & 4D tomographic SAR imaging. Meanwhile, a pixel by pixel orbit baseline estimation method is developed to address the research gaps of baseline estimation for 3D & 4D spaceborne SAR tomography imaging. Moreover, a SAR tomography imaging algorithm and a differential tomography four-dimensional SAR imaging algorithm based on compressive sensing, SAR interferometry phase (InSAR) calibration reference to DEM with DEM error correction, a new phase error calibration and compensation algorithm, based on PS, SVD, PGA, weighted least squares and minimum entropy, are developed to obtain accurate 3D & 4D tomographic SAR imaging results. The new baseline estimation method and consequent TomoSAR processing results showed that an accurate baseline estimation is essential to build up the TomoSAR model. After baseline estimation, phase calibration experiments (via FFT and Capon method) indicate that a phase calibration step is indispensable for TomoSAR imaging, which eventually influences the inversion results. A super-resolution reconstruction CS based study demonstrates X band data with the CS method does not fit for forest reconstruction but works for reconstruction of large civil engineering structures such as dams and urban buildings. Meanwhile, the L band data with FFT, Capon and the CS method are shown to work for the reconstruction of large manmade structures (such as bridges) and urban buildings