Disaster debris estimation using high-resolution polarimetric stereo-SAR

AbstractThis paper addresses the problem of debris estimation which is one of the most important initial challenges in the wake of a disaster like the Great East Japan Earthquake and Tsunami. Reasonable estimates of the debris have to be made available to decision makers as quickly as possible. Current approaches to obtain this information are far from being optimal as they usually rely on manual interpretation of optical imagery. We have developed a novel approach for the estimation of tsunami debris pile heights and volumes for improved emergency response. The method is based on a stereo-synthetic aperture radar (stereo-SAR) approach for very high-resolution polarimetric SAR. An advanced gradient-based optical-flow estimation technique is applied for optimal image coregistration of the low-coherence non-interferometric data resulting from the illumination from opposite directions and in different polarizations. By applying model based decomposition of the coherency matrix, only the odd bounce scattering contributions are used to optimize echo time computation. The method exclusively considers the relative height differences from the top of the piles to their base to achieve a very fine resolution in height estimation. To define the base, a reference point on non-debris-covered ground surface is located adjacent to the debris pile targets by exploiting the polarimetric scattering information. The proposed technique is validated using in situ data of real tsunami debris taken on a temporary debris management site in the tsunami affected area near Sendai city, Japan. The estimated height error is smaller than 0.6m RMSE. The good quality of derived pile heights allows for a voxel-based estimation of debris volumes with a RMSE of 1099m3. Advantages of the proposed method are fast computation time, and robust height and volume estimation of debris piles without the need for pre-event data or auxiliary information like DEM, topographic maps or GCPs

Electromagnetic modeling for SAR polarimetry and interferometry

Investigation of the globe remotely from hundreds of kilometers altitude, and fast growing of environmental and civil problems, triggered the necessity of development of new and more advanced techniques. Electromagnetic modeling of polarimetry and interferometry has always been a key driver in remote sensing research, ever since of the First pioneering sensors were launched. Polarimetric and interferometric SAR (Synthetic Aperture Radar) surveillance and mapping of the Earth surface has been attracting lots of interest since 1970s. This thesis covers two SAR's main techniques: (1) space-borne Interferometric Synthetic Aperture Radar (InSAR), which has been used to measure the Earth's surface deformation widely, and (2) SAR Polarimetry, which has been used to retrieve soil and vegetation physical parameters in wide areas. Time-series InSAR methodologies such as PSI (Permanent Scatterer Interferometry) are designed to estimate the temporal characteristics of the Earth's deformation rates from multiple InSAR images acquired over time. These techniques also enable us to overcome the limitations that conventional InSAR suffer, with a very high accuracy and precision. In this thesis, InSAR time-series analysis and modeling basis, as well as a case study in the Campania region (Italy), have been addressed. The Campania region is characterized by intense urbanization, active volcanoes, complicated fault systems, landslides, subsidence, and hydrological instability; therefore, the stability of public transportation structures is highly concerned. Here Differential Interferometric Synthetic Aperture Radar (DInSAR), and PSI techniques have been applied to a stack of 25 X-band radar images of Cosmo-SkyMed (CSK) satellites collected over an area in Campania (Italy), in order to monitor the railways' stability. The study area was already under investigation with older, low-resolution sensors like ERS1&2 and ENVISAT-ASAR before, but the number of obtained persistent scatterers (PSs) was too limited to get useful results. With regard to SAR polarimetry, in this thesis a fully polarimetirc SAR simulator has been presented, which is based on the use of sound direct electromagnetic models and it is able to provide as output the simulated raw data of all the three polarization channels in such a way as to obtain the correct covariance or coherence matrixes on the final focused polarimetic radar images. A fast Fourier-domain approach is used for the generation of raw signals. Presentation of theory is supplemented by meaningful experimental results, including a comparison of simulations with real polarimetric scattering data

Compaction of C-band synthetic aperture radar based sea ice information for navigation in the Baltic Sea

In this work operational sea ice synthetic aperture radar (SAR) data products were improved and developed. A SAR instrument is transmitting electromagnetic radiation at certain wavelengths and measures the radiation which is scattered back towards the instrument from the target, in our case sea and sea ice. The measured backscattering is converted to an image describing the target area through complex signal processing. The images, however, differ from optical images, i.e. photographs, and their visual interpretation is not straightforward. The main idea in this work has been to deliver the essential SAR-based sea ice information to end-users (typically on ships) in a compact and user-friendly format. The operational systems at Finnish Institute of Marine Research (FIMR) are currently based on the data received from a Canadian SAR-satellite, Radarsat-1. The operational sea ice classification, developed by the author with colleagues, has been further developed. One problem with the SAR data is typically that the backscattering varies depending on the incidence angle. The incidence angle is the angle in which the transmitted electromagnetic wave meets the target surface and it varies within each SAR image and between different SAR images depending on the measuring geometry. To improve this situation, an incidence angle correction algorithm to normalize the backscattering over the SAR incidence angle range for Baltic Sea ice has been developed as part of this work. The algorithm is based on SAR backscattering statistics over the Baltic Sea. To locate different sea ice areas in SAR images, a SAR segmentation algorithm based on pulse-coupled neural networks has been developed and tested. The parameters have been tuned suitable for the operational data in use at FIMR. The sea ice classification is based on this segmentation and the classification is segment-wise rather than pixel-wise. To improve SAR-based distinguishing between sea ice and open water an open water detection algorithm based on segmentation and local autocorrelation has been developed. Also ice type classification based on higher-order statistics and independent component analysis have been studied to get an improved SAR-based ice type classification. A compression algorithm for compressing sea ice SAR data for visual use has been developed. This algorithm is based on the wavelet decomposition, zero-tree structure and arithmetic coding. Also some properties of the human visual system were utilized. This algorithm was developed to produce smaller compressed SAR images, with a reasonable visual quality. The transmission of the compressed images to ships with low-speed data connections in reasonable time is then possible. One of the navigationally most important sea ice parameters is the ice thickness. SAR-based ice thickness estimation has been developed and evaluated as part of this work. This ice thickness estimation method uses the ice thickness history derived from digitized ice charts, made daily at the Finnish Ice Service, as its input, and updates this chart based on the novel SAR data. The result is an ice thickness chart representing the ice situation at the SAR acquisition time in higher resolution than in the manually made ice thickness charts. For the evaluation of the results a helicopter-borne ice thickness measuring instrument, based on electromagnetic induction and laser altimeter, was used.reviewe

Satellite measurement of ocean turbulence

Turbulence and mixing in the surface layer of the ocean is a significant element in the combined ocean-atmosphere system, and plays a considerable role in the transfer of heat, gas and momentum across the air-sea boundary. Furthermore, improving knowledge of the evolution of energy within the ocean system, both globally and locally, holds importance for improving our understanding of the dynamics of the ocean at large- and small-scales. As such, insight into turbulence and turbulent flows at the ocean surface is becoming increasingly important for its role in ocean-atmosphere exchange and, from a wider perspective, climate change.A research project was initiated to understand the role that spacecraft remote-sensing may play in improving observation of “turbulence” (in a broad sense) in the ocean, and for identifying how steps towards such observation may be made. An initial, exploratory study identified the potential benefit of Synthetic Aperture Radar in “bridging the gap” between in-situ and remote observations o