Polarimetric Interferometric Synthetic Aperture Radar Research Articles

A Multi-Baseline Forest Height Estimation Method Combining Analytic and Geometric Expression of the RVoG Model

As an important parameter of forest biomass, forest height is of great significance for the calculation of forest carbon stock and the study of the carbon cycle in large-scale regions. The main idea of the current forest height inversion methods using multi-baseline P-band polarimetric interferometric synthetic aperture radar (PolInSAR) data is to select the best baseline for forest height inversion. However, the approach of selecting the optimal baseline for forest height inversion results in the process of forest height inversion being unable to fully utilize the abundant observation data. In this paper, to solve the problem, we propose a multi-baseline forest height inversion method combining analytic and geometric expression of the random volume over ground (RVoG) model, which takes into account the advantages of the selection of the optimal observation baseline and the utilization of multi-baseline information. In this approach, for any related pixel, an optimal baseline is selected according to the geometric structure of the coherence region shape and the functional model for forest height inversion is established by the RVoG model’s analytic expression. In this way, the other baseline observations are transformed into a constraint condition according to the RVoG model’s geometric expression and are also involved in the forest height inversion. PolInSAR data were used to validate the proposed multi-baseline forest height inversion method. The results show that the accuracy of the forest height inversion with the algorithm proposed in this paper in a coniferous forest area and tropical rainforest area was improved by 17% and 39%, respectively. The method proposed in this paper provides a multi-baseline PolInSAR forest height inversion scheme for exploring regional high-precision forest height distribution. The scheme is an applicable method for large-scale, high-precision forest height inversion tasks.

Volume Structure Retrieval Using Drone-Based SAR Interferometry with Wide Fractional Bandwidth

Synthetic aperture radar (SAR) interferometry (InSAR) is a well-established remote sensing technique capable of providing accurate topographic information of an area. Multiple scattering occurring at different heights within a single resolution cell, however, cannot be resolved using a single baseline and results in a degradation of the height accuracy. Techniques such as polarimetric InSAR and SAR tomography tackle this problem using additional measurements to obtain the three-dimensional structure of the volume. However, polarimetry-based methods assume orthogonal and deterministic scattering mechanisms and tomography requires numerous observations. This paper presents a novel approach to retrieving the three-dimensional structure of semitransparent media using wide fractional bandwidth InSAR. The frequency dependency of the InSAR coherence is exploited to retrieve the three-dimensional structure, assuming that signals with wide fractional bandwidth are available, as is the case for unmanned aerial vehicles (UAVs) or drones. As the extinction is, in general, frequency-dependent, a third observation may improve the inversion results. Simulations of different scenarios show that the frequency profile can be inverted to obtain information about the three-dimensional scattering structure. The acquired data also enable retrieval of the phase center depth as a function of frequency. This method allows for frequent and accurate monitoring of semitransparent targets, e.g., forests and vegetation, using a swarm of two or three drones.

Underlying Topography Estimation over Forest Using Maximum a Posteriori Inversion with Spaceborne Polarimetric SAR Interferometry

This paper presents a method for extracting the digital elevation model (DEM) of forested areas from polarimetric interferometric synthetic aperture radar (PolInSAR) data. The method models the ground phase as a Von Mises distribution, with a mean of the topographic phase computed from an external DEM. By combining the prior distribution of the ground phase with the complex Wishart distribution of the observation covariance matrix, we derive the maximum a posterior (MAP) inversion method based on the RVoG model and analyze its Cramer–Rao Lower Bound (CRLB). Furthermore, considering the characteristics of the objective function, this paper introduces a Four-Step Optimization (FSO) method based on gradient optimization, which solves the inefficiency problem caused by exhaustive search in solving ground phase using the MAP method. The method is validated using spaceborne L-band repeat-pass SAOCOM data from a test forest area. The test results for FSO indicate that it is approximately 5.6 times faster than traditional methods without compromising accuracy. Simultaneously, the experimental results demonstrate that the method effectively solves the problem of elevation jumps in DEM inversion when modeling the ground phase with the Gaussian distribution. ICESAT-2 data are used to evaluate the accuracy of the inverted DEM, revealing that our method improves the root mean square error (RMSE) by about 23.6% compared to the traditional methods.

Rice Height Estimation with Multi-Baseline PolInSAR Data and Optimal Detection Baseline Combination Analysis

Rice is a primary food source, and height is a crucial parameter affecting its growth status. Consequently, high-precision, real-time monitoring of quantitative changes in crop height are required for improved crop production. Polarimetric interferometric SAR (PolInSAR) has both polarization and interferometric observation capabilities. Due to the short height of crops and rapid growth changes, the large spatial and short temporal baselines of PolInSAR data are essential for effective crop height inversion; however, relevant satellite-borne SAR and airborne SAR data are currently limited. This study presents a PolInSAR rice height inversion algorithm that uses the oriented volume over ground (OVoG) mode with PolInSAR 0-time and controllable spatial baseline data from a LAMP microwave anechoic chamber. Exploiting the advantages of microwave anechoic chamber measurement data, which includes continuous frequency bands, multiple baselines, and varied incidence angles, the influences of incident angles, baseline length, number of baselines, and baseline combinations are assessed. The highest accuracy rice plant height inversion has a root mean square deviation (RMSE) of 0.1093 m and is achieved with an incidence angle of 35–55°, baseline length of 0.25°, and three to four baselines. Furthermore, an imaging geometric equivalence analysis provides reliable foundation data to guide research into new earth observation SAR systems. The results indicate that, under simulated observations from the GF3 satellite at an altitude of 755 km, the optimal spatial baseline ranges for various frequency bands are as follows: L-band: 10.93–42.59 km; S-band: 4.10–15.97 km; C-band: 2.48–9.64 km; X-band: 1.36–5.32 km; Ku-band: 0.87–3.40 km. Notably, the measurement modes corresponding to the C, X, and Ku bands are ultimately identified as the most suitable for PolInSAR rice height inversion.

A Fourier–Legendre Polynomial Forest Height Inversion Model Based on a Single-Baseline Configuration

In this article, we propose a Fourier–Legendre (FL) polynomial forest height estimation algorithm based on low-frequency single-baseline polarimetric interferometric synthetic aperture radar (PolInSAR) data. The algorithm can obtain forest height with a single-baseline PolInSAR configuration while capturing a high-resolution vertical profile for the forest volume. This is based on the consideration that the forest height remains constant within neighboring pixels. Meanwhile, we also assume that the coefficients of the FL polynomials remain unchanged within neighboring pixels, except for the last polynomial coefficient. The idea of using neighboring pixels to increase the observations provides us with the possibility to obtain high-order FL polynomials. With this approach, it is possible to obtain a high-resolution vertical profile that is suitable for forest height estimation without losing too much spatial resolution. P-band PolInSAR data acquired in Mabounie in Gabon and Krycklan in Sweden were selected for testing the proposed algorithm. The results show that the algorithm outperforms the random volume over ground (RVoG) model by 18% and 16.7% in forest height estimation for the Mabounie and Krycklan study sites, respectively.

A new snow depth retrieval method by improved hybrid DEM differencing and coherence amplitude algorithm for PolInSAR

Snow depth is a critical parameter for measuring the variations of the snowpack, and its measurement also influences the accuracy of hydrological and climate models. Polarimetric SAR interferometry (PolInSAR) technology combines the capabilities of Polarimetric SAR (PolSAR) and Interferometric SAR (InSAR), and has significant potential for retrieving snow depth. The hybrid DEM differencing and coherence amplitude algorithm (HDCA) is a robust algorithm for PolInSAR data, but encounters such problems as uncertainty in finding the ground scattering phase center, impure volume decorrelation, and the lack of consideration of the dense medium characteristic of snow. In this study, the HDCA is improved by optimizing the phase and decorrelation parameters and by considering the dense medium characteristic of snow. First, the Freeman-Duren decomposition method is put forward to extract the ground and volume scattering phase centers from the snow. Second, the coherence region boundary estimation method is proposed to optimize the volume decorrelation. Following this, the parameter assumed by the original method in free space is adjusted by considering the dense medium characteristic of snow. Finally, the Ku band UAV SAR data were applied to validate the new method, and in-situ data were used to assess its accuracy. The correlation coefficient (R), root mean square error (RMSE), and mean absolute error (MAE) of the proposed method are 0.93, 3.59 cm, and 2.98 cm, respectively, and were significantly superior to those of the original HDCA method of 0.74, 7.44 cm, and 6.04 cm, respectively.

Pol-InSAR-Island - A benchmark dataset for multi-frequency Pol-InSAR data land cover classification

This paper presents Pol-InSAR-Island, the first publicly available multi-frequency Polarimetric Interferometric Synthetic Aperture Radar (Pol-InSAR) dataset labeled with detailed land cover classes, which serves as a challenging benchmark dataset for land cover classification. In recent years, machine learning has become a powerful tool for remote sensing image analysis. While there are numerous large-scale benchmark datasets for training and evaluating machine learning models for the analysis of optical data, the availability of labeled SAR or, more specifically, Pol-InSAR data is very limited. The lack of labeled data for training, as well as for testing and comparing different approaches, hinders the rapid development of machine learning algorithms for Pol-InSAR image analysis. The Pol-InSAR-Island benchmark dataset presented in this paper aims to fill this gap. The dataset consists of Pol-InSAR data acquired in S- and L-band by DLR's airborne F-SAR system over the East Frisian island Baltrum. The interferometric image pairs are the result of a repeat-pass measurement with a time offset of several minutes. The image data are given as 6 × 6 coherency matrices in ground range on a 1 m × 1m grid. Pixel-accurate class labels, consisting of 12 different land cover classes, are generated in a semi-automatic process based on an existing biotope type map and visual interpretation of SAR and optical images. Fixed training and test subsets are defined to ensure the comparability of different approaches trained and tested prospectively on the Pol-InSAR-Island dataset. In addition to the dataset, results of supervised Wishart and Random Forest classifiers that achieve mean Intersection-over-Union scores between 24% and 67% are provided to serve as a baseline for future work. The dataset is provided via KITopenData: https://doi.org/10.35097/1700.

A novel framework for combining polarimetric Sentinel-1 InSAR time series in subsidence monitoring - A case study of Sydney

The rapid growth of the city of Sydney, Australia over the last decades, has led to significant development of residential and transportation infrastructure. Land subsidence associated with the urban development can lead to serious issues which should be thoroughly understood and carefully managed. To address this challenge, an enhanced polarisation time-series InSAR (Pol-TS-InSAR) processing framework was developed, using the dual polarisation (DP) Sentinel-1 data to integrate information from different polarimetric channels with different weighting during the TS-InSAR deformation analysis. Ninety DP Sentinel-1 images acquired between 2019 and 2022 are analysed using Pol-TS-InSAR to map the land subsidence in Sydney, with the assistance of the GPS measurements. Improvement of measurement points density from Pol-TS-InSAR is observed compared to the single polarimetric TS-InSAR counterpart for all land use types (ranging between 68% and 208%). The comparison between the Pol-TS-InSAR measurements and GPS measurements shows an absolute mean difference and RMS difference of 0.75 mm/yr and 0.95 mm/yr, respectively, in vertical direction. The results of the ground subsidence analysis revealed that the main subsidence factors in Sydney are related to groundwater extraction, mining activities, underground tunnel construction and landfill. The latter two factors were less well-known prior to this study. In additional to these factors, land subsidence related to high-rise building construction has also been observed, even though the impact seems to be less significant than other factors.

Patch‐based nonlocal estimation of scattering mechanism vector in PolInSAR

AbstractNowadays, polarimetric interferometric synthetic aperture radar (PolInSAR) has attracted increasing attention for the simultaneous acquisition of scattering and terrain information. The scattering mechanism vectors corresponding to the optimal coherences have shown great effect in land cover classification and forest mapping. To extract scattering mechanism vectors, different optimization algorithms are proposed but always accompanied with spatial discontinuities. This letter proposes a method to optimize scattering mechanism vectors based on nonlocal estimations. A pixel‐by‐pixel optimization model is established by calculating patch‐based nonlocal similarity, which effectively reduces the uncertainty of scattering mechanism vector estimation while maintaining high resolution. The phase recovery experiment is conducted on simulated data and the root mean square errors (RMSEs) of three algorithms are compared. The proposed method performs better with a RMSE of 0.22 in the whole image, against 0.25 and 1.27 for other algorithms. The recovery results in the target region (RMSE 0.11 against 0.21 and 0.99) and linear edge region (RMSE 0.26 against 0.49 and 1.52) also indicate strong ability in detail preservation. Then the observed ESAR data is utilized for further assessment. It demonstrates that the proposed method has a significant improvement to reduce singularities and distinguish artificial buildings.

A Dual-Baseline PolInSAR Method for Forest Height and Vertical Profile Function Inversion Based on the Polarization Coherence Tomography Technique

Forest height and vertical structure profile functions can be estimated using polarimetric interferometric synthetic aperture radar (PolInSAR) data based on the random volume over ground (RVoG) model and polarization coherence tomography (PCT) theory, respectively. For each resolution cell, considering different forest vertical scattering structure functions to solve the corresponding forest height, the accuracy of PolInSAR forest height inversion will be improved. In this study, a forest vertical structure profile function and forest height inversion algorithm based on PCT technology was developed by using dual-baseline PolInSAR data. Then the deviation of forest height was corrected according to the inverted forest vertical structure. Finally, the LiDAR and PolInSAR data were employed to verify the proposed method. The experimental results show that the accuracy of the proposed method (tropical forest: RMSE = 5.96 m, boreal forest: RMSE = 3.11 m) is 25.5% and 30.43% higher than that of the dual-baseline RVoG model algorithm (tropical forest: RMSE = 8 m, boreal forest: RMSE = 4.47 m).

Improved ESPO Method Based on Spatial Variation of Scattering Mechanisms for MT-PolInSAR

The existing exhaustive search polarimetric optimization (ESPO) method exploiting polarimetric diversity has been widely used to improve the phase quality in polarimetric interferometric synthetic aperture radar (PolInSAR). However, the optimization ceiling of the ESPO method based on the coherence metric is largely restricted by the variation of scattering mechanisms in the spatial domain, especially for high-resolution SAR data. To this end, this letter proposes an improved ESPO (ImESPO) method that considers the changes of the scattering mechanisms in local windows and applies it to the time-series phase optimization framework. Both simulated and real experiments demonstrated the effectiveness of the proposed method, which shows that the changes of the scattering mechanisms are more serious in the large window or for high-resolution SAR data. In addition, the effects of the filtering window size and the number of interferograms on optimization were also analyzed in detail.

Bayesian analysis for uncertainty of forest height inversed by polarimetric interferometric SAR data

Bayesian analysis for uncertainty of forest height inversed by polarimetric interferometric SAR data

Distributed Scatterer Processing Based on Binary Partition Trees with Multi-Baseline PolInSAR Data

Distributed scatterers (DSs) are necessary to increase point density in multi-temporal InSAR (MT-InSAR) monitoring. The identification of homogeneous pixels (HPs) is the first and key step for DS processing to overcome the low signal-to-noise ratio condition. Since multi-polarization data are good at describing geometrical structures and dielectric properties of ground objects, they have been applied for HP identification. However, polarimetric information is not enough for identifying areas with similar ground objects but different deformation. We propose a novel DS preprocessing algorithm based on polarimetric interferometric homogeneous pixel (PIHP) identification. Firstly, a novel Polarimetric InSAR (PolInSAR) similarity that combines polarimetric intensity, interferometric coherence, and phase is proposed, which is readily available in multi-baseline and multi-polarization data and flexible by controlling weighting factors. Secondly, based on the binary partition tree (BPT) framework, object-orientated multi-scale PIHP identification is achieved, which is suitable for complex deformation scenes. Tested with simulated quad-polarization data, our method shows improvement in phase quality and point density, especially in the deformed areas, compared with the traditional HP identification method based on the polarimetric homogeneity (PolHom) test and the method with ground object type map. Tested with 30 quad-polarization Radarsat-2 images over Kilauea Volcano, the point density of our method is three times higher than that of the PolHom test in vegetation areas. Our method is proven to be more sensitive and mechanically more advanced to homogeneous pixels identification than the traditional ones, which is helpful for phase optimization, spatial enlargement of monitoring points, and stability of the MT-InSAR algorithm.

Rice Crop Height Inversion from TanDEM-X PolInSAR Data Using the RVoG Model Combined with the Logistic Growth Equation

The random volume over ground (RVoG) model has been widely used in the field of vegetation height retrieval based on polarimetric interferometric synthetic aperture radar (PolInSAR) data. However, to date, its application in a time-series framework has not been considered. In this study, the logistic growth equation was introduced into the PolInSAR method for the first time to assist in estimating crop height, and an improved inversion scheme for the corresponding RVoG model parameters combined with the logistic growth equation was proposed. This retrieval scheme was tested using a time series of single-pass HH-VV bistatic TanDEM-X data and reference data obtained over rice fields. The effectiveness of the time-series RVoG model based on the logistic growth equation and the convenience of using equation parameters to evaluate vegetation growth status were analyzed at three test plots. The results show that the improved method can effectively monitor the height variation of crops throughout the whole growth cycle and the rice height estimation achieved an accuracy better than when single dates were considered. This proved that the proposed method can reduce the dependence on interferometric sensitivity and can achieve the goal of monitoring the whole process of rice height evolution with only a few PolInSAR observations.

Forest Height Retrieval Based on the Dual PolInSAR Images

A new algorithm for forest height estimation based on dual polarimetric interferometric SAR data is presented in this study. The main objective is to consider the efficiency of the dual-polarization data compared to the full polarimetric images with respect to forest height retrieval. Accordingly, the forest height estimation based on the random volume over the ground model is examined using a geometrical procedure named the three-stage method. An exhaustive search polarization optimization technique is also applied to improve the results by employing the efficiency of all the polarization bases based on the four-dimensional lexicographic PolInSAR vector. The repeat-pass experimental SAR (ESAR) images, which include both L- and P-band full polarimetric data, are employed for the accuracy assessment of the dual PolInSAR data and the newly proposed method for forest height estimation. The experimental results on the L-band PolInSAR data show the ability of the dual PolInSAR data for forest height estimation with an average root mean square error (RMSE) of 4.97 m against Lidar data based on the conventional three-stage method. Additionally, the proposed method results in an accuracy of 2.95 m for forest height estimation, indicating its high potential for tree height retrieval.

A Polarization Stacking Method for Optimizing Time-Series Interferometric Phases of Distributed Scatterers

For time-series interferometric phases optimization of distributed scatterers (DSs), the SqueeSAR technology used the phase linking (PL) to extract the equivalent single-master (ESM) interferometric phases from the multilooking time-series coherence matrix. The Cramer–Rao lower bound (CRLB) for the PL describes the highest achievable estimation accuracy of the ESM phases, which depends on the number of looks and the time-series coherence magnitude matrix. With the abundance of time-series polarimetric SAR data, many scholars have studied the coherence magnitude-based polarimetric optimization methods for optimizing the DS’s time-series interferometric phases, for example, the widely-used exhaustive search polarimetric optimization (ESPO) algorithm. However, the traditional polarimetric optimization methods select the boundary extremums of the coherence region (CR) as the optimized complex coherence, which is usually biased from the free-noise one. Currently, in the polarimetric InSAR (PolInSAR) technology, Shen et al. innovatively considered polarimetric information as statistical samples and proposed the total power (TP) coherency matrix construction method for increasing the number of looks and reducing the interferometric phase noise. Therefore, to optimize the time-series interferometric phases for DS, this paper proposes performing a polarization stacking and extending the PolInSAR TP construction to the time-series PolInSAR (TSPolInSAR) data configuration, called the time-series TP (TSTP) method. Simulated and real experiments prove that the new TSTP construction method has better performance and higher efficiency than the single polarimetric and the traditional ESPO algorithms.

Evaluation of Multilooking Size on Single-Baseline PolInSAR Forest Height Inversion

In polarimetric interferometric SAR (PolInSAR) technology, the random volume over ground (RVoG) model establishes the mapping relationship between polarimetric complex coherence and forest biophysical parameters (e.g., forest height). However, due to speckle noise and the finite multilooking effect, the real observed coherence region in the complex unitary circle (CUC) is an ellipse, which is biased by the ideal noise-free coherence region represented as a straight line by the RVoG model. Multilooking processing can reduce speckle noise at the cost of resolution loss. Therefore, this paper analyzes the influence of different multilooking sizes on forest height inversion. Experimental results show that the accuracy of forest height inversion first increases and then decreases with the increase in multilooking size, which means there exists an optimal size for PolInSAR forest estimation. From statistical analysis of the forest height estimation error, inversion accuracy mainly depends on estimation bias rather than estimation variance. This is mainly because, in a homogeneous forest area, a large multilooking size helps to reduce the statistical bias effect; in the textured area, the inversion accuracy benefits from a small multilooking size for avoiding the mixing of multiple types of ground targets.

The Iterative Extraction of the Boundary of Coherence Region and Iterative Look-Up Table for Forest Height Estimation Using Polarimetric Interferometric Synthetic Aperture Radar Data

In this paper, we introduce a refined three-stage inversion algorithm (TSIA) for forest height estimation using polarimetric interferometric synthetic aperture radar (PolInSAR). Specifically, the iterative extraction of the boundary of the coherence region (IEBCR) and iterative look-up table (ILUT) are proposed to improve the efficiency of traditional TSIA. A class of refined TSIA utilizes the boundary of the coherence region (BCR) to alleviate the underestimation phenomenon in forest height estimation. Given many eigendecompositions in the extraction of BCR (EBCR), we analyze the relationship of eigenvectors between the adjacent points on the BCR and propose the IEBCR utilizing the power methods. In the final inversion stage of TSIA, the look-up table (LUT) uses the exhaustive search method to minimize the loss function in the 2-D grid with defined step sizes and thus costs high computational complexity. To alleviate the deficiency, we define the random volume over ground (RVoG) function based on the RVoG model and prove its monotonicity and convergence from the analytical and numerical points of view. After analyzing the relationship between the RVoG function and the loss function, we propose the ILUT for the inversion stage. The simulation and experiments based on the BioSAR 2008 campaign data illustrate that the IEBCR and ILUT greatly improve the computational efficiency almost without compromising on accuracy.

Multi-Rotor UAV-Borne PolInSAR Data Processing and Preliminary Analysis of Height Inversion in Urban Area

Polarimetric Interferometric Synthetic Aperture Radar (PolInSAR) has many useful applications, especially in forest areas. With the development of SAR miniaturization technology, researchers can install PolInSAR on small unmanned aerial vehicles (UAV), which can reduce flight costs. Limited by size and power, UAV-borne SAR usually works in a high-frequency band, which restricts its application to such things as vegetation height inversion. While on the other hand, the high resolution acquired under a short wavelength promises its application in urban areas. However, there are fewer studies on the application of PolInSAR in urban areas compared with that in forest areas. In this paper, we propose a processing method for a Ku-band multi-rotor-UAV-borne PolInSAR and provide a preliminary analysis of height inversion results on its data from the Fudan campus in Shanghai. We obtain the digital surface model (DSM) of different polarization modes and the DSM of polarimetric interferometry optimal decomposition in this area, whose RMSE is 2.88 m. On this basis, the elevation inversion results of targets such as buildings, lampposts, and trees are compared and analyzed. We preliminarily explore and analyze the reasons for the different results of different targets. To this end, we propose a mathematical derivation of the relationship between the interferometric phase between PolInSAR and InSAR of Pauli decomposition. We also perform a simulation to analyze the relationship between the phase center height of Pauli decomposition and PolInSAR under different cases. It provides a reference for the application of small UAV-borne PolInSAR in urban areas.

A Forest Height Joint Inversion Method Using Multibaseline PolInSAR Data

Estimating vegetation height from polarimetric interferometric synthetic aperture radar (PolInSAR) data using the random volume over ground (RVoG) model has long been used. Most of these methods propose models and apply them to real airborne data to demonstrate their potential. The single-baseline PolInSAR forest height estimation based on the RVoG model lacks sufficient observation information. For this reason, multibaseline data are introduced to address this. This paper fits the relationship of model parameters in multibaseline observation scenarios, and focuses the forest height inversion on the calculation of pure volume decorrelation. Subsequently, a multibaseline forest height joint inversion method based on the least squares principle is adopted. Finally, we use airborne PolInSAR data from the Lope and Mondah sites collected by UAVSAR and F-SAR systems during AfriSAR 2016 to verify the proposed method. The experimental results show that the accuracy of the proposed method (Lope: root mean square error (RMSE) = 5.8 m, Mondah: RMSE = 5.12 m) is 38.1% and 34.53% higher than the coherence separation product (Lope: RMSE = 9.37 m, Mondah: RMSE = 7.82 m).