Study of the speckle noise effects over the eigen decomposition of polarimetric SAR data: a review

This paper is focused on considering the effects of speckle noise on the eigen decomposition of the co- herency matrix. Based on a perturbation analysis of the matrix, it is possible to obtain an analytical expression for the mean value of the eigenvalues and the eigenvectors, as well as for the Entropy, the Anisotroopy and the dif- ferent a angles. The analytical expressions are compared against simulated polarimetric SAR data, demonstrating the correctness of the different expressions.Peer ReviewedPostprint (published version

The BIOMASS level 2 prototype processor : design and experimental results of above-ground biomass estimation

BIOMASS is ESA’s seventh Earth Explorer mission, scheduled for launch in 2022. The satellite will be the first P-band SAR sensor in space and will be operated in fully polarimetric interferometric and tomographic modes. The mission aim is to map forest above-ground biomass (AGB), forest height (FH) and severe forest disturbance (FD) globally with a particular focus on tropical forests. This paper presents the algorithms developed to estimate these biophysical parameters from the BIOMASS level 1 SAR measurements and their implementation in the BIOMASS level 2 prototype processor with a focus on the AGB product. The AGB product retrieval uses a physically-based inversion model, using ground-canceled level 1 data as input. The FH product retrieval applies a classical PolInSAR inversion, based on the Random Volume over Ground Model (RVOG). The FD product will provide an indication of where significant changes occurred within the forest, based on the statistical properties of SAR data. We test the AGB retrieval using modified airborne P-Band data from the AfriSAR and TropiSAR campaigns together with reference data from LiDAR-based AGB maps and plot-based ground measurements. For AGB estimation based on data from a single heading, comparison with reference data yields relative Root Mean Square Difference (RMSD) values mostly between 20% and 30%. Combining different headings in the estimation process significantly improves the AGB retrieval to slightly less than 20%. The experimental results indicate that the implemented retrieval scheme provides robust results that are within mission requirements

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

Wetland mapping and monitoring using polarimetric and interferometric synthetic aperture radar (SAR) data and tools

Wetlands are home to a great variety of flora and fauna species and provide several unique environmental functions, such as controlling floods, improving water-quality, supporting wildlife habitat, and shoreline stabilization. Detailed information on spatial distribution of wetland classes is crucial for sustainable management and resource assessment. Furthermore, hydrological monitoring of wetlands is also important for maintaining and preserving the habitat of various plant and animal species. This thesis investigates the existing knowledge and technological challenges associated with wetland mapping and monitoring and evaluates the limitations of the methodologies that have been developed to date. The study also proposes new methods to improve the characterization of these productive ecosystems using advanced remote sensing (RS) tools and data. Specifically, a comprehensive literature review on wetland monitoring using Synthetic Aperture Radar (SAR) and Interferometric SAR (InSAR) techniques is provided. The application of the InSAR technique for wetland mapping provides the following advantages: (i) the high sensitivity of interferometric coherence to land cover changes is taken into account and (ii) the exploitation of interferometric coherence for wetland classification further enhances the discrimination between similar wetland classes. A statistical analysis of the interferometric coherence and SAR backscattering variation of Canadian wetlands, which are ignored in the literature, is carried out using multi-temporal, multi-frequency, and multi-polarization SAR data. The study also examines the capability of compact polarimetry (CP) SAR data, which will be collected by the upcoming RADARSAT Constellation Mission (RCM) and will constitute the main source of SAR observation in Canada, for wetland mapping. The research in this dissertation proposes a methodology for wetland classification using the synergistic use of intensity, polarimetry, and interferometry features using a novel classification framework. Finally, this work introduces a novel model based on the deep convolutional neural network (CNN) for wetland classification that can be trained in an end-to-end scheme and is specifically designed for the classification of wetland complexes using polarimetric SAR (PolSAR) imagery. The results of the proposed methods are promising and will significantly contribute to the ongoing efforts of conservation strategies for wetlands and monitoring changes. The approaches presented in this thesis serve as frameworks, progressing towards an operational methodology for mapping wetland complexes in Canada, as well as other wetlands worldwide with similar ecological characteristics

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

TomoSAR Mapping of 3D Forest Structure: Contributions of L-Band Configurations

Synthetic Aperture Radar (SAR) measurements are unique for mapping forest 3D structure and its changes in time. Tomographic SAR (TomoSAR) configurations exploit this potential by reconstructing the 3D radar reflectivity. The frequency of the SAR measurements is one of the main parameters determining the information content of the reconstructed reflectivity in terms of penetration and sensitivity to the individual vegetation elements. This paper attempts to review and characterize the structural information content of L-band TomoSAR reflectivity reconstructions, and their potential to forest structure mapping. First, the challenges in the accurate TomoSAR reflectivity reconstruction of volume scatterers (which are expected to dominate at L-band) and to extract physical structure information from the reconstructed reflectivity is addressed. Then, the L-band penetration capability is directly evaluated by means of the estimation performance of the sub-canopy ground topography. The information content of the reconstructed reflectivity is then evaluated in terms of complementary structure indices. Finally, the dependency of the TomoSAR reconstruction and of its structural information to both the TomoSAR acquisition geometry and the temporal change of the reflectivity that may occur in the time between the TomoSAR measurements in repeat-pass or bistatic configurations is evaluated. The analysis is supported by experimental results obtained by processing airborne acquisitions performed over temperate forest sites close to the city of Traunstein in the south of Germany

Modeling of Subsurface Scattering from Ice Sheets for Pol-InSAR Applications

Remote sensing is a fundamental tool to measure the dynamics of ice sheets and provides valuable information for ice sheet projections under a changing climate. There is, however, the potential to further reduce the uncertainties in these projections by developing innovative remote sensing methods. One of these remote sensing techniques, the polarimetric synthetic aperture radar interferometry (Pol-InSAR), is known since decades to have the potential to assess the geophysical properties below the surface of ice sheets, because of the penetration of microwave signals into dry snow, firn, and ice. Despite this, only very few studies have addressed this topic and the development of robust Pol-InSAR applications is at an early stage. Two potential Pol-InSAR applications are identified as the motivation for this thesis. First, the estimation and compensation of the penetration bias in digital elevation models derived with SAR interferometry. This bias can lead to errors of several meters or even tens of meters in surface elevation measurements. Second, the estimation of geophysical properties of the subsurface of glaciers and ice sheets using Pol-InSAR techniques. There is indeed potential to derive information about melt-refreeze processes within the firn, which are related to density and affect the mass balance. Such Pol-InSAR applications can be a valuable information source with the potential for monthly ice sheet wide coverage and high spatial resolution provided by the next generation of SAR satellites. However, the required models to link the Pol-InSAR measurements to the subsurface properties are not yet established. The aim of this thesis is to improve the modeling of the vertical backscattering distribution in the subsurface of ice sheets and its effect on polarimetric interferometric SAR measurements at different frequencies. In order to achieve this, polarimetric interferometric multi-baseline SAR data at different frequencies and from two different test sites on the Greenland ice sheet are investigated. This thesis contributes with three concepts to a better understanding and to a more accurate modeling of the vertical backscattering distribution in the subsurface of ice sheets. First, the integration of scattering from distinct subsurface layers. These are formed by refrozen melt water in the upper percolation zone and cause an interesting coherence undulation pattern, which cannot be explained with previously existing models. This represents a first link between Pol-InSAR data and geophysical subsurface properties. The second step is the improved modeling of the general vertical backscattering distribution of the subsurface volume. The advantages of more flexible volume models are demonstrated, but interestingly, the simple modification of a previously existing model with a vertical shift parameter lead to the best agreement between model and data. The third contribution is the model based compensation of the penetration bias, which is experimentally validated. At the investigated test sites, it becomes evident that the model based estimates of the surface elevations are more accurate than the interferometric phase center locations, which are conventionally used to derive surface elevations of ice sheets. This thesis therefore improves the state of the art of subsurface scattering modeling for Pol-InSAR applications, demonstrates the model-based penetration bias compensation, and makes a further research step towards the retrieval of geophysical subsurface information with Pol-InSAR

Spatial and temporal statistics of SAR and InSAR observations for providing indicators of tropical forest structural changes due to forest disturbance

Tropical forests are extremely important ecosystems which play a substantial role in the global carbon budget and are increasingly dominated by anthropogenic disturbance through deforestation and forest degradation, contributing to emissions of greenhouse gases to the atmosphere. There is an urgent need for forest monitoring over extensive and inaccessible tropical forest which can be best accomplished using spaceborne satellite data. Currently, two key processes are extremely challenging to monitor: forest degradation and post-disturbance re-growth. The thesis work focuses on these key processes by considering change indicators derived from radar remote sensing signal that arise from changes in forest structure. The problem is tackled by exploiting spaceborne Synthetic Aperture Radar (SAR) and Interferometric SAR (InSAR) observations, which can provide forest structural information while simultaneously being able to collect data independently of cloud cover, haze and daylight conditions which is a great advantage over the tropics. The main principle of the work is that a connection can be established between the forest structure distribution in space and signal variation (spatial statistics) within backscatter and Digital Surface Models (DSMs) provided by SAR. In turn, forest structure spatial characteristics and changes are used to map forest condition (intact or degraded) or disturbance. The innovative approach focuses on looking for textural patterns (and their changes) in radar observations, then connecting these patterns to the forest state through supporting evidence from expert knowledge and auxiliary remote sensing observations (e.g. high resolution optical, aerial photography or LiDAR). These patterns are descriptors of the forest structural characteristics in a statistical sense, but are not estimates of physical properties, such as above-ground biomass or canopy height. The thesis tests and develops methods using novel remote sensing technology (e.g. single-pass spaceborne InSAR) and modern image statistical analysis methods (wavelet-based space-scale analysis). The work is developed on an experimental basis and articulated in three test cases, each addressing a particular observational setting, analytical method and thematic context. The first paper deals with textural backscatter patterns (C-band ENVISAT ASAR and L-band ALOS PALSAR) in semi-deciduous closed forest in Cameroon. Analysis concludes that intact forest and degraded forest (arising from selective logging) are significantly different based on canopy structural properties when measured by wavelet based space-scale analysis. In this case, C-band data are more effective than longer wavelength L-band data. Such a result could be explained by the lower wave penetration into the forest volume at shorter wavelength, with the mechanism driving the differences between the two forest states arising from upper canopy heterogeneity. In the second paper, wavelet based space-scale analysis is also used to provide information on upper canopy structure. A DSM derived from TanDEM-X acquired in 2014 was used to discriminate primary lowland Dipterocarp forest, secondary forest, mixed-scrub and grassland in the Sungai Wain Protection Forest (East Kalimantan, Indonesian Borneo) which was affected by the 1997/1998 El Niño Southern Oscillation (ENSO). The Jeffries- Matusita separability of wavelet spectral measures of InSAR DSMs between primary and secondary forest was in some cases comparable to results achieved by high resolution LiDAR data. The third test case introduces a temporal component, with change detection aimed at detecting forest structure changes provided by differencing TanDEM-X DSMs acquired at two dates separated by one year (2012-2013) in the Republic of Congo. The method enables cancelling out the component due to terrain elevation which is constant between the two dates, and therefore the signal related to the forest structure change is provided. Object-based change detection successfully mapped a gradient of forest volume loss (deforestation/forest degradation) and forest volume gain (post-disturbance re-growth). Results indicate that the combination of InSAR observations and wavelet based space-scale analysis is the most promising way to measure differences in forest structure arising from forest fires. Equally, the process of forest degradation due to shifting cultivation and post-disturbance re-growth can be best detected using multiple InSAR observations. From the experiments conducted, single-pass InSAR appears to be the most promising remote sensing technology to detect forest structure changes, as it provides three-dimensional information and with no temporal decorrelation. This type of information is not available in optical remote sensing and only partially available (through a 2D mapping) in SAR backscatter. It is advised that future research or operational endeavours aimed at mapping and monitoring forest degradation/regrowth should take advantage of the only currently available high resolution spaceborne single-pass InSAR mission (TanDEM-X). Moreover, the results contribute to increase knowledge related to the role of SAR and InSAR for monitoring degraded forest and tracking the process of forest degradation which is a priority but still highly challenging to detect. In the future the techniques developed in the thesis work could be used to some extent to support REDD+ initiatives

Growing stock volume estimation in temperate forsted areas using a fusion approach with SAR Satellites Imagery

Forest monitoring plays a central role in the context of global warming mitigation and in the assessment of forest resources. To meet these challenges, significant efforts have been made by scientists to develop new feasible remote sensing techniques for the retrieval of forest parameters. However, much work remains to be done in this area, in particular in establishing global assessments of forest biomass. In this context, this Ph.D. Thesis presents a complete methodology for estimating Growing Stock Volume (GSV) in temperate forested areas using a fusion approach based on Synthetic-Aperture Radar (SAR) satellite imagery. The investigations which were performed focused on the Thuringian Forest, which is located in Central Germany. The satellite data used are composed of an extensive set of L-band (ALOS PALSAR) and X-band (TerraSAR-X, TanDEM-X, Cosmo-SkyMed) images, which were acquired in various sensor configurations (acquisition modes, polarisations, incidence angles). The available ground data consists of a forest inventory delivered by the local forest offices. Weather measurements and a LiDAR DEM complete the datasets. The research showed that together with the topography, the forest structure and weather conditions generally limited the sensitivity of the SAR signal to GSV. The best correlations were obtained with ALOS PALSAR (R2 = 0.61) and TanDEM-X (R2 = 0.72) interferometric coherences. These datasets were chosen for the retrieval of GSV in the Thuringian Forest and led with regressions to an root-mean-square error (RMSE) in the range of 100─200 m3ha-1. As a final achievement of this thesis, a methodology for combining the SAR information was developed. Assuming that there are sufficient and adequate remote sensing data, the proposed fusion approach may increase the biomass maps accuracy, their spatial extension and their updated frequency. These characteristics are essential for the future derivation of accurate, global and robust forest biomass maps

Crop development monitoring from Synthetic Aperture Radar (SAR) imagery

Satellite remote sensing plays a vital role in providing large-scale and timely data to stakeholders of the agricultural supply chain. This allows for informed decision-making that promotes sustainable and cost-effective crop management practices. In particular, data derived from satellite-based Synthetic Aperture Radar (SAR) systems, provide opportunities for continuous crop monitoring, taking advantage of its ability to acquire images during day or night and under almost all weather conditions. Moreover, an abundance of SAR data can be anticipated in the next 5 years with the launch of several international SAR missions. However, research on crop development monitoring with data from SAR satellites has not been as widely studied as with data derived from passive multi-spectral satellites and contributions can be made to the current state-of-the-art techniques. This thesis aims at improving the current knowledge on the use of satellite-based SAR imagery for crop development monitoring. This is approached by developing novel methodologies and detailed interpretations of multitemporal SAR and Polarimetric SAR (PolSAR) responses to crop growth in three different test sites. Chapter two presents a detailed analysis of the Sentinel-1 SAR satellite response to asparagus crop development in Peru, investigating the capabilities of the sensor to capture seasonality effects as well as providing an interpretation of the temporal backscatter signature. This is complemented with a case study where a multiple-output random forest regression algorithm is used to successfully retrieve crop growth stage from Sentinel-1 data and temperature measurements. Following the limitations identified with this approach, a methodology that builds upon ideas of Bayesian Filtering Frameworks (BFFs) for crop monitoring is proposed in chapter three. It incorporates Gaussian processes to model crop dynamics as well as to model the remote sensing response to the crop state. Using this approach, it is possible to derive daily predictions with the associated uncertainties, to combine in near-real-time data from active and passive satellites as well as to estimate past and future crop key events that are of strategic importance for different stakeholders. The final section of this thesis looks at the new developments of the SAR technology considering that future open access missions will provide Quad Polarimetric SAR data. An algorithm based on multitemporal PolSAR change detection is introduced in chapter four. It defines a Change Matrix to encode an interpretable representation of the crop dynamics as captured by the evolution of the scattering mechanisms over time. We use rice fields in Spain and multiple cereal crops in Canada to test the use of the algorithm for crop monitoring. A supervised learning-based crop type classification methodology is then proposed with the same method by using the encoded scattering mechanisms as input for a neural-network-based classifier, achieving comparable performances to state-of-the-art classifiers. The results obtained in this thesis represent novel additions to the literature that contribute to our understanding and successful use of SAR imagery for agricultural monitoring. For the first time, a detailed analysis of asparagus crops is presented. It is a key crop for agricultural exports of Peru, the largest exporter of asparagus in the world. Secondly, two key contributions to the state of the art BFFs for crop monitoring are presented: a) A better exploitation of the SAR temporal dimension and an application with freely available data and b) given that it is a learning-based approach, it overcomes current limitations of transferability among crop types and regions. Finally, the PolSAR change detection approach presented in the last thesis chapter, provides a novel and easy-to-interpret tool for both crop monitoring and crop type mapping applications