Novel clustering schemes for full and compact polarimetric SAR data: An application for rice phenology characterization

Information on rice phenological stages from Synthetic Aperture Radar (SAR) images is of prime interest for in-season monitoring. Often, prior in-situ measurements of phenology are not available. In such situations, unsupervised clustering of SAR images might help in discriminating phenological stages of a crop throughout its growing period. Among the existing unsupervised clustering techniques using full-polarimetric (FP) SAR images, the eigenvalue-eigenvector based roll-invariant scattering-type parameter, and the scattering entropy parameter are widely used in the literature. In this study, we utilize a unique target scattering-type parameter, which jointly uses the Barakat degree of polarization and the elements of the polarimetric coherency matrix. Likewise, we also utilize an equivalent parameter proposed for compact-polarimetric (CP) SAR data. These scattering-type parameters are analogous to the Cloude-Pottier’s parameter for FP SAR data and the ellipticity parameter for CP SAR data. Besides this, we also introduce new clustering schemes for both FP and CP SAR data for segmenting diverse scattering mechanisms across the phenological stages of rice. In this study, we use the RADARSAT-2 FP and simulated CP SAR data acquired over the Indian test site of Vijayawada under the Joint Experiment for Crop Assessment and Monitoring (JECAM) initiative. The temporal analysis of the scattering-type parameters and the new clustering schemes help us to investigate detailed scattering characteristics from rice across its phenological stages.This work was supported in part by the Spanish Ministry of Science, Innovation and Universities, the State Agency of Research (AEI), and the European Funds for Regional Development (EFRD) under Project TEC 2017-85244-C 2-1-P. The work of Dipankar Mandal was supported by the Ministry of Human Resource Development, Government of India (New Delhi, India) towards his Ph.D. assistantship through grant no. RSPHD0210

Study of the speckle noise effects over the eigen decomposition of polarimetric SAR data: a review

This paper is focused on considering the effects of speckle noise on the eigen decomposition of the co- herency matrix. Based on a perturbation analysis of the matrix, it is possible to obtain an analytical expression for the mean value of the eigenvalues and the eigenvectors, as well as for the Entropy, the Anisotroopy and the dif- ferent a angles. The analytical expressions are compared against simulated polarimetric SAR data, demonstrating the correctness of the different expressions.Peer ReviewedPostprint (published version

Polarimetric Synthetic Aperture Radar

This open access book focuses on the practical application of electromagnetic polarimetry principles in Earth remote sensing with an educational purpose. In the last decade, the operations from fully polarimetric synthetic aperture radar such as the Japanese ALOS/PalSAR, the Canadian Radarsat-2 and the German TerraSAR-X and their easy data access for scientific use have developed further the research and data applications at L,C and X band. As a consequence, the wider distribution of polarimetric data sets across the remote sensing community boosted activity and development in polarimetric SAR applications, also in view of future missions. Numerous experiments with real data from spaceborne platforms are shown, with the aim of giving an up-to-date and complete treatment of the unique benefits of fully polarimetric synthetic aperture radar data in five different domains: forest, agriculture, cryosphere, urban and oceans

Electromagnetic backscatter modelling of icebergs at c-band in an ocean environment

This thesis outlines the development of an electromagnetic (EM) backscatter model of icebergs. It is a necessary first step for the generation of in-house synthetic aperture radar (SAR) data of icebergs to support optimum iceberg/ship classifier design. The EM modelling was developed in three stages. At first, an EM backscatter model was developed to generate simulated SAR data chips of iceberg targets at small incidence angles. The model parameters were set to mimic a dual polarized dataset collected at C-Band with the Sentinel-1A satellite. The simulated SAR data chips were compared with signatures and radiometric properties of the satellite data, including total radar cross section (TRCS). A second EM model was developed to mimic the parameters of a second SAR data collection with RADARSAT-2; this second data collection was at larger incidence angles and was fully polarimetric (four channels and interchannel phase). The full polarimetric SAR data allowed for a comparison of modelled TRCS and polarimetric decompositions. Finally, the EM backscatter models were tested in the context of iceberg/ship classification by comparing the performance of various computer vision classifiers using both simulated and real SAR image data of iceberg and vessel targets. This step is critical to check the compatibility of simulated data with the real data, and the ability to mix real and simulated SAR imagery for the generation of skilled classifiers. An EM backscatter modelling tool called GRECOSAR was used for the modelling work. GRECOSAR includes the ability to generate small scenes of the ocean using Pierson-Moskowitz spectral parameters. It also allows the placement of a 3D target shape into that ocean scene. Therefore, GRECOSAR is very useful for simulating SAR targets, however it can only model single layer scattering from the targets. This was found to be limiting in that EM scattering throughout volume of the iceberg could not be generated. This resulted in EM models that included only surface scattering of the iceberg. In order to generate realistic SAR scenes of icebergs on the ocean, 3D models of icebergs were captured in a series of field programs off the coast of Newfoundland and Labrador, Canada. The 3D models of the icebergs were obtained using a light detection and ranging (LiDAR) and multi-beam sonar data from a specially equipped vessel by a team of C-CORE. While profiling the iceberg targets, SAR images from satellites were captured for comparison with the simulated SAR images. The analysis of the real and simulated SAR imagery included comparisons of TRCS, SAR signature morphology and polarimetric decompositions of the targets. In general, these comparisons showed a good consistency between the simulated and real SAR scene. Simulations were also performed with varying target orientation and sea conditions (i.e., wind speed and direction). A wide variability of the TRCS and SAR signature morphology was observed with varying scene parameters. Icebergs were modelled using a high dielectric constant to mimic melting iceberg surfaces as seen during field work. Given that GRECOSAR could only generate surface backscatter, a mathematical model was developed to quantify the effect of melt water on the amount of surface and volume backscatter that might be expected from the icebergs. It was found that the icebergs in a high state of melt should produce predominantly surface scatter, thus validating the use of GRECOSAR for icebergs in this condition. Once the simulated SAR targets were validated against the real SAR data collections, a large dataset of simulated SAR chips of ships and icebergs were created specifically for the purpose of target classification. SAR chips were generated at varying imaging parameters and target sizes and passed on to an iceberg/ship classifier. Real and simulated SAR chips were combined in varying quantities (or targets) resulting in a series of different classifiers of varying skill. A good agreement between the classifier’s performance was found. This indicates the compatibility of the simulated SAR imagery with this application and provides an indication that the simulated data set captures all the necessary physical properties of icebergs for ship and iceberg classification

Monitoring Soil Moisture and Freeze/Thaw State Using C-band Imaging Radar

Soil moisture is an important state variable in many hydrological and meteorological applications. This thesis explores the use of the C-band synthetic aperture radar (SAR) parameters to monitor soil moisture and freeze/thaw state in a cold-season hydrologic environment. The circular-linear compact polarimetric (CP) configuration is considered as a possible alternative of the quad polarimetric (QP) system because it acquires images with wider swath and reduced complexity, cost and energy requirement of the radar system while maintaining the information content of the acquired imagery. In this study, 15 RADARSAT-2 QP images were acquired from October 2013 to June 2014 and CP images were simulated from each RADARSAT-2 QP imagery acquired. Field measurements of soil properties were collected along with the radar imagery acquisitions. The backscattering coefficients in all polarizations were able to discriminate frozen and unfrozen soils. But their correlations with soil moisture content were weak if examining frozen or unfrozen soils separately. The Oh et al. (1992) model was implemented in this study to compare with acquired RADARSAT-2 data. A good agreement was found between the linear polarimetric backscattering coefficients simulated by the Oh model and the RADARSAT-2 data, indicating that the study site even with 10 cm tall standing hay was consistent with a bare soil site at C-band and the Oh model can be applied to frozen soils. With respect to CP parameters, the first and fourth Stokes parameters and m-δ surface and volume scattering components can detect soil freeze/thaw state and have potential for frozen/unfrozen soils mapping. The influence of vegetation on selected CP parameters was also evident in this study. Results demonstrated the utility of C-band radar in detecting soil freeze/thaw state rather than monitoring the changes of soil moisture content. More image acquisitions during the freezing and thawing periods, continuous field measurements of soil moisture and state, and ground measurements collected over wider study area can help further develop understanding of the CP parameters and facilitate future use of the CP mode. The contribution of this thesis is to provide better understanding of the CP parameters at a specific site and to demonstrate that CP parameters can replicate QP SAR variables to detect surface soil conditions

Oil Spills and Slicks imaged by Synthetic Aperture Radar

Oil spills and slicks occur in the ocean around the world due to natural seeps, oil extraction, transportation, and consumption. Satellite synthetic aperture radar (SAR) has proven to be an efficient tool for identifying and classifying oil on the sea surface. This information can be used to monitor areas for potential illegal marine discharge or to respond to an oil spill incident. When used to monitor shipping lanes or drilling platforms, timely analysis can identify offending parties and lead to prosecution. Following an oil spill such as that from the Deepwater Horizon rig in the Gulf of Mexico in 2010, SAR can be used to direct response activities and optimize available resources

Microwave Satellite Measurements for Coastal Area and Extreme Weather Monitoring

In this project report, the main outcomes relevant to the Sino-European Dragon-4 cooperation project ID 32235 “Microwave satellite measurements for coastal area and extreme weather monitoring” are reported. The project aimed at strengthening the Sino-European research cooperation in the exploitation of European Space Agency, Chinese and third-party mission Earth Observation (EO) microwave satellite data. The latter were exploited to perform an effective monitoring of coastal areas, even under extreme weather conditions. An integrated multifrequency/polarization approach based on complementary microwave sensors (e.g., Synthetic Aperture Radar, scatterometer, radiometer), together with ancillary information coming from independent sources, i.e., optical imagery, numerical simulations and ground measurements, was designed. In this framework, several tasks were addressed including marine target detection, sea pollution, sea surface wind estimation and coastline extraction/classification. The main outcomes are both theoretical (i.e., new models and algorithms were developed) and applicative (i.e., user-friendly maps were provided to the end-user community of coastal area management according to smart processing of remotely sensed data). The scientific relevance consists in the development of new algorithms, the effectiveness and robustness of which were verified on actual microwave measurements, and the improvement of existing methodologies to deal with challenging test cases

Remote sensing based assessment of land cover and soil moisture in the Kilombero floodplain in Tanzania

Wetlands provide important ecological, biological, and social-economic services that are critical for human existence. The increasing demand for food, arable land shortage and changing climate conditions in East Africa have created a paradigm shift from upland cultivation to wetland use due to their year-round soil water availability. However, there is need to control and manage the activities within the wetlands to ensure sustainable use while negating any negative effects caused by these activities. This is implemented through the decisions made by the land managers within the wetlands. Providing the users of the wetlands with scientific knowledge acts as a support tool for policy-making geared towards the sustainable use of the wetlands. The overall research contains two main components: First, the need for timely land cover maps at a reasonable scale, and secondly, the assessment of soil moisture as a major contributor to agricultural production. The objectives of the study were to generate land cover maps from multi-sensor optical datasets and to assess the performance of single-polarized Sentinel-1 Gray Level Co-occurrence Matrix (GLCM) texture and Principal Component Analysis (PCA) features by applying multiple classification algorithms in a floodplain in the Kilombero catchment. Furthermore, soil moisture spatial-temporal patterns over three hydrological zones was assessed, estimation of soil moisture from radar data and generation of soil moisture products from global products was investigated. The correlation of the merged products to Normalized Difference Vegetation Index (NDVI) measures was also investigated. RapidEye, Sentinel-2 and Landsat images were used in determining the areal extents of four major land cover classes namely vegetated, bare, water and built up. The acquisition period of the images ranges from August 2013 to June 2015 for the RapidEye images, December 2015 to August 2016 for the Sentinel-2 images and 2013 to 2016 Landsat-8 images were included in the land cover time series dynamic study. However, the major challenge arising was cloud coverage and hence Sentinel-1 images were tested in the application of Synthetic Aperture Radar (SAR) in wetland mapping. Variograms were used in spatial-temporal assessment of soil moisture data collected from three hydrological zones, riparian, middle and fringe. A roughness parameter was derived from a semi-empirical model. Soil moisture was retrieved from TerraSAR-X and RadarSAT-2 with the retrieved roughness parameter as an input in a linear regression equation. Triple collocation was applied in error assessment of the global soil moisture products prior to development of a merged product. Cross-correlation was applied in relating NDVI to soil moisture. Optical data (RapidEye, Landsat-8, and Sentinel-2) generated land cover maps used in assessing the land cover dynamics over time. The land cover ratios were related to depth to groundwater. As the depth to groundwater reduced in June the bare land coverage was 45-57% while that of vegetation was 34-47%. In December when the depth to groundwater was highest, bare land coverage was 62-69% while that of the vegetated area was 27-25%. This indicates that depth of groundwater and vegetation coverage responds to seasonality. During the dry season, 68-81% of the total vegetation class is within the riparian zone. In the classification of the SAR images, the overall accuracies for the single polarized VV images ranged from 54-76%, 60-81% and 61-80% for Random Forest (RF), Neural Network (NN) and Support Vector Machine (SVM) respectively. GLCM features had overall accuracies of 64-86%, 65-88% and 65-86% for RF, NN, and SVM respectively. PCA derived images had similar overall accuracies of 68-92% for NN, RF, and SVM respectively. The PCA images had the highest overall accuracy for the entire time series indicating that reduction in the number of texture features to layers containing the maximum variance improves the accuracy. The standard deviation of soil moisture was noted to increase with increasing soil moisture. Soil texture plays a key role in soil moisture retention. The riparian fields had a high water content explained by the high clay and organic matter content. A roughness parameter was derived and utilized in the retrieval of soil moisture from SAR resulting to R2 of 0.88- 0.92 between observed and simulated soil moisture values from co-polarized RadarSAT-2 HH and TerraSAR-X HH and VV. Merged soil moisture product from FEWSNET Land Data Assimilation System_NOAH (FLDAS_NOAH), ECMWF Re-Analysis Interim (ERA-Interim) and Soil Moisture and Ocean Salinity (SMOS) and FLDAS_Variable Infiltration Capacity (VIC), ERA-Interim and SMOS had similar patterns attributed to FLDAS_NOAH and FLDAS_VIC forced by the same precipitation product (RFE). Cross-correlation of Moderate-resolution Imaging Spectrometer (MODIS) NDVI and the merged soil moisture products revealed a 2-month lag of NDVI. Hence, the relationship is useful in determining the Start of Season from soil moisture products. In conclusion, the successful land cover mapping of the study area demonstrated the use of satellite imagery for wetland characterization. The vast coverage and frequent acquisitions of optical and microwave remotely sensed data additionally make the approaches transferable to other locations and allow for mapping at larger scales. Soil moisture assessment from point data revealed varied soil moisture patterns whereas global remotely sensed and modeled products rather provide complementary information about growing conditions, and hence a situational assessment tool of potential of physical availability dimension of food security. This study forms a baseline upon which additional monitoring and assessment of the Kilombero wetland ecosystem can be performed with the current results marked as a reference. Moreover, the study serves as a demonstration case of remote sensing based approaches for land cover and soil moisture mapping, whose results are useful to stakeholders to aid in the implementation of adapted production techniques for yield optimization while minimizing the unsustainable use of the natural resources.Feuchtgebiete erbringen wichtige ökologische, biologische und sozial-ökonomische Dienstleistungen, welche entscheidend für das menschliche Dasein sind. Der steigende Bedarf an Nahrung, der Mangel an landwirtschaftlichen Nutzflächen und die Veränderung der klimatischen Bedingungen in Ostafrika haben zu einem Paradigmenwechsel vom Anbau im Hochland hin zur Nutzung von Feuchtgebieten geführt. Allerdings sind Kontrolle und Management der Aktivitäten in Feuchtgebieten notwendig, um die nachhaltige Nutzung zu sichern und negative Effekte dieser Aktivitäten zu vermeiden. Die Implementierung erfolgt durch die Landverwalter in den Feuchtgebieten. Den Nutzern von Feuchtgebieten wissenschaftliche Erkenntnisse bereitzustellen dient als Hilfsmittel zur politischen Entscheidungsfindung für die nachhaltige Feuchtgebietsnutzung. Die Forschung im Rahmen der Dissertation beinhaltet zwei Hauptkomponenten: erstens den Bedarf an aktuellen Landbedeckungskarten auf einer angemessenen Skalenebene und zweitens die Erfassung der Bodenfeuchte als wichtiger Einflussfaktor auf die landwirtschaftliche Produktion. Das Ziel der Untersuchung war, Landbedeckungskarten auf Grundlage von multisensorischen optischen Daten zu erstellen und die Eignung der Textur der einfach polarisierten Sentinel-1 Grauwertmatrix (GLCM) sowie der einer Hauptkomponentenanalyse (PCA) bei Anwendung unterschiedlicher Klassifikationsalgorithmen zu beurteilen. Des Weiteren wurden raum-zeitliche Bodenfeuchtemuster über drei hydrologische Zonen hinweg modelliert, die Bodenfeuchte aus Radardaten abgeleitet sowie die Erstellung von Bodenfeuchteprodukten auf Basis von globalen Produkten untersucht. Die Korrelation der Bodenfeuchteprodukte mit dem Normalisierten Differenzierten Vegetationsindex (NDVI) wurde ebenfalls analysiert. RapidEye, Sentinel-2 und Landsat Bilder wurden genutzt um die räumliche Ausdehnung der vier Hauptklassen (Vegetation, freiliegender Boden, Wasser und Bebauung) der Landbedeckung zu ermitteln. Für die Zeitreihenanalyse der der Landbedeckungsdynamik wurden RapidEye-Daten von August 2013 bis Juni 2015, Sentinel-2-Bilder von Dezember 2015 bis August 2016 und Landsat-8-Bilder von 2013 bis 2016 verwendet. Die größte Herausforderung war jedoch die Wolkenbedeckung, weshalb die Anwendung von Synthetic Aperture Radar (SAR) für die Feuchtgebietskartierung getestet wurde. Die gemessene Bodenfeuchte wurde mittels Variogrammen für die drei hydrologischen Zonen (Uferzone, Mitte und Randgebiete) raum-zeitlich interpoliert. Ein Rauhigkeitsparameter wurde aus einem semi-empirischen Modell hergeleitet. Die Bodenfeuchte wurde aus TerraSAR-X und RadarSAT-2- Bildern unter Verwendung des Rauhigkeitsparameters als Eingangsgröße in einer linearen Regression abgeleitet. Vor der Zusammenführung der Produkte wurde das globale Bodenfeuchteprodukt mithilfe von dreifacher Kollokation auf Fehler überprüft. Die Kreuzkorrelation zwischen NDVI und Bodenfeuchte wurde berechnet. Optische Daten (RapidEye, Landsat-8 und Sentinel-2) wurden genutzt, um die zeitliche Dynamik der Landbedeckung zu bestimmen. Die Landbedeckungsverhältnisse wurde mit der Höhe des Grundwasserspiegels korreliert. Ein hoher Grundwasserstand im Juni resultierte in 45-57% unbedecktem Boden, während der Anteil der Vegetation 34-47% betrug. Im Dezember, als der Grundwasserspiegel seinen Tiefststand hatte, erhöhte sich der Anteil des freiliegenden Bodens auf 62-69% und der Anteil der Vegetation verringerte sich auf 27-25%. Das zeigt, dass Grundwasserspiegel und Vegetation saisonalen Schwankungen unterworfen sind. Während der Trockenzeit liegen 68-81% der gesamten als Vegetation klassifizierten Fläche innerhalb der Uferzone. In der Klassifikation der SAR-Bilder liegt die Gesamtgenauigkeit der einfach polarisierten VV-Bilder im Rahmen von 54-76%, 60-81% und 61-80%, entsprechend für Random Forest (RF), Neuronale Netze (NN) und Support Vector Machine (SVM). Die GLCM ergab eine Gesamtgenauigkeit von 64-86%, 65-88% und 65-86% für RF, NN und SVM. Die über eine PCA abgeleiteten Bilder erreichten eine ähnliche Genauigkeit von 68-92% für NN, RF und SVM. Die PCA-Bilder weisen die höchste Gesamtgenauigkeit der gesamten Zeitreihe auf, was darauf hinweist, dass eine Reduktion von Textureigenschaften auf Layer der maximalen Varianz enthalten, die Genauigkeit erhöht. Die Standardabweichung der Bodenfeuchte stieg mit zunehmender Bodenfeuchte. Die Bodentextur spielt dabei eine Schlüsselrolle für das Wasserhaltevermögen des Bodens. Die Uferzone wies einen hohen Wassergehalt auf, was durch den hohen Anteil von Ton und Humus zu erklären ist. Die beobachteten und simulierten Bodenfeuchtewerte von co-polarisierten RadarSAT-2 HH, TerraSAR-X HH und VV Daten korrelieren mit einem R2 von 0.88 - 0.92. Die zusammengesetzten globalen Bodenfeuchteprodukte von FLDAS_NOAH, ERA-Interim sowie SMOS und FLDAS_VIC, ERA-Interim und SMOS zeigen ähnliche Muster wie FLDAS_NOAH und FLDAS_VIC, was über die Verwendung desselben Niederschlagsproduktes (RFE) zu erklären ist. Die Kreuzkorrelation von MODIS NDVI und den zusammengeführten Bodenfeuchteprodukten ergab eine zeitliche Verzögerung des NDVI von zwei Monaten. Dieser Zusammenhang kann daher bei der Bestimmung des Saisonbeginns aus Bodenfeuchtigkeitsprodukten nützlich sein. Zusammengefasst hat die Studie gezeigt, wie Satellitenbilder zur Charakterisierung von Wetlands genutzt werden können. Die große Abdeckung und häufige Aufnahme der optischen und Mikrowellen-Fernerkundungsdaten ermöglichen darüber hinaus die Übertragung der Ansätze auf weitere Gebiete und Kartierung auf größeren Skalen. Die Punktmessungen zeigen kleinräumige Muster der Bodenfeuchte, während globale Fernerkundungsprodukte und Modelle Informationen über die Wachstumsbedingungen liefern und somit ein Bewertungsinstrument der Ernährungssicherheit darstellen können. Weiterhin bildet die Studie eine Basis, auf der ein weitergehendes Monitoring und eine Bewertung des Feuchtgebietsökosystems durchgeführt werden kann. Sie ist ein Beispiel für fernerkundungsbasierte Ansätze zur Landbedeckungs- und Bodenfeuchtekartierung; ihre Ergebnisse sind nützlich, um Akteuren bei der Implementierung von Produktionstechniken zu unterstützen, welche die Erträge maximieren und gleichzeitig die nicht nachhaltige Nutzung der natürlichen Ressourcen minimieren