# Polarimetric SAR time-series for identification of winter land use

J. Denize, Laurence Hubert-Moy, E. Pottier

<figure><img src="image_1.png" /></figure>

To cite this version:

J. Denize, Laurence Hubert-Moy, E. Pottier. Polarimetric SAR time-series for identification of winter land use. Sensors, 2019, 19 (24), pp.5574. 10.3390/s19245574. hal-02442758

HAL Id: hal-02442758

https://univ-rennes.hal.science/hal-02442758

Submitted on 16 Jan 2020

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_3.png" /></figure>

Article

# Polarimetric SAR Time-Series for Identification of Winter Land Use

Julien Denize<sup>1,2,*</sup>, Laurence Hubert-Moy<sup>2</sup> and Eric Pottier<sup>1</sup>

Received: 3 October 2019; Accepted: 11 December 2019; Published: 17 December 2019

<figure><img src="image_4.png" /></figure>

check for
updates

**Abstract:** In the past decade, high spatial resolution Synthetic Aperture Radar (SAR) sensors have provided information that contributed significantly to cropland monitoring. However, the specific configurations of SAR sensors (e.g., band frequency, polarization mode) used to identify land-use types remains underexplored. This study investigates the contribution of C/L-Band frequency, dual/quad polarization and the density of image time-series to winter land-use identification in an agricultural area of approximately 130 km<fnref n="2" /> located in northwestern France. First, SAR parameters were derived from RADARSAT-2, Sentinel-1 and Advanced Land Observing Satellite 2 (ALOS-2) time-series, and one quad-pol and six dual-pol datasets with different spatial resolutions and densities were calculated. Then, land use was classified using the Random Forest algorithm with each of these seven SAR datasets to determine the most suitable SAR configuration for identifying winter land-use. Results highlighted that (i) the C-Band (F1-score 0.70) outperformed the L-Band (F1-score 0.57), (ii) quad polarization (F1-score 0.69) outperformed dual polarization (F1-score 0.59) and (iii) a dense Sentinel-1 time-series (F1-score 0.70) outperformed RADARSAT-2 and ALOS-2 time-series (F1-score 0.69 and 0.29, respectively). In addition, Shannon Entropy and SPAN were the SAR parameters most important for discriminating winter land-use. Thus, the results of this study emphasize the interest of using Sentinel-1 time-series data for identifying winter land-use.

**Keywords:** crops; RADARSAT-2; sentinel-1; ALOS-2; C-band frequency; L-band frequency; dual-polarization; quad polarization; random forest classification

## 1. Introduction

The importance of vegetation cover during winter to preserve soil quality and water resources is now well-recognized by scientists, decision makers and citizens, and land-use mapping is considered a relevant input into decision-making to implement appropriate policy responses [1]. However, although identifying land use in agricultural areas is a major environmental and scientific issue [2], it remains challenging due to its high spatio-temporal dynamics [3]. In this context, remotely sensed time-series data are a valuable tool to identify land use by providing precise and timely information about the phenological status and development of vegetation at different scales, from local to global extents [4,5]. In the past few decades, progress has been made with the development of high and very high spatial and temporal resolution optical (e.g., Satellite pour l'observation de la Terre (SPOT-6/7), Sentinel-2) and Synthetic Aperture Radar (SAR) (e.g., TerraSAR-X, RADARSAT-2, Advanced Land Observing Satellite 2 (ALOS-2), Sentinel-1) sensors [6–8]. However, using optical time-series to identify land use in winter is limited by cloud cover and/or low solar irradiance [9], and late winter is a critical period during which vegetation begins to grow [3]. Conversely, SAR time-series provide a reliable solution to address the limitations of optical images because they are not sensitive to atmospheric conditions and can operate day and night [10].

<footnote n="1">¹ University of Rennes & IETR UMR 6164, 35 000, Rennes, France; eric.pottier@univ-rennes1.fr</footnote>

<footnote n="2">² University of Rennes & LETG UMR 6554, 35 000, Rennes, France; laurence.moy@univ-rennes2.fr</footnote>

<footnote>* Correspondence: julien.denize@univ-rennes1.fr</footnote>
These advantages, along with a sensitivity of microwave scattering to soil and vegetation characteristics [11,12], have led scientists to evaluate the potential of using SAR sensors to monitor agriculture [13–16]. The potential of using SAR data to identify land use is based on the sensitivity of the radar signal to the dielectric constant of the objects and to their structure (i.e., the distribution of size shape and orientation of the scatterers) [11,17]. Thus, some studies have demonstrated the ability of SAR backscattering-coefficient and polarimetric data to classify land use using the dielectric properties of soil, surface roughness and vegetation canopy structure [18–21]. Most of these studies used the C-Band frequency rather than other frequencies (e.g., X- or L-Bands). Since microwave penetration depends on wavelength and incident angle, L-Band (~20 cm) wavelengths penetrate further into crop canopies than those of the C-Band (~5.5 cm). Consequently, C-Band wavelengths interact more with the upper canopy, while L-Band wavelength responses result from greater interaction with the soil/canopy and scattering directly from the soil [11,12]. Thus, McNairn and Brisco [13] demonstrated the potential of vertically (V) and horizontally (H) polarized C-Band microwaves to identify land use. Indeed, since V microwaves respond to predominantly vertical structure, they penetrate the canopy less. Conversely, H microwaves tend to penetrate the canopy more than V microwaves. Similarly, others [17,18] illustrated the ability of the L-Band and the combined use of C- and L-Bands to classify land use. Many studies have also shown the potential of radar data to map crop residues [13,22] or crop types during summer [23,24].

However, only a few studies have demonstrated the potential of SAR time-series to identify land-use types in winter. Haldar et al. [16] illustrated the potential of polarimetric C-Band SAR time-series data to derive useful information, such as biophysical parameters, from mustard (*Sinapis alba* L.) and wheat (*Triticum aestivum* L.) crops. Similarly, few studies have shown Sentinel-1's potential to identify and characterize land-cover and land-use dynamics. For example, some studies [5,25–27] have shown benefits of using Sentinel-1 time-series to understand crop behavior and dynamics and classify land use. Likewise, Minh et al. [28] used Sentinel-1 time-series to produce a winter vegetation quality map with five classes ("bare soil" to "high quality") based on a deep-learning method, with an overall accuracy (OA) >98%. However, the specific configurations of SAR sensors (e.g., band frequency, polarization mode) used to identify land-use types remain underexplored, and even unexplored for winter land-use.

This study aimed to evaluate and compare the value of multi-temporal ALOS-2 (ALS-2), RADARSAT-2 (RST-2) and Sentinel-1 (S-1) data for monitoring winter land-use in agricultural areas. Specifically, we addressed the main question: What is the most appropriate SAR configuration (polarization, frequency, density of time-series) for characterizing winter land-use? To this end, we first calculated SAR parameters from quad-polarization (pol) C-Band, dual-pol C-Band and dual-pol L-Band time-series to generate several parameter datasets. We then applied a classification procedure using the Random Forest (RF) algorithm to determine the most discriminating SAR configuration for land-use monitoring in winter. Finally, we identified advantages and disadvantages of the method and its results.

## 2. Materials and Methods

### 2.1. *Study Site*

The study site is a relatively flat area located to the south of the Bay of Mont-Saint-Michel (8° 31′ 0″ N, 1° 31′ 30″ W) in France and covers an area of ca. 130 km² (Figure 1). This site has been included in long-term ecological research (LTER) networks at the European (LTER-Europe) and international (ILTER) levels since 1993 to assess relationships between changes in farming activities, landscape dynamics and ecological processes related to biodiversity, water quality and climate [29]. This study site is also referenced in the "Kalidéos program", which is coordinated and managed by the CNES (French Space Agency) to promote and demonstrate the use of spatial data by supporting research and development, prototyping and user demonstration activities. The site is composed of ca. 7 000
agricultural fields ranging from 0.1–65 ha and is characterized by a fragmented agricultural landscape (from a hedgerow network to open fields). This temperate climate area is exposed to mean annual precipitation of 600–700 mm and mean average temperature >12 °C. In summer, farming systems are based on one main crop per field: maize (*Zea mays* L.), wheat (*Triticum aestivum* L.), rapeseed (*Brassica napus* L.) or barley (*Hordeum vulgare* L.). In winter, in addition to grasslands, which play a major role in regulating water flows and nutrient cycling, catch crops are sown to decrease nitrogen leaching, as required by the European Union's "Nitrates Directive" [30].

<figure><img src="image_5.png" /><figcaption>**Figure 1.** Study site location, ground surveys (RGB composite image constructed from Shannon Entropy extracted from Advanced Land Observing Satellite 2 (ALOS-2) data for three dates: 03-09-2017, 04-15-2017 and 05-13-2017. ©Kalidéos data 2017 and JAXA data).</figcaption></figure>

## 2.2. Field Data

The mainland-use types encountered and investigated in winter in the study area (Figure 2 and Table 1) are the following:

* Winter crops, which cover ca. 40% of the UAA (utilized agricultural area) and include three main annual crops: winter wheat, winter barley and rapeseed

* Grasslands, which cover ca. 30% of the UAA and can be mown or grazed

* Catch crops, which are sown after harvest of the main crop from August to October, cover ca. 25% of the UAA and include a wide variety of crops

* Crop residues, which cover ca. 5% of the UAA and correspond to maize stalks that are left in fields when the maize is harvested after 1 November [3]

Field observations were conducted in 231 crop fields monthly from November 2016 to February 2017 to calibrate and validate classification of remote sensing data (Figure 1). Samples were randomly distributed throughout the study site, with half of the fields (116) being used for training and the other half for validation. The fields inventoried ranged in size from 0.1–65 ha. The number of training samples of each winter land-use class inventoried was 54 for winter crops, 14 for grasslands, 42 for catch crops and 6 for crop residues.
<figure><img src="image_6.png" /><figcaption>**Figure 2.** The mainland-use types encountered in winter in the study area: (a) winter crops (winter barley), (b) catch crops (mustard), (c) grasslands and (d) crop residues (maize stalks).</figcaption></figure>

<table><caption>**Table 1.** Land-use classification.</caption><thead><tr><th>Winter Land Use Type</th><th>Main Crops</th></tr></thead><tbody><tr><td>Winter crops</td><td>Winter wheat<br>Winter barley<br>Rapeseed</td></tr><tr><td>Grasslands</td><td>Mown grasslands<br>Grazed grasslands</td></tr><tr><td>Catch crops</td><td>Oat<br>Fodder cabbage<br>Ryegrass and clover<br>Phacelia<br>Phacelia and mustard<br>Phacelia and oat</td></tr><tr><td>Crop residues</td><td>Maize stalks</td></tr></tbody></table>

## 2.3. Satellite Data

### 2.3.1. RADARSAT-2 Time-Series

A series of 10 RST-2 SAR images from October 2016 to May 2017 were acquired by MacDonald, Dettwiler and Associates and provided by the VIGISAT project managed by Collecte Localisation Satellites in the framework of the GIS Bretel ("Groupement Bretagne Télédétection", Bretagne, France) (Table 2). RST-2 images were acquired in Single Look Complex (SLC) mode (delivered in quad-pol mode: HH, HV, VH and VV polarization states) with an incidence angle of 35°. The range and azimuth spatial resolutions were 8.2 and 4.7 m, respectively.

### 2.3.2. Sentinel-1 Time-Series

Two series of eight and twenty S-1 SAR images, respectively, from August 2016 to May 2017 were acquired by the European Space Agency and provided by its data hub [31] (Table 2). S-1 images were acquired in SLC mode (delivered in dual-pol mode: VH and VV) with an incidence angle of 31° to 46° and an angle of 40° on the study area. The range and azimuth spatial resolutions were 2.3 and 13.9 m, respectively.

### 2.3.3. ALOS-2 Time-Series

A series of six ALS-2 SAR images from January–June 2017 were acquired by the Japan Aerospace Exploration Agency and provided by the Kalidéos program [32] (Table 2). ALS-2 images were acquired in SLC mode (delivered in dual-pol mode: HH and HV) with an incidence angle of 40° and range and azimuth spatial resolutions of 1.4 and 1.9 m, respectively.
<table><caption>**Table 2.** Characteristics of the RADARSAT-2, Sentinel-1 and Advanced Land Observing Satellite 2 (ALOS-2) images used in the study. Dates in bold text for Sentinel-1 images indicate those with sparse time-series, while asterisks indicate those with dense time-series.</caption><thead><tr><th></th><th>RADARSAT-2</th><th>Sentinel-1</th><th>ALOS-2</th></tr></thead><tbody><tr><td></td><td></td><td><strong>08-25-2016</strong><br>09-18-2016*<br>09-30-2016*<br>10-12-2016*<br>10-24-2016*<br><strong>11-05-2016*</strong><br>11-17-2016*<br>11-29-2016*<br>12-11-2016*<br><strong>12-23-2016*</strong><br>01-04-2017*<br><strong>01-16-2017*</strong><br>01-28-2017*<br>02-09-2017*<br><strong>02-21-2017*</strong><br>03-05-2017*<br>03-17-2017*<br><strong>03-29-2017*</strong><br><strong>04-10-2017*</strong><br>04-22-2017*<br>05-04-2017*<br><strong>05-16-2017</strong></td><td></td></tr><tr><td><strong>Dates</strong><br>(M-D-Y)</td><td>10-23-2016<br>11-16-2016<br>12-10-2016<br>01-03-2017<br>01-27-2017<br>02-20-2017<br>03-16-2017<br>04-09-2017<br>05-03-2017<br>05-27-2017</td><td></td><td>01-04-2017<br>02-04-2017<br>03-06-2017<br>04-15-2017<br>05-13-2017<br>06-10-2017</td></tr><tr><td><strong>Ground Resolution</strong></td><td>8.2 m</td><td>2.3 m</td><td>1.4 m</td></tr><tr><td><strong>Azimuth Resolution</strong></td><td>4.7 m</td><td>13.9 m</td><td>1.9 m</td></tr><tr><td><strong>Polarization</strong></td><td>Quad (HH-VV-HV-VH)</td><td>Dual (VV - VH)</td><td>Dual (HH-HV)</td></tr><tr><td><strong>Frequency</strong></td><td>C-Band</td><td>C-Band</td><td>L-Band</td></tr><tr><td><strong>Mode</strong></td><td>Fine Quad Polarization (SLC)</td><td>Interferometric wide (SLC)</td><td>Spotlight (SLC)</td></tr><tr><td><strong>Incidence Angle</strong></td><td>35° (right descending)</td><td>31° to 46° (right descending)</td><td>40° (left ascending)</td></tr><tr><td><strong>Coverage</strong></td><td>18 km × 25 km</td><td>>250 km × 100 km</td><td>25 km × 25 km</td></tr></tbody></table>

## 2.4. Extraction of SAR Parameters

SAR quad-pol and dual-pol parameters from the three sensor datasets (RST-2, S-1 and ALS-2) were calculated using Sentinel Application Platform (SNAP) v6.0 (Develop by Esa, Paris, France) and PolSARpro v5.1.3 software [33] (Develop by IETR, Rennes, France).

### 2.4.1. Quad-Pol Time-Series

* Backscattering coefficients ($\sigma^0HH$, $\sigma^0HV$, $\sigma^0VH$, $\sigma^0VV$) were calculated from the radiometrically calibrated RST-2 time-series using SNAP according to Equation (1) [34]:

$$ \sigma_j^o = \beta_j^o + 10 \times \log_{10}(\sin I_j) \quad (1) $$

where $\beta$ is the radar brightness and $I_j$ is the incidence angle at the $j^{th}$ range pixel.

Equation (1) assumes that the Earth is a smooth ellipsoid at sea level. A Lee Sigma filter [35] was applied with a window of 7 × 7 pixels and a sigma value of 0.8. RST-2 images were then geocoded at an 8-m resolution using Shuttle Radar Topography Mission 3s data to correct topographic deformations. The accuracy of geometric correction was less than 8 m per pixel. Next,
two backscattering ratios were calculated ($\sigma^\circ HH: \sigma^\circ VV$, $\sigma^\circ HH: \sigma^\circ HV$) that highlight scattering mechanisms of each target.

* Polarimetric parameters were calculated from SLC RST-2 time-series. First, a 3 × 3 coherency matrix $T_3$ was extracted from the scattering matrix (S) of each image using PolSARpro. Next, a Lee Sigma filter was applied with a window of 7 × 7 pixels and a sigma value of 0.8. The elements of the matrix, which are independent of the polarimetric absolute phase [36], were then geocoded directly using SNAP with an 8-m resolution.

Second, Cloude-Pottier decomposition [37] was then calculated based on the $T_3$ matrix. From the eigenvalues extracted, we calculated three independent parameters: (i) entropy (H), which expresses the randomness of the scatter; (ii) alpha angle (α), which describes the dominant scattering mechanism and (iii) anisotropy (A), which represents the relative power of the dominant mechanism. In addition, four parameters based on the Cloude-Pottier decomposition ($H \times A$; $H \times (1-A)$; $(1-H) \times A$; $(1-H)(1-A)$) were calculated because they can provide the number of scattering mechanisms in each resolution cell.

Third, Freeman-Durden decomposition [38] was used to model the 3 × 3 covariance matrix ($C_3$) as the contribution of three scattering mechanisms for each pixel: volume, double-bounce and surface/single-bounce. Fourth, SPAN (total scattered power) and Shannon Entropy (SE), which equals the sum of two parameters related to the intensity ($SE_i$) and degree of polarization ($SE_p$) [36], were calculated from the $T_3$ matrix. SE measures the disorder encountered in polarimetric SAR images.

Finally, two polarimetric parameters were extracted from the $T_3$ matrix: pedestal height and the Radar Vegetation Index (RVI). Pedestal height is the ratio of the maximum received intensity to the minimum received intensity; it indicates the presence of unpolarized scattering and thus the degree of polarization of a scattered wave [39]. RVI is a function of incidence angle, since the path length through the vegetation canopy will increase as the incidence angle increases [40]. Ranging from 0 to 1, RVI measures the randomness of the scatter according to Equation (2):

$$RVI = \frac{8\sigma^0HV}{\sigma^0HH + \sigma^0VV + 2\sigma^0HV} \qquad (2)$$

A total of 25 quad-pol parameters were calculated for each of the 10 quad-pol 8-m RST-2 images, yielding a dataset with 250 variables.

### 2.4.2. Dual-Pol Time-Series

The dual-polarization preprocessing step included (i) conversion of RST-2 time-series images from quad- to dual-polarization mode, (ii) resampling of RST-2 and ALS-2 time-series images at the spatial resolution of S-1 images and (iii) extraction of dual-pol C-Band and dual-pol L-Band SAR parameters.

* Converting the RST-2 time-series polarization mode: Each 3 × 3 coherency matrix $T_3$ extracted from RST-2 quad-polarization images was converted to a 2 × 2 covariance matrix $C_2$ using PolSARpro. The converted RST-2 images had the same polarizations (HH and HV) as ALS-2 images.

* Calculating S-1 and ALS-2 covariance matrices: A 2 × 2 covariance matrix ($C_2$) was extracted from the two polarizations of each ALS-2 2-m image and each S-1 image (for S-1 dense and sparse time-series).

* Resampling RST-2 and ALS-2 time-series: A multi-looking function was applied using a 2 × 1 pixel window for RST-2 images (i.e., 2 × 4.7 m and 1 × 8.2 m) and a 5 × 6 pixel window for ALS-2 images (i.e., 5 × 1.9 m and 6 × 1.4 m) using PolSARpro. Resampled RST-2 images had a resolution of 9.4 × 8.2 m, which was close to that of resampled ALS-2 images (9.5 × 8.4 m) and 10-m corrected S-1 images.

* Calibrating backscattering coefficients: Backscattering coefficients $\sigma^\circ HH$ and $\sigma^\circ HV$ were simultaneously calibrated radiometrically from dual-pol converted and resampled RST-2
images and original and resampled ALS-2 images using SNAP. Backscattering coefficients σ°VV and σ°VH were simultaneously calibrated radiometrically from dual-pol S-1 images. Then, a Lee Sigma filter [35] with a window of 7 × 7 pixels and a sigma value of 0.8 was applied to all images to attenuate speckle noise. Next, geometric correction was performed using the Shuttle radar topographic mission (SRTM) for each time-series dataset, with a 2-m resolution for original ALS-2 images; 8-m resolution for original RST-2 images and 10-m resolution for resampled ALS-2, resampled RST-2 and the two S-1 images. Finally, the σ°HH : σ°HV ratio and σ°HH-σ°HV difference were calculated from ALS-2 and RST-2 backscattering coefficients, and the σ°VH : σ°VV ratio and σ°VH-σ°VV difference were calculated from S-1 backscattering coefficients.

*   Extracting polarimetric parameters: Dual-polarimetric parameters were simultaneously extracted from dual-pol S-1 images, converted and resampled RST-2 images and original and resampled ALS-2 images. To this end, the same Lee Sigma filter was applied to the C₂ matrices to filter out speckle noise. Then, geometric corrections were applied to all polarimetric parameters using the SRTM with the previously used 2-m, 8-m and 10-m resolutions. SPAN, SE, SEi, SEp, normalized SE, normalized SEi and normalized SEp were also extracted.

A total of 11 dual-pol parameters (2 backscattering coefficients, 1 ratio and 1 difference based on the backscattering coefficients, 1 SPAN, 2 SE, 2 SEi and 2 SEp) were calculated for each of the 10 RST-2, 6 ALS-2 and 8 or 20 S-1 images. Finally, six dual-pol datasets were created:

1.  original RST-2 8-m dataset with 110 variables (11 parameters × 10 dates)

2.  resampled RST-2 10-m dataset with 110 variables (11 parameters × 10 dates)

3.  original ALS-2 2-m dataset with 66 variables (11 parameters × 6 dates)

4.  resampled ALS-2 10-m dataset with 66 variables (11 parameters × 6 dates)

5.  sparse S-1 10-m dataset with 88 variables (11 parameters × 8 dates)

6.  dense S-1 10-m dataset with 220 variables (11 parameters × 20 dates)

2.5. Classification of SAR Parameter Datasets

A two-step approach was performed to identify winter land-use using the seven SAR parameter
datasets (one quad-pol and six dual-pol). First, parameter importance was analyzed for each SAR
parameter dataset using the RF importance function, The RF importance function was performed to
rank the features in order of importance based on the mean decrease in the Gini index (Breiman et al.,
2001) repeated 100 times by changing the training sample, to identify the parameters most important for
identifying winter land-use types. A measure of variable importance was provided for each candidate
predictor and each classification using the heuristic method based on the Gini Index [41–43].

Second, RF, as a supervised classification algorithm, was used to classify land-use during winter 2016–2017. The RF was chosen for its high performance and accurate classification of land use [44]. RF is an ensemble classifier that uses classification and regression trees to make predictions [42]. The trees are created by drawing a subset of training samples through replacement (a bagging approach). In this way, some samples may be selected several times, while other samples may not be selected at all. The “randomForest” package (v.4.6–14) developed by [45] and implemented in R software (v.3.5.1) [46] (Develop by Bell Laboratories, Murray Hill, New Jersey, United-Sates) was used to perform classifications. Two RF parameters were tuned using the “name” function in this package. The first, the number of trees randomly created using the training dataset, was set to 1000, since the number of errors decreases little with more than 1000 trees [47]. The second parameter, the number of variables randomly sampled as candidates at each split node (mtry), was defined starting with a mtry equals the square root of the number of input variables, then searching for the optimal value to improve the quality of the model.

From each image (one per date), 1200 samples (pixels) were randomly selected. Half (600 "in-bag"
samples) were used to train the trees, while the rest (600 "out-of-the-bag" samples) were used to
estimate the accuracy of the RF model [48]. This process being repeated 100 times with replacement of samples, we can consider that bagging was used during the overall process [3]. Classification accuracy was assessed using the F1-score and the Kappa index. The F1-score is a standard measure of classification accuracy defined as the weighted average of precision and recall [49], while the Kappa index expresses the proportional decrease in error generated by the classification compared to the error of a completely random classification [50].

A four-step classification procedure was performed to evaluate the potential of SAR time-series to identify winter land-use types:

1. RST-2 quad-pol, S-1 dense and ALS-2 2-m time-series datasets were classified to demonstrate the full potential of these SAR sensors

2. RST-2 10-m dual-pol, ALS-2 10-m and S-1 sparse time-series datasets were classified to identify the best band frequency

3. RST-2 quad-pol and RST-2 8-m dual-pol time-series datasets were classified to identify the best polarization mode

4. S-1 dense and sparse time-series datasets were classified to evaluate the influence of the number of images

## 3. Results

### 3.1. Importance of SAR Parameters for Discriminating Winter Land-Use

The most important SAR parameters for discriminating winter land-use depended on the polarization mode and band frequency (Figures 3 and 4). Concerning polarization mode, the most important parameters in the quad-pol data were SPAN (7.8%), SE (7.7%), normalized SE (7.5%) and, to a lesser extent, SEᵢ (5.9%) and the normalized SEᵢ (5.8%) (Figure 3). Conversely, the least important parameters were the Freeman-Durden double-bounce (2.1%) and anisotropy (2%). Parameters related to backscattering coefficients were less important (4.3% and 2.4% for VV backscattering and the HV:HH ratio, respectively). For the dual-pol data, the most important parameters were SE and normalized SE, with an importance of 11.5% for both for the ALS-2 L-Band; 13.1% and 12.5%, respectively, for RST-2 and 11.7% and 11.6%, respectively, for S-1. In addition, for the ALS-2 L-Band, SPAN was the most important parameter (11.7%). Conversely, the least important parameter was the HH:HV ratio for the ALS-2 L-Band and RST-2 C-Band (5.6% and 6.4%, respectively) and the difference VV – VH for the S-1 C-Band (6.7%).

<figure><img src="image_7.png" /><figcaption>**Figure 3.** Importance (in %) of quad-pol SAR parameters based on 100 random forest classifications. Parameters related to backscattering coefficients are in black, while polarimetric parameters are in gray. SE: Shannon Entropy.</figcaption></figure>
Concerning band frequency, in the dual-pol L-Band (Dual-L) configuration, parameters related to backscattering coefficients were less important than polarimetric parameters (mean importance of 6.6% and 10.5%, respectively) (Figure 4). In the dual-pol C-Band (Dual-C) configuration, the ranks of the parameters related to backscattering coefficients were similar for RST-2 and S-1 classification models. For each of them, a backscattering coefficient was important: the HV coefficient for RST-2 and VV coefficient for S-1 (10.4% and 10.7%, respectively). Similarly, the HH:HV and VV:VH ratios were less important (6.4% and 6.7%, respectively). The importance of polarimetric parameters was similar for Dual-C and Dual-L configurations, with SE and SEᵢ being the most important parameters. However, SPAN was the most important parameter for Dual-L but one of the least important for Dual-C.

<figure><img src="image_8.png" /><figcaption>**Figure 4.** Importance (in %) of dual-pol SAR parameters based on 100 random forest classifications using (A) ALOS-2 parameters, (B) RADARSAT-2 parameters and (C) Sentinel-1 parameters. Parameters related to backscattering coefficients are in black, while polarimetric parameters are in gray. SE: Shannon Entropy.</figcaption></figure>
## 3.2. Contribution of Polarization Mode to Accuracy of Winter Land-Use Classification

The F1-score of land-use classification was higher using the quad-pol dataset (median F1-score 0.69, standard deviation (SD) 0.05; median Kappa 0.68, SD 0.02;) than the dual-pol datasets (median F1-score 0.59, SD 0.05; median Kappa 0.54, SD 0.02) (Figure 5). The superiority of quad-pol mode over dual-pol mode was observed for all winter land-use classes (Figure 5). The largest differences in median F1-score between quad-pol and dual-pol mode concerned the oat (0.64 and 0.50, respectively), ryegrass and clover (0.60 and 0.46, respectively) and phacelia (0.75 and 0.66, respectively) classes. Conversely, the rapeseed class was discriminated slightly better using the quad-pol rather than dual-pol dataset (median F1-score 0.87 and 0.85, respectively). Moreover, Standard Deviation (SD) in F1-score accuracy was slightly lower using the quad-pol rather than dual-pol dataset for fodder cabbage (0.05 and 0.06, respectively) and crop residue (0.03 and 0.04, respectively) classes, but the opposite is true for other classes such as ryegrass and clover (0.05 and 0.04, respectively) (Figure 5).

<figure><img src="image_9.png" /><figcaption>**Figure 5.** Comparison of classification accuracy of each land-use class between dual and quad polarization (pol) modes. Box-and-whisker plots represent the variation in random forest classification accuracy based on 100 iterations. Whiskers indicate 1.5 times the interquartile range.</figcaption></figure>
### 3.3. Contribution of Band Frequency to Accuracy of Winter Land-Use Classification

The F1-score of land-use classification was slightly higher for the RST-2 C-Band frequency (median F1-score 0.64, SD 0.05; median Kappa 0.64, SD 0.02) than S-1 C-Band frequency (median F1-score 0.61, SD 0.05; median Kappa 0.59, SD 0.02) or the ALS-2 L-Band frequency (median F1-score 0.57, SD 0.05; median Kappa 0.60, SD 0.02) (Figure 6). More precisely, classification accuracy depended on the band frequency and land-use class (Figure 6). F1-score was highest for the C-Band frequency, except for the crop residue, phacelia, fodder cabbage and mown grassland classes, for which the ALS-2 L-Band frequency had the highest accuracy (median F1-score 0.89, 0.80, 0.77 and 0.49, respectively). Classification accuracy using the S-1 parameter dataset was highest for grazed grassland, winter wheat and oat classes (median F1-score 0.83, 0.65 and 0.57, respectively). For the other classes (phacelia and mustard, ryegrass and clover, phacelia and oat, rapeseed and winter barley), F1-score was highest using the RSR-2 parameter dataset. Regardless of the band frequency, the F1-score of the rapeseed, grazed grassland and crop residue classes was higher than those of the other classes (Figure 6). The Standard Deviation in classification accuracy of land-use classes was similar among the ALS-2, RST-2 and S-1 models (0.05, 0.049 and 0.05, respectively).

<figure><img src="image_10.png" /><figcaption>**Figure 6.** Comparison of classification accuracy of each land-use class among band frequencies. Box-and-whisker plots represent the variation in random forest classification accuracy based on 100 iterations. Whiskers indicate 1.5 times the interquartile range.</figcaption></figure>
### 3.4. Contribution of Time-Series Density to Accuracy of Winter Land-Use Classification

The F1-score of land-use classification was higher for the dense S-1 time-series (median F1-score 0.70, SD 0.05; median Kappa 0.69, SD 0.02) than for the sparse S-1 time-series (median F1-score 0.61, SD 0.05; median Kappa 0.60, SD 0.02) (Figure 7). The largest differences between dense and sparse time-series concerned the phacelia and mustard (median F1-score 0.59 and 0.36, respectively), phacelia and oat (median F1-score 0.60 and 0.39, respectively) and phacelia (median F1-score 0.80 and 0.66, respectively) classes. Moreover, the variance in classification accuracy was slightly lower for the dense rather than sparse S-1 time-series of the phacelia and oat (SD 0.06 and 0.07, respectively) and grazed grassland (SD 0.03 and 0.04, respectively) classes (Figure 8). Conversely, the variance in classification accuracy in sparse and dense time-series of some winter land-use classes, such as winter wheat (SD 0.04 and 0.05, respectively) and winter barley (SD 0.05 and 0.06, respectively), was higher using the dense time-series.

<figure><img src="image_11.png" /><figcaption>**Figure 7.** Comparison of classification accuracy of each land-use class between sparse and dense Sentinel-1 time-series. Box-and-whisker plots represent the variation in random forest classification accuracy based on 100 iterations. Whiskers indicate 1.5 times the interquartile range.</figcaption></figure>
<figure><img src="image_12.png" /><figcaption>**Figure 8.** Comparison of classification accuracy of each land-use class among SAR sensors. Box-and-whisker plots represent the variation in RF classification accuracy based on 100 iterations. Whiskers indicate 1.5 times the interquartile range.</figcaption></figure>

### 3.5. *Definition of the Best SAR Configuration*

The F1-score of land-use classification was slightly higher using the S-1 dense time-series (median F1-score 0.70, SD 0.00; median Kappa 0.69, SD 0.02) than using the RST-2 quad-pol (median F1-score 0.69, SD 0.05; median Kappa 0.68, SD 0.02) or ALS-2 original dual-pol (median F1-score 0.29, SD 0.05; median Kappa 0.32, SD 0.02) time-series (Figure 8). More precisely, classification accuracy depended on the sensor and land-use class (Figure 8). The F1-score was higher with the S-1 model, except for the winter wheat, winter barley, rapeseed, grazed grassland and crop residue classes, for which the RST-2 model had higher accuracy. The accuracy of the S-1 model was slightly higher (by 0.01 in the F1-score 1%) than that of RST-2 for the phacelia and mustard and mown grassland classes, while that of the ALS-2 model was always lower than those of S-1 and RST-2. Finally, the variance in classification accuracy was similar among the S-1, RST-2 and ALS-2 models (SD 0.05, 0.05 and 0.05, respectively).
## 3.6. Spatial Distribution of Winter Land-Use Classes

Winter land-use classes were mapped at the 1:100 000 scale using the best classification model: the RF classification using the SAR parameter dataset extracted from S-1 dense time-series images. While grasslands were classified well (F1-score, 0.71), some misclassification errors and artifacts occurred, mainly between catch-crop classes, such as oat and phacelia and oat classes (Table 3). In general, catch crops and winter crops were located on the largest fields, while grasslands were located on the smallest ones (Figure 9).

<table><caption>**Table 3.** Confusion matrix of winter land-use obtained from RF classification using the SAR parameter dataset extracted from the S-1 dense time-series dataset.</caption><thead><tr><th>Land Use</th><th>1</th><th>2</th><th>3</th><th>4</th><th>5</th><th>6</th><th>7</th><th>8</th><th>9</th><th>10</th><th>11</th><th>12</th><th>Commission Error (%)</th></tr></thead><tbody><tr><td>1: Winter wheat</td><td>32</td><td>1</td><td>4</td><td>1</td><td>0</td><td>0</td><td>2</td><td>0</td><td>1</td><td>0</td><td>2</td><td>0</td><td>74.4</td></tr><tr><td>2: Winter barley</td><td>11</td><td>22</td><td>0</td><td>6</td><td>2</td><td>2</td><td>4</td><td>0</td><td>1</td><td>0</td><td>5</td><td>0</td><td>41.5</td></tr><tr><td>3: Rapeseed</td><td>1</td><td>3</td><td>49</td><td>2</td><td>1</td><td>3</td><td>2</td><td>2</td><td>1</td><td>1</td><td>0</td><td>1</td><td>74.2</td></tr><tr><td>4: Mown grasslands</td><td>0</td><td>3</td><td>2</td><td>28</td><td>7</td><td>0</td><td>1</td><td>0</td><td>0</td><td>2</td><td>1</td><td>2</td><td>60.9</td></tr><tr><td>5: Grazed grasslands</td><td>0</td><td>0</td><td>0</td><td>0</td><td>44</td><td>0</td><td>1</td><td>2</td><td>1</td><td>0</td><td>0</td><td>0</td><td>91.7</td></tr><tr><td>6: Oat</td><td>0</td><td>2</td><td>0</td><td>5</td><td>0</td><td>33</td><td>8</td><td>0</td><td>1</td><td>4</td><td>2</td><td>1</td><td>58.9</td></tr><tr><td>7: Phacelia and oat</td><td>0</td><td>1</td><td>0</td><td>2</td><td>0</td><td>4</td><td>26</td><td>0</td><td>5</td><td>3</td><td>3</td><td>0</td><td>59.1</td></tr><tr><td>8: Fodder cabbage</td><td>1</td><td>0</td><td>1</td><td>3</td><td>0</td><td>1</td><td>1</td><td>41</td><td>2</td><td>0</td><td>2</td><td>1</td><td>77.4</td></tr><tr><td>9: Ryegrass and clover</td><td>1</td><td>1</td><td>1</td><td>1</td><td>2</td><td>2</td><td>0</td><td>4</td><td>36</td><td>0</td><td>1</td><td>0</td><td>73.5</td></tr><tr><td>10: Phacelia</td><td>0</td><td>0</td><td>1</td><td>0</td><td>0</td><td>1</td><td>5</td><td>0</td><td>0</td><td>41</td><td>4</td><td>0</td><td>78.9</td></tr><tr><td>11: Phacelia and mustard</td><td>1</td><td>0</td><td>0</td><td>0</td><td>0</td><td>1</td><td>0</td><td>0</td><td>0</td><td>0</td><td>34</td><td>1</td><td>91.9</td></tr><tr><td>12: Crop residues</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>3</td><td>0</td><td>0</td><td>2</td><td>0</td><td>4</td><td>44</td><td>83.0</td></tr><tr><td>Omission error (%)</td><td>68.1</td><td>66.7</td><td>84.5</td><td>58.3</td><td>78.6</td><td>66.0</td><td>52.0</td><td>83.7</td><td>72.0</td><td>80.4</td><td>58.6</td><td>88.0</td><td>71.7</td></tr></tbody></table>

<figure><img src="image_13.png" /><figcaption>**Figure 9.** Map of winter land-use classes obtained using a parameter dataset derived from the Sentinel-1 dense time-series. Classification was performed using the random forest algorithm.</figcaption></figure>
# 4. Discussion

This study aimed to evaluate and compare the value of multi-temporal ALS-2, RST-2 and S-1 data for monitoring winter land-use in agricultural areas. Winter land-use classes, which depend on vegetation characteristics and phenology, can be discriminated using the specific properties of SAR sensors to identify crops [5,51,52]. A two-step method was applied to the study site. First, SAR parameters were derived from RST-2, S-1, and ALS-2 time-series, and one quad-pol and six dual-pol parameter datasets with different spatial resolutions and densities were calculated from these time-series. Then, land use was classified using the RF algorithm with each of these seven SAR parameter datasets to determine the most suitable SAR configuration for identifying winter land-use.

## 4.1. Which SAR Configuration for Mapping Winter Land-Use?

Our study took advantage of the potential of SAR sensor characteristics to identify winter land-use types. The physical properties and specific resolutions (spatial, temporal and frequency) of ALS-2, RST-2 and S-1 sensors allowed us to discriminate winter land-use classes in our study area, including areas where the agricultural landscape was fragmented. However, some SAR configurations appear to be more effective than others.

Results highlight that the C-Band outperformed the L-Band, the quad-polarization mode outperformed the dual-polarization mode and S-1 dense time-series outperformed RST-2 and ALOS-2 time-series. Overall, the dense S-1 time-series was the most suitable SAR configuration data for winter land-use identification, with classification accuracy slightly higher than those of the others. This result is consistent with other studies that have reported the added value of S-1 data for mapping cropland [44,53] or identifying and characterizing crop phenology [27]. Nevertheless, the quad-pol C-Band RST-2 dataset also showed potential to identify winter land-use, with accuracies similar to those obtained using the dense dual-pol C-Band S-1 dataset. This result is consistent with the studies of [54,55], which also showed the potential of polarimetric RST-2 data for monitoring and classifying crops. Conversely, the dual-pol L-Band ALS-2 dataset appeared to be of little use for identifying and mapping land-use in winter, except for certain classes such as crop residues or phacelia and mustard. These results add information to the research to date, which has demonstrated the ability of the L-Band to classify crops in summer [17,56]. However, these results should be interpreted with caution because, as shown by the results obtained with S-1 images, time-series density has a significant impact on the identification of winter land-use classes, as shown by [57]. Thus, if the classification score obtained with the ALS-2 time-series is lower than that obtained with the RST-2 and S-1 sparse time-series (F1-score 0.57 vs 0.61 and 0.64, respectively), it should be noted that the ALS-2 time-series had the least number of images (6 against 10 and 8 for RST-2 and S-1, respectively). In addition, the difference between the incidence angles of the images, although small (35° for RTS-2 and 40° for S-1 and ALS-2), can have an impact on the accuracy of the results, as noted by [58,59].

Results indicate that SE and SPAN were the most important parameters, while difference and ratio parameters had low importance. These results are consistent with the research conducted to date. [60] demonstrated the importance of SE in studying land use and land cover. Some studies also demonstrated the ability of SE to characterize the canopy of vegetated areas [61]. Likewise, many studies (e.g., [62], [63]) have described the ability of SPAN parameters to classify crops.

The best accuracy, obtained using the dense dual-pol C-Band S-1 dataset (OA 72%), demonstrated the limits of using one SAR dataset alone to identify winter land-use accurately. Nevertheless, our results demonstrate the ability of other band frequencies to discriminate land-use classes. For example, ALS-2 data appear to be useful for classifying crop residue and phacelia and mustard classes. Similarly, RST-2 data appear to be useful for classifying rapeseed and grazed grassland classes. These results highlight the utility of combining SAR parameters. For example, [17] demonstrated the utility of combining C-Band and L-Band SAR data for classifying crops. Hence, use of combined ALS-2, RST-2 and S-1 SAR datasets should be considered to obtain a denser time-series to discriminate between land-use types better.
## 4.2. *Advantages and Disadvantages of the Classification Approach*

Results demonstrated the potential of this classification approach based on the RF algorithm, although the highest F1-score achieved was 0.70. Some studies based on SAR data classified using the RF algorithm also highlighted this complexity in identifying winter crops, unlike for annual crops such as maize [5,27]. Moreover, the results also confirm the ability of SAR data to classify crop residues, as some previous studies have shown [19,64].

Several other approaches can be considered. In this study, the term time-series refers to a set of images acquired during a given period, and each acquisition was processed as an independent image without capturing properties of time-series data. New approaches that take time explicitly into account such as TWDTW (Time-Weighted Dynamic Time Warping) algorithm have proven to be very effective in solving complex classification problems. Although few studies of winter land-use classification have been performed, [28] for example showed the utility of deep learning for classifying vegetation quality (cover density) in winter using S-1 data, with preliminary results achieving a Kappa index of 0.98. However, these approaches require high-performance calculations and need to run hyperparameter searches, adjustments and tests, which need further investigations that will be done in future work. Furthermore, as mentioned, the results highlighted the value of specific radar parameters in the L- and C-Bands for discriminating winter land-use classes. These results highlight the utility of approaches that combine data to classify crops, as demonstrated in certain studies [17,18].

## 5. Conclusions

This study evaluated advantages of using S-1, RST-2 and/or ALS-2 time-series to determine the best SAR configuration to identify winter land-use classes accurately using the RF algorithm. Several SAR configurations were tested to discriminate land-use types during winter using the RF algorithm; to our knowledge, this is the first time such a study has been undertaken. Results show that the best SAR configuration was the dense dual-pol C-Band S-1 time-series, although RST-2 and ALS-2 time-series provided useful information about vegetation cover. Finally, our results demonstrated the limits of using one SAR dataset alone to identify winter land use accurately, the highest F1-score reaching only 0.70. Thus, future research could study the utility of combining SAR parameters or using new classification approaches based on deep learning to improve the accuracy of classifying land-use types in winter. Better understanding of SAR signal behaviors of agricultural practices and environmental conditions could also help to identify winter land use, which has important implications for developing sustainable agriculture.

**Author Contributions:** Conceptualization, J.D., E.P. and L.H.-M.; methodology, J.D., E.P. and L.H.-M.; software, E.P. and J.D.; validation, J.D.; resources, J.D.; data curation, J.D., writing—original draft preparation, J.D.; writing—review and editing, J.D., L.H.-M. and E.P.; visualization, J.D., E.P. and L.H.-M.; supervision, E.P. and L.H.-M.; project administration, E.P. and L.H.-M.; funding acquisition, E.P.

**Funding:** This research was funded by the French Ministry of Higher Education and Research, Ph.D. grant 2016

**Acknowledgments:** This study was supported by the Kalideos project, funded by the CNES (https://www.kalideos.fr/drupal/) and the Zone Atelier Armorique project. This study was also supported by public funds (Ministère de l'Éducation Nationale, de l'Enseignement Supérieur et de la Recherche, FEDER, Région Bretagne, Conseil Général du Finistère, Brest Métropole), VIGISAT project and the GIS Bretel ("Groupement Bretagne Télédétection"). We thank Napo N'Bohn and Marianne Balaresque for their help collecting field data.

**Conflicts of Interest:** The authors declare no conflict of interest.

## References

1. Fasona, M.J.; Omojola, A.S. Climate change, human security and communal clashes in Nigeria. In Proceedings of the International Workshop on Human Security and Climate Change, Oslo, Norway, 21–23 June 2005; pp. 21–23.

2. Corgne, S. Hiérarchisation des facteurs structurant les dynamiques pluriannuelles des sols nus hivernaux. Application au bassin versant du Yar (Bretagne). *Norois Environ. Aménage. Société* **2004**, *193*, 17–29. [CrossRef]
3. Denize, J.; Hubert-Moy, L.; Betbeder, J.; Corgne, S.; Baudry, J.; Pottier, E. Evaluation of using sentinel-1 and-2 time-series to identify winter land use in agricultural landscapes. *Remote Sens.* **2019**, *11*, 37. [CrossRef]

4. Duchemin, B.; Fieuzal, R.; Rivera, M.A.; Ezzahar, J.; Jarlan, L.; Rodriguez, J.C.; Hagolle, O.; Watts, C. Impact of sowing date on yield and water use efficiency of wheat analyzed through spatial modeling and FORMOSAT-2 images. *Remote Sens.* **2015**, *7*, 5951–5979. [CrossRef]

5. Veloso, A.; Mermoz, S.; Bouvet, A.; Le Toan, T.; Planells, M.; Dejoux, J.-F.; Ceschia, E. Understanding the temporal behavior of crops using Sentinel-1 and Sentinel-2-like data for agricultural applications. *Remote Sens. Environ.* **2017**, *199*, 415–426. [CrossRef]

6. Morena, L.C.; James, K.V.; Beck, J. An introduction to the RADARSAT-2 mission. *Can. J. Remote Sens.* **2004**, *30*, 221–234. [CrossRef]

7. Potin, P.; Bargellini, P.; Laur, H.; Rosich, B.; Schmuck, S. Sentinel-1 mission operations concept. In *Proceedings of the 2012 IEEE International Geoscience and Remote Sensing Symposium, Munich, Germany*, 22–27 July 2012; pp. 1745–1748.

8. Drusch, M.; Del Bello, U.; Carlier, S.; Colin, O.; Fernandez, V.; Gascon, F.; Hoersch, B.; Isola, C.; Laberinti, P.; Martimort, P. Sentinel-2: ESA's optical high-resolution mission for GMES operational services. *Remote Sens. Environ.* **2012**, *120*, 25–36. [CrossRef]

9. Huang, H.; Chen, Y.; Clinton, N.; Wang, J.; Wang, X.; Liu, C.; Gong, P.; Yang, J.; Bai, Y.; Zheng, Y. Mapping major land cover dynamics in Beijing using all Landsat images in Google Earth Engine. *Remote Sens. Environ.* **2017**, *202*, 166–176. [CrossRef]

10. Smith, L.C. Satellite remote sensing of river inundation area, stage, and discharge: a review. *Hydrol. Process.* **1997**, *11*, 1427–1439. [CrossRef]

11. Ulaby, F.T.; Moore, R.K.; Fung, A.K. Microwave Remote Sensing Active and Passive-Volume III: From Theory to Applications. 1986. Available online: https://ntrs.nasa.gov/search.jsp?R=19860041708 (accessed on 1 January 1986).

12. Hosseini, M.; McNairn, H.; Merzouki, A.; Pacheco, A. Estimation of Leaf Area Index (LAI) in corn and soybeans using multi-polarization C-and L-band radar data. *Remote Sens. Environ.* **2015**, *170*, 77–89. [CrossRef]

13. McNairn, H.; Brisco, B. The application of C-band polarimetric SAR for agriculture: a review. *Can. J. Remote Sens.* **2004**, *30*, 525–542. [CrossRef]

14. Mascolo, L. Polarimetric SAR for the Monitoring of Agricultural Crops. Ph.D. Thesis, Universita’degli Studi di Cagliari, Cagliari, Italy, 2015.

15. Hütt, C.; Koppe, W.; Miao, Y.; Bareth, G. Best accuracy land use/land cover (LULC) classification to derive crop types using multitemporal, multisensor, and multi-polarization SAR satellite images. *Remote Sens.* **2016**, *8*, 684. [CrossRef]

16. Haldar, D.; Rana, P.; Yadav, M.; Hooda, R.S.; Chakraborty, M. Time series analysis of co-polarization phase difference (PPD) for winter field crops using polarimetric C-band SAR data. *Int. J. Remote Sens.* **2016**, *37*, 3753–3770. [CrossRef]

17. Skriver, H. Crop classification by multitemporal C-and L-band single-and dual-polarization and fully polarimetric SAR. *IEEE Trans. Geosci. Remote Sens.* **2012**, *50*, 2138–2149. [CrossRef]

18. Skriver, H.; Svendsen, M.T.; Thomsen, A.G. Multitemporal C-and L-band polarimetric signatures of crops. *IEEE Trans. Geosci. Remote Sens.* **1999**, *37*, 2413–2429. [CrossRef]

19. McNairn, H.; Duguay, C.; Boisvert, J.; Huffman, E.; Brisco, B. Defining the sensitivity of multi-frequency and multi-polarized radar backscatter to post-harvest crop residue. *Can. J. Remote Sens.* **2001**, *27*, 247–263. [CrossRef]

20. Jiao, X.; McNairn, H.; Shang, J.; Liu, J. The sensitivity of multi-frequency (X, C and L-band) radar backscatter signatures to bio-physical variables (LAI) over corn and soybean fields. In *Proceedings of the ISPRS TC VII Symposium—100 Years ISPRS*, Vienna, Austria, 5–7 July 2010; pp. 317–325.

21. Nurtyawan, R.; Saepuloh, A.; Budiharto, A.; Wikantika, K. Modeling Surface Roughness to Estimate Surface Moisture Using Radarsat-2 Quad Polarimetric SAR Data. In *Proceedings of the Journal of Physics: Conference Series, Bandung, Indonesia*, 19–20 August 2016; Volume 739, p. 012105.

22. Smith, A.M.; Major, D.J. Radar backscatter and crop residues. *Can. J. Remote Sens.* **1996**, *22*, 243–247. [CrossRef]
23. Paris, J.F. Radar backscattering properties of corn and soybeans at frequencies of 1.6, 4.75, and 13.3 Ghz. *IEEE Trans. Geosci. Remote Sens.* **1983**, *3*, 392–400. [CrossRef]

24. Lee, J.-S.; Grunes, M.R.; Pottier, E. Quantitative comparison of classification capability: Fully polarimetric versus dual and single-polarization SAR. *IEEE Trans. Geosci. Remote Sens.* **2001**, *39*, 2343–2351.

25. Abdikan, S.; Sanli, F.B.; Ustuner, M.; Calò, F. Land cover mapping using sentinel-1 SAR data. *Int. Arch. Photogramm. Remote Sens. Spat. Inf. Sci.* **2016**, *41*, 757. [CrossRef]

26. Dimov, D.; Löw, F.; Ibrakhimov, M.; Stulina, G.; Conrad, C. SAR and optical time series for crop classification. In Proceedings of the 2017 IEEE International Geoscience and Remote Sensing Symposium (IGARSS), Fort Worth, TX, USA, 23–28 July 2017; pp. 811–814.

27. Bargiel, D. A new method for crop classification combining time series of radar images and crop phenology information. *Remote Sens. Environ.* **2017**, *198*, 369–383. [CrossRef]

28. Minh, D.H.T.; Ienco, D.; Gaetano, R.; Lalande, N.; Ndikumana, E.; Osman, F.; Maurel, P. Deep recurrent neural networks for winter vegetation quality mapping via multitemporal SAR Sentinel-1. *IEEE Geosci. Remote Sens. Lett.* **2018**, *15*, 464–468. [CrossRef]

29. ZA Armorique. Available online: https://osur.univ-rennes1.fr/za-armorique/ (accessed on 24 May 2019).

30. Nitrates—Water Pollution—Environment—European Commission. Available online: http://ec.europa.eu/environment/water/water-nitrates/index_en.html (accessed on 24 May 2019).

31. Open Access Hub. Available online: https://scihub.copernicus.eu/ (accessed on 10 September 2019).

32. Kalideos. Available online: https://www.kalideos.fr/drupal/fr (accessed on 24 May 2019).

33. Pottier, E.; Ferro-Famil, L.; Fitrzyk, M.; Desnos, Y.-L. PolSARpro-BIO: The new Scientific Toolbox for ESA & third party fully Polarimetric SAR Missions. In Proceedings of the 12th European Conference on Synthetic Aperture Radar, EUSAR 2018, Aachen, Germany, 4–7 June 2018; pp. 1–4.

34. Radarsat 1 Products Image Resolution Geomatics. Available online: https://www.scribd.com/document/73559622/Radarsat-1-Products (accessed on 20 November 2017).

35. Lee, J.-S.; Wen, J.-H.; Ainsworth, T.L.; Chen, K.-S.; Chen, A.J. Improved sigma filter for speckle filtering of SAR imagery. *IEEE Trans. Geosci. Remote Sens.* **2009**, *47*, 202–213.

36. Lee, J.-S.; Pottier, E. *Polarimetric Radar Imaging: From Basics to Applications*; CRC Press: Boca Raton, FL, USA, 2009.

37. Cloude, S.R.; Pottier, E. A review of target decomposition theorems in radar polarimetry. *IEEE Trans. Geosci. Remote Sens.* **1996**, *34*, 498–518. [CrossRef]

38. Freeman, A.; Durden, S.L. A three-component scattering model for polarimetric SAR data. *IEEE Trans. Geosci. Remote Sens.* **1998**, *36*, 963–973. [CrossRef]

39. Kim, Y.; Jackson, T.; Bindlish, R.; Lee, H.; Hong, S. Radar vegetation index for estimating the vegetation water content of rice and soybean. *IEEE Geosci. Remote Sens. Lett.* **2012**, *9*, 564–568.

40. McNairn, H.; Hochheim, K.; Rabe, N. Applying polarimetric radar imagery for mapping the productivity of wheat crops. *Can. J. Remote Sens.* **2004**, *30*, 517–524. [CrossRef]

41. Kostelich, E.J.; Schreiber, T. Noise reduction in chaotic time-series data: A survey of common methods. *Phys. Rev. E* **1993**, *48*, 1752. [CrossRef]

42. Breiman, L. Random Forests. *Mach. Learn.* **2001**, *45*, 5–32. [CrossRef]

43. Pal, M. Random forest classifier for remote sensing classification. *Int. J. Remote Sens.* **2005**, *26*, 217–222. [CrossRef]

44. Belgiu, M.; Csillik, O. Sentinel-2 cropland mapping using pixel-based and object-based time-weighted dynamic time warping analysis. *Remote Sens. Environ.* **2018**, *204*, 509–523. [CrossRef]

45. Liaw, A.; Wiener, M. Classification and regression by randomForest. *R News* **2002**, *2*, 18–22.

46. R: The R Project for Statistical Computing. Available online: https://www.r-project.org/ (accessed on 10 September 2019).

47. Lawrence, R.L.; Wood, S.D.; Sheley, R.L. Mapping invasive plants using hyperspectral imagery and Breiman Cutler classifications (RandomForest). *Remote Sens. Environ.* **2006**, *100*, 356–362. [CrossRef]

48. Belgiu, M.; Drăguț, L. Random forest in remote sensing: A review of applications and future directions. *ISPRS J. Photogramm. Remote Sens.* **2016**, *114*, 24–31. [CrossRef]

49. Audebert, N.; Le Saux, B.; Lefèvre, S. Beyond RGB: Very high resolution urban remote sensing with multimodal deep networks. *ISPRS J. Photogramm. Remote Sens.* **2018**, *140*, 20–32. [CrossRef]
50. Congalton, R.G. A review of assessing the accuracy of classifications of remotely sensed data. *Remote Sens. Environ.* **1991**, *37*, 35–46. [CrossRef]

51. Jiao, X.; Kovacs, J.M.; Shang, J.; McNairn, H.; Walters, D.; Ma, B.; Geng, X. Object-oriented crop mapping and monitoring using multi-temporal polarimetric RADARSAT-2 data. *ISPRS J. Photogramm. Remote Sens.* **2014**, *96*, 38–46. [CrossRef]

52. Yusoff, N.M.; Muharam, F.M.; Takeuchi, W.; Darmawan, S.; Abd Razak, M.H. Phenology and classification of abandoned agricultural land based on ALOS-1 and 2 PALSAR multi-temporal measurements. *Int. J. Digit. Earth* **2017**, *10*, 155–174. [CrossRef]

53. Immitzer, M.; Vuolo, F.; Atzberger, C. First experience with Sentinel-2 data for crop and tree species classifications in central Europe. *Remote Sens.* **2016**, *8*, 166. [CrossRef]

54. McNairn, H.; Shang, J.; Champagne, C.; Jiao, X. TerraSAR-X and RADARSAT-2 for crop classification and acreage estimation. In Proceedings of the 2009 IEEE International Geoscience and Remote Sensing Symposium, Cape Town, South Africa, 12–17 July 2009; Volume 2, pp. 2–898.

55. Liu, C.; Shang, J.; Vachon, P.W.; McNairn, H. Multiyear crop monitoring using polarimetric RADARSAT-2 data. *IEEE Trans. Geosci. Remote Sens.* **2012**, *51*, 2227–2240. [CrossRef]

56. Haldar, D.; Das, A.; Mohan, S.; Pal, O.; Hooda, R.S.; Chakraborty, M. Assessment of L-band SAR data at different polarization combinations for crop and other landuse classification. *Prog. Electromagn. Res.* **2012**, *36*, 303–321. [CrossRef]

57. Blaes, X.; Vanhalle, L.; Defourny, P. Efficiency of crop identification based on optical and SAR image time series. *Remote Sens. Environ.* **2005**, *96*, 352–365. [CrossRef]

58. Baghdadadi, N.; Boyer, N.; Todoroff, P.; El Hajj, M.; Bégué, A. Potential of SAR sensors TerraSAR-X, ASAR/ENVISAT and PALSAR/ALOS for monitoring sugarcane crops on Reunion Island. *Remote Sens. Environ.* **2009**, *113*, 1724–1738. [CrossRef]

59. Inoue, Y.; Kurosu, T.; Maeno, H.; Uratsuka, S.; Kozu, T.; Dabrowska-Zielinska, K.; Qi, J. Season-long daily measurements of multifrequency (Ka, Ku, X, C, and L) and full-polarization backscatter signatures over paddy rice field and their relationship with biological variables. *Remote Sens. Environ.* **2002**, *81*, 194–204. [CrossRef]

60. Loosvelt, L.; Peters, J.; Skriver, H.; De Baets, B.; Verhoest, N.E. Impact of reducing polarimetric SAR input on the uncertainty of crop classifications based on the random forests algorithm. *IEEE Trans. Geosci. Remote Sens.* **2012**, *50*, 4185–4200. [CrossRef]

61. Dufour, S.; Bernez, I.; Betbeder, J.; Corgne, S.; Hubert-Moy, L.; Nabucet, J.; Rapinel, S.; Sawtschuk, J.; Trollé, C. Monitoring restored riparian vegetation: How can recent developments in remote sensing sciences help? *Knowl. Manag. Aquat. Ecosyst.* **2013**, *410*, 10. [CrossRef]

62. Hong, G.; Wang, S.; Li, J.; Huang, J. Fully polarimetric synthetic aperture radar (SAR) processing for crop type identification. *Photogramm. Eng. Remote Sens.* **2015**, *81*, 109–117. [CrossRef]

63. Zhang, Y.; Wu, L. Crop classification by forward neural network with adaptive chaotic particle swarm optimization. *Sensors* **2011**, *11*, 4721–4743. [CrossRef]

64. Adams, J.R.; Berg, A.A.; McNairn, H.; Merzouki, A. Sensitivity of C-band SAR polarimetric variables to unvegetated agricultural fields. *Can. J. Remote Sens.* **2013**, *39*, 1–16. [CrossRef]

<figure><img src="image_14.png" /></figure>