On the use of the l(2)-norm for texture analysis of polarimetric SAR data

In this paper, the use of the l2-norm, or Span, of the scattering vectors is suggested for texture analysis of polarimetric synthetic aperture radar (SAR) data, with the benefits that we need neither an analysis of the polarimetric channels separately nor a filtering of the data to analyze the statistics. Based on the product model, the distribution of the l2-norm is studied. Closed expressions of the probability density functions under the assumptions of several texture distributions are provided. To utilize the statistical properties of the l2-norm, quantities including normalized moments and log-cumulants are derived, along with corresponding estimators and estimation variances. Results on both simulated and real SAR data show that the use of statistics based on the l2-norm brings advantages in several aspects with respect to the normalized intensity moments and matrix variate log-cumulants.Peer ReviewedPostprint (published version

CFAR Ship Detection in Polarimetric Synthetic Aperture Radar Images Based on Whitening Filter

Polarimetric whitening filter (PWF) can be used to filter polarimetric synthetic aperture radar (PolSAR) images to improve the contrast between ships and sea clutter background. For this reason, the output of the filter can be used to detect ships. This paper deals with the setting of the threshold over PolSAR images filtered by the PWF. Two parameter-constant false alarm rate (2P-CFAR) is a common detection method used on whitened polarimetric images. It assumes that the probability density function (PDF) of the filtered image intensity is characterized by a log-normal distribution. However, this assumption does not always hold. In this paper, we propose a systemic analytical framework for CFAR algorithms based on PWF or multi-look PWF (MPWF). The framework covers the entire log-cumulants space in terms of the textural distributions in the product model, including the constant, gamma, inverse gamma, Fisher, beta, inverse beta, and generalized gamma distributions (GΓDs). We derive the analytical forms of the PDF for each of the textural distributions and the probability of false alarm (PFA). Finally, the threshold is derived by fixing the false alarm rate (FAR). Experimental results using both the simulated and real data demonstrate that the derived expressions and CFAR algorithms are valid and robust

Recent Advancement of Synthetic Aperture Radar (SAR) Systems and Their Applications to Crop Growth Monitoring

Synthetic aperture radars (SARs) propagate and measure the scattering of energy at microwave frequencies. These wavelengths are sensitive to the dielectric properties and structural characteristics of targets, and less affected by weather conditions than sensors that operate in optical wavelengths. Given these advantages, SARs are appealing for use in operational crop growth monitoring. Engineering advancements in SAR technologies, new processing algorithms, and the availability of open-access SAR data, have led to the recent acceleration in the uptake of this technology to map and monitor Earth systems. The exploitation of SAR is now demonstrated in a wide range of operational land applications, including the mapping and monitoring of agricultural ecosystems. This chapter provides an overview of—(1) recent advancements in SAR systems; (2) a summary of SAR information sources, followed by the applications in crop monitoring including crop classification, crop parameter estimation, and change detection; and (3) summary and perspectives for future application development

Circular polarization measurement in millimeter-wavelength spectral-line VLBI observations

This paper considers the problem of accurate measurement of circular polarization in imaging spectral-line VLBI observations in the lambda=7 mm and lambda=3 mm wavelength bands. This capability is especially valuable for the full observational study of compact, polarized SiO maser components in the near-circumstellar environment of late-type, evolved stars. Circular VLBI polarimetry provides important constraints on SiO maser astrophysics, including the theory of polarized maser emission transport, and on the strength and distribution of the stellar magnetic field and its dynamical role in this critical circumstellar region. We perform an analysis here of the data model containing the instrumental factors that limit the accuracy of circular polarization measurements in such observations, and present a corresponding data reduction algorithm for their correction. The algorithm is an enhancement of existing spectral line VLBI polarimetry methods using autocorrelation data for calibration, but with innovations in bandpass determination, autocorrelation polarization self-calibration, and general optimizations for the case of low SNR, as applicable at these wavelengths. We present an example data reduction at $\lambda=7$ mm and derive an estimate of the predicted accuracy of the method of m_c < 0.5% or better at lambda=7 mm and m_c < 0.5-1% or better at lambda=3 mm. Both the strengths and weaknesses of the proposed algorithm are discussed, along with suggestions for future work.Comment: 23 pages, 13 figure

Robust CFAR Detector Based on Truncated Statistics for Polarimetric Synthetic Aperture Radar

Constant false alarm rate (CFAR) algorithms using a local training window are widely used for ship detection with synthetic aperture radar (SAR) imagery. However, when the density of the targets is high, such as in busy shipping lines and crowded harbors, the background statistics may be contaminated by the presence of nearby targets in the training window. Recently, a robust CFAR detector based on truncated statistics (TS) was proposed. However, the truncation of data in the format of polarimetric covariance matrices is much more complicated with respect to the truncation of intensity (single polarization) data. In this article, a polarimetric whitening filter TS CFAR (PWF-TS-CFAR) is proposed to estimate the background parameters accurately in the contaminated sea clutter for PolSAR imagery. The CFAR detector uses a polarimetric whitening filter (PWF) to turn the multidimensional problem to a 1-D case. It uses truncation to exclude possible statistically interfering outliers and uses TS to model the remaining background samples. The algorithm does not require prior knowledge of the interfering targets, and it is performed iteratively and adaptively to derive better estimates of the polarimetric covariance matrix (although this is computationally expensive). The PWF-TS-CFAR detector provides accurate background clutter modeling, a stable false alarm property, and improves the detection performance in high-target-density situations. RadarSat2 data are used to verify our derivations, and the results are in line with the theory

Method for landslides detection with semi-automatic procedures: The case in the zone center-east of Cauca department, Colombia

Landslides are a common natural hazard that causes human casualties, but also infrastructure damage and land-use degradation. Therefore, a quantitative assessment of their presence is required by means of detecting and recognizing the potentially unstable areas. This research aims to develop a method supported on semiautomatic methods to detect potential mass movements at a regional scale. Five techniques were studied: Morphometry, SAR interferometry (InSAR), Persistent Scatterer InSAR (PS-InSAR), SAR polarimetry (PolSAR) and NDVI composites of Landsat 5, Landsat 7, and Landsat 8. The case study was chosen within the mid-eastern area of the Cauca state, which is characterised by its mountainous terrain and the presence of slope instabilities, officially registered in the CGS-SIMMA landslide inventory. This inventory revealed that the type `slide' occurred with 77.4% from the entire registries, `fall' with 16.5%, followed by `creeps' with 3%, flows with 2.6%, and `lateral spread' with 0.43%. As a result, we obtained the morphometric variables: slope, CONVI, TWI, landform, which were highly associated with landslides. The effect of a DEM in the processing flow of the InSAR method was similar for the InSAR coherence variable using the DEMs ASTER, PALSAR RTC, Topo-map, and SRTM. Then, a multiInSAR analysis gave displacement velocities in the LOS direction between -10 and 10 mm/year. With the dual-PolSAR analysis (Sentinel-1), VH and VV C-band polarised radar energy emitted median values of backscatters, for landslides, about of -14.5 dB for VH polarisation and -8.5 dB for VV polarisation. Also, L-band fully polarimetric NASA-UAVSAR data allowed to nd the mechanism of dispersion of CGS landslide inventory: 39% for surface scattering, 46.4% for volume dispersion, and 14.6% for double-bounce scattering. The optical remote sensing provided NDVI composites derived from Landsat series between 2012 and 2016, showing that NDVI values between 0.40 and 0.70 had a high correlation to landslides. In summary, we found the highest categories related to landslides by Weight of Evidence method (WofE) for each spaceborne technique applied. Finally, these results were merged to generate the landslide detection model by using the supervised machine learning method of Random Forest. By taking training and test samples, the precision of the detection model was of about 70% for the rotational and translational types.Los deslizamientos son una amenaza natural que causa pérdidas humanas, daños a la infraestructura y degradación del suelo. Una evaluación cuantitativa de su presencia se requiere mediante la detección y el reconocimiento de potenciales áreas inestables. Esta investigación tuvo como alcance desarrollar un método soportado en métodos semi-automáticos para detectar potenciales movimientos en masa a escala regional. Cinco técnicas fueron estudiadas: Morfometría, Interferometría radar, Interferometría con Persistent Scatterers, Polarimetría radar y composiciones del NDVI con los satélites Landsat 5, Landsat 7 y Landsat 8. El caso de estudio se seleccionó dentro de la región intermedia al este del departamento del Cauca, la cual se caracteriza por terreno montañoso y la presencia de inestabilidades de la pendiente oficialmente registrados en el servicio SIMMA del Servicio Geológico Colombiano. Este inventario reveló que el tipo de movimiento deslizamiento ocurrió con una frecuencia relativa de 77.4%, caidos con el 16.5% de los casos y reptaciones con 3%, flujos con 2.6% y propagación lateral con 0.43%. Como resultado, se obtuvo las variables morfométricas: pendiente, convergencia, índice topográfico de humedad y forma del terreno altamente asociados con los deslizamientos. El efecto de un DEM en el procesamiento del método InSAR fue similar para la variable coherencia usando los DEMs: ASTER, PAlSAR RTC, Topo-map y SRTM. Un análisis Multi-InSAR estimó velocidades de desplazamiento en dirección de vista del radar entre -10 y 10 mm/año. El análisis de polarimetría dual del Sentinel-1 arrojó valores de retrodispersión promedio de -14.5 dB en la banda VH y -8.5dB en la banda VV. Las cuatro polarimetrías del sensor aéreo UAVSAR permitió caracterizar el mecanismo de dispersión del Inventario de Deslizamiento así: 39% en el mecanismo de superficie, 46.4% en el mecanismo de volumen y 14.6% en el mecanismo de doble rebote. La información generada en el rango óptico permitió obtener composiciones de NDVI derivados de la plataforma Landsat entre los años 2012 y 2016, mostrando que el rango entre 0.4 y 0.7 tuvieron una alta asociación con los deslizamientos. En esta investigación se determinaron las categorías de las variables de Teledetección más altamente relacionadas con los movimientos en masa mediante el método de Pesos de Evidencias (WofE). Finalmente, estos resultados se fusionaron para generar el modelo de detección de deslizamientos usando el método supervisado de aprendizaje de máquina Random Forest. Tomando muestras aleatorias para entrenar y validar el modelo en una proporción 70:30, el modelo de detección, especialmente los movimientos de tipo rotacional y traslacional fueron clasificados con una tasa general de éxito del 70%.Ministerio de CienciasConvocatoria 647 de 2014Research line: Geotechnics and Geoenvironmental HazardDoctorad

The BIOMASS level 2 prototype processor : design and experimental results of above-ground biomass estimation

BIOMASS is ESA’s seventh Earth Explorer mission, scheduled for launch in 2022. The satellite will be the first P-band SAR sensor in space and will be operated in fully polarimetric interferometric and tomographic modes. The mission aim is to map forest above-ground biomass (AGB), forest height (FH) and severe forest disturbance (FD) globally with a particular focus on tropical forests. This paper presents the algorithms developed to estimate these biophysical parameters from the BIOMASS level 1 SAR measurements and their implementation in the BIOMASS level 2 prototype processor with a focus on the AGB product. The AGB product retrieval uses a physically-based inversion model, using ground-canceled level 1 data as input. The FH product retrieval applies a classical PolInSAR inversion, based on the Random Volume over Ground Model (RVOG). The FD product will provide an indication of where significant changes occurred within the forest, based on the statistical properties of SAR data. We test the AGB retrieval using modified airborne P-Band data from the AfriSAR and TropiSAR campaigns together with reference data from LiDAR-based AGB maps and plot-based ground measurements. For AGB estimation based on data from a single heading, comparison with reference data yields relative Root Mean Square Difference (RMSD) values mostly between 20% and 30%. Combining different headings in the estimation process significantly improves the AGB retrieval to slightly less than 20%. The experimental results indicate that the implemented retrieval scheme provides robust results that are within mission requirements

Wetland mapping and monitoring using polarimetric and interferometric synthetic aperture radar (SAR) data and tools

Wetlands are home to a great variety of flora and fauna species and provide several unique environmental functions, such as controlling floods, improving water-quality, supporting wildlife habitat, and shoreline stabilization. Detailed information on spatial distribution of wetland classes is crucial for sustainable management and resource assessment. Furthermore, hydrological monitoring of wetlands is also important for maintaining and preserving the habitat of various plant and animal species. This thesis investigates the existing knowledge and technological challenges associated with wetland mapping and monitoring and evaluates the limitations of the methodologies that have been developed to date. The study also proposes new methods to improve the characterization of these productive ecosystems using advanced remote sensing (RS) tools and data. Specifically, a comprehensive literature review on wetland monitoring using Synthetic Aperture Radar (SAR) and Interferometric SAR (InSAR) techniques is provided. The application of the InSAR technique for wetland mapping provides the following advantages: (i) the high sensitivity of interferometric coherence to land cover changes is taken into account and (ii) the exploitation of interferometric coherence for wetland classification further enhances the discrimination between similar wetland classes. A statistical analysis of the interferometric coherence and SAR backscattering variation of Canadian wetlands, which are ignored in the literature, is carried out using multi-temporal, multi-frequency, and multi-polarization SAR data. The study also examines the capability of compact polarimetry (CP) SAR data, which will be collected by the upcoming RADARSAT Constellation Mission (RCM) and will constitute the main source of SAR observation in Canada, for wetland mapping. The research in this dissertation proposes a methodology for wetland classification using the synergistic use of intensity, polarimetry, and interferometry features using a novel classification framework. Finally, this work introduces a novel model based on the deep convolutional neural network (CNN) for wetland classification that can be trained in an end-to-end scheme and is specifically designed for the classification of wetland complexes using polarimetric SAR (PolSAR) imagery. The results of the proposed methods are promising and will significantly contribute to the ongoing efforts of conservation strategies for wetlands and monitoring changes. The approaches presented in this thesis serve as frameworks, progressing towards an operational methodology for mapping wetland complexes in Canada, as well as other wetlands worldwide with similar ecological characteristics

Advanced machine learning algorithms for Canadian wetland mapping using polarimetric synthetic aperture radar (PolSAR) and optical imagery

Wetlands are complex land cover ecosystems that represent a wide range of biophysical conditions. They are one of the most productive ecosystems and provide several important environmental functionalities. As such, wetland mapping and monitoring using cost- and time-efficient approaches are of great interest for sustainable management and resource assessment. In this regard, satellite remote sensing data are greatly beneficial, as they capture a synoptic and multi-temporal view of landscapes. The ability to extract useful information from satellite imagery greatly affects the accuracy and reliability of the final products. This is of particular concern for mapping complex land cover ecosystems, such as wetlands, where complex, heterogeneous, and fragmented landscape results in similar backscatter/spectral signatures of land cover classes in satellite images. Accordingly, the overarching purpose of this thesis is to contribute to existing methodologies of wetland classification by proposing and developing several new techniques based on advanced remote sensing tools and optical and Synthetic Aperture Radar (SAR) imagery. Specifically, the importance of employing an efficient speckle reduction method for polarimetric SAR (PolSAR) image processing is discussed and a new speckle reduction technique is proposed. Two novel techniques are also introduced for improving the accuracy of wetland classification. In particular, a new hierarchical classification algorithm using multi-frequency SAR data is proposed that discriminates wetland classes in three steps depending on their complexity and similarity. The experimental results reveal that the proposed method is advantageous for mapping complex land cover ecosystems compared to single stream classification approaches, which have been extensively used in the literature. Furthermore, a new feature weighting approach is proposed based on the statistical and physical characteristics of PolSAR data to improve the discrimination capability of input features prior to incorporating them into the classification scheme. This study also demonstrates the transferability of existing classification algorithms, which have been developed based on RADARSAT-2 imagery, to compact polarimetry SAR data that will be collected by the upcoming RADARSAT Constellation Mission (RCM). The capability of several well-known deep Convolutional Neural Network (CNN) architectures currently employed in computer vision is first introduced in this thesis for classification of wetland complexes using multispectral remote sensing data. Finally, this research results in the first provincial-scale wetland inventory maps of Newfoundland and Labrador using the Google Earth Engine (GEE) cloud computing resources and open access Earth Observation (EO) collected by the Copernicus Sentinel missions. Overall, the methodologies proposed in this thesis address fundamental limitations/challenges of wetland mapping using remote sensing data, which have been ignored in the literature. These challenges include the backscattering/spectrally similar signature of wetland classes, insufficient classification accuracy of wetland classes, and limitations of wetland mapping on large scales. In addition to the capabilities of the proposed methods for mapping wetland complexes, the use of these developed techniques for classifying other complex land cover types beyond wetlands, such as sea ice and crop ecosystems, offers a potential avenue for further research

A Novel Phenology Based Feature Subset Selection Technique Using Random Forest for Multitemporal PolSAR Crop Classification

Feature selection techniques intent to select a subset of features that minimizes redundancy and maximizes relevancy for classification problems in machine learning. Standard methods for feature selection in machine learning seldom take into account the domain knowledge associated with the data. Multitemporal crop classification studies with full-polarimetric synthetic aperture radar (PolSAR) data ought to consider the changes in the scattering mechanisms with their phenological growth stages. Hence, it is desirable to incorporate these changes while determining a feature subset for classification. In this study, a random forest (RF) based feature selection technique is proposed that takes into account the changes in the physical scattering mechanism with crop phenological stages for multitemporal PolSAR classification. The partial probability plot, which is an attribute of RF, provides information about the marginal effect of a polarimetric parameter on the desired crop class. Moreover, it is used to identify the specific range of a parameter where the probability of the presence of a particular crop class is high. The proposed technique identifies features that change significantly with crop phenology. The selected features are the ones whose ranges show maximum separation amongst crop classes. Additionally, the feature subset is refined by eliminating correlated features. The E-SAR L-band dataset of the AgriSAR-2006 campaign over the Demmin test site in Germany is used in this study. The classification accuracy using the novel feature selection technique is 99.12%. This is nominally better than using the features obtained from a standard feature selection method used in RF, such as mean decrease Gini (98.73%) and mean decrease accuracy (98.68%) that do not take into account the information based on crop phenology.This work was supported in part by the Spanish Ministry of Economy, Industry and Competitiveness, in part by the State Agency of Research (AEI), and in part by the European Funds for Regional Development under Projects TIN2014-55413-C2-2-P and TEC2017-85244-C2-1-P