Analysis of Combining SAR and Optical Optimal Parameters to Classify Typhoon-Invasion Lodged Rice: A Case Study Using the Random Forest Method

Lodging, a commonly occurring rice crop disaster, seriously reduces rice quality and production. Monitoring rice lodging after a typhoon event is essential for evaluating yield loss and formulating suitable remedial policies. The availability of Sentinel-1 and Sentinel-2 open-access remote sensing data provides large-scale information with a short revisit time to be freely accessed. Data from these sources have been previously shown to identify lodged crops. In this study, therefore, Sentinel-1 and Sentinel-2 data after a typhoon event were combined to enable monitoring of lodging rice to be quickly undertaken. In this context, the sensitivity of synthetic aperture radar (SAR) features (SF) and spectral indices (SI) extracted from Sentinel-1 and Sentinel-2 to lodged rice were analyzed, and a model was constructed for selecting optimal sensitive parameters for lodging rice (OSPL). OSPL has high sensitivity to lodged rice and strong ability to distinguish lodged rice from healthy rice. After screening, Band 11 (SWIR-1) and Band 12 (SWIR-2) were identified as optimal spectral indices (OSI), and VV, VV + VH and Shannon Entropy were optimal SAR features (OSF). Three classification results of lodging rice were acquired using the Random Forest classification (RFC) method based on OSI, OSF and integrated OSI–OSF stack images, respectively. Results indicate that an overall level of accuracy of 91.29% was achieved with the combination of SAR and optical optimal parameters. The result was 2.91% and 6.05% better than solely using optical or SAR processes, respectively.


Introduction
Rice is a major food source for nearly half of the world's population. It is widely planted in China [1], making China the largest global rice producer [2]. Rice lodging is a common disaster which is mainly caused by typhoons or hurricanes. During these events, strong winds bend the rice stems and heavy rain causes flooding of the crop fields, resulting in rice grains and roots to rot after prolonged soaking, having a serious effect on paddy rice production [3]. The magnitude of rice production and grain quality will consequently decrease as photosynthesis processes are disrupted due to self-shading, and translocation of nutrients and carbon from roots to grain is interpreted after lodging [4,5]. It is, therefore, important to quickly and accurately determine the range and location of lodged rice fields for disaster assessment, emergency response and resource allocation.
The traditional method of lodging rice monitoring uses on-site visual observation and assessment, methods which are inefficient and expensive. The rapid development of remote sensing technology has enabled the development of a cost-effective and multi-scale method to monitor lodging crops [6]. As rice lodging is mainly caused by typhoons, which are often accompanied by storms and heavy rain, monitoring studies of lodging rice using remote sensing technology has mainly focused on optical and infrared sensors on unmanned aerial vehicles (UAVs) that operate below cloud cover. Liu et al. [7] used principal component analysis (PCA) and probabilistic neural network (PNN) to demonstrate that healthy and lodged rice can be distinguished using hyperspectral data. Yang et al. [8] combined spectral indices, texture information and the digital surface model (DSM) extracted from UAV images to acquire high classification accuracy of lodging rice using a decision tree method. Liu et al. [9] constructed a model to recognize lodged rice using support vector machine (SVM) and particle swarm optimization (PSO) by analyzing temperature, color and texture information of lodged rice collected from UAV visible and thermal infrared images. Although UAV plays an important role in monitoring rice lodging, the imaging width of UAV is limited to the local or experimental scale; therefore, it cannot meet the requirement of rapid regional-scale monitoring [10].
As synthetic aperture radar (SAR) microwaves can penetrate clouds, haze and light rain, SAR backscattering signals obtained from radar satellites are not generally influenced by weather conditions. SAR data are rapidly collected and can cover areas. SAR signals are also sensitive to the morphological structure of rice canopies [11,12]. Radarsat-2 and Sentinel-1 data have been explored to monitor lodging crops. Previous investigations have shown that multiple polarimetric features have the ability to distinguish between lodging wheat and sugarcane from healthy crops based on Randarsat-2 fully polarimetric data [6,13,14]. By analyzing the sensitivity of polarimetric features extracted from Sentinel-1 to lodged corn, Han and Shu et al. [15,16] constructed a model to classify lodged corn with different lodging angles. Chauhan et al. [17,18] explored the ability of Radarsat-2 and Sentinel-1 data with multiple incidence angles to classify healthy wheat and lodged wheat with different crop angles of inclination (CAI). They found that images with low incidence had good performance for high CAI wheat identification, and vice versa.
Although satellite-based optical images are often contaminated by cloud shadows [19], it has been demonstrated that multi-spectral indices of optical data are sensitive to lodging crops [20]; optical and SAR data provide more information for distinguishing and monitoring lodging crops. Although the information synergy of SAR features and optical indices has been widely studied in crop variety identification and crop mapping [1,19,21], few applications for lodging crop mapping have been undertaken. Having the characteristics of a fast-revisiting period, a wide-coverage area and free-access, Sentinel-1/2 are suited to efficiently monitor lodging rice after a typhoon or hurricane [22]. Wu et al. [23] proposed a workflow for mapping lodging paddy rice area based on Sentinel-1 and Sentinel-2 data in the Google Earth Engine (GEE) platform; however, limitation parameters (VV, VH, NDVI and EVI) were used and analyzed, and the single parameter image was used in the lodging rice classification step.
In this study, 10 Radar features of Sentinel-1 and 10 spectral indices of Sentinel-2 were evaluated for lodging rice sensitivity, and a selection model was developed for acquiring the optimal sensitivity parameters for lodging rice (OSPL). Optimal SAR features (OSF) and optimal spectral indices (OSI) were then selected based on the OSPL model. Three classification results were acquired using the Random Forest (RF) classification method with the input of OSF, OSI and integrated OSF-OSI. The accuracy of the three classification results was evaluated using the validation data and confusion matrix.

Study Area
The study area is situated within Qindeli farm (Figure 1), Tongjiang City, Heilongjiang Province, China (47 • 51 31" N, 133 • 7 32" E). The study area is located in a humid monsoon climate zone, in the third accumulated temperature zone in Heilongjiang Province, being one of the important commodity grain production bases in China. The total cultivated area of Qindeli Farm is about 52,800 hectares, of which rice planting accounts for more than 80% of the total cultivated area. Single season rice is cultivated in the test area, having a growing period from early May to late September.
On 7 September 2019, typhoon "Lingling" passed over this area, resulting in heavy rainfall that lasted for two days. During this event, a large area of rice plants (the majority of which were in the mature stage) were lodged, as can be seen in Figure 2a

Experimental Data
In order to explore the sensitivity of different SAR features and spectral indices to lodged rice, three Sentinel-1B and two Sentinel-2(A/B) images over the study area were downloaded from the Copernicus Open Access Hub of European Space Agency (ESA) before and after Typhoon Lingling; the main parameters of Sentinel-1/2 data used in this study were detailed in (Table 1). A single look complex (SLC) format of Sentinel-1 data was used to extract backscatter coefficients of different polarizations and H/A/Alpha polarization decomposition parameters.
Two scenes of Sentinel-2 Level-2A optical images in the UTM/WGS84 projection were cloudless over the study area, and their reflectance was extracted, except for bands 1, 9 and 10, which are irrelevant for crop studies. At the same time, wind speed and rainfall hourly data were collected from the China Meteorological Data Service Center (CMDC) to aid interpretation of results.
Due to the large scale of lodged rice, ground data provided by local farm worker could not provide field scale verification. Therefore, a Pleiades satellite 0.5 panchromatic high spatial resolution image captured on 14 September (within one week of the typhoon event) was introduced to downscale ground data from a local scale to a field scale. In addition, the SRTM (Shuttle Radar Topography Mission) 30 m resolution digital elevation model (DEM) was used for geocoding of Sentinel-1 SAR imagery.

Pre-Processing of Remote Sensing Images
Backscatter coefficients from SLC format data of Sentinel-1 were extracted using SNAP 6.0. After applying the orbit file correction, the thermal noise removal and radiometric calibration operator were conducted so that pixel values were directly related to the radar backscatter of the scene. The refined Lee speckle filter was then applied with a 3 × 3 window to eliminate the inherent speckles of the SAR image, and the terrain correction operator was performed so that geometric representation of the image was as close as possible to the real world. Finally, the value of resulting images was converted to dB from linear (σ 0 VV and σ 0 VH ). σ 0 VV+VH ,σ 0 VV−VH and σ 0 VH/VV were also calculated using the band math operator.
In relation to polarimetric decomposition parameter extraction, IW1 subswath was initially extracted using PolSARpro 6.0. Multi-looking (Azimuth: 1, Range: 4) was then performed followed by C2 matrix extraction. SNAP software assisted in completing part of the data preprocessing, including TOPSAR spilt, Apply Orbit File, Complex Calibration and Deburst. After geocoding and refined Lee filter, a total of five polarimetric features were acquired using H/A/Alpha polarimetric decomposition, including Alpha, Entropy Anisotropy, Span and Shannon Entropy. Different from fully polarimetric data, H/A/Alpha decomposition parameters of dual-pol data were generated from C2 matrix, as: The process of Sentinel-2 L2A products mainly includes bottom-of-atmosphere (BOA) corrected reflectance, and images were subset to reduce processing time. Finally, the Sen2Res tool was used to improve the spatial resolution of Sentinel-2 products to 10 m resolution whilst simultaneously preserving their reflectance.

Statistical Analysis
Ten lodged (L1-L10) and ten healthy (H1-H10) rice sample plots were selected in the study area based on downscaling field data ( Figure 1). Inhomogeneity in each plot was eliminated by extracting mean values of SAR features (SF) and spectral indices (SI) (see Section 3.1 for data pre-processing). Sensitivity differences of SF and SI to lodged rice and the reasons for these differences were then analyzed.

OSPL Screening Model
After undertaking sensitivity analysis of SF and SI to lodged rice, an OSPL screening model was developed. Parameters selected by the OSPL screening model were highly sensitive to lodged rice, having a strong ability to distinguish lodged rice from healthy rice. Selection of the OSPL screening model included three important steps: (i) A sensitivity criterion to rice lodging, γ, was proposed to select parameters sensitive to lodged rice, as: where LB av (LB av = (L1 + L2 + . . . + L10)/10) is the average value of ten lodged rice plots before the typhoon event; LA av (LA av = (L1 + L2 + . . . + L10)/10) is the average value of ten lodged rice plots after the typhoon event; HB av (HB av = (H 1 + H 2 + . . . + H 10 )/10) is the average value of ten healthy rice plots before the typhoon event and HA av (HA av = (H 1 + H 2 + . . . + H 10 )/10) is the average value of ten healthy rice plots after the typhoon event.
When γ is greater than 0, the variation of the corresponding parameter value of the lodged rice plot before and after typhoon is much greater than that of the healthy rice plot. Therefore, parameters meeting the criterion that γ is greater than 0 are sensitive to lodged rice.
It is worth noting that each element in the formula of sensitivity criterion γ is the average value of ten rice plots. As accuracy and reliability of γ value determination will be reduced by the presence of extreme noise values, it was therefore important to evaluate the value of parameters screened in the first step for each plot in the following steps.
(ii) Two parameters (β and ε) were constructed to analyze the variation tendency of SF and SI of all the lodged rice sample plots after lodging, as: where, ε i is the difference between parameter values LB i (lodged sample plot before the typhoon) and LA i (lodged sample plot after the typhoon); β is the number of the sign of ε i with a higher frequency (ε has six positive signs and four negative signs, then β = 6). When β is ≥9, the corresponding parameter will be selected. After this step, sensitivity of the selected parameters to lodged rice has good consistency.
(iii) Parameters that fulfill the following criteria were used to distinguish lodged rice from healthy rice. (

Classification of Lodging and Healthy Rice
Before lodged rice and healthy rice fields were classified, interference due to buildings, water bodies, woods or other crops to lodged rice mapping were eliminated using the maximum likelihood classification (MLC) method [24] to extract the rice region as a rice mask shapefile using Sentinel-2 images. The extracted rice mask shapefile was then applied to all OSPL images to remove non-rice regions, and a Random Forest (RF) classifier was introduced to lodged rice mapping, having three inputs of OSF, OSI and integrated OSF-OSI stack images.
RF, a supervised machine learning method, is built on Classification and Regression Trees (CARTs) using the GINI index as the calculation principle. This method was developed by Breiam [25] and it has been widely applied to classify crops at different scales [21,[26][27][28]. The classification result is generated from ensemble voting results of multi trees for the most popular class. Each tree is trained using sub-sample dataset grabbed from training dataset by bootstrap sampling, and the proportional randomly selected features from the feature dataset are used at each tree node. By excluding the 10 + 10 sample fields previously mentioned, a total of 4303 additional sample point pixels were randomly selected on the basis of downscaling field data and high-resolution Google Earth imagery, of which 2876 pixels were assigned for training, and 1427 pixels were set as a validation dataset. After being fine-tuned numerous times, the RF optimizer had 100 trees, two maximum features and a 0.2 min impurity. Finally, the classification accuracies were evaluated using the validation dataset and confusion matric. The flowchart of the detailed processes for lodged rice mapping is shown in Figure 3.

SAR Parameter Analysis
The backscatter behavior of rice and other typical features (buildings, woods and barren land) were analyzed before, during and after the typhoon. Backscatter coefficients of VV and VH polarizations (Figure 4a,b, respectively) before (26 August), during (7 September) and after the typhoon (19 September) indicate that rice and other typical features had similar backscatter behavior. This behavior initially increased (26 August to 7 September) before decreasing (19 September). The significant increase in backscatter coefficients on 7 September was mainly caused by heavy rain (63 mm/d, Figure 4c), similar to findings from previous studies [29,30].
The significant increase in backscatter results recorded on 7 September associated to heavy rain, resulted in rice lodging to be overwhelmed. Therefore, images on 26 August and 19 September were only used to explore the capability of SAR in rice lodging mapping. Backscatter results (Figure 5a,b) indicate that backscatter coefficients of lodged rice sample plots (LS) significantly increased from 26 August to 19 September at VV and VH polarizations, and slight difference was recorded from 26 August to 19 September in healthy rice sample plots (HS).
Backscatter coefficients of σ 0 VV and σ 0 VH in LS increased by 1.15 dB and 1.61 dB after the typhoon, respectively ( Table 2). The increase of σ 0 VV backscatter coefficients was probably due to a decrease of signal attenuation, associated to the destruction of the vertical structure of the rice canopy during rice lodging [31]; the increase of σ 0 VH may be the enhancement of multiple scatterings [6]. In addition, the increment of σ 0 VH was greater than σ 0 VV for LS (Table 2). Two probable factors can account for this phenomenon: (i) When VV polarized wave interacts with rice plants, double scatter is the main component of the total backscatter signals [32], however, as the angle between rice stems and the ground surface decreases after lodging, backscatter signals weaken. (ii) Volume scatter, the main component of total backscatter signals for σ 0 VH , increases due to enhancement of multi scattering after lodging [6].   The polarimetric decomposition parameters, such as entropy (H), anisotropy (A) and average scattering angle (α) of healthy and lodged rice, were also analyzed. H, A and α recorded a low sensitivity to rice lodging, with L σ change values of 0.09, 0.02 and 0.02 after rice lodging, respectively ( Table 2).
H/A/Alpha, reflecting target scatter mechanisms, have been shown to be strongly sensitive to lodged crops when generated from fully polarimetric data [6,17]. However, the same parameters from dual-pol data were reported by Chauhan et al. [17] to have a lack of sensitivity to lodged crops. This is probably the reason why H/A/Alpha were not sensitive to rice lodging in this study. The polarimetric decomposition parameters, such as entropy (H), anisotropy (A) and average scattering angle (α) of healthy and lodged rice, were also analyzed. H, A and α recorded a low sensitivity to rice lodging, with σ L change values of 0.09, 0.02 and 0.02 after rice lodging, respectively ( Table 2).
H/A/Alpha, reflecting target scatter mechanisms, have been shown to be strongly sensitive to lodged crops when generated from fully polarimetric data [6,17]. However, the same parameters from dual-pol data were reported by Chauhan et al. [17] to have a lack of sensitivity to lodged crops. This is probably the reason why H/A/Alpha were not sensitive to rice lodging in this study.

Reflectance Analysis
LS and HS spectral behavior were analyzed before (2 September) and after (17 September) the typhoon event (Figure 5a,b). Although rice lodging had not yet occurred, LS and HS reflectance differences in the majority of spectral bands were recorded before the typhoon occurred (Figure 5a), with variety differences between non-resistant and resistant lodging probably accounting for this situation. Due to canopy architecture and biophysical factors (such as chlorophyll and LAI (Leaf area index)) of rice changing after rice lodging, LS reflectance significantly increased in all spectral bands after the typhoon event ( Figure 5b) [33,34]. In addition, different spectral band reflectance values for each sample plot before and after the typhoon were analyzed (Figure 6a,b). Results indicate that HS reflectance values in almost all spectral bands also increased after the typhoon (Figure 6a), probably indicating a stronger illumination in the study area on 17 September compared to 2 September. Finally, by combining the above conditions of canopy architecture, biophysical changes and illumination difference, the biggest reflectance difference was recorded on 17 September (Figure 6b), recording the greatest contrast in LS reflectance before and after the typhoon.
Results from this analysis indicate that the increase in the magnitude of reflectance values within the VIS-NIR-SWIR region for LS after the typhoon was evident (Figure 6b); reasons for this result have been previously investigated by Chauhan et al. [20]. The reduction in rice chlorophyll content due to disrupted photosynthesis is a major factor increasing reflectance values in the VIS region [5]. Reflectance increases in NIR and SWIR are associated with canopy structure and plant water content changes.

Selecting of OSPL
OSPL screening model incorporating γ and β are shown in Table 3 (for SF) and Table 4 (for SI). According to the criterion of γ > 0, VV, VH, VV + VH, VV − VH, Shannon Entropy, and Span were selected as SAR features with parameters sensitive to lodged rice. As the magnitude of β was highly correlated with positive and negative γ values (Table 3), the ability of γ to determine whether SF (SAR feature) is sensitive to lodged rice was also confirmed.
Correlation scatterplots of SF screened by γ were constructed and each scatterplot included ten points ( Figure 7). As the X and Y coordinate values of each point denotes SF values of one lodged sample plot before and after the typhoon, respectively, the distribution of these points reflects the consistency of SF changes after rice lodging. Here, VV, VH, VV + VH, VV − VH, Shannon Entropy and Span values changed after rice lodging, recording good consistency (scatter points are concentrated on the onside of line 1:1) (Figure 7). The value of β also indicated the consistency of parameter changes after rice lodging. With regard to SI, based on the criterion of γ > 0 (Table 4), all spectral bands apart for Band 4 (Red) were selected as sensitive parameters to rice lodging, indicating that the majority of the band indices of Sentinel-2 are sensitive to lodged rice. A similar conclusion was recorded in a study mapping lodged wheat [20]. However, different from SF, β and γ for SI did not have a strong consistency for judging rice lodging sensitivity as, except for Band 3 (Green), β values for all other spectral bands reached 10 (including Band4). The higher β values of SI, therefore, were not credible because differences in illumination between 2 and 17 September interfered with results (see Section 4.2). A total of 14 parameters were selected using the first and second steps of the OSPL screening model, including six SF (VV, VH, VV + VH, VV − VH, Shannon Entropy and Span) and eight SI (VIS-Blue (Band 2), Vegetation Red Edge (Bands 5, 6, 7 and 8A), NIR (Band 8) and SWIR (Bands 11 and 12)). Parameters with good performance for distinguishing healthy and lodged rice (detailed information is provided in Section 3.2.2) were obtained using the third step. Finally, VH, VV + VH, Shannon Entropy, SWIR1 (Band 11) and SWIR2 (Band 12) were selected as optimal sensitive parameters to rice lodging using the OSPL screening model.
Discrimination scatter plots of all sample plots (S1-S10) for SF and SI (Figures 8 and 9) were selected using screening steps 1 and 2, respectively. Results indicate that all selected OSPLs had almost no overlap between healthy and lodged sample rice plots, indicating that OSPL has sufficient ability to distinguish lodged rice from healthy rice.
In this study, σ 0 VV and σ 0 VH were found to be sensitive to lodged rice (see Section 4.1). Similar conclusions were recorded when the C-band spaceborne SAR was used to monitor wheat and corn lodging [16,18,20]. Han et al. also used VH and VV + VH as optimal polarization parameters to classify lodged corn with a different lodging degree [15].

Mapping Healthy and Lodged Rice Based on RF and OSPL
Rice region extraction was conducted using optical images and MLC (Figure 10a). The study area contained a total of 2,187,630 pixels with 10 m space resolution, of which 1,754,260 pixels (80.19%) were paddy rice and the remaining 433,370 pixels (19.81%) were non-rice region (Figure 10a). The recognition ability of MLC to various and small features was identified using the detailed classification features of bare land, water body, roads and buildings (Figure 10a). The shape of the paddy rice region was used to eliminate the non-rice region of OSPL image, and lodged rice mapping was undertaken using a fine-tuned RF classifier (Section 3.3) according to OSFs and OSIs selected by the OSPL selection model. Three cases were designed to evaluate the capability of SAR and optical data in lodged rice mapping. Case I: Optimal SAR features were used for classification; Case II: Optimal optical spectral indices were used as the inputs of classification; Case III: The integration of optimal SAR and optical parameters were used as the inputs of classification. Accuracy results of lodged rice mapping using the three cases are shown in Table 5.  Through the observation and comparison for intersection pixel distribution of ground truth (GT) and classification category (CC) and values of PA and UA, we found that commission error of lodge class and omission error of healthy class were larger in Case I, while commission error of healthy class and omission error of lodge class were larger in Case II. Finally, the commission and omission errors for healthy and lodge classes were reduced by combining SAR and optical information with mutually complement characteristics, and the values of OA and kappa coefficient also increased in Case III compared to Case I and II. The integrated OSF-OSI classification result image (Figure 10b) indicates lodging rice (red region), healthy rice (green regions) and non-rice regions (white regions). In Case III, lodging rice accounted for 34.63% of the overall paddy rice region, being concentrated in the Northeast, Northwest and Central regions in study area (Figure 10b).
This method used free-to-download Sentinel-1/2 spaceborne data and a whole set of processing flow paths executed in available software, and RF classification was employed based on fewer adjusted parameters. However, interferences of cloud and mist to Sentinel-2 optical data and heavy rain to Sentinel-1 SAR data limit the applicability of this method.

Conclusions
In this study, multiple SAR and optical parameters extracted from Sentinel-1/2 were analyzed for lodged rice sensitivity. An OSPL selection model was developed to select optimal parameters for lodged rice mapping, and five parameters were selected based on OSPL screening model, of which three parameters were OSF (VH, VV + VH and Shannon Entropy) and the remaining were OSI (Band 11 and Band 12). Finally, a higher classification accuracy was achieved using ORF-OSI combination images based on fine-tuned RF classifier. Some of the findings are as follows: (i) The VV and VH polarization signals of Sentinel-1 C-band data might be significantly enhanced under heavy storms in farmland regions; (ii) alpha (α), entropy (H) and anisotropy (A) of H/A/Alpha polarimetric decomposition based on dual-polarized data have a low level of sensitivity to lodged rice; (iii) VV and VH of Sentinel-1 are largely sensitive to lodging rice due to structural changes of canopy and stalks. The majority of spectral bands were sensitive to lodging rice mainly due to biophysical/biochemical and rice canopy architecture changes.
The limitations of this work and areas for improvement in future investigations have the following points: (i) In this study, the time interval between the typhoon transit time and acquisition time of Sentinel-1 and Sentinel-2 image data was long; available data recorded closer to typhoon events should be collected and analyzed. (ii) Large scale lodged rice mapping was not undertaken as lack of the appropriate field data and Google Earth imagery; the classification effects of OSF, OSI and OSF-OSI using RFs on regional-scale can be executed further. (iii) The dual polarized information (VV and VH) used in this study has limitations for exploring the radar polarimetric response of lodging rice; fully polarimetric spaceborne data should be explored in the future as it not only has richer polarization information, it can also provide more variety of polarimetric decomposition information.
Author Contributions: J.W. and K.L. proposed and performed the experiments and wrote the paper. Y.S., F.Z. and Z.W. checked the paper and provided suggestions. X.G., Y.Q. and X.L. contributed to the data analysis,

Conclusions
In this study, multiple SAR and optical parameters extracted from Sentinel-1/2 were analyzed for lodged rice sensitivity. An OSPL selection model was developed to select optimal parameters for lodged rice mapping, and five parameters were selected based on OSPL screening model, of which three parameters were OSF (VH, VV + VH and Shannon Entropy) and the remaining were OSI (Band 11 and Band 12). Finally, a higher classification accuracy was achieved using ORF-OSI combination images based on fine-tuned RF classifier. Some of the findings are as follows: (i) The VV and VH polarization signals of Sentinel-1 C-band data might be significantly enhanced under heavy storms in farmland regions; (ii) alpha (α), entropy (H) and anisotropy (A) of H/A/Alpha polarimetric decomposition based on dual-polarized data have a low level of sensitivity to lodged rice; (iii) VV and VH of Sentinel-1 are largely sensitive to lodging rice due to structural changes of canopy and stalks. The majority of spectral bands were sensitive to lodging rice mainly due to biophysical/biochemical and rice canopy architecture changes.
The limitations of this work and areas for improvement in future investigations have the following points: (i) In this study, the time interval between the typhoon transit time and acquisition time of Sentinel-1 and Sentinel-2 image data was long; available data recorded closer to typhoon events should be collected and analyzed. (ii) Large scale lodged rice mapping was not undertaken as lack of the appropriate field data and Google Earth imagery; the classification effects of OSF, OSI and OSF-OSI using RFs on regional-scale can be executed further. (iii) The dual polarized information (VV and VH) used in this study has limitations for exploring the radar polarimetric response of lodging rice; fully polarimetric spaceborne data should be explored in the future as it not only has richer polarization information, it can also provide more variety of polarimetric decomposition information.
Author Contributions: J.W. and K.L. proposed and performed the experiments and wrote the paper. Y.S., F.Z. and Z.W. checked the paper and provided suggestions. X.G., Y.Q. and X.L. contributed to the data analysis, making charts and the paper modification, respectively. The order of the authors reflects their level of contribution. All authors have read and agreed to the published version of the manuscript.

Funding:
The work described in this paper is funded by the National Natural Science Foundation of China (41871272, 41431174) and Project supported by the Special Director Foundation of the Institute of Remote Sensing and Digital Earth, Chinese Academy of Sciences (Y6SJ1400CX).