Next Article in Journal
An Occlusion-Robust Feature Selection Framework in Pedestrian Detection
Next Article in Special Issue
Evaluation of Atmospheric Effects on Interferograms Using DEM Errors of Fixed Ground Points
Previous Article in Journal
Evaluating Geometric Measurement Accuracy Based on 3D Reconstruction of Automated Imagery in a Greenhouse
Previous Article in Special Issue
Deceptive Jamming Detection for SAR Based on Cross-Track Interferometry
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Inversion of Rice Biophysical Parameters Using Simulated Compact Polarimetric SAR C-Band Data

1
College of Geomatics, Shandong University of Science and Technology, Qingdao 266590, China
2
Laboratory of Radar Remote Sensing Application Technology, Institute of Remote Sensing and Digital Earth (RADI), Chinese Academy of Sciences, Datun Road, Chaoyang District, Beijing 100101, China
3
Laboratory of Target Microwave Properties (LAMP), Zhongke Academy of Satellite Application in Deqing (DASA), Deqing 313200, China
4
School of Earth Science and Resources, China University of Geosciences (Beijing), Beijing 100083, China
5
Institute of Transmission and Transformation Engineering, China Electric Power Research Institute, Beijing 100055, China
*
Authors to whom correspondence should be addressed.
Sensors 2018, 18(7), 2271; https://doi.org/10.3390/s18072271
Submission received: 24 May 2018 / Revised: 10 July 2018 / Accepted: 11 July 2018 / Published: 13 July 2018

Abstract

:
Timely and accurate estimation of rice parameters plays a significant role in rice monitoring and yield forecasting for ensuring food security. Compact-polarimetric (CP) synthetic aperture radar (SAR), a good compromise between the dual- and quad-polarized SARs, is an important part of the new generation of Earth observation systems. In this paper, the ability of CP SAR data to retrieve rice biophysical parameters was explored using a modified water cloud model. The results showed that S1 was superior to other CP variables in rice height inversion with a coefficient of determination (R2) of 0.92 and a root-mean-square error (RMSE) of 5.81 cm. RL was the most suitable for inverting the volumetric water content of the rice canopy, with an R2 of 0.95 and a RMSE of 0.31 kg/m3. The m-χ decomposition produced the highest accuracies for the ear biomass: R2 was 0.89 and RMSE was 0.17 kg/m2. The highest accuracy of leaf area index (LAI) retrieval was obtained for RH (right circular transmit and horizontal linear receive) with an R2 of 0.79 and a RMSE of 0.33. This study illustrated the capability of CP SAR data with respect to retrieval of rice biophysical parameters, especially for height, volumetric water content of the rice canopy, and ear biomass, and this mode may offer the best option for rice-monitoring applications because of swath coverage.

1. Introduction

Rice, one of the world’s leading food crops, occupies more than 11 percent of the world’s total arable land and provides food for more than half of the world’s population [1]. Rice biophysical parameters, such as the biomass, the crop height, the leaf area index (LAI), and the water content, are vital indicators of rice growth condition. For example, the dry biomass of rice ears is a direct manifestation of the rice yield. Rice crop height, an important index of rice phenology, is a comprehensive reflection of environmental factors, such as the climate, hydrology, and soil. The LAI is a key parameter in plant population and community growth analysis, and volumetric water content (VWC) can fully characterize any vegetation water shortage, both of which are closely related to vegetation growth. Therefore, timely and accurate estimation of rice biophysical properties is of great significance and consequently plays a role in forecasting rice productivity for ensuring regional or national food security.
Synthetic aperture radar (SAR) has been proven to have a great potential for rice monitoring [2,3,4] and inversion of rice biophysical parameters [5,6] because of its all-weather day–night imaging capability, its ability to penetrate microwave signals, and the sensitivity of its signals to target properties. Erten et al. [5] dealt with the retrieval of agricultural rice height from space by using multi-polarization SAR images, and coherent and incoherent crop height estimation methods were discussed for the first time. Tan et al. [6] built a biomass inversion algorithm and inverted rice biomass, leaf area index, canopy height and related eco-physiological canopy variables, and the inversion results showed that the multitemporal SAR images at the C-band can be used to monitor the growth of rice. Although the application effect is very good, the swath coverage is small, and it is difficult to meet a wide range of application requirements. Compact-polarimetric (CP) SAR transmits on only one polarization and receives on two orthogonal polarizations, retaining their relative phase. The major motivation for compact polarimetry is to strive for quantitative applications of the same level as those from a fully polarized system, while avoiding the principal disadvantages (mass, power, and limited coverage) associated with a quad-polarized (quad-pol) SAR [7]. As a good compromise between the dual-polarized and quad-pol SARs, CP SAR is an important part of the new generation of Earth observation systems (EOS). RISAT-1, the first radar satellite with CP measurement capability, was successfully launched in April 2012 [8]. ALOS-2, launched in 2014, also includes CP as in the experimental imaging model in 2014 [9]. Moreover, in near future the RADARSAT Constellation Mission (RCM) of Canada will launch a C-band CP SAR satellite as well [10]. With the emergence of CP SARs of EOS, several studies have reported on the applications of CP SAR in rice monitoring such as rice mapping [11,12,13] and rice phenology [14,15,16].
However, the potential of CP SAR for the inversion of rice biophysical parameters has not been fully explored. Kumar et al. [17] found that the CP parameters are sensitive to the growth changes of wheat and corn, but did not look at the relationship between CP parameters and biophysical parameters. Zhang et al. [18] not only studied the response of rape to C-band CP SAR, they also carried out the inversion of rape biophysical parameters based on the random forest method, thus demonstrating the potential of CP SAR for the inversion of crop biophysical parameters. However, the physical mechanism of the crop was not taken into account by the inversion strategy. H. Yang et al. [19] devised a scattering index based on the m-χ decomposition components of simulated CP SAR data and constructed a regression model to obtain inversion results for oilseed rape biomass. The above studies indicated the capability of CP SAR in crop biophysical parameters inversion. However, the potential of CP SAR for the inversion of rice biophysical parameters has not been fully explored. Therefore, in this study, five temporal RADARSAT-2 fully polarimetric single look complex (SLC) C-band datasets were used to simulate CP SAR data, and 14 CP parameters, such as backscattering coefficients, Stokes vectors, and decomposition parameters, were extracted in each acquisition. Then, the capability of CP SAR data for the inversion of rice biophysical parameters was fully investigated by introducing a classic semi-empirical model, the water cloud model (WCM) [20], and a modified WCM (MWCM) in which the heterogeneity of the rice canopy and its phenological changes, as well as the double-bounce scattering between the rice canopy and the underlying surface were considered.

2. Study Site, Data Source, and Ground Measurement Campaign

The study site was located at Jinhu, Xuyi, and Hongze, Jiangsu province (118°41′34′′~119°16′27′′ E, 33°17′05′′~33°56′39′′ N). Lacustrine plain is the main terrain type there, and so, the terrain is flat with an altitude of 5–9 m. The climate type is humid subtropical monsoon with abundant sunshine and rainfall. The site is one of the major grain-producing areas in the east of China, with rice being the main crop. The farmland plots in this area are large and produce one rice crop per year. The rice is at the growing stage from mid-June to late October or early November. The main types of rice in the region are hybrid and japonica, and the planting method used is transplanting and sowing. In this study, hybrid rice in a transplanted rice paddy was used. The phenological changes of indica rice are shown in Figure 1.
As the variations of biophysical parameters change with the development of the rice plants, five temporal RADARSAT-2 fully polarimetric single look complex (SLC) C-band datasets (Table 1), selected from 21 July to 15 October 2012, were used to simulate CP SAR data.
During the passage of the RADARSAT-2 satellite, field work was carried out. Thirty transplanted hybrid rice fields were selected as the sample plots for measuring the rice parameters. Each rice was recorded with an accurate global positioning system (GPS). Secondly, the rice parameters were measured from three rice plants, which were randomly selected in each rice field. The sizes of the rice plants (such as plant height, stem length and diameter, leaf length and width, ear length and width, etc.), the planting density, and other structural parameters involved in MWCM were measured with a steel ruler and a Vernier caliper. The LAI was measured using a SunScan Canopy analysis system [22]. The biomass and water content were obtained by destructive measurements of the sampled rice plants based on normal gravimetric processes.
If volume scattering is the predominant mechanism responsible for the backscatter from vegetation, it seems appropriate to model a vegetation canopy as a cloud containing a volumetric water content (VWC) mv (kg/m3), which is defined as the water content per unit volume of the rice canopy [23]. Here, considering the heterogeneity of the rice canopy, the VWC of each layer in rice canopy, such as ear, leaf, stem, and so on, were proposed and calculated according to the definition of VWC:
m v _ e =   ( F e - D e ) N / h e m v _ s =   ( F s - D s ) N / h s m v _ l =   ( F l - D l ) N / h l m v =   ( F - D ) N / h }  
where N (plants·m−2) is the density of the rice plants; mv_e, mv_s, mv_l, and mv are the VWC (ear, stem, leaf, and total, respectively); F (kg) is the fresh weight of the total plant, and Fe, Fs, and Fl (kg) are the fresh weight of the ear, stem, and leaf of each rice plant, respectively; D (kg) is the dry weight of the total plant, and De, Ds, and Dl (kg) are the dry weight of the ear, stem, and leaf of each rice plant, respectively; the difference between Fi and Di (i = e, s, l) indicates the water content corresponding to each part of a rice plant, and he, hs, hl, and h (m) are the height (ear, stem, leaf, and total).

3. Methodology

The flow chart of the methodology adopted in this study is shown in Figure 2. This methodology consisted mainly of the CP SAR data simulation, CP parameter extraction, data preprocessing, the inversion model building using the WCM/MWCM model, the genetic algorithm (GA), and finally, validation and analysis.

3.1. CP SAR Data Simulation

The CP SAR data we used were simulated using the RADARSAT-2 full polarimetric (FP) SAR data as follows [24]:
(1) The sensor emits right-hand circularly polarized microwave radiation, denoted R. Based on the Sinclair matrix for the FP SAR data [ Γ ], the electric vector, EB, of the target scattering field is constructed as follows:
E B = [ Γ ] R , R = ( 1 2 ) [ 1 , j ] T ,
E B = [ Γ ] R = ( 1 2 ) [ S H H j S H V , S H V j S V V ] T ,
where SHH, SHV, and SVV are the elements of the Sinclair matrix for the FP SAR data.
(2) The EH and EV components are calculated in the case where the antenna sends and receives linear polarization (H and V) signal as follows:
E H = [ 1 0 ] E B = ( 1 2 ) ( S H H j S H V ) ,
E V = [ 0 1 ] E B = ( 1 2 ) ( S H V j S V V ) .
(3) The four elements of the coherence matrix J are calculated as follows:
2 J H H = | S H H | 2 + | S H V | 2 + j S H H S H V j S H V S H H ,
2 J H V = S H H S H V S H V S V V j | S H V | 2 j S H H S V V ,
J V H = J H V * ,
2 J V V = | S V V | 2 + | S H V | 2 j S V V S H V + j S H V S V V .
(4) Based on the elements of matrix J, the Stokes vector, S, is calculated as follows:
S 1 = J H H + J V V ,
S 2 = J H H J V V ,
S 3 = Re { S H H S H V + S H V S V V } Im S H H S V V ,
S 4 = Im { S H H S H V S H V S V V } Re S H H S V V + | S H V | 2 ,
where S1, S2, S3, and S4 are the four elements of the Stokes vector.
(5) According to the relation between the Sinclair matrix [ Γ ] and the covariance matrix, C3, a relation between the matrix elements of the FP data and the Stokes vector for the CP mode is established as follows:
S 1 = 1 / 2 C 11 + 1 / 2 C 22 + 1 / 2 C 33 + ( 1 / 2 ) Im C 12 + ( 1 / 2 ) Im C 23 ,
S 2 = 1 / 2 C 11 1 / 2 C 33 + ( 1 / 2 ) Im C 12 ( 1 / 2 ) Im C 23 ,
S 3 = ( 1 / 2 ) Re C 12 + ( 1 / 2 ) Re C 23 Im C 13 ,
S 4 = Re C 13 1 / 2 C 22 ( 1 / 2 ) Im C 12 + ( 1 / 2 ) Re C 23 .
The simulated CP SAR data were transmitted as right-circular polarization (R) and received as horizontal (H) and vertical (V) polarizations, with a resolution of 30 m and a noise floor of −25 dB. Figure 3 shows the color composite images of the CP backscattering coefficients generated from the simulated CP SAR data.

3.2. Extraction of CP SAR Parameters and Data Preprocessing

Following the process described in Section 3.1, the simulated CP SAR data was stored in Stokes vector S = {S1, S2, S3, S4}. A total of 14 CP parameters, including Stokes vectors, four backscattering coefficients, three m-δ decomposition parameters, and three m–χ decomposition parameters, were extracted in each acquisition, as described in Section 2. A 5 × 5 frost filter was applied to all the CP parameters, and the resulting images were georectified. Then, image subsets were taken within the boundaries of the study site, and the areas designated as urban, water, road, forest, and other land cover types were removed. Only the areas classified as hybrid rice were extracted for the inversion of the rice biophysical parameters.
The four Stokes parameters represent the intensity and polarization state of the radar echoes. S1 and S2 are closely related to the overall backscatter. Figure 4 shows the temporal behaviors of the four Stokes parameters from the elongation stage to the mature stage. It can be seen that the energy is mainly concentrated in the S1 and S2 components, whereas the contributions of S3 and S4 are very small. Although there were significant changes in the rice biophysical parameters during rice growth, the energy values of S3 and S4 were close to zero and did not change much during each rice phenology stage. It can be seen that S3 and S4 are not sensitive to the changes in rice growth and biophysical parameters. Therefore, S3 and S4 are not suitable for the inversion of the rice biophysical parameters. For S1 and S2, from the elongation stage to the booting stage, the rice height increased, and the water content of the rice also increased. In addition, the leaves widened and enlarged. Moreover, the values of S1 and S2 also closely reached the maximum values. However, from the booting stage to the heading stage, the dough stage, and the mature stage, the appearance of the rice ear and the change in leaf color and the underlay surface of the rice resulted in the decrease of the value of S1 and S2. Therefore, S1 and S2 deserve to be used to invert the rice biophysical parameters.
CP backscattering parameters have been proven to have a great potential in crop parameters’ inversion [18,19]. Four CP backscattering coefficients were simulated respectively in the RH (right circular transmit and horizontal linear receive), RV (right circular transmit and horizontal linear receive), RR (right circular transmit and right circular receive), and RL (right circular transmit and left circular receive) from the fully-polarized Stokes vector S. The calculation details are described by Equations (18)–(21) as follows:
σ RH 0 = [ 1 , 1 , 0 , 0 ] × S ,
σ RV 0 = [ 1 , 1 , 0 , 0 ] × S ,
σ RR 0 = [ 1 , 0 , 0 , 1 ] × S ,
σ RL 0 = [ 1 , 0 , 0 , 1 ] × S .
Figure 5 shows the temporal behaviors of the CP backscattering coefficients in the RH, RV, RR, and RL polarizations from the elongation stage to the mature stage. All the backscatter powers were high values at the elongation stage. The backscattering coefficients of RV dropped sharply in the booting stage, and the other three backscatter values showed little difference. The increase in rice height resulted in high attenuation of vertically-polarized return. From the booting stage to the heading, dough, and mature stages, the backscattering coefficients of RV had declined slowly. This was because the rice had grown to a certain height so that degree of attenuation of vertically-polarized return had not changed much. From the booting stage to the heading and dough stages, the backscattering coefficients of RH\RR\RV declined dramatically, and the increase in the rice leave layer and the appearance of the rice ear layer resulted in high attenuation of horizontal-polarized return. Moreover, as the leaves withered and turned yellow gradually, the volume scattering and double-bounce scattering were decreased, and the backscattering power of RH\RR\RV decreased gradually as well. With rice leaves continuing to fall off and wither, although the volume scattering decreased, the double-bounce scattering surface scattering from the underlying surface increased. Therefore, from the dough to the mature stage, all the backscatter powers increased slowly with a gentle slope. All the backscattering coefficients were excellent indicators for different phenological stages of rice and also had great potential for the inversion of the rice biophysical parameters.
The three scattering components of the m-χ and m-δ decomposition, the physical meanings of which were similar to that of the Freeman–Durden decomposition for full polarimetric SAR, were generated from the simulated CP SAR data as described by Equations (22)–(26) [25] as follows:
m = ( S 2 2 + S 3 2 + S 4 2 ) / S 1 ,
sin 2 χ = S 4 / ( m S 1 ) ,
δ = arg ( S 3 + j S 4 ) ,
where m is the degree of polarization, δ is the relative phase, and χ is the ellipticity angle of the Poincaré sphere. They are all the child parameters of the Stokes vector. The degree of polarization m is negatively correlated with the scattering entropy; the greater the scattering entropy, the smaller the m. The phase angle δ, which is the relative phase angle between the RH and RV polarizations in CP SAR data, includes the phase information of the backscatter of the targets. The sign of the ellipticity angle χ is an unambiguous indicator of even versus odd bounce backscatter, and sin2χ is known formally as the degree of circularity [26].
The Stokes parameter S1, which is the total backscattered energy received by the CP systems, in combination with the CP parameters m, δ, and χ were used to generate the m-δ (as shown in Equation (25)) and m-χ decompositions (as shown in Equation (26)), respectively. Pd, Pv, and Ps correspond to the double-bounce, the randomly polarized constituent (volume), and single-bounce (and Bragg) backscattering.
[ P d P v p s ] m - δ = [ m S 1 ( 1 + sin δ ) / 2 S 1 ( 1 m ) m S 1 ( 1 sin δ ) / 2 ] ,
[ P d P v p s ] m - χ = [ m S 1 ( 1 sin 2 χ ) / 2 S 1 ( 1 m ) m S 1 ( 1 + sin 2 χ ) / 2 ] .
Figure 6 shows the temporal behaviors of the three scattering components of the m-χ and m-δ decomposition from the elongation stage to the mature stage. It can be seen that the main scattering mechanism for the reflected radar signal is multiple scattering (volume scattering), which is much higher than the even scattering (double-bounce scattering) and single scattering (surface scattering). As the rice crop was in the elongation stage on 21 July, it was short, and the underlay surface was mostly water, which prompted the three kinds of scattering to be weak. With the growth of the rice, the rice biophysical parameters, such as rice height, biomass, and water content, increased on 4 August (booting stage), which resulted in the increase of the three kinds of scattering. However, from heading stage to the dough and the mature stages, with the emergence, growth, and maturity of the rice ear, the leaves turned yellow gradually, the rice water content decreased, and the rice biomass increased, which made the three kinds of scattering energy decline. It can be seen that with the growth of the rice and the change of the rice biophysical parameters, the three scattering energies of the rice also changed to varying degrees. Therefore, the scattering parameters are very sensitive to the rice biophysical parameters in the rice growth, which provides a possibility for building a model relationship between the rice biophysical parameters and the CP decomposition parameters.

3.3. The Inversion Model of Rice Parameters Building with CP SAR

The CP parameters, including the backscattering coefficients, the Stokes parameters, and the three scattering components of the m-χ and m-δ decomposition, was investigated for the inversion of the rice biophysical parameters by introducing a classic semi-empirical model, the water cloud model (WCM) [20,23], and a modified WCM (MWCM) in which the heterogeneity of the rice canopy and its phenological changes, as well as the double-bounce scattering between the rice canopy and the underlying surface, were considered [27].

3.3.1. The Inversion Model of Rice Parameters Using CP Backscatter and WCM

The WCM is one of the most popular semi-empirical models. It was first proposed by Attema [23] and Ulaby et al. [20] and was extensively applied to the inversion problems of radar remote sensing [28,29,30]. It treats the canopy as a water cloud, consisting of a collection of identical water particles, characterized by a uniform scattering phase function [23]. Ignoring second-order contributions resulting from multiple scattering between the canopy particles and the soil surface, the backscattering coefficient of the canopy is given by the following equation:
σ c a n 0 ( θ ) = σ v e g 0 ( θ ) + σ s 0 ( θ ) ,
where σ v e g 0 ( θ ) is the contribution of the vegetation volume, σ s 0 ( θ ) is the backscattering contribution of the soil surface in the presence of vegetation cover, and θ is the angle of incidence relative to nadir [20].
With respect to the rice canopy, as there is a large difference of volumetric water content between the leaf and stem layer in the rice canopy, it was assumed to consist of two layers: an upper layer of height h1, dominated by leaves, and a lower layer of height h2, dominated by stem. We established the scattering energy from each layer of the rice canopy in the framework of the WCM to express the rice biophysical parameters:
σ c a n 0 ( θ ) = σ l 0 ( θ ) + σ s t 0 ( θ ) + σ s 0 ( θ ) ,
σ l 0 ( θ ) = A l · ( 1 exp ( B l L / h 1 ) ) cos ( θ ) · ( 1 γ l 2 ( θ ) ) ,
σ s t 0 ( θ ) = A s t ( θ ) m v · h 2 · γ l 2 ( θ ) ,
σ s 0 ( θ ) = C S ( θ ) m s · γ s t 2 ( θ ) γ l 2 ( θ ) ,
γ l 2 ( θ ) = exp ( 2 α l L sec ( θ ) ) ,
γ s t 2 ( θ ) = exp ( α s t m v h 2 ) ,
where θ is the incidence angle of the CP SAR data. The parameters with subscript l are the parameters related to the leaves layer of the rice, and the parameters with subscript st are the parameters related to the stem layer of the rice. Moreover, the parameters with subscript s are the parameters related to the soil layer of rice. Because h1 and h2 are proportional to the total plant height h, and h1 and h2 were replaced with h in Equations (29), (30), and (33). L is the LAI, and γ is the transmission coefficient. However, the WCM treats the canopy as a water cloud, consisting of a collection of identical water particles, characterized by a uniform scattering phase function. So, the VWC of the rice canopy (mv) is identical to that of the stem layer (mv_s). It can be seen that Equation (28) contains six model coefficients (Ast, Al, Bl Cs, αst, and αl) and four rice biophysical parameters (L, h, ms and mv).
As the CP backscattering coefficients and Stokes parameters are the power of the backscatter wave, they can replace σ c a n 0 and be expressed as a function of the rice biophysical parameters by substituting Equations (28)–(33) into Equation (27). For each equation for σ c a n 0 ( θ ) , there is an ill-posed problem. Therefore, it is obviously impossible to solve for four rice biophysical parameters in one equation. Therefore, to solve for one of the parameters, the values of the other three parameters must be assumed to ensure that one equation can be used to calculate the value of a rice biophysical parameter. Taking the seedling stage to the heading stage as an example, we take h as the parameter to be solved, and the other three parameters must be assumed. From the seedling stage to the heading stage, the underlying surface is calm water, so ms is assumed to be 1. The quantity mvh1 for the top layer is related to the wet and dry biomasses of the leaves, which in turn are related to the green LAI [20]. So, we can substitute mvh1 for the LAI in Equations (29) and (32). In each phenological stage, we measured the mv of the sample plots, and the entire image was interpolated by the nearest neighbor interpolation algorithm using the measured mv so that each pixel had a value of mv. So, when every pixel in the image was used Equation (28) to calculate h, the mv was replaced by corresponding value. Therefore, an unknown parameter (h) has only one equation to solve for h. When solving for other parameters, such as the LAI, we assume that the other three rice parameters are specific values according to the same method so as to effectively solve for the LAI.

3.3.2. Construction of the MWCM

In order to further investigate the capability of CP SAR data in the rice parameters inversion, especially the CP decomposition parameters, such as surface, volume, and double-bounce scattering, a modified WCM (MWCM) was introduced in which the heterogeneity of the rice canopy and its phenological changes, as well as the double-bounce scattering between the rice canopy and the underlying surface were considered [27]. In the MWCM, the rice canopy was assumed to consist of many scattering cells in the horizontal direction, containing the rice and space part. A volume fraction coefficient F and a water content fraction coefficient n are defined to quantitatively describe the heterogeneity of the rice canopy in the horizontal direction. While in the vertical direction, the rice canopy was divided into leaves, stems, and ears. The multi-layers, parts, and their volumetric water content changes with the phenological stages.
In the frame work of the MWCM, ten kinds of scattering mechanisms of the rice canopy were considered during the whole growing season, in which four periods were included (Table 2). In the first period, the seedling stage, because the rice stem is short and the leaves are small, the rice part can be regarded as a whole. However, no leaves were distributed in the space part, and the water content in the space part is negligible. Therefore, there are four kinds of scattering, including the underlying surface scattering through the rice part (Sg_r), the underlying surface scattering through the space part (Sg_s), the volume scattering from the rice part (Vf_r), and the double-bounce scattering between the underlying surface and the rice part (Dg_f). From the tillering stage to booting stage, the leaves and the stems grew longer, the water content of the space part increased, and the water contents of the stem layer and the leaf layer were different. So, there are three additional scattering mechanisms, including the volume scattering from the leaf layer in the space part (Vf_s), the surface scattering from the stem layer (St), and the double-bounce scattering between the underlying surface and the stem (Dg_t). From the heading stage to the flowering stage, due to the growth of the rice ear, the rice growth reached its peak and the underlying surface of the rice was covered by the plant, and so, the double-bounce scattering between the leaves layer and the underlying surface is negligible. Moreover, the volume scattering from the ear layer in the rice part (Ve_r) and the double-bounce scattering between the underlying surface and the ear layer (De_s) is considered. From the dough stage to the mature stage, some of the rice ear moves to the space part due to increase of the biomass of the rice ear. Therefore, it is then necessary to increase the volume scattering (Ve_s) of the space part. Table 2 shows the rice scattering mechanisms for different the phenological stages. All the 10 scattering mechanisms can be expressed as a function of the rice biophysical parameters with 16 model coefficients, including F, n1, n2, Ae1(θ), Ae2(θ), Af1, Bf1, Af2, Bf2, At1, At2, Cg1(θ), Cg2(θ), αf, αt, and αe, the physical meanings of which were described in [27].

3.3.3. Coupling the CP Parameters from the m-χ and m-δ Decomposition with the MWCM.

With the CP decomposition, the surface scattering Ps, the double-bounce scattering Pd, and the volume scattering Pv are calculated. Based on the MWCM, three equations relating the three scattering components to the rice variables could be built as follows:
[ P d p v P s ] m - χ = [ D g _ e + D g _ t + D g _ f V e _ r + V e _ s + V f _ r + V f _ s S g _ r + S g _ s + S t ] ,
[ P d p v P s ] m - δ = [ D g _ e + D g _ t + D g _ f V e _ r + V e _ s + V f _ r + V f _ s S g _ r + S g _ s + S t ] ,
where the parameters of V, D, and S with subscripts can be expressed as a function of the rice parameters (i.e., De, LAI, h, mv_s, and ms) by introducing 16 model coefficients [27].
However, only three rice variables could be estimated simultaneously from Equation (34) or (35), which leads to the ill-posed problem of inversion. Therefore, as for the MWCM, the different phenological stages were assumed. Before the heading stage, no rice grains have been formed, and so, the rice biomass De is 0. The underlying surface throughout these stages is water, and so, the moisture content of the underlying surface, ms, is 1. Therefore, there are only three rice parameters, that LAI, h, and mv_s, which could be calculated easily. From the heading to the mature stage, the rice particles have been presented, and the underlying surface of the rice field is wet soil. Two hypotheses were applied as follows: (1) assume that the backscattering from the underlying surface is constant so that ms can be eliminated; and (2) the stem layer is ignored considering that the attenuation from the ear and the leaf layers are large after the heading stage, which results in the variable mv being eliminated. Three unknown rice parameters (the LAI, h, and De) can thus be calculated using Equations (34) or (35).

3.3.4. Calculating the Model Parameters Based on the GA of Rice Characteristics

A genetic algorithm [31] was used to solve Equations (28), (33), and (34) for estimating the rice biophysical parameters in five phenological stages. A total of 10 rice fields, distributed uniformly in the study area, were used as training data for the model training, and the remaining twenty rice fields were used for validation. The specific steps were as follows.
(1)
Initialization and Genetic Representation
As mentioned in Section 3.3.1 and Section 3.3.2, there are six and 16 model coefficients in the inversion models of the WCM (Equation (28)) and the MWCM (Equations (33) or (34)), respectively. In the genetic representation of the GA, each coefficient of Equations (28) and (33)/(34) was encoded by using binary coding as the gene substring.
In the WCM, the initial values and ranges of four model coefficients were initialized by solving the WCM with a set of training data. For example, the attenuation coefficients of the leaf layer and the stem layer, αl and αst, respectively, were initialized to [0, 4] and [0, 2], respectively.
As to the MWCM, the initial values and ranges of the volume fraction coefficient F and VWC fraction coefficient ni were determined by the ground data measurements and empirical assumption. In terms of the other coefficients, their initial values and ranges were initialized by solving the MWCM with a set of training data. For example, the attenuation coefficients αf and αe were initialized to [0, 1] and [0, 4], respectively.
Based on their ranges, the gene substring lengths of the coefficients were determined when the desired solution accuracy to four decimal places was considered. Then, the length of the chromosome (total of the gene substring) was determined (i.e., for the inversion model of the WCM, 66 from the elongation (21 July) to the booting stage (4 August), 51 in the heading stage (28 August), and 55 from the dough (21 September) to the mature stages (15 October); In the inversion model using the MWCM, 175 from the elongation (21 July) to the booting stage (4 August), 135 in the heading stage (28 August), and 149 from the dough (21 September) to the mature stages (15 October)). The population of the coefficients was generated using a random generator, and the population size of 100 was initialized.
(2)
Reproduction
The chromosomes generated in the initial population were then chosen for participation in the reproduction process based on their fitness values [32]. In this study, the fitness value was calculated as Equations (36)–(38) and proportionate fitness selection was used. A chromosome with a higher fitness value had a higher probability of being copied into the cross-over pool.
f = 1 1 + S S E ,
F o r   W C M ,   S S E = i = 1 m [ ( S i , e s t S i , o b s ) 2 ]   o r   i = 1 m [ ( σ i , e s t 0 σ i , o b s 0 ) 2 ] ,  
F o r   M W C M , S S E = i = 1 m [ ( P v , o b s , i P v , e s t , i ) 2 + ( P d , o b s , i P d , e s t , i ) 2 + ( P s , o b s , i P s , e s t , i ) 2 ] m χ   o r   m δ ,
where σ   i , e s t   0 is the estimated backscattering coefficient for the ith training sample, σ   i , o b s   0 is the backscattering coefficient for the observed value of the ith training sample, S i , e s t is the estimated Stokes scattering parameters for the ith training sample, and S i , o b s is the observed value of the Stokes parameters of the ith training sample. P v , e s t , i , P d , e s t , i , and P s , e s t , i represent the decomposition of the three scattering components for the ith training sample estimates of the scattering component values; P v , o b s , i , P d , o b s , i and P s , o b s , i represent the decomposition of the three scattering components to observe the scattering of the ith training sample component values. m is the number of training samples.
(3)
Cross-over
Cross-over is a recombinant operator that selects two chromosomes from the cross-over pool at random and cuts them into bits at a randomly chosen position [31]. The number of strings participating in mating depended on the cross-over probability. In this study, the cross-over probability was assumed to 0.9, and one-point cross-over was selected.
(4)
Mutation
Mutation is an important process that permits new genetic material to be introduced to a population. A mutation probability is specified that permits random mutations to be made to individual genes (e.g., changing 1 to 0, and vice versa, for binary GAs). The mutation probability of 0.01 was selected in this study.
(5)
Termination Criteria
Finally, the termination criterion of the GA process was determined. The GA process could be stopped when the fitness criterion was satisfied or the maximum number of generations was exceeded. In this study, the maximum number of generations was 2000.

4. Results and Discussions

Based on the schemes of Section 3, the rice biophysical parameters, including the rice height (h), the LAI, the VWC of the rice canopy or the stem layer (mv/mv_s), and the biomass of the rice ear (De), were estimated using the simulated CP SAR data. A total of 10 rice fields were used as training data for the GA, and the remaining twenty rice fields were used for validation. The estimated results of the rice biophysical parameters in different phenological stages were shown in Figure 7, Figure 8, Figure 9, Figure 10 and Figure 11.

4.1. Inversion of Rice Biophysical Parameters Using CP Backscattering Coefficients

Figure 7 shows the estimated results of the rice biophysical parameters using CP backscattering coefficients and the validation with the ground truth data. As the WCM treats the canopy as a water cloud, consisting of a collection of identical water particles, characterized by a uniform scattering phase function, the VWC of the rice canopy (mv) is identical to that of the stem layer (mv_s). In general, the inversion results of h and mv were superior to those of De and the LAI, when the four CP backscattering coefficients (RH/RV/RR/RL) were used. The coefficient of determinations (R2) for h and mv were both above 0.81, and the root-mean-square errors (RMSEs) were less than 10 cm and 0.39 kg/m3, respectively. This indicated that the backscattering coefficients of the four compact polarizations were more sensitive to the rice height and the VWC of the rice canopy compared with the LAI and De. Among the four compact polarizations (RH/RV/RR/RL), RV had the best performance in the inversion of h and mv. The preferable results with the RV polarization were mainly due to its sensitivity to the volume scattering of the rice canopy. For the inversion of De, RH and RR got the accuracies of R2 above 0.85 and a RMSE of less than 0.19 kg/m3, a little higher than that of RV but much better than that of RL. This suggested that RL was not as sensitive to the ear biomass as other compact polarizations. The reason is that the volume scattering was dominant in the ear layer, while RL was more related to the surface scattering. Compared with h, mv, and De, the LAI had the poorest inversion results when using the backscattering coefficients of the four compact polarizations, of which RH was the best for the retrieval of the LAI, with a R2 of 0.79 and an RMSE of 0.33.

4.2. Inversion of Rice Biophysical Parameters Using Stokes Parameters

Figure 8 shows the estimated results of the rice biophysical parameters using the Stokes parameters (S1 and S2) and the validation with the ground truth data. In general, similar to CP backscattering coefficients, the inversion results of h and mv were superior to those of De and the LAI, when the Stokes parameters (S1 and S2) were used. The R2 for h and mv were both above 0.85, and the RMSEs were less than 8 cm and 0.49 kg/m3, respectively. This indicated that the Stokes parameters (S1 and S2) were more sensitive to the rice height and the VWC of the rice canopy compared with the LAI and De. In addition, based on the equation of the water cloud model and the idea of establishing the water cloud model, h and mv contributed a lot to the model, which also led to the high precision of the inversion of h and mv. Among the Stokes parameters (S1 and S2), S1 had the best performance in the inversion of h and mv. As can be seen from the above Equations (14) and (15) and Figure 4, S1 contained multi-polarization characteristics, which contained more energy than S2. The preferable results with S1 were mainly due to its sensitivity to the total scattering of the rice canopy. For the inversion of De, S1 got the accuracies of R2 of 0.79 and RMSE of 0.21 kg/m3, much better than that of S2 with a R2 of 0.41 and a RMSE of 0.27. The main scattering energy of the rice ear came from its volume scattering, indicating that S2 was not as good as S1 in the characterization of volume scattering characteristics. Especially on 21 September and 15 October, for the inversion of De, the points of inversion results deviated from the line “y = x” significantly, which meant that the accuracy of S2 was far less than that of S1 in the dough stage and the mature stage. Compared with h, mv, and De, the LAI had the poorest inversion results when using the Stokes parameters (S1 and S2), of which S1 was the best for the retrieval of the LAI, with a R2 of 0.75 and a RMSE of 0.45.

4.3. Inversion of Rice Biophysical Parameters Using m-χ and m-δ Decomposition Parameters

Figure 9 shows the estimated results of the rice biophysical parameters using the m-χ and m-δ decomposition parameters and the validation with the ground truth data. In MWCM, the heterogeneity of the rice canopy and its phenological changes, as well as the double-bounce scattering between the rice canopy and the underlying surface were considered. Therefore, the VWC of the rice canopy (mv) is not identical to that of the stem layer (mv_s). In general, the inversion results of h and De were superior to those of mv_s and the LAI, when the m-χ and m-δ decomposition parameters were used. The R2 for h and De were both above 0.81, and the RMSEs were less than 10 cm and 0.22 kg/m2, respectively. This indicated that the m-χ and m-δ decomposition parameters were more sensitive to the rice height and De compared with the LAI and the VWC of the rice canopy. Moreover, there were close relations between h and Vf_r, Vf_s, St, Sg_r, Dg_f, and Dg_t in the MWCM. The three scattering types are all related to h and so make a significant contribution to the model, producing highly accurate results for the inversion of the rice height. Comparing the inversion results for h based on the two kinds of decomposition, the inversion accuracy based on the m-χ decomposition was higher than that based on the m-δ decomposition—the R2 was 0.07 larger, and the RMSE was 1.7 cm smaller.
According to Figure 6, the dominant scattering type for both the m-χ and m-δ decomposition is the volume scattering. Moreover, the scattering characteristics of the rice ear were mainly reflected in the volume scattering. Therefore, the De inversion using the decomposition parameters had a higher accuracy, with an R2 above 0.8 and a RMSE below 0.22 kg/m2. Comparing the inversion results for De based on the two kinds of decomposition, the inversion accuracy based on the m-χ decomposition was higher than that for the m-δ decomposition: the R2 was 0.059 bigger, and the RMSE was 0.05 kg/m2 smaller. From 28 August to 15 October, m-χ-db was bigger than m-δ-db (see Figure 6) which is related to the effect of Dg_e in the MWCM and indicates that m-χ-db is more sensitive than m-δ-db to De inversion.
For the inversion of mv_s, the m-δ decomposition parameters got the accuracies of R2 of about 0.83 and RMSE of 0.63 kg/m3, much better than that of the m-χ decomposition parameters. From Figure 9, it can be seen that both sets of inversion results produced overestimates on 21 July, with the overestimation by the m-χ decomposition being more obvious. This may be related to the three kinds of scattering components; also, at this stage, the double-bounce scattering and surface scattering values for m-χ were higher than those for the m-χ decomposition.
Compared with h, mv_s, and De, the LAI had the poorest inversion results when using the m-χ and m-δ decomposition parameters. The m-δ decomposition parameters were better than the m-χ decomposition parameters for the retrieval of the LAI, with a R2 of 0.7 and a RMSE of 0.48, but the difference in the accuracy was small: the R2 was 0.05 higher, and the RMSE was 0.14 smaller.
Figure 10 shows the inversion results for the rice biophysical parameters h, mv_s, and the LAI using the m-χ decomposition parameters. It can be seen from Figure 10 that the inversion results of h, mv_s, and the LAI have low values in the elongation stage (21 July), and the inversion results increase from the elongation stage (21 July) to the booting stage (4 August). The rice height was about 60–80 cm, the LAI is around 2–3, and the value of mv_s is about 0.5–3 in the elongation stage and the booting stage. During the heading stage (28 August) and the dough stage (21 September), the three biophysical parameters reached their peak values. The peak height of rice was then about 120 cm, the peak LAI is around 3–5, and the peak value of mv_s is about 3–4 kg/m3. These values then decreased slightly during the mature stage (15 October). This was because the rice as mature, the plant turned yellow, the stem contained less water, the size of the ear layer increased, and the rice got shorter. The inversion results for the other CP parameters show similar patterns to those seen in Figure 10.
Figure 11 shows the inversion results for De using the m-χ decomposition parameters. Because the ears began to form during the heading stage (28 August), De then is low. De reached its peak during the dough stage (21 September). The rice ears on September 21 were nearly mature, and so, the value De remained basically unchanged on 15 October. The results for inverting De using the other CP parameters show similar patterns to those seen in Figure 11.

5. Conclusions

The objective of this study was to investigate the capability of CP SAR data in the inversion of rice biophysical parameters during the growing season. Five temporal RADARSAT-2 fully polarimetric SAR datasets were used to simulate CP SAR data, and the CP parameters, such as backscattering coefficients, Stokes vectors and decomposition parameters were extracted in each acquisition. Then, the inversion models of rice biophysical parameters with CP SAR data were developed by introducing a classic semi-empirical model, water cloud model, and a modified WCM in which the heterogeneity of the rice canopy and its phenological changes as well as the double-bounce scattering between the rice canopy and the underlying surface were considered. Finally, rice biophysical parameters, including rice height (h), LAI, the VWC of the rice canopy or the stem layer (mv/mv_s), and biomass of the ear (De), were obtained, and the validation was conducted using the ground truth data. The detailed conclusions are as follows:
(1)
Using the four CP backscattering coefficients (RH/RV/RR/RL), the inversion results of h and mv_s were superior to those of LAI and De, with the R2 above 0.81 and the RMSE less than 10 cm and 0.39 kg/m3, respectively. RV had the best performance in the inversion of h and mv_s. For De, RH and RR performed better than RV and RL, giving an R2 above 0.85 and an RMSE of less than 0.18 kg/m2. Again, the inversion results for the LAI were the poorest, with RV and RL giving better results than RH and RR. For RL, the R2 was 0.73, and the RMSE was less than 0.55.
(2)
Using the Stokes parameters, the inversion results for the four rice biophysical parameters were good. For h, the R2 had a value of about 0.92, and the RMSE was less than 5.8 cm. For mv_s, the R2 was about 0.93, and the RMSE was less than 0.36 kg/m3. For De, the R2 was about 0.80, and the RMSE was less than 0.21 kg/m2. The inversion results for the LAI were the poorest, with an R2 of 0.75 and an RMSE of up to 0.45.
(3)
Using the m-χ and m-δ decomposition parameters, the inversion accuracy for h and De was higher than for mv_s and the LAI. For h, the coefficient of determination, R2, was about 0.88, and the RMSE was less than 8.2 cm. For De, the R2 was 0.89, and the RMSE was less than 0.17 kg/m2. For mv_s, the coefficient of determination was 0.83, and the RMSE was less than 0.6 kg/m3. As before, the inversion results for the LAI had the lowest accuracy, with an R2 of 0.73 and an RMSE up to 0.55. Overall, the h and De inversion results obtained using the m-χ decomposition parameters were better than those obtained using the m-δ decomposition. For LAI and mv_s, the accuracy obtained using the m-δ decomposition was slightly higher than that found using the m-χ decomposition. In the early stages of rice growth, mv_s was slightly overestimated; in the case of the m-χ decomposition parameters, this overestimation was even greater.
(4)
In general, S1, RV, RL, and the m-χ decomposition parameters proved to be more suitable for the inversion of h, and high inversion accuracy was obtained. S1, RH, RV, RR, and RL are more suitable for inverting mv. RH, RR, and the decomposition parameters are more suitable for inverting De.
In this paper, the capability of CP SAR data in the inversion of rice biophysical parameters has been proven, especially for the inversion of h, De, and mv_s/mv, with an R2 of higher than 0.89 and an RMSE of less than 10 cm, 0.68 kg/m3, and 0.27 kg/m2. Further validation and improvement will be conducted with the subsequent dataset of this site and the dataset of other regions.

Author Contributions

X.G. and K.L. proposed and performed the experiments and wrote the paper; Y.S. and Z.W. checked the paper and provided suggestions. H.L., Z.Y., L.L., and S.W. contributed to the data analysis, making charts and the paper modification, respectively. The order of the authors reflects their level of contribution.

Acknowledgments

This research was supported by funding from the National Natural Science Foundation of China (61471358, 41431174, Y6SJ1400CX, 41601375) and the GEO international cooperation project (Joint Experiment of Crop Assessment and Monitoring, JECAM).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Maclean, J.L.; Dawe, D.C.; Hardy, B.; Patel, G.P. Rice Almanac: Source Book for the Most Important Economic Activity on Earth; International Rice Research Institute: Los Banos, Philippines, 2002; p. 253. ISBN 971-22-0172-4. [Google Scholar]
  2. Toan, T.L.; Ribbes, F.; Wang, L.F.; Floury, N.; Ding, K.H.; Kong, J.A.; Fujita, M.; Kurosu, T. Rice crop mapping and monitoring using ERS-1 data based on experiment and modeling results. IEEE Trans. Geosci. Remote Sens. 1997, 35, 41–56. [Google Scholar] [CrossRef]
  3. Bouvet, A.; Toan, T.L.; Lam-Dao, N. Monitoring of the Rice Cropping System in the Mekong Delta Using ENVISAT/ASAR Dual Polarization Data. IEEE Trans. Geosci. Remote Sens. 2009, 47, 517–526. [Google Scholar] [CrossRef] [Green Version]
  4. Shao, Y.; Fan, X.; Liu, H.; Xiao, J.; Ross, S.; Brisco, B.; Brown, R.; Staples, G. Rice monitoring and production estimation using multitemporal RADARSAT. Remote Sens. Environ. 2001, 76, 310–325. [Google Scholar] [CrossRef]
  5. Erten, E.; Lopez-Sanchez, J.M.; Yuzugullu, O.; Hajnsek, I. Retrieval of agricultural crop height from space: A comparison of SARtechniques. Remote Sens. Environ. 2016, 187, 130–144. [Google Scholar] [CrossRef]
  6. Tan, L.; Chen, Y.; Jia, M.; Tong, L.; Li, X.; He, L. Rice biomass retrieval from advanced synthetic aperture radar image based on radar backscattering measurement. J. Appl. Remote Sens. 2015, 9, 097091. [Google Scholar] [CrossRef] [Green Version]
  7. Raney, R.K.; Spudis, P.D.; Bussey, B.; Crusan, J.; Jensen, J.R.; Marinelli, W.; Mckerracher, P.; Neish, C.; Palsetia, M.; Schulze, R. The Lunar Mini-RF Radars: Hybrid Polarimetric Architecture and Initial Results. Proc. IEEE 2011, 99, 808–823. [Google Scholar] [CrossRef]
  8. Kunhikrishnan, P.; Sowmianarayanan, L.; Nampoothiri, M.V. PSLV-C19/RISAT-1 mission: The launcher aspects. Curr. Sci. 2013, 104, 472–476. [Google Scholar]
  9. Arikawa, Y.; Suzuki, S. ALOS-2 launch and initial checkout result. In Proceedings of the Sensors, Systems, and Next-Generation Satellites XVIII, Amsterdam, The Netherlands, 11 November 2014. [Google Scholar]
  10. Sanden, J.J.V.D.; Geldsetzer, T. Compact Polarimetry in Support of Lake Ice Breakup Monitoring: Anticipating the RADARSAT Constellation Mission. Can. J. Remote Sens. 2015, 41, 440–457. [Google Scholar] [CrossRef]
  11. Brisco, B.; Li, K.; Tedford, B.; Charbonneau, F.; Yun, S.; Murnaghan, K. Compact polarimetry assessment for rice and wetland mapping. Int. J. Remote Sens. 2013, 34, 1949–1964. [Google Scholar] [CrossRef]
  12. Uppala, D.; Kothapalli, R.V.; Poloju, S.; Mullapudi, S.S.V.R.; Dadhwal, V.K. Rice Crop Discrimination Using Single Date RISAT1 Hybrid (RH, RV) Polarimetric Data. Photogramm. Eng. Remote Sens. 2015, 81, 557–563. [Google Scholar] [CrossRef]
  13. Li, K.; Brisco, B.; Shao, Y.; Touzi, R. Polarimetric decomposition with RADARSAT-2 for rice mapping and monitoring. Can. J. Remote Sens. 2012, 38, 169–179. [Google Scholar] [CrossRef]
  14. Watanabe, T.; Izumi, Y.; Tetuko, J.; Sumantyo, S. Comparison of three compact-SAR mode through rice phenology monitoring. In Proceedings of the International Symposium on Antennas and Propagation (ISAP), Phuket, Thailand, 30 October–2 November 2017. [Google Scholar]
  15. Yang, Z.; Li, K.; Liu, L.; Shao, Y.; Brisco, B.; Li, W. Rice growth monitoring using simulated compact polarimetric C band SAR. Radio Sci. 2014, 49, 1300–1315. [Google Scholar] [CrossRef] [Green Version]
  16. Lopez-Sanchez, J.M.; Vicente-Guijalba, F.; Ballester-Berman, J.D.; Cloude, S.R. Polarimetric Response of Rice Fields at C-Band: Analysis and Phenology Retrieval. IEEE Trans. Geosci. Remote Sens. 2014, 52, 2977–2993. [Google Scholar] [CrossRef] [Green Version]
  17. Kumar, V.; Mcnairn, H.; Bhattacharya, A.; Rao, Y.S. Temporal Response of Scattering From Crops for Transmitted Ellipticity Variation in Simulated Compact-Pol SAR Data. IEEE. Sel. Top. Appl. Earth Obs. Remote Sens. 2017, 10, 5163–5174. [Google Scholar] [CrossRef]
  18. Zhang, W.; Li, Z.; Chen, E.; Zhang, Y.; Yang, H.; Zhao, L.; Ji, Y. Compact Polarimetric Response of Rape (Brassica napus L.) at C-Band: Analysis and Growth Parameters Inversion. Remote Sens. 2017, 9, 591. [Google Scholar]
  19. Yang, H.; Xie, L.; Chen, E.X.; Zhang, H.; Yang, G.J.; Li, Z.H.; Gu, X.H. Biomass estimation of oilseed rape using simulated compact polarimtric SAR imagery. In Proceedings of the 2016 IEEE International Geoscience and Remote Sensing Symposium (IGARSS), Beijing, China, 10–15 July 2016. [Google Scholar]
  20. Ulaby, F.; Allen, C.; Eger, G.; Kanemasu, E. Relating the microwave backscattering coefficient to leaf area index. Remote Sens. Environ. 1984, 14, 113–133. [Google Scholar] [CrossRef]
  21. Thompson, A.A.; Mcleod, I.H. The RADARSAT-2 SAR processor. Can. J. Remote Sens. 2004, 30, 336–344. [Google Scholar] [CrossRef]
  22. Shi, Z.Y.; Gao, X.F.; Xie, Y. The application of SUNSCAN canopy analysis system in the measurement of field ecosystem. Agric. Res. Arid Areas 2005, 23, 78–82. [Google Scholar]
  23. Attema, E.; Ulaby, F.T. Vegetation modeled as a water cloud. Radio Sci. 1978, 13, 357–364. [Google Scholar] [CrossRef]
  24. Raney, R.K. Hybrid-Polarity SAR Architecture. IEEE Trans. Geosci. Remote Sens. 2007, 45, 3397–3404. [Google Scholar] [CrossRef]
  25. Yin, J.; Yang, J. Ship detection by using the M-Chi and M-Delta decompositions. In Proceedings of the 2016 IEEE International Geoscience and Remote Sensing Symposium (IGARSS), Quebec, QC, Canada, 13–18 July 2014. [Google Scholar]
  26. Raney, R.K.; Cahill, J.S.T.S.; Patterson, G.W.; Bussey, D.B.J. The m-chi decomposition of hybrid dual-polarimetric radar data. In Proceedings of the 2012 IEEE International Geoscience and Remote Sensing Symposium, Munich, Germany, 22–27 July 2012. [Google Scholar]
  27. Yang, Z.; Li, K.; Shao, Y.; Brisco, B.; Liu, L. Estimation of Paddy Rice Variables with a Modified Water Cloud Model and Improved Polarimetric Decomposition Using Multi-Temporal RADARSAT-2 Images. Remote Sens. 2016, 8, 878. [Google Scholar] [CrossRef]
  28. Beriaux, E.; Lucau-Danila, C.; Auquiere, E.; Defourny, P. Multiyear independent validation of the water cloud model for retrieving maize leaf area index from SAR time series. Int. J. Remote Sens. 2013, 34, 4156–4181. [Google Scholar] [CrossRef]
  29. Lievens, H.; Verhoest, N.E.C. On the Retrieval of Soil Moisture in Wheat Fields from L-Band SAR Based on Water Cloud Modeling, the IEM, and Effective Roughness Parameters. IEEE Geosci. Remote Sens. Lett. 2011, 8, 740–744. [Google Scholar] [CrossRef]
  30. Bériaux, E.; Waldner, F.; Collienne, F.; Bogaert, P.; Defourny, P. Maize Leaf Area Index Retrieval from Synthetic Quad Pol SAR Time Series Using the Water Cloud Model. Remote Sens. 2015, 7, 16204–16225. [Google Scholar] [CrossRef] [Green Version]
  31. Golberg, D.E. Genetic Algorithms in Search, Optimization, and Machine Learning; Addison Wesley: Boston, MA, USA, 1989. [Google Scholar]
  32. Kumar, K.; Hari Prasad, K.S.; Arora, M.K. Estimation of water cloud model vegetation parameters using a genetic algorithm. Hydrol. Sci. J. 2012, 57, 776–789. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Photographs showing the eight phenological stages of indica rice in the field.
Figure 1. Photographs showing the eight phenological stages of indica rice in the field.
Sensors 18 02271 g001
Figure 2. Methodology flow chart for the inversion of the rice biophysical parameters with the water cloud model (WCM) and modified WCM (MWCM) and compact-polarimetric (CP) synthetic aperture radar (SAR) parameters.
Figure 2. Methodology flow chart for the inversion of the rice biophysical parameters with the water cloud model (WCM) and modified WCM (MWCM) and compact-polarimetric (CP) synthetic aperture radar (SAR) parameters.
Sensors 18 02271 g002
Figure 3. The simulated CP SAR data of 21 July 2012—color composite images of CP backscattering coefficients (R = RR (right circular transmit and right circular receive), G = RV (right circular transmit and horizontal linear receive), and B = RH (right circular transmit and horizontal linear receive)). The red circles in the figure represents 30 typical rice fields as the experimental area.
Figure 3. The simulated CP SAR data of 21 July 2012—color composite images of CP backscattering coefficients (R = RR (right circular transmit and right circular receive), G = RV (right circular transmit and horizontal linear receive), and B = RH (right circular transmit and horizontal linear receive)). The red circles in the figure represents 30 typical rice fields as the experimental area.
Sensors 18 02271 g003
Figure 4. The temporal behaviors of the four Stokes parameters from the elongation stage to the mature stage.
Figure 4. The temporal behaviors of the four Stokes parameters from the elongation stage to the mature stage.
Sensors 18 02271 g004
Figure 5. The temporal behaviors of the CP backscattering coefficients in the RH, RV, RR and RL polarizations from elongation stage to mature stage.
Figure 5. The temporal behaviors of the CP backscattering coefficients in the RH, RV, RR and RL polarizations from elongation stage to mature stage.
Sensors 18 02271 g005
Figure 6. The temporal behaviors of the three scattering components of the m-χ and m-δ decomposition from the elongation stage to the mature stage, db = double bounce, vol = volume, and s = surface.
Figure 6. The temporal behaviors of the three scattering components of the m-χ and m-δ decomposition from the elongation stage to the mature stage, db = double bounce, vol = volume, and s = surface.
Sensors 18 02271 g006
Figure 7. Results for the rice biophysical parameters using the CP backscattering coefficients for the different phenological stages of rice.
Figure 7. Results for the rice biophysical parameters using the CP backscattering coefficients for the different phenological stages of rice.
Sensors 18 02271 g007aSensors 18 02271 g007b
Figure 8. Results for the rice biophysical parameters using the Stokes parameters for the different phenological stages of rice.
Figure 8. Results for the rice biophysical parameters using the Stokes parameters for the different phenological stages of rice.
Sensors 18 02271 g008aSensors 18 02271 g008b
Figure 9. Results for the rice biophysical parameters using m-χ and m-δ decomposition for the different phenological stages of rice.
Figure 9. Results for the rice biophysical parameters using m-χ and m-δ decomposition for the different phenological stages of rice.
Sensors 18 02271 g009
Figure 10. The inversion results for the rice biophysical parameters at five different phenological stages using the m-χ decomposition parameters. (The five images in the first row are the inversion results for h, the five images in the second row are the results for the leaf area index (LAI), and the five images in the third row are the results for mv_s).
Figure 10. The inversion results for the rice biophysical parameters at five different phenological stages using the m-χ decomposition parameters. (The five images in the first row are the inversion results for h, the five images in the second row are the results for the leaf area index (LAI), and the five images in the third row are the results for mv_s).
Sensors 18 02271 g010
Figure 11. The inversion results for the rice biophysical parameter De at three phenological stages using the m-χ decomposition parameters.
Figure 11. The inversion results for the rice biophysical parameter De at three phenological stages using the m-χ decomposition parameters.
Sensors 18 02271 g011
Table 1. Details of the full polarimetric RADARSAT-2 SAR data.
Table 1. Details of the full polarimetric RADARSAT-2 SAR data.
DateModeProductIncidence Angle (°)Pixel Spacing (A × R, m)Phenology
Stage
21 July 2012FQ20W 1SLC 238–415.2 × 7.6Elongation
4 August 2012FQ9WSLC27–305.2 × 7.6Booting
28 August 2012FQ9WSLC27–305.2 × 7.6Heading
21 September 2012FQ9WSLC27–305.2 × 7.6Dough
15 October 2012FQ9WSLC27–305.2 × 7.6Mature
1 FQW = fine quad-polarimetry wide, 20 or 9 is the number of the beam position, which is related to the incidence angles; 2 SLC = single look complex [21].
Table 2. Rice scattering mechanisms for different phenological stages.
Table 2. Rice scattering mechanisms for different phenological stages.
Phenological StagesUnderlying SurfaceRice Part ScatteringSpace Part Scattering
SeedlingWater surfaceSg_rVf_r Sg_sDg_f
Tillering to BootingWater surfaceSg_rVf_rStDg_t Sg_sDg_fVf_s
Heading to FloweringMoist soilSg_rVf_rStDg_tVe_rSg_s Vf_sDg_e
Dough to MatureMoist soilSg_rVf_rStDg_tVe_rSg_s Vf_sDg_eVe_s

Share and Cite

MDPI and ACS Style

Guo, X.; Li, K.; Shao, Y.; Wang, Z.; Li, H.; Yang, Z.; Liu, L.; Wang, S. Inversion of Rice Biophysical Parameters Using Simulated Compact Polarimetric SAR C-Band Data. Sensors 2018, 18, 2271. https://doi.org/10.3390/s18072271

AMA Style

Guo X, Li K, Shao Y, Wang Z, Li H, Yang Z, Liu L, Wang S. Inversion of Rice Biophysical Parameters Using Simulated Compact Polarimetric SAR C-Band Data. Sensors. 2018; 18(7):2271. https://doi.org/10.3390/s18072271

Chicago/Turabian Style

Guo, Xianyu, Kun Li, Yun Shao, Zhiyong Wang, Hongyu Li, Zhi Yang, Long Liu, and Shuli Wang. 2018. "Inversion of Rice Biophysical Parameters Using Simulated Compact Polarimetric SAR C-Band Data" Sensors 18, no. 7: 2271. https://doi.org/10.3390/s18072271

APA Style

Guo, X., Li, K., Shao, Y., Wang, Z., Li, H., Yang, Z., Liu, L., & Wang, S. (2018). Inversion of Rice Biophysical Parameters Using Simulated Compact Polarimetric SAR C-Band Data. Sensors, 18(7), 2271. https://doi.org/10.3390/s18072271

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop