Accurate flood mapping is important for both planning activities during emergencies and as a support for the successive assessment of damaged areas. A valuable information source for such a procedure can be remote sensing synthetic aperture radar (SAR) imagery. However, flood scenarios are typical examples of complex situations in which different factors have to be considered to provide accurate and robust interpretation of the situation on the ground. For this reason, a data fusion approach of remote sensing data with ancillary information can be particularly useful. In this paper, a Bayesian network is proposed to integrate remotely sensed data, such as multitemporal SAR intensity images and interferometric-SAR coherence data, with geomorphic and other ground information. The methodology is tested on a case study regarding a flood that occurred in the Basilicata region (Italy) on December 2013, monitored using a time series of COSMO-SkyMed data. It is shown that the synergetic use of different information layers can help to detect more precisely the areas affected by the flood, reducing false alarms and missed identifications which may affect algorithms based on data from a single source. The produced flood maps are compared to data obtained independently from the analysis of optical images; the comparison indicates that the proposed methodology is able to reliably follow the temporal evolution of the phenomenon, assigning high probability to areas most likely to be flooded, in spite of their heterogeneous temporal SAR/InSAR signatures, reaching accuracies of up to 89%.

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.

Flood area detection from multipass Synthetic Aperture Radar (SAR) data can be performed via amplitude change detection techniques. These methods allow flooded zones to be discriminated only when they are flooded at the time of the second passage, and not at the time of the first one. Coherence derived from multipass SAR interferometry can be used instead, as an indicator of changes in the electromagnetic scattering behaviour of the surface, thus potentially revealing all the areas affected by the flood event at any time between the two passes. The paper presents a prototype application of such techniques, that is, a flood map obtained from ERS-1/2 data taken over Beziers (Southern France), through proper thresholding of a combination of amplitude and coherence information. Produced in the framework of an ESA project, the map consists of a DXF vector file which can be imported directly into most commercial GIS software.

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.