Advances in Radar Remote Sensing of Agricultural Crops: A Review

There are enormous advantages of a review article in the field of emerging technology like radar remote sensing applications in agriculture. This paper aims to report select recent advancements in the field of Synthetic Aperture Radar (SAR) remote sensing of crops. In order to make the paper comprehensive and more meaningful for the readers, an attempt has also been made to include discussion on various technologies of SAR sensors used for remote sensing of agricultural crops viz. basic SAR sensor, SAR interferometry (InSAR), SAR polarimetry (PolSAR) and polarimetric interferometry SAR (PolInSAR). The paper covers all the methodologies used for various agricultural applications like empirically based models, machine learning based models and radiative transfer theorem based models. A thorough literature review of more than 100 research papers indicates that SAR polarimetry can be used effectively for crop inventory and biophysical parameters estimation such are leaf area index, plant water content, and biomass but shown less sensitivity towards plant height as compared to SAR interferometry. Polarimetric SAR Interferometry is preferable for taking advantage of both SAR polarimetry and SAR interferometry. Numerous studies based upon multi-parametric SAR indicate that optimum selection of SAR sensor parameters enhances SAR sensitivity as a whole for various agricultural applications. It has been observed that researchers are widely using three models such are empirical, machine learning and radiative transfer theorem based models. Machine learning based models are identified as a better approach for crop monitoring using radar remote sensing data. It is expected that the review article will not only generate interest amongst the readers to explore and exploit radar remote sensing for various agricultural applications but also provide a ready reference to the researchers working in this field

Dual polarimetric radar vegetation index for crop growth monitoring using sentinel-1 SAR data

Sentinel-1 Synthetic Aperture Radar (SAR) data have provided an unprecedented opportunity for crop monitoring due to its high revisit frequency and wide spatial coverage. The dual-pol (VV-VH) Sentinel-1 SAR data are being utilized for the European Common Agricultural Policy (CAP) as well as for other national projects, which are providing Sentinel derived information to support crop monitoring networks. Among the Earth observation products identified for agriculture monitoring, indicators of vegetation status are deemed critical by end-user communities. In literature, several experiments usually utilize the backscatter intensities to characterize crops. In this study, we have jointly utilized the scattering information in terms of the degree of polarization and the eigenvalue spectrum to derive a new vegetation index from dual-pol (DpRVI) SAR data. We assess the utility of this index as an indicator of plant growth dynamics for canola, soybean, and wheat, over a test site in Canada. A temporal analysis of DpRVI with crop biophysical variables (viz., Plant Area Index (PAI), Vegetation Water Content (VWC), and dry biomass (DB)) at different phenological stages confirms its trend with plant growth dynamics. For each crop type, the DpRVI is compared with the cross and co-pol ratio (σVH0/σVV0) and dual-pol Radar Vegetation Index (RVI = 4σVH0/(σVV0 + σVH0)), Polarimetric Radar Vegetation Index (PRVI), and the Dual Polarization SAR Vegetation Index (DPSVI). Statistical analysis with biophysical variables shows that the DpRVI outperformed the other four vegetation indices, yielding significant correlations for all three crops. Correlations between DpRVI and biophysical variables are highest for canola, with coefficients of determination (R2) of 0.79 (PAI), 0.82 (VWC), and 0.75 (DB). DpRVI had a moderate correlation (R2≳ 0.6) with the biophysical parameters of wheat and soybean. Good retrieval accuracies of crop biophysical parameters are also observed for all three crops.This work was supported by the Spanish Ministry of Science, Innovation and Universities, the State Agency of Research (AEI) and the European Funds for Regional Development (EFRD) under Project TEC2017-85244-C2-1-P

Crop monitoring and yield estimation using polarimetric SAR and optical satellite data in southwestern Ontario

Optical satellite data have been proven as an efficient source to extract crop information and monitor crop growth conditions over large areas. In local- to subfield-scale crop monitoring studies, both high spatial resolution and high temporal resolution of the image data are important. However, the acquisition of optical data is limited by the constant contamination of clouds in cloudy areas. This thesis explores the potential of polarimetric Synthetic Aperture Radar (SAR) satellite data and the spatio-temporal data fusion approach in crop monitoring and yield estimation applications in southwestern Ontario. Firstly, the sensitivity of 16 parameters derived from C-band Radarsat-2 polarimetric SAR data to crop height and fractional vegetation cover (FVC) was investigated. The results show that the SAR backscatters are affected by many factors unrelated to the crop canopy such as the incidence angle and the soil background and the degree of sensitivity varies with the crop types, growing stages, and the polarimetric SAR parameters. Secondly, the Minimum Noise Fraction (MNF) transformation, for the first time, was applied to multitemporal Radarsat-2 polarimetric SAR data in cropland area mapping based on the random forest classifier. An overall classification accuracy of 95.89% was achieved using the MNF transformation of the multi-temporal coherency matrix acquired from July to November. Then, a spatio-temporal data fusion method was developed to generate Normalized Difference Vegetation Index (NDVI) time series with both high spatial and high temporal resolution in heterogeneous regions using Landsat and MODIS imagery. The proposed method outperforms two other widely used methods. Finally, an improved crop phenology detection method was proposed, and the phenology information was then forced into the Simple Algorithm for Yield Estimation (SAFY) model to estimate crop biomass and yield. Compared with the SAFY model without forcing the remotely sensed phenology and a simple light use efficiency (LUE) model, the SAFY incorporating the remotely sensed phenology can improve the accuracy of biomass estimation by about 4% in relative Root Mean Square Error (RRMSE). The studies in this thesis improve the ability to monitor crop growth status and production at subfield scale

Crop biophysical parameter retrieval from Sentinel-1 SAR data with a multi-target inversion of Water Cloud Model

Estimation of bio-and geophysical parameters from Earth observation (EO) data is essential for developing applications on crop growth monitoring. High spatio-temporal resolution and wide spatial coverage provided by EO satellite data are key inputs for operational crop monitoring. In Synthetic Aperture Radar (SAR) applications, a semi-empirical model (viz., Water Cloud Model (WCM)) is often used to estimate vegetation descriptors individually. However, a simultaneous estimation of these vegetation descriptors would be logical given their inherent correlation, which is seldom preserved in the estimation of individual descriptors by separate inversion models. This functional relationship between biophysical parameters is essential for crop yield models, given that their variations often follow different distribution throughout crop development stages. However, estimating individual parameters with independent inversion models presume a simple relationship (potentially linear) between the biophysical parameters. Alternatively, a multi-target inversion approach would be more effective for this aspect of model inversion compared to an individual estimation approach. In the present research, the multi-output support vector regression (MSVR) technique is used for inversion of the WCM from C-band dual-pol Sentinel-1 SAR data. Plant Area Index (PAI, m2 m−2) and wet biomass (W, kg m−2) are used as the vegetation descriptors in the WCM. The performance of the inversion approach is evaluated with in-situ measurements collected over the test site in Manitoba (Canada), which is a super-site in the Joint Experiment for Crop Assessment and Monitoring (JECAM) SAR inter-comparison experiment network. The validation results indicate a good correlation with acceptable error estimates (normalized root mean square error–nRMSE and mean absolute error–MAE) for both PAI and wet biomass for the MSVR approach and a better estimation with MSVR than single-target models (support vector regression–SVR). Furthermore, the correlation between PAI and wet biomass is assessed using the MSVR and SVR model. Contrary to the single output SVR, the correlation between biophysical parameters is adequately taken into account in MSVR based simultaneous inversion technique. Finally, the spatio-temporal maps for PAI and W at different growth stages indicate their variability with crop development over the test site.This research was supported in part by Shastri Indo-Candian Institute, New Delhi, India and the Spanish Ministry of Economy, Industry and Competitiveness, in part by the State Agency of Research (AEI), in part by the European Funds for Regional Development under project TEC2017-85244-C2-1-P

Generative feature extraction from sentinel 1 and 2 data for prediction of forest aboveground biomass in the Italian Alps

—Aboveground biomass (AGB) is an important forest attribute directly linked to the forest carbon pool. The use of satellite remote sensing (RS) data has increased for AGB prediction due to their large footprint and low-cost availability. However, they have been limited due to saturation effect that leads to low prediction precision. In this article, we propose an innovative and dynamic architecture based on generative neural network that extracts target oriented generative features for predicting forest AGB using satellite RS data. These features are more resilient to mixed forest types and geographical conditions as compared to the traditional features and models. The effectiveness of the proposed features was assessed by experiments performed using multispectral, synthetic aperture radar, and combined dual-source datasets. The proposed model achieved best performance in terms of precision, model agreement, and overfitting as compared to the other conventional models for all analyzed datasets. The t-distributed stochastic neighbor embedding scatterplots of the generative features clearly show one dimension of the feature space associated with the target AGB. Feature importance analysis indicated that the produced generative features were more significant than the conventional analytical features. Also, the model provided a robust framework for homogeneous fusion of multisensor features from satellite RS data for predicting AG

Spatial Adaptive Speckle Filtering Driven by Temporal Polarimetric Statistics and Its Application to PSI

Persistent scatterer (PS) interferometry (PSI) techniques are designed to measure ground deformations using satellite synthetic aperture radar (SAR) data. They rely on the identification of pixels not severely affected by spatial or temporal decorrelation, which, in general, correspond to pointlike PSs commonly found in urban areas. However, in urban areas, we can find not only PSs but also distributed scatterers (DSs) whose phase information may be exploited for PSI applications. Estimation of DS parameters requires speckle filtering to be applied to the complex SAR data, but conventional speckle filtering approaches tend to mask PS information due to spatial averaging. In the context of single-polarization PSI, adaptive speckle filtering strategies based on the exploitation of amplitude temporal statistics have been proposed, which seek to avoid spatial filtering on nonhomogeneous areas. Given the growing interest on polarimetric PSI techniques, i.e., those using polarimetric diversity to increase performance over conventional single-polarization PSI, in this paper, we propose an adaptive spatial filter driven by polarimetric temporal statistics, rather than single-polarization amplitudes. The proposed approach is able to filter DS while preserving PS information. In addition, a new methodology for the joint processing of PS and DS in the context of PSI is introduced. The technique has been tested for two different urban data sets: 41 dual-polarization TerraSAR-X images of Murcia (Spain) and 31 full-polarization Radarsat-2 images of Barcelona (Spain). Results show an important improvement in terms of number of pixels with valid deformation information, hence denser area coverage.This work was supported in part by the Spanish Ministerio de Economía y Competitividad and in part by the European Union FEDER funds under Project TEC2011-28201-C02-02

Monitoring Rice Phenology Based on Backscattering Characteristics of Multi-Temporal RADARSAT-2 Datasets

Accurate estimation and monitoring of rice phenology is necessary for the management and yield prediction of rice. The radar backscattering coefficient, one of the most direct and accessible parameters has been proved to be capable of retrieving rice growth parameters. This paper aims to investigate the possibility of monitoring the rice phenology (i.e., transplanting, vegetative, reproductive, and maturity) using the backscattering coefficients or their simple combinations of multi-temporal RADARSAT-2 datasets only. Four RADARSAT-2 datasets were analyzed at 30 sample plots in Meishan City, Sichuan Province, China. By exploiting the relationships of the backscattering coefficients and their combinations versus the phenology of rice, HH/VV, VV/VH, and HH/VH ratios were found to have the greatest potential for phenology monitoring. A decision tree classifier was applied to distinguish the four phenological phases, and the classifier was effective. The validation of the classifier indicated an overall accuracy level of 86.2%. Most of the errors occurred in the vegetative and reproductive phases. The corresponding errors were 21.4% and 16.7%, respectively

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.This research was funded in part by the National Natural Science Foundation of China (grant nos. 41820104005, 42030112, 41904004) and in part by the and the Spanish Ministry of Science and Innovation (grant no. PID2020-117303GB-C22)

SAR Amplitude Probability Density Function Estimation Based on a Generalized Gaussian Model

International audienceIn the context of remotely sensed data analysis, an important problem is the development of accurate models for the statistics of the pixel intensities. Focusing on synthetic aperture radar (SAR) data, this modeling process turns out to be a crucial task, for instance, for classification or for denoising purposes. In this paper, an innovative parametric estimation methodology for SAR amplitude data is proposed that adopts a generalized Gaussian (GG) model for the complex SAR backscattered signal. A closed-form expression for the corresponding amplitude probability density function (PDF) is derived and a specific parameter estimation algorithm is developed in order to deal with the proposed model. Specifically, the recently proposed “method-of-log-cumulants” (MoLC) is applied, which stems from the adoption of the Mellin transform (instead of the usual Fourier transform) in the computation of characteristic functions and from the corresponding generalization of the concepts of moment and cumulant. For the developed GG-based amplitude model, the resulting MoLC estimates turn out to be numerically feasible and are also analytically proved to be consistent. The proposed parametric approach was validated by using several real ERS-1, XSAR, E-SAR, and NASA/JPL airborne SAR images, and the experimental results prove that the method models the amplitude PDF better than several previously proposed parametric models for backscattering phenomena