Improvement of PolSAR Decomposition Scattering Powers Using a Relative Decorrelation Measure

In this letter, a methodology is proposed to improve the scattering powers obtained from model-based decomposition using Polarimetric Synthetic Aperture Radar (PolSAR) data. The novelty of this approach lies in utilizing the intrinsic information in the off-diagonal elements of the 3$\times$3 coherency matrix $\mathbf{T}$ represented in the form of complex correlation coefficients. Two complex correlation coefficients are computed between co-polarization and cross-polarization components of the Pauli scattering vector. The difference between modulus of complex correlation coefficients corresponding to $\mathbf{T}^{\mathrm{opt}}$ (i.e. the degree of polarization (DOP) optimized coherency matrix), and $\mathbf{T}$ (original) matrices is obtained. Then a suitable scaling is performed using fractions \emph{i.e.,} $(T_{ii}^{\mathrm{opt}}/\sum\limits_{i=1}^{3}T_{ii}^{\mathrm{opt}})$ obtained from the diagonal elements of the $\mathbf{T}^{\mathrm{opt}}$ matrix. Thereafter, these new quantities are used in modifying the Yamaguchi 4-component scattering powers obtained from $\mathbf{T}^{\mathrm{opt}}$. To corroborate the fact that these quantities have physical relevance, a quantitative analysis of these for the L-band AIRSAR San Francisco and the L-band Kyoto images is illustrated. Finally, the scattering powers obtained from the proposed methodology are compared with the corresponding powers obtained from the Yamaguchi \emph{et. al.,} 4-component (Y4O) decomposition and the Yamaguchi \emph{et. al.,} 4-component Rotated (Y4R) decomposition for the same data sets. The proportion of negative power pixels is also computed. The results show an improvement on all these attributes by using the proposed methodology.Comment: Accepted for publication in Remote Sensing Letter

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

Forest attributes mapping with SAR data in the romanian South-Eastern Carpathians requirements and outcomes

Esta tesis doctoral se centra en la estimación de variables forestales en la zona Sureste de los Cárpatos Rumanos a partir de imágenes de radar de apertura sintética. La investigación abarca parte del preprocesado de las imágenes, métodos de generación de mosaicos y la extracción de la cobertura de bosque, sus subtipos o su biomasa. La tesis se desarrolló en el Instituto Nacional de Investigación y Desarrollo Forestal Marín Dracea (INCDS) y la Universidad de Alcalá (UAH) gracias a varios proyectos: el proyecto EO-ROFORMON del INCDS (Prototyping an Earth-Observation based monitoring and forecasting system for the Romanian forests), y el proyecto EMAFOR de la UAH (Synthetic Aperture Radar (SAR) enabled Analysis Ready Data (ARD) cubes for efficient monitoring of agricultural and forested landscapes). El proyecto EO-ROFORMON fue financiado por la Autoridad Nacional para la Investigación Científica de Rumania y el Fondo Europeo de Desarrollo Regional. El proyecto EMAFOR fue financiado por la Comunidad Autónoma de Madrid (España). El objetivo de esta tesis es el desarrollo de algoritmos para la extracción de variables forestales de uso general como la cobertura, el tipo o la biomasa del bosque a partir de imagen de radar de apertura sintética. Para alcanzar dicho propósito se analizaron posibles fuentes de sesgo sistemático que podrían aparecer en zonas de montaña (ej., normalización topográfica, generación de mosaicos), y se aplicaron técnicas de aprendizaje de máquina para tareas de clasificación y regresión. La tesis contiene ocho secciones: una introducción, cinco publicaciones en revistas o actas de congresos indexados, una pendiente de publicación (quinto capítulo) y las conclusiones. La introducción contextualiza la importancia del bosque, cómo se recoge la información sobre su estado (ej., inventario forestal) y las iniciativas o marcos legislativos que requieren dicha información. A continuación, se describe cómo la teledetección puede complementar la información de inventario forestal, detallando el contexto histórico de las distintas tecnologías, su funcionamiento, y cómo pueden ser aplicadas para la extracción de información forestal. Por último, se describe la problemática y el monitoreo del bosque en Rumanía, detallando el objetivo de la tesis y su estructura. El primer capítulo analiza la influencia del modelo digital de elevaciones (MDE) en la calidad de la normalización topográfica, analizando tres MDE globales (SRTM, AW3D y TanDEM-X DEM) y uno nacional (PNOA-LiDAR). Los experimentos se basan en la comparación entre órbitas, con un MDE de referencia, y la variación del acierto en la clasificación dependiendo del MDE empleado para la normalización. Los resultados muestran una menor diferencia ente órbitas al utilizar un MDE con una mejor resolución (ej. TanDEM-X, PNOA-LIDAR), especialmente en el caso de zonas con fuertes pendientes o formas del terreno complejas, como pueden ser los valles. En zonas de alta montaña las imágenes de radar de apertura sintética (SAR) sufren frecuentes distorsiones. Estas distorsiones dependen de la geometría de adquisición, por lo que es posible combinar imágenes adquiridas desde varias órbitas para que la cobertura sea lo más completa posible. El segundo capítulo evalúa dos metodologías para la clasificación de usos del suelo utilizando datos de Sentinel-1 adquiridos desde varias órbitas. El primer método crea clasificaciones por órbita y las combina, mientras que el segundo genera un mosaico con datos de múltiples órbitas y lo clasifica. El acierto obtenido mediante combinación de clasificaciones es ligeramente mayor, mientras que la clasificación de mosaicos tiene importantes omisiones de las zonas boscosas debido a problemas en la normalización topográfica y a los efectos direccionales. El tercer capítulo se enfoca en separar la cobertura forestal de otras coberturas del suelo (urbano, vegetación baja, agua) analizando la utilidad de las variables basadas en la coherencia interferométrica. En él se realizan tres clasificaciones de máquina vector-soporte basadas en un conjunto concreto de variables. El primer conjunto contiene las estadísticas anuales de la retrodispersión (media y desviación típica anual), el segundo añade la coherencia a largo plazo (separación temporal mayor a un año), el tercero incluye las estadísticas de la coherencia a corto plazo (mínima separación temporal). Utilizar variables basadas en la coherencia aumenta el acierto de la clasificación hasta un 5% y reduce los errores de omisión de la cobertura forestal. El cuarto capítulo evalúa la posibilidad de detectar talas selectivas utilizando datos de Sentinel-1 y Sentinel-2. Sus resultados muestran que la detección resulta muy difícil debido a la saturación de los sensores y la confusión introducida por el efecto de la fenología. El quinto capítulo se centra en la clasificación de tipos de bosque basado en una serie temporal de datos Sentinel-1. Se basa en la creación de un conjunto de modelos que describen la relación entre la retrodispersión y el ángulo local de incidencia para un determinado tipo de bosque y fecha concreta. Para cada píxel se calcula el residuo respecto al modelo de cada uno de los tipos de bosque, acumulando dichos residuos a lo largo de la serie temporal. Hecho esto, cada píxel es asignado al tipo de bosque que acumula un menor residuo. Los resultados son prometedores, mostrando que frondosas y coníferas tienen un comportamiento distintivo, y que es posible separar ambos tipos de bosque con un alto grado de acierto. El sexto capítulo está dedicado a la estimación de biomasa utilizando datos Sentinel-1, ALOS PALSAR y regresión Random Forest. Se obtiene un error similar para ambos sensores a pesar de utilizar una banda diferente (band-C vs. -L), con poca reducción en el error cuando ambas bandas se utilizan conjuntamente. Sin embargo, el ajuste de un estimador adaptado a las condiciones locales de Rumanía sí ofreció una reducción de del error al ser comparado con las estimaciones globales de biomasa

Evaluation of the soil moisture prediction accuracy of a space radar using simulation techniques

Image simulation techniques were employed to generate synthetic aperture radar images of a 17.7 km x 19.3 km test site located east of Lawrence, Kansas. The simulations were performed for a space SAR at an orbital altitude of 600 km, with the following sensor parameters: frequency = 4.75 GHz, polarization = HH, and angle of incidence range = 7 deg to 22 deg from nadir. Three sets of images were produced corresponding to three different spatial resolutions; 20 m x 20 m with 12 looks, 100 m x 100 m with 23 looks, and 1 km x 1 km with 1000 looks. Each set consisted of images for four different soil moisture distributions across the test site. Results indicate that, for the agricultural portion of the test site, the soil moisture in about 90% of the pixels can be predicted with an accuracy of = + or - 20% of field capacity. Among the three spatial resolutions, the 1 km x 1 km resolution gave the best results for most cases, however, for very dry soil conditions, the 100 m x 100 m resolution was slightly superior

Applications of ground-based radar to mine slope monitoring

"Slope failure accidents were responsible for about 12% of U.S. surface mine fatalities from 1995 to 2003. Small surface movements on a mine highwall may be precursors of failure that, if detected, could provide sufficient warning to enable workers and machinery to be withdrawn to safety. Radar interferometry offers the necessary precision to detect these movements. Radar has some advantages over other methods in its ability to cover large surface areas for true two-dimensional monitoring day and night under almost any weather condition. Radar's active transmit/receive mode of operation provides for more direct sampling than passive optical methods that depend on solar illumination. Improvements in microprocessor speeds and capacities have led to the development of a number of small, portable, ground-based systems. Such systems are now being deployed at several locations around the world. As part of an ongoing study of monitoring technologies, researchers from NIOSH and Brigham Young University, Provo, UT, cooperated to assess the feasibility of using interferometric radar to monitor mine slope stability. Field tests of a device incorporating prototype equipment were successful in that small, centimeter-scale displacements on rock slopes were detected." - p. NIOSHTIC-

Comparing synthetic aperture radar and LiDAR for above-ground biomass estimation in Glen Affric, Scotland

Quantifying above-ground biomass (AGB) and carbon sequestration has been a significant focus of attention within the UNFCCC and Kyoto Protocol for improvement of national carbon accounting systems (IPCC, 2007; UNFCCC, 2011). A multitude of research has been carried out in relatively flat and homogeneous forests (Ranson & Sun, 1994; Beaudoin et al.,1994; Kurvonen et al., 1999; Austin et al., 2003; Dimitris et al., 2005), yet forests in the highlands, which generally form heterogeneous forest cover and sparse woodlands with mountainous terrain have been largely neglected in AGB studies (Cloude et al., 2001; 2008; Lumsdon et al., 2005; 2008; Erxue et al., 2009, Tan et al., 2010; 2011a; 2011b; 2011c; 2011d). Since mountain forests constitute approximately 28% of the total global forest area (Price and Butt, 2000), a better understanding of the slope effects is of primary importance in AGB estimation. The main objective of this research is to estimate AGB in the aforementioned forest in Glen Affric, Scotland using both SAR and LiDAR data. Two types of Synthetic Aperture Radar (SAR) data were used in this research: TerraSAR-X, operating at X-band and ALOS PALSAR, operating at L-band, both are fully polarimetric. The former data was acquired on 13 April 2010 and of the latter, two scenes were acquired on 17 April 2007 and 08 June 2009. Airborne LiDAR data were acquired on 09 June 2007. Two field measurement campaigns were carried out, one of which was done from winter 2006 to spring 2007 where physical parameters of trees in 170 circular plots were measured by the Forestry Commission team. Another intensive fieldwork was organised by myself with the help of my fellow colleagues and it comprised of tree measurement in two transects of 200m x 50m at a relatively flat and dense plantation forest and 400m x 50m at hilly and sparse semi-natural forest. AGB is estimated for both the transects to investigate the effectiveness of the proposed method at plot-level. This thesis evaluates the capability of polarimetric Synthetic Aperture Radar data for AGB estimation by investigating the relationship between the SAR backscattering coefficient and AGB and also the relationship between the decomposed scattering mechanisms and AGB. Due to the terrain and heterogeneous nature of the forests, the result from the backscatter-AGB analysis show that these forests present a challenge for simple AGB estimation. As an alternative, polarimetric techniques were applied to the problem by decomposing the backscattering information into scattering mechanisms based on the approach by Yamaguchi (2005; 2006), which are then regressed to the field measured AGB. Of the two data sets, ALOS PALSAR demonstrates a better estimation capacity for AGB estimation than TerraSAR-X. The AGB estimated results from SAR data are compared with AGB derived from LiDAR data. Since tree height is often correlated with AGB (Onge et al., 2008; Gang et al., 2010), the effectiveness of the tree height retrieval from LiDAR is evaluated as an indicator of AGB. Tree delineation was performed before AGB of individual trees were calculated allometrically. Results were validated by comparison to the fieldwork data. The amount of overestimation varies across the different canopy conditions. These results give some indication of when to use LiDAR or SAR to retrieve forest AGB. LiDAR is able to estimate AGB with good accuracy and the R2 value obtained is 0.97 with RMSE of 14.81 ton/ha. The R2 and RMSE obtained for TerraSAR-X are 0.41 and 28.5 ton/ha, respectively while for ALOS PALSAR data are 0.70 and 23.6 ton/ha, respectively. While airborne LiDAR data with very accurate height measurement and consequent three-dimensional (3D) stand profiles which allows investigation into the relationship between height, number density and AGB, it's limited to small coverage area, or large areas but at large cost. ALOS PALSAR, on the other hand, can cover big coverage area but it provide a lower resolution, hence, lower estimation accuracy

Shuttle imaging radar-C science plan

The Shuttle Imaging Radar-C (SIR-C) mission will yield new and advanced scientific studies of the Earth. SIR-C will be the first instrument to simultaneously acquire images at L-band and C-band with HH, VV, HV, or VH polarizations, as well as images of the phase difference between HH and VV polarizations. These data will be digitally encoded and recorded using onboard high-density digital tape recorders and will later be digitally processed into images using the JPL Advanced Digital SAR Processor. SIR-C geologic studies include cold-region geomorphology, fluvial geomorphology, rock weathering and erosional processes, tectonics and geologic boundaries, geobotany, and radar stereogrammetry. Hydrology investigations cover arid, humid, wetland, snow-covered, and high-latitude regions. Additionally, SIR-C will provide the data to identify and map vegetation types, interpret landscape patterns and processes, assess the biophysical properties of plant canopies, and determine the degree of radar penetration of plant canopies. In oceanography, SIR-C will provide the information necessary to: forecast ocean directional wave spectra; better understand internal wave-current interactions; study the relationship of ocean-bottom features to surface expressions and the correlation of wind signatures to radar backscatter; and detect current-system boundaries, oceanic fronts, and mesoscale eddies. And, as the first spaceborne SAR with multi-frequency, multipolarization imaging capabilities, whole new areas of glaciology will be opened for study when SIR-C is flown in a polar orbit

Structure-From-Motion Photogrammetry of Antarctic Historical Aerial Photographs in Conjunction with Ground Control Derived from Satellite Data

A longer temporal scale of Antarctic observations is vital to better understanding glacier dynamics and improving ice sheet model projections. One underutilized data source that expands the temporal scale is aerial photography, specifically imagery collected prior to 1990. However, processing Antarctic historical aerial imagery using modern photogrammetry software is difficult, as it requires precise information about the data collection process and extensive in situ ground control is required. Often, the necessary orientation metadata for older aerial imagery is lost and in situ data collection in regions like Antarctica is extremely difficult to obtain, limiting the use of traditional photogrammetric methods. Here, we test an alternative methodology to generate elevations from historical Antarctic aerial imagery. Instead of relying on pre-existing ground control, we use structure-from-motion photogrammetry techniques to process the imagery with manually derived ground control from high-resolution satellite imagery. This case study is based on vertical aerial image sets collected over Byrd Glacier, East Antarctica in December 1978 and January 1979. Our results are the oldest, highest resolution digital elevation models (DEMs) ever generated for an Antarctic glacier. We use these DEMs to estimate glacier dynamics and show that surface elevation of Byrd Glacier has been constant for the past ∼40 years