Davide Oscar Nitti

Impact of DEM-Assisted Coregistration on High-Resolution SAR Interferometry

Image alignment is a crucial step in synthetic aperture radar (SAR) interferometry. Interferogram formation requires images to be coregistered with an accuracy of better than a few tenths of a resolution cell to avoid significant loss of phase coherence. In conventional interferometric precise coregistration methods for full-resolution SAR data, a 2-D polynomial of low degree is usually chosen as warp function, and the polynomial parameters are estimated through least squares fit from the shifts measured on image windows. In case of rough topography or long baselines, the polynomial approximation may become inaccurate, leading to local misregistrations. These effects increase with spatial resolution of the sensor. An improved elevation-assisted image-coregistration procedure can be adopted to provide better prediction of the offset vectors. This approach computes pixel by pixel the correspondence between master and slave acquisitions by using the orbital data and a reference digital elevation model (DEM). This paper aims to assess the performance of this procedure w.r.t. the “standard” one based on polynomial approximation. Analytical relationships and simulations are used to evaluate the improvement of the DEM-assisted procedure w.r.t. the polynomial approximation as well as the impact of the finite vertical accuracy of the DEM on the final coregistration precision for different resolutions and baselines. The two approaches are then evaluated experimentally by processing high-resolution SAR data provided by the COnstellation of small Satellites for the Mediterranean basin Observation (COSMO/SkyMed) and TerraSAR-X missions, acquired over mountainous areas in Italy and Tanzania, respectively. Residual-range pixel offsets and interferometric coherence are used as quality figure.

Using C/X-band SAR interferometry and GNSS measurements for the Assisi landslide analysis

This work presents an analysis of the applicability of synthetic aperture radar (SAR) interferometry to landslide monitoring. This analysis was carried out by using different interferometric approaches, different spaceborne SAR data (both in the C-band and in the X-band), and in situ global navigation satellite system (GNSS) measurements. In particular, we investigated both the reliability of displacement monitoring and the issues of the cross-comparison and validation of the interferometric synthetic aperture radar (InSAR) results. The work was focused on the slow-moving landslide that affects a relevant part of the urban area of the historical town of Assisi (Italy). A C-band ENVISAT advanced synthetic aperture radar (ENVISAT ASAR) dataset acquired between 2003 and 2010 was processed by using two different interferometric techniques, to allow cross-comparison of the obtained displacement maps. Good correspondence between the results was found, and a deeper analysis of the movement field was possible. Results were further compared to a set of GNSS measurements with a 7 year overlap with SAR data. A comparison was made for each GNSS marker with the surrounding SAR scatterers, trying to take into account local topological effects, when possible. Further, the high-resolution X-band acquired on both ascending and descending tracks by the COSMO-SkyMed (CSK) constellation was processed. The resultant displacement fields show good agreement with C-band and GNSS measurements and a sensible increase in the density of measurements.

TerraSAR-X InSAR multipass analysis on Venice (Italy)

The TerraSAR-X (copyright) mission, launched in 2007, carries a new X-band Synthetic Aperture Radar (SAR) sensor optimally suited for SAR interferometry (InSAR), thus allowing very promising application of InSAR techniques for the risk assessment on areas with hydrogeological instability and especially for multi-temporal analysis, such as Persistent Scatterer Interferometry (PSI) techniques, originally developed at Politecnico di Milano. The SPINUA (Stable Point INterferometry over Unurbanised Areas) technique is a PSI processing methodology which has originally been developed with the aim of detection and monitoring of coherent PS targets in non or scarcely-urbanized areas. The main goal of the present work is to describe successful applications of the SPINUA PSI technique in processing X-band data. Venice has been selected as test site since it is in favorable settings for PSI investigations (urban area containing many potential coherent targets such as buildings) and in view of the availability of a long temporal series of TerraSAR-X stripmap acquisitions (27 scenes in all). The Venice Lagoon is affected by land sinking phenomena, whose origins are both natural and man-induced. The subsidence of Venice has been intensively studied for decades by determining land displacements through traditional monitoring techniques (leveling and GPS) and, recently, by processing stacks of ERS/ENVISAT SAR data. The present work is focused on an independent assessment of application of PSI techniques to TerraSAR-X stripmap data for monitoring the stability of the Venice area. Thanks to its orbital repeat cycle of only 11 days, less than a third of ERS/ENVISAT C-band missions, the maximum displacement rate that can be unambiguously detected along the Line-of-Sight (LOS) with TerraSAR-X SAR data through PSI techniques is expected to be about twice the corresponding value of ESA C-band missions, being directly proportional to the sensor wavelength and inversely proportional to the revisit time. When monitoring displacement phenomena which are known to be within the C-band rate limits, the increased repeat cycle of TerraSAR-X offers the opportunity to decimate the stack of TerraSAR-X data, e.g. by doubling the temporal baseline between subsequent acquisitions. This strategy can be adopted for reducing both economic and computational processing costs. In the present work, the displacement rate maps obtained through SPINUA with and without decimation of the number of Single Look Complex (SLC) acquisitions are compared. In particular, it is shown that with high spatial resolution SAR data, reliable displacement maps could be estimated through PSI techniques with a number of SLCs much lower than in C-band.

On the use of COSMO/SkyMed data and Weather Models for interferometric DEM generation

Abstract This work experiments the potentialities of COSMO/SkyMed (CSK) data in providing interferometric Digital Elevation Model (DEM). We processed a stack of CSK data for measuring with meter accuracy the ground elevation on the available coherent targets, and used these values to check the accuracy of DEMs derived from 5 tandem-like CSK pairs. In order to suppress the atmospheric signal we experimented a classical spatial filtering of the differential phase as well as the use of numerical weather prediction (NWP) model RAMS. Tandem-like pairs with normal baselines higher than 300 m allows to derive DEMs fulfilling the HRTI Level 3 specifications on the relative vertical accuracy, while the use of NWP models still seems unfeasible especially for X-band.

Feasibility of Using Synthetic Aperture Radar to Aid UAV Navigation

This study explores the potential of Synthetic Aperture Radar (SAR) to aid Unmanned Aerial Vehicle (UAV) navigation when Inertial Navigation System (INS) measurements are not accurate enough to eliminate drifts from a planned trajectory. This problem can affect medium-altitude long-endurance (MALE) UAV class, which permits heavy and wide payloads (as required by SAR) and flights for thousands of kilometres accumulating large drifts. The basic idea is to infer position and attitude of an aerial platform by inspecting both amplitude and phase of SAR images acquired onboard. For the amplitude-based approach, the system navigation corrections are obtained by matching the actual coordinates of ground landmarks with those automatically extracted from the SAR image. When the use of SAR amplitude is unfeasible, the phase content can be exploited through SAR interferometry by using a reference Digital Terrain Model (DTM). A feasibility analysis was carried out to derive system requirements by exploring both radiometric and geometric parameters of the acquisition setting. We showed that MALE UAV, specific commercial navigation sensors and SAR systems, typical landmark position accuracy and classes, and available DTMs lead to estimate UAV coordinates with errors bounded within ±12 m, thus making feasible the proposed SAR-based backup system.

C/X-Band SAR Interferometry Used to Monitor Slope Instability in Daunia, Italy

We apply multi-temporal Persistent Scatterer Interferometry (PSI) analysis to investigate slope instability in the Daunia region in the Southern Apennine Mountains. Daunia includes many small hill-top towns affected by landslides and is of particular interest for the Civil Protection – Regione Puglia Authority, one of the end users of the PSI deformation maps. The SPINUA algorithm is used to perform interferometric analysis and detect, with mm precision, the presence of ground surface movements. Consistent results on very slow displacements are obtained using the radar imagery acquired between 2002 and 2010 by the ENVISAT ESA satellite (C-band, medium spatial resolution sensor) and the images acquired between 2010 and 2011 by the X-band high resolution sensor onboard the TerraSAR-X satellite. Thanks to the finer spatial resolution the X-band PSI applications are very promising for monitoring single man-made structures and slope/ground instability in areas where C-band PS density is low.

High Resolution PSI for Mapping Ground Deformations and Infrastructure Instability

Persistent scatterers interferometry (PSI) based on space-borne radar data can provide thousands km2 coverage and precise (mm resolution), spatially dense measurements (from hundreds to over thousands points/km2) on ground deformations and infrastructure instability. Furthermore, the practical applicability of PSI is now improved thanks to the increased data availability and the better capabilities of the new space radar sensors (Cosmo-SkyMed, TerraSAR-X) in terms of resolution (from ~3 to 1 m) and revisit time (from 11 to 4 days). We compare results from medium and high resolution PSI investigations of subsidence, slope and associated infrastructure instability in two areas with different geo-environmental characteristics to illustrate (i) the potential in providing site-specific information, (ii) the advantages of high resolution radar data. We also focus on technical and interpretation issues in PSI applications and offer specific user guidelines, with emphasis on the benefits resulting from the exploitation of new generation radar sensors.

A Bayesian network for flood detection

We apply a Bayesian Network (BN) paradigm to the problem of monitoring flood events through synthetic aperture radar (SAR) and interferometric SAR (InSAR) data. BNs are well-founded statistical tools which help formalizing the information coming from heterogeneous sources, such as remotely sensed images, LiDAR data, and topography. The approach is tested on the fluvial floodplains of the Basilicata region (southern Italy), which have been subject to recurrent flooding events in the last years. Results show maps efficiently representing the different scattering/coherence classes with high accuracy, and also allowing separating the multitemporal dimension of the data, where available. The BN approach proves thus helpful to gain insight into the complex phenomena related to floods, possibly also with respect to comparisons with modeling data.

Corner reflectors and multi-temporal SAR inteferometry for landslide monitoring

The application of Persistent Scatterer Interferometry (PSI) to slope instability monitoring poses challenges related to the complex kinematics of the phenomenon, as well as to the unfavourable settings of the area affected by landslides, often occurring on sites of limited extension, characterized by steep topography and variable vegetation cover. New-generation SAR sensors, such as TerraSAR-X (TSX) thanks to their higher spatial resolution, make PSI applications very promising for monitoring areas with low density man-made. Nevertheless, the application of techniques still remains problematic or impossible in rural and mountainous areas. This is the case, for instance, for the Municipality of Carlantino, in Southern Italy. Both C-band medium resolution SAR data from ESA satellites, and X-band high resolution SAR data from the TSX satellite, were processed through the PSI algorithm SPINUA. Despite the higher spatial density of PS from TSX, the landslide body is lacking coherent targets, due to vegetation and variable land cover. To allow stability monitoring, a network of six CRs was designed and deployed over the landslide test site. Twenty-six TSX stripmap images were processed by using both PSI and an ad hoc procedure based on double-difference analysis of DInSAR phase values on the CR pixels, constrained by the accurate CR height measurements provided by DGPS. Despite the residual noise due to the sub-optimal CR network and the strong atmospheric signal, displacement estimation on the CRs allows to propagate the PSI results downslope, proving the stability of the landslide area subjected to consolidation works.

High resolution satellite multi-temporal interferometry for monitoring infrastructure instability hazards

Advanced remote sensing techniques are now capable of delivering more rapidly high quality information that is sufficiently detailed (and cost-effective) for many engineering applications. Here we focus on synthetic aperture radar (SAR), multi-temporal interferometry (MTI). With radar satellites periodically re-visiting the same area, MTI provides information on distance changes between the on-board radar sensor and the targets on the ground (e.g., human-made structures such as buildings, roads and other infrastructure). The detected distance changes are thus interpreted as evidence of ground and/or structure instability. In settings with limited vegetation cover, MTI can deliver very precise (mm resolution), spatially dense information (from hundreds to thousands measurement points/km2) on slow (mm-cm/year) deformations affecting the ground and engineering structures. Radar satellites offer wide-area coverage (thousands km2) and, with the sensors that actively emit electromagnetic radiation and thus can “see” through the clouds, one can obtain deformation measurements even under bad weather conditions. We illustrate the potential of high resolution MTI and explain what this technique can deliver to assist in infrastructure instability hazard assessment. This is done by presenting selected examples of MTI applied to monitor post-construction behavior of engineering structures. The examples are from Italy and include: an earthfill dam, an off-shore vertical breakwater built to protect an oil terminal, city buildings and a highway. We also stress that the current approach to the assessment of instability hazard can be transformed by capitalizing more on the presently underexploited advantage of the MTI technique, i.e., the capability to provide regularly spatially dense quantitative information for large areas where engineering infrastructure may currently be unaffected by instability, but where the terrain and infrastructure history (e.g., aging) may indicate potential for future failures.