Polarimetric Interferometric SAR Data Research Articles

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.

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 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).

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

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.

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.

Machine-Learning Inversion of Forest Vertical Structure Based on 2-D-SGVBVoG Model for P-Band Pol-InSAR

The polarimetric interferometric synthetic aperture radar (Pol-InSAR) model under P-band observations exhibits vertical structure diversity. Compared with the exponential-based random volume over ground (RVoG) model, the Gaussian vertical backscatter volume over ground (GVBVoG) model expresses a more complex forest vertical structure via introducing more parameters. On account of the influence of topographic fluctuation on the model, this article establishes the sloped Gaussian vertical backscatter volume over ground (SGVBVoG) model by drawing into the terrain slope. Based upon the SGVBVoG model, this article develops the 2-D SGVBVoG (2-D-SGVBVoG) model by defining the structure factor, which effectively reduces the model complexity from three to two dimensions. In the 2-D-SGVBVoG model inversion, in view of the diversity of forest species, age, shape, density, etc., in the natural scene and the variation of specific radar systems, a structure factor prediction scheme relying on machine learning is proposed. In the machine-learning model training, the radar incidence angle and the PDHigh coherence acquired by coherence optimization with terrain phase removal are utilized as the variables for characterizing the structure factor. Ultimately, in the case of fixing structure factor, a geometric inversion process on the complex plane is put forward to extract the forest height. The BIOSAR 2008 P-band Pol-InSAR data validation shows that the proposed method achieves an RMSE of 3.07 m, which is 24.0% better than the three-baseline SRVoG inversion.

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).

The analysis on uncertainty resulting from method and wavelength selecting in forest height inversion using simulated polarimetric interferometric SAR data

The analysis on uncertainty resulting from method and wavelength selecting in forest height inversion using simulated polarimetric interferometric SAR data

EVALUATING UNIFORM MANIFOLD APPROXIMATION AND PROJECTION FOR DIMENSION REDUCTION AND VISUALIZATION OF POLINSAR FEATURES

Abstract. In this paper, the nonlinear dimension reduction algorithm Uniform Manifold Approximation and Projection (UMAP) is investigated to visualize information contained in high dimensional feature representations of Polarimetric Interferometric Synthetic Aperture Radar (PolInSAR) data. Based on polarimetric parameters, target decomposition methods and interferometric coherences a wide range of features is extracted that spans the high dimensional feature space. UMAP is applied to determine a representation of the data in 2D and 3D euclidean space, preserving local and global structures of the data and still suited for classification. The performance of UMAP in terms of generating expressive visualizations is evaluated on PolInSAR data acquired by the F-SAR sensor and compared to that of Principal Component Analysis (PCA), Laplacian Eigenmaps (LE) and t-distributed Stochastic Neighbor embedding (t-SNE). For this purpose, a visual analysis of 2D embeddings is performed. In addition, a quantitative analysis is provided for evaluating the preservation of information in low dimensional representations with respect to separability of different land cover classes. The results show that UMAP exceeds the capability of PCA and LE in these regards and is competitive with t-SNE.

A Novel Four-Stage Method for Vegetation Height Estimation with Repeat-Pass PolInSAR Data via Temporal Decorrelation Adaptive Estimation and Distance Transformation

Vegetation height estimation plays a pivotal role in forest mapping, which significantly promotes the study of environment and climate. This paper develops a general forest structure model for vegetation height estimation using polarimetric interferometric synthetic aperture radar (PolInSAR) data. In simple terms, the temporal decorrelation factor of the random volume over ground model with volumetric temporal decorrelation (RVoG-vtd) is first modeled by random motions of forest scatterers to solve the problem of ambiguity. Then, a novel four-stage algorithm is proposed to improve accuracy in forest height estimation. In particular, to compensate for the temporal decorrelation mainly caused by changes between multiple observations, one procedure of temporal decorrelation adaptive estimation via Expectation-Maximum (EM) algorithm is added into the novel method. On the other hand, to extract the features of amplitude and phase more effectively, in the proposed method, we also convert Euclidean distance to a generalized distance for the first time. Assessments of different algorithms are given based on the repeat-pass PolInSAR data of Gabon Lope Park acquired in AfriSAR campaign of German Aerospace Center (DLR). The experimental results show that the proposed method presents a significant improvement of vegetation height estimation accuracy with a root mean square error (RMSE) of 6.23 m and a bias of 1.28 m against LiDAR heights, compared to the results of the three-stage method (RMSE: 8.69 m, bias: 4.81 m) and the previous four-stage method (RMSE: 7.72 m, bias: −2.87 m).

PolInSAR Coherence and Entropy‐Based Hybrid Decomposition Model

AbstractTarget characterization is an essential aspect of polarimetric decomposition. This technique is capable of categorizing polarimetric signatures for different types of targets based on the scattering mechanisms they follow, enabling straightforward physical interpretation of the targets. The geometric anomalies associated with human‐made targets escalate the degree of randomness in the scattering process, which causes scattering ambiguity for such targets. The second‐order model descriptors do not relate to the actual physical structure and yield predominant volume scattering power. Such urban targets are decomposed as natural targets leading to irrelevant decomposition results. The methods developed to curb the problem are unable to maintain the consistency in the decomposition modeling as they underestimate volume scattering powers for natural landcover. A hybrid decomposition model is proposed herein to solve the problem of predominant volume scattering observed from urban targets by preserving volume scattering powers for natural targets. The model uses eigenvalue‐based decomposition parameters and polarimetric interferometric synthetic aperture radar (PolInSAR) coherence to decompose ambiguous targets. The proposed model has been tested on NISAR UAVSAR PolInSAR data acquired over the Greenville region, MS, USA. The proposed model has increased the double‐bounce scattering from the urban targets and enhanced the volume scattering from natural landcover as well. By comparing the results with existing decomposition models, it is observed that the proposed model gives a more robust representation of the landcover than the compared decomposition models.

Corrections to “A Machine-Learning Approach to PolInSAR and LiDAR Data Fusion for Improved Tropical Forest Canopy Height Estimation Using NASA AfriSAR Campaign Data” [Oct 19 3453-3463

This paper investigates the benefits of integrating multibaseline polarimetric interferometric SAR (PolInSAR) data with LiDAR measurements using a machine-learning approach in order to obtain improved forest canopy height estimates. Multiple interferometric baselines are required to ensure consistent height retrieval performance across a broad range of tree heights. Previous studies have proposed multibaseline merging strategies using metrics extracted from PolInSAR measurements. Here, we introduce the multibaseline merging using a support vector machine trained by sparse LiDAR samples. The novelty of this method lies in the new way of combining the two datasets. Its advantage is that it does not require a complete LiDAR coverage, but only sparse LiDAR samples distributed over the PolInSAR image. LiDAR samples are not used to obtain the best height among a set of height stacks, but rather to train the retrieval algorithm in selecting the best height using the variables derived through PolInSAR processing. This enables more accurate height estimation for a wider scene covered by the SAR with only partial LiDAR coverage. We test our approach on NASA AfriSAR data acquired over tropical forests by the L-band UAVSAR and the LVIS LiDAR instruments. The estimated height from this approach has a higher accuracy (r2 = 0.81, RMSE = 7.1 m) than previously introduced multibaseline merging approach (r2 = 0.67, RMSE = 9.2 m). This method is beneficial to future spaceborne missions, such as GEDI and BIOMASS, which will provide a wealth of near-contemporaneous LiDAR samples and PolInSAR measurements for mapping forest structure at global scale.

Unsupervised classification using polarimetric interferometric decomposition and Wishart classifier

ABSTRACTIn this paper, a decomposition scheme of the coherency matrix is presented to parse the information of polarimetric interferometric synthetic aperture radar (PolInSAR) images in detail. First, the decomposition method is improved by the polarimetric interferometric similarity parameter (PISP) to relief the overestimation occurred in the traditional four-component decomposition method. Second, after using the improved four-component decomposition results as the original inputs, the decomposition method is applied to retrieve scattering mechanisms or identify scatters, with the image separated into seven subsections. Finally, based on the modified decomposition results, the basic classification results are regarded as the feature training sets, and the Wishart classifier is then used as an optimized classification process. The applications of the decomposition and classification scheme are shown with typical representative L-band E-SAR images, which are used to show the robustness of the method, as well as with the first published airborne C-band PolInSAR data collected by the Institute of Electronics, Chinese Academy of Sciences, in November 2017. Experimental results demonstrate that the obtained decomposition and classification results are in good agreement with the actual physical scattering mechanisms.

Robust, Model-Based External Calibration of Multi-Channel Airborne SAR Sensors Using Range Compressed Raw Data

The paper describes a method for the accurate calibration of multi-channel SAR instruments, such as those required to support SAR polarimetry, single-pass interferometry and digital beam-forming (DBF), on the basis of dedicated SAR acquisitions containing reference targets with known properties. Unlike conventional approaches, the method is based entirely on the analysis of range-compressed raw data. It leverages the pulse-by-pulse analysis of amplitude, phase and delay variations observed within the range histories of reference targets to fully characterize and correct propagation direction dependent calibration issues such as those related to antenna pointing or antenna phase center positions. The fact that the approach does not require SAR image focusing in azimuth is especially relevant in the context of DBF, where individual channels need to be calibrated but are, by themselves, under-sampled. The calibration techniques presented are illustrated and validated using multi-channel polarimetric and single-pass interferometric SAR data acquired by DLR’s airborne F-SAR and DBFSAR instruments.

Assessing the performance of indicators resulting from three-component Freeman–Durden polarimetric SAR interferometry decomposition at P-and L-band in estimating tropical forest aboveground biomass

ABSTRACT Aboveground Biomass (AGB) estimation performance in a dense tropical forest using Polarimetric Interferometric Synthetic Aperture Radar (PolInSAR) data, at P-and L-band, is evaluated. At high levels of biomass, the backscatter saturation effect leads to a low sensitivity of the backscattered intensity for biomass. This study tries to overcome this problem based on decomposition of PolInSAR data into ground, ground-trunk, and volume components to retrieve the vertical structure information of the forest and polarimetric characteristics of layers. Some sensitive parameters, which extracted from Polarimetric Synthetic Aperture Radar (PolSAR) and PolInSAR data are chosen. Then, many sets of these features are used for assessing biomass estimation using Linear Regression (LR) and Support Vector Machine (SVM) regression models. The data analyzed in this paper are from the P-and L-band airborne dataset acquired by Office National d’Études et de Recherches Aérospatiales (ONERA) over French Guiana in 2009, in the frame of the European Space Agency (ESA) campaign, TropiSAR. The average forest biomass is 374.54 t ha−1 and goes up to 503 t ha−1 for the in-situ plots. The results display that the indicator of double-bounce (C DB) contribution matrix, extracted from PolInSAR data decomposition, is the highest correlated feature to AGB, at both P-and L-band. The combination of PolInSAR indicators improved the Root-Mean-Squared Error (RMSE) values up to 24.76 t ha−1 at P-band and 18.02 t ha−1 at L-band. At the best state, AGB estimation RMSE reduces to 40.95 t ha−1 (12.73%) in simultaneous use of the features derived from PolSAR and PolInSAR data, at P-band. SVM produces a more realistic biomass value than LR. However, cross-validated AGB RMSE of two regression models is close to each other.

Polarimetric Interferometric SAR Change Detection Discrimination

A coherent change detection (CCD) image, computed from a geometrically matched, temporally separated pair of complex-valued synthetic aperture radar (SAR) image sets, conveys the pixel-level equivalence between the two observations. Low-coherence values in a CCD image are typically due to either some physical change in the corresponding pixels or a low signal-to-noise observation. A CCD image does not directly convey the nature of the change that occurred to cause low coherence. In this paper, we introduce a mathematical framework for discriminating between different types of change within a CCD image. We utilize the extra degrees of freedom and information from polarimetric interferometric SAR (PolInSAR) data and PolInSAR processing techniques to define a 29-dimensional feature vector that contains information capable of discriminating between different types of change in a scene. We also propose two change-type discrimination functions that can be trained with feature vector training data and demonstrate change-type discrimination on an example image set for three different types of change. Furthermore, we also describe and characterize the performance of the two proposed change-type discrimination functions by way of receiver operating characteristic curves, confusion matrices, and pass matrices.

Vegetation Height Estimate in Rice Fields Using Single Polarization TanDEM-X Science Phase Data

This study presents the retrieval of rice paddy height using single polarization (single-pol) interferometric SAR (InSAR) data by means of model-based height inversion without an external DEM. A total of eight TanDEM-X (TDX) scenes were used and the TDX images were; acquired using a large cross-track baseline configuration during the TDX Science Phase (June–August 2015). A single-pol inversion approach for a flooded rice field is proposed and evaluated over the Buan test site in South Korea. A novel approach is adopted for the estimation of the ground (i.e., water level) interferometric phase within a flooded rice paddy from TDX data acquired during a period of early rice growth. It is consequently possible to apply the model-based inversion algorithm to the single-pol InSAR data. Rice height maps during the rice growth cycle presented and validated by field measurements. The results demonstrated the high performance of the inversion with a correlation coefficient of 0.78 and an RMSE of 0.10 m. The proposed methodology will be useful to monitor rice plants and to predict a gross rice yield, along with dual and fully polarimetric interferometric SAR data.