Raffaele Nutricato

SAR and InSAR for Flood Monitoring: Examples With COSMO-SkyMed Data

We apply high-resolution, X-band, stripmap COSMO-SkyMed data to the monitoring of flood events in the Basilicata region (Southern Italy), where multitemporal datasets are available with short spatial and temporal baselines, allowing interferometric (InSAR) processing. We show how the use of the interferometric coherence information can help to detect more precisely the areas affected by the flood, reducing false alarms and missed identifications which affect algorithms based on SAR intensity alone. The effectiveness of using the additional InSAR information layer is illustrated by RGB composites of various combinations of intensity and coherence data. Analysis of multitemporal SAR intensity and coherence trends reveals complex behavior of various field types, which we interpret through a Bayesian inference approach, based on a manual identification of representative scattering and coherence signatures of selected homogeneous fields. The approach allows to integrate external, ancillary information to derive a posteriori probabilistic maps of flood inundation accounting for different scattering responses to the presence of water. First results of this semiautomated methodology, using simple assumptions for the SAR signatures and a priori information based on the distance from river courses, show encouraging results, and open a path to improvement through use of more complex hydrologic and topo-hydrographic information.

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.

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.

Coastal Observation through Cosmo-SkyMed High-Resolution SAR Images

ABSTRACT Bruno, M.F.; Molfetta M.G.; Mossa, M., Nutricato, R., Morea, A., and Chiaradia, M.T., 2016. Coastal observation through Cosmo-SkyMed high-resolution SAR images. In: Vila-Concejo, A.; Bruce, E.; Kennedy, D.M., and McCarroll, R.J. (eds.), Proceedings of the 14th International Coastal Symposium (Sydney, Australia). Journal of Coastal Research, Special Issue, No. 75, pp. 795–799. Coconut Creek (Florida), ISSN 0749-0208. The study deals with the application and further improvement of an advanced Earth Observation system, named COSMO-Beach, developed for semi-automatic shoreline extraction and coastal morphology identification. The system exploits SAR Single-Look-Complex data acquired by the COSMO-SkyMed constellation, which is able to provide X-band images with a short revisiting time. The implemented procedures have been tested over a very popular beach in Apulia Region (Italy), affected by erosion problems induced by human activities. The outcomes of the COSMO-Beach system are presented and discussed.

Persistent Scatterers Interferometry Provides Insight on Slope Deformations and Landslide Activity in the Mountains of Zhouqu, Gansu, China

We present the results of Persistent Scatterers Interferometry (PSI) applied to investigate slope instabilities in a remote high mountain region of Southern Gansu, known to be prone to large magnitude (M7-8) earthquakes and catastrophic slope failures. The PSI processing of high resolution (~3 m) COSMO/SkyMed (CSK) satellite images produced spatially dense information (more than 1,000 PS/km2) on ground surface displacements in the area of Zhouqu, a town located in the Bailong River valley. A substantial portion of the radar targets showed significant displacements (from few to over 50 mm/year), denoting widespread occurrence of slope instabilities. In particular, the PSI results provided valuable information on the activity of some very large, apparently slow landslides that represent a persistent hazard to the local population and infrastructure. Monitoring movements of large long-lived landslides is important especially when, as in the case of the Bailong valley, they are known to undergo periods of increased activity resulting in river damming and disastrous flooding. Given the general lack of monitoring data on large landslides at Zhouqu and on other similar major failures that are common in Southern Gansu, the PSI-derived displacements offer unique information, which, following expert judgment, can be used for preliminary wide-area assessments of hazards linked to landslide activity. Furthermore, this study shows that with the high resolution CSK data resulting in high radar target density, PSI can also assist in slope/landslide-specific assessments.

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.

Optimum interpolation and resampling for PSC identification

The first basic step in the permanent scatterers analysis technique is the location of the stable reflectors. For this purpose, an inter-image amplitude analysis, on a pixel-by-pixel basis, is performed, and the scatterers showing stable amplitude response are named Permanent Scatterer Candidates (PSC). The effects of the interpolation of the SLC (Single Look Complex) images are considered, with particular attention to the oversampling factor and to the interpolation kernels. Finally, experimental results are shown to confirm the PSC theory and the interpolation analysis.

Automated calibration of multi-temporal ERS SAR data

Quantitative analysis of multi-temporal SAR datasets requires accurate radiometric calibration. This can usually be achieved by considering a number of multiplicative factors, to be applied to statistically-homogeneous areas on each image to be calibrated. We propose an automated procedure to easily obtain a relative radiometric calibration of a stack of an arbitrary number of ERS SAR images, all coregistered to a unique master. The procedure relies upon the theoretical invariance of both system and scene parameters for targets whose amplitude values are statistically correlated with each other. In practice, with the above-mentioned assumptions, the only varying quantity to take into account is the temporal variation of the calibration factor. We show, based on experiments on two multi-temporal datasets, that the ERS-1 calibration constant exhibits practically no variations except for the expected statistical fluctuations of about /spl plusmn/1 dB. This observation is consistent with the reports on the operations of the ERS-2 satellite, published by ESA. The proposed procedure then consists in applying the full absolute calibration procedure only to the master image, and then relatively correct all other images in the stack by the method described.

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.