Flood mapping in vegetated areas using an unsupervised clustering approach on Sentinel-1 and-2 imagery

The European Space Agency's Sentinel-1 constellation provides timely and freely available dual-polarized C-band Synthetic Aperture Radar (SAR) imagery. The launch of these and other SAR sensors has boosted the field of SAR-based flood mapping. However, flood mapping in vegetated areas remains a topic under investigation, as backscatter is the result of a complex mixture of backscattering mechanisms and strongly depends on the wave and vegetation characteristics. In this paper, we present an unsupervised object-based clustering framework capable of mapping flooding in the presence and absence of flooded vegetation based on freely and globally available data only. Based on a SAR image pair, the region of interest is segmented into objects, which are converted to a SAR-optical feature space and clustered using K-means. These clusters are then classified based on automatically determined thresholds, and the resulting classification is refined by means of several region growing post-processing steps. The final outcome discriminates between dry land, permanent water, open flooding, and flooded vegetation. Forested areas, which might hide flooding, are indicated as well. The framework is presented based on four case studies, of which two contain flooded vegetation. For the optimal parameter combination, three-class F1 scores between 0.76 and 0.91 are obtained depending on the case, and the pixel- and object-based thresholding benchmarks are outperformed. Furthermore, this framework allows an easy integration of additional data sources when these become available

Evaluación de la degradación de la tierra usando la entropía de shannon sobre imágenes polarimétricas en desiertos costeros Patagónicos

En esta investigación se focalizó en la Entropía de Shannon (ES) para la caracterización de imágenes polarimétricas de apertura sintética. Este parámetro analiza la contribución de la información por pixeles individuales para toda la imagen en la evaluación de la degradación de la tierra en imágenes ALOS PALSAR. Escenas de polarización dual y cuádruple fueron adquiridas bajo el proyecto SAOCOM (Satélite Argentino de Observación con Microondas) en 2010 y 2011, del desierto costero noreste patagónico, Argentina. Los mapas fueron verificados con información de alta verosimilitud para la misma área de estudio. Los resultados muestran que la ES puede describir y precisar las características de las imágenes de manera obvia, de tal manera que representa un valor de referencia para la detección de la degradación de la tierra y la extracción de las características de los diferentes estados y transiciones.We focus on Shannon Entropy (SE) for the characterization of polarimetric Synthetic Aperture Radar (PolSAR) images. This approach analyzes the information contribution made by individual pixels to the whole image for assessment of land degradation in the information content of ALOS PALSAR images. Additionally, the performance of other polarization parameters, and polarization decomposition is illustrated and discussed. Dual-Pol and Quad-Pol scenes have been acquired under the SAOCOM (Satélite Argentino de Observación con Microondas, Spanish for Argentine Microwaves Observation Satellite) project in 2010 and 2011, from northeastern Patagonian coastal desert, Argentina. The accuracy of the SE map was assessed using a set of ground observations based on remotely sensed data that have higher accuracy. The results show that the SE can describe and determine the image features more obviously in the study area, so that it represents an important reference value for land degradation detection and land status characteristics extraction .Fil: del Valle, Hector Francisco. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Nacional Patagonico; ArgentinaFil: Hardtke, Leonardo Andrés. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Nacional Patagonico; ArgentinaFil: Blanco, Paula Daniela. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Nacional Patagonico; ArgentinaFil: Sione, Walter Fabian. Universidad Autónoma de Entre Ríos. Fac de Ciencia y Tecnologia. Centro Regional de Geomatica; Argentina. Universidad Nacional de Luján; Argentin

Detection of temporarily flooded vegetation using time series of dual polarised C-band synthetic aperture radar data

The intense research of the last decades in the field of flood monitoring has shown that microwave sensors provide valuable information about the spatial and temporal flood extent. The new generation of satellites, such as the Sentinel-1 (S-1) constellation, provide a unique, temporally high-resolution detection of the earth's surface and its environmental changes. This opens up new possibilities for accurate and rapid flood monitoring that can support operational applications. Due to the observation of the earth's surface from space, large-scale flood events and their spatiotemporal changes can be monitored. This requires the adaptation of existing or the development of new algorithms, which on the one hand enable precise and computationally efficient flood detection and on the other hand can process a large amounts of data. In order to capture the entire extent of the flood area, it is essential to detect temporary flooded vegetation (TFV) areas in addition to the open water areas. The disregard of temporary flooded vegetation areas can lead to severe underestimation of the extent and volume of the flood. Under certain system and environmental conditions, Synthetic Aperture Radar (SAR) can be utilized to extract information from under the vegetation cover. Due to multiple backscattering of the SAR signal between the water surface and the vegetation, the flooded vegetation areas are mostly characterized by increased backscatter values. Using this information in combination with a continuous monitoring of the earth's surface by the S-1 satellites, characteristic time series-based patterns for temporary flooded vegetation can be identified. This combination of information provides the foundation for the time series approach presented here. This work provides a comprehensive overview of the relevant sensor and environmental parameters and their impact on the SAR signal regarding temporary open water (TOW) and TFV areas. In addition, existing methods for the derivation of flooded vegetation are reviewed and their benefits, limitations, methodological trends and potential research needs for this area are identified and assessed. The focus of the work lies in the development of a SAR and time series-based approach for the improved extraction of flooded areas by the supplementation of TFV and on the provision of a precise and rapid method for the detection of the entire flood extent. The approach developed in this thesis allows for the precise extraction of large-scale flood areas using dual-polarized C-band time series data and additional information such as topography and urban areas. The time series features include the characteristic variations (decrease and/or increase of backscatter values) on the flood date for the flood-related classes compared to the whole time series. These features are generated individually for each available polarization (VV, VH) and their ratios (VV/VH, VV-VH, VV+VV). The generation of the time series features was performed by Z-transform for each image element, taking into account the backscatter values on the flood date and the mean value and standard deviation of the backscatter values from the nonflood dates. This allowed the comparison of backscatter intensity changes between the image elements. The time series features constitute the foundation for the hierarchical threshold method for deriving flood-related classes. Using the Random Forest algorithm, the importance of the time series data for the individual flood-related classes was analyzed and evaluated. The results showed that the dual-polarized time series features are particularly relevant for the derivation of TFV. However, this may differ depending on the vegetation type and other environmental conditions. The analyses based on S-1 data in Namibia, Greece/Turkey and China during large-scale floods show the effectiveness of the method presented here in terms of classification accuracy. Theiv supplementary integration of temporary flooded vegetation areas and the use of additional information resulted in a significant improvement in the detection of the entire flood extent. It could be shown that a comparably high classification accuracy (~ 80%) was achieved for the flood extent in each of study areas. The transferability of the approach due to the application of a single time series feature regarding the derivation of open water areas could be confirmed for all study areas. Considering the seasonal component by using time series data, the seasonal variability of the backscatter signal for vegetation can be detected. This allows for an improved differentiation between flooded and non-flooded vegetation areas. Simultaneously, changes in the backscatter signal can be assigned to changes in the environmental conditions, since on the one hand a time series of the same image element is considered and on the other hand the sensor parameters do not change due to the same acquisition geometry. Overall, the proposed time series approach allows for a considerable improvement in the derivation of the entire flood extent by supplementing the TOW areas with the TFV areas

Integration of L-band SAR data into Land Surface Models

Abstract—Land surface process modelling might be limited due to lack of reliable model input data. Key surface variables as land cover information or soil moisture conditions have been proven to be observable by remote sensing systems. The integration of remote sensing data into land surface process models might therefore help to improve their simulations results. Longer wavelength SAR data has a higher sensitivity to soil moisture content than higher frequency systems. Recent (ALOS) and planed (e.g. TerraSAR-L) SAR systems are therefore expected to provide valuable information about soil moisture dynamics. The present study investigates the potential to retrieve land cover information and geophysical parameters from L-band SAR data. The retrieval results are assimilated into a state-of-the-art land surface model to evaluate the merit of L-band SAR data assimilation

Modifying the Yamaguchi Four-Component Decomposition Scattering Powers Using a Stochastic Distance

Model-based decompositions have gained considerable attention after the initial work of Freeman and Durden. This decomposition which assumes the target to be reflection symmetric was later relaxed in the Yamaguchi et al. decomposition with the addition of the helix parameter. Since then many decomposition have been proposed where either the scattering model was modified to fit the data or the coherency matrix representing the second order statistics of the full polarimetric data is rotated to fit the scattering model. In this paper we propose to modify the Yamaguchi four-component decomposition (Y4O) scattering powers using the concept of statistical information theory for matrices. In order to achieve this modification we propose a method to estimate the polarization orientation angle (OA) from full-polarimetric SAR images using the Hellinger distance. In this method, the OA is estimated by maximizing the Hellinger distance between the un-rotated and the rotated $T_{33}$ and the $T_{22}$ components of the coherency matrix $\mathbf{[T]}$. Then, the powers of the Yamaguchi four-component model-based decomposition (Y4O) are modified using the maximum relative stochastic distance between the $T_{33}$ and the $T_{22}$ components of the coherency matrix at the estimated OA. The results show that the overall double-bounce powers over rotated urban areas have significantly improved with the reduction of volume powers. The percentage of pixels with negative powers have also decreased from the Y4O decomposition. The proposed method is both qualitatively and quantitatively compared with the results obtained from the Y4O and the Y4R decompositions for a Radarsat-2 C-band San-Francisco dataset and an UAVSAR L-band Hayward dataset.Comment: Accepted for publication in IEEE J-STARS (IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing

SAR data and field surveys combination to update rainfall-induced shallow landslide inventory

The Campania region has been recurrently hit by severe landslides in volcanoclastic deposits. The city of Naples, and in particular the Camaldoli and Agnano hills (Phlegraean Fields), also suffered several landslide crises in weathered volcanoclastic rocks as a consequence of intense rainfalls or wildfires. To identify slope failures phenomena occurred in the winter season 2019–2020 an innovative procedure has been proposed. The purpose of this procedure is to highlight areas where major land cover changes occurred within our area of study, which can be potentially related to mass movements. The amplitude of spaceborne SAR images has been exploited for the change detection analysis and the output derived from the segmentation procedure has been compared with field observations. The amplitude-based method has been already applied in the detection of landslides, but never on the event with limited extensions, such as for this application. The achieved outcomes allowed the mapping of 62 new landslides that have been used to update the current landslide inventory database. This type of information is expected to help decision-makers with land planning and risk assessment

Fusion of VNIR Optical and C-Band Polarimetric SAR Satellite Data for Accurate Detection of Temporal Changes in Vegetated Areas

In this paper, we propose a processing chain jointly employing Sentinel-1 and Sentinel-2 data, aiming to monitor changes in the status of the vegetation cover by integrating the four 10 m visible and near-infrared (VNIR) bands with the three red-edge (RE) bands of Sentinel-2. The latter approximately span the gap between red and NIR bands (700 nm–800 nm), with bandwidths of 15/20 nm and 20 m pixel spacing. The RE bands are sharpened to 10 m, following the hypersharpening protocol, which holds, unlike pansharpening, when the sharpening band is not unique. The resulting 10 m fusion product may be integrated with polarimetric features calculated from the Interferometric Wide (IW) Ground Range Detected (GRD) product of Sentinel-1, available at 10 m pixel spacing, before the fused data are analyzed for change detection. A key point of the proposed scheme is that the fusion of optical and synthetic aperture radar (SAR) data is accomplished at level of change, through modulation of the optical change feature, namely the difference in normalized area over (reflectance) curve (NAOC), calculated from the sharpened RE bands, by the polarimetric SAR change feature, achieved as the temporal ratio of polarimetric features, where the latter is the pixel ratio between the co-polar and the cross-polar channels. Hyper-sharpening of Sentinel-2 RE bands, calculation of NAOC and modulation-based integration of Sentinel-1 polarimetric change features are applied to multitemporal datasets acquired before and after a fire event, over Mount Serra, in Italy. The optical change feature captures variations in the content of chlorophyll. The polarimetric SAR temporal change feature describes depolarization effects and changes in volumetric scattering of canopies. Their fusion shows an increased ability to highlight changes in vegetation status. In a performance comparison achieved by means of receiver operating characteristic (ROC) curves, the proposed change feature-based fusion approach surpasses a traditional area-based approach and the normalized burned ratio (NBR) index, which is widespread in the detection of burnt vegetation