Next Article in Journal
Research Trends in the Remote Sensing of Phytoplankton Blooms: Results from Bibliometrics
Previous Article in Journal
Attention Enhanced U-Net for Building Extraction from Farmland Based on Google and WorldView-2 Remote Sensing Images
Previous Article in Special Issue
Sentinel-1 Time Series for Crop Identification in the Framework of the Future CAP Monitoring
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Unsupervised Classification of Crop Growth Stages with Scattering Parameters from Dual-Pol Sentinel-1 SAR Data

1
Microwave Remote Sensing Lab, Centre of Studies in Resources Engineering, Indian Institute of Technology Bombay, Mumbai 400076, India
2
Centre Eau Terre Environnement, INRS, 490 Couronne St, Quebec City, QC G1K 9A9, Canada
3
Ottawa Research and Development Centre, Agriculture and Agri-Food Canada, 1341 Baseline Road, Ottawa, ON K1A 0C5, Canada
*
Author to whom correspondence should be addressed.
Remote Sens. 2021, 13(21), 4412; https://0-doi-org.brum.beds.ac.uk/10.3390/rs13214412
Submission received: 18 September 2021 / Revised: 21 October 2021 / Accepted: 25 October 2021 / Published: 2 November 2021

Abstract

:
Global crop mapping and monitoring requires high-resolution spatio-temporal information. In this regard, dual polarimetric Synthetic Aperture Radar (SAR) sensors provide high temporal and high spatial resolutions with large swath width. Generally, crop phenological development studies utilized SAR backscatter intensity-based descriptors. However, these descriptors are derived either from the covariance matrix elements or from the eigendecomposition. Therefore, this approach fails to utilize the complete polarization information of the scattered wave. In this study, we propose a target characterization parameter, θ xP that utilizes the 2D Barakat degree of polarization and the elements of the covariance matrix. We also propose an unsupervised clustering scheme using θ xP and the scattering entropy, H xP . We utilize time-series Sentinel-1 data of canola and wheat fields over a Canadian test site to show the sensitivity of θ xP to the development of crop morphology at different phenological stages. During the initial growth stages, θ xP values are low due to the low vegetation density. In contrast, at advanced phenological stages, we observe decreased values of θ xP due to the appearance of complex canopy structure. Similarly, the effectiveness of the unsupervised H xP / θ xP clustering plane is also evident from the temporal clustering plots. This innovative clustering framework is beneficial for the operational use of Sentinel-1 SAR data for agricultural applications.

Graphical Abstract

1. Introduction

Identification and monitoring of crop phenological stages are essential factors in agriculture for estimating crop production. In this context, Synthetic Aperture Radar (SAR) data have been extensively used, especially for crop classification [1,2], yield estimation, and biophysical parameter retrieval [3,4,5], formulation of vegetation indices [6,7] and descriptors [8]. Moreover, SAR sensors have high spatial resolution and all-weather imaging capability. The scattered SAR signal is significantly affected by the geometry of crop canopy, underlying soil roughness, and dielectric properties of both crop and soil. In addition to the acquisition wavelength, polarization and angle of incidence have a significant influence on the received information [9,10]. Over agricultural fields, crop biophysical parameters impact SAR backscatter, particularly at advanced vegetative stages. Soil roughness and moisture remain the governing factors for SAR backscatter response during initial growth periods [11].
In particular, Sentinel-1 dual-pol images provide opportunities for crop classification and biophysical retrieval-related applications. Sentinel-1 also provides large swaths and lower data volumes due to reduced polarimetric dimensionality [12]. In the literature, dual-pol intensity parameters are widely used for crop type identification [1,13,14,15,16,17,18,19]. The existing studies also extend to retrieval of crop biophysical parameters [4,20,21] and crop phenology monitoring [22,23].
For dual polarimetric HH-HV or VV-VH data, Cloude et al. [24] proposed an eigen-based decomposition technique. In this technique, the 2 × 2 covariance matrix is decomposed into two orthogonal eigenvectors that are further utilized to derive a target characterization parameter. This average scattering angle α ¯ is obtained from the two orthogonal polarization states weighted by their respective pseudo probabilities obtained from the eigenvalues. The target scattering entropy H xP is obtained from the pseudo probabilities. This target characterization parameter and H xP are further utilized to propose an unsupervised clustering scheme. Similarly, Ainsworth et al. [25] proposed another characterization parameter, θ , utilizing an eigen-based technique. The parameter is related to the cross, σ XY and co-pol, σ XX ratios. In addition, θ and H xP are used to propose a clustering technique for scattering target identification. In their study, the clustering plane consists of eight zones that represent different scattering characteristics.
Sugimoto et al. [26] provided a comparison between polarimetric parameters obtained from dual and full polarimetric SAR data. They reported a high 2D correlation between the parameters derived from the full-pol H / α decomposition and the dual-pol (HH+VV) H / α decomposition. For the full polarimetric data, the analysis was performed using model-based three [27], and four [28] component decompositions. Several other descriptors were introduced for diverse applications: clustering [29], maritime applications [30,31], etc.
To monitor the growth of vegetation, different vegetation indices, such as the Radar Vegetation Index (RVI) [32], Dual-Pol SAR Vegetation Index (DPSVI) [33], and Dual-pol Radar Vegetation Index (DpRVI) [34], can be derived. In general, these descriptors capture the randomness from the complex crop canopy to describe their phenological stages. Apart from these, machine learning models have shown promising results for crop biophysical parameter estimation [22,35,36]. It was observed from these studies that the cross-pol ratio is the most important parameter for vegetation studies. Nasrallah et al. [37] utilized the Gaussian function on the time series SAR data to find the date of significant phenology stages for wheat. Wali et al. [38] analyzed the sensitivity of the temporal backscatter for rice crops using the line regression technique. Similarly, the interferometric phase for identifying wheat phenological stages is observed in a study by Schlund et al. [39].
Usually, polarimetric parameters have been suitably attributed to the physical properties of the crop canopy [40,41,42], and have therefore helped monitor crop phenology. Recently, Dey et al. [43,44] proposed a model-free target characterization parameter using full, compact and dual co-pol data [22]. This parameter utilizes the Barakat degree of polarization [45] and elements of the coherency matrix to characterize diverse target types. In their study, it was shown that this parameter enhanced the target characterization capability as compared to the α ¯ parameter.
Motivated by this work, we derive a target characterization parameter for dualpol SAR data in this study. This parameter jointly utilizes the 2D Barakat degree of polarization and the elements of the covariance matrix. Therefore, this parameter, θ xP , captures the additional scattering information which the other existing descriptors might miss. Alongside this, we also propose an unsupervised clustering technique using θ xP and entropy, H xP parameters. We utilize θ xP and the clustering scheme to analyze the temporal variation of canola and wheat over a Canadian test site derived from Sentinel-1 C-band SAR data.

2. Methodology

2.1. Study Area and Dataset

The test site in Canada is located in southern Manitoba (49°40′32″N, 97°59′57″S). This is one of the Canadian Joint Experiment for Crop Assessment and Monitoring (JECAM) test sites. Nominal field sizes range from 20–30 hectares to 50–60 hectares, with wheat, soybean, canola, and corn as major annual crops. The land under permanent grassland is < 5 . The in-situ measurements were conducted during Soil Moisture Active Passive Validation Experiment 2016 Manitoba (SMAPVEX16-MB) campaign. Figure 1 shows the location of the test site and the distribution of canola and wheat fields.
The in-situ measurement campaign was carried out from June to July 2016 over 50 different fields. The measurement period mainly consisted of two temporal windows: 8 June to 22 June and 8 July to 22 July 2016. Within this temporal extent, most crops advanced from early to advanced phenological stages achieving peak biomass. The nominal size of each field is approximately 800   m × 800   m . Vegetation sampling, which includes biomass, Plant Area Index (PAI), and plant height, was performed in three points for each field. Crop biomass was collected via destructive sampling. A 0.5 m × 0.5 m square was placed over the canopy for canola and wheat. All above-ground biomass was collected by cutting all vegetation at the soil level within the square. The PAI was measured with hemispherical digital photos. In this technique, a camera with a fisheye lens captures photos of the crop canopy with the camera positioned at least 50 cm above or below the canopy. Details of the in-situ biophysical parameters (viz. PAI, Dry biomass and VWC) of Canola and Wheat are provided in Appendix A.3. One can find further details on the sampling schema in the SMPAVEX16-MB field report [46].

2.2. SAR Data Pre-Processing

We acquired the Terrain Observation with Progressive Scans SAR (TOPSAR) mode Level-1 Sentinel-1 Interferometric Wide (IW) swath Single Look Complex (SLC) data over the test site. To this SLC configuration, the swath length is ≈250 km and the spatial resolution is 5   m × 20   m . The whole IW swath is divided into three sub-swaths (IW1, IW2, and IW3) with nine bursts in the azimuth direction. Hence, prior to application, these images are pre-processed with standard correction steps [47].
The Sentinel-1 images are ingested in the SNAP 8.0 platform. Following this, the sub-swaths and bursts are selected based on the test site location, and the state vectors are updated with the Sentinel-1 precise orbit file. The images are then radiometrically calibrated. As we are interested in generating the complex covariance C 2 matrix information, and hence, we need to extract both amplitude and phase values. Due to this reason, the calibration output is set to complex output values.
All Sentinel-1 images are then back-geocoded using the Sentinel-1 Back Geocoding operator. This operation essentially co-registers the temporal images with sub-pixel accuracy. We utilize the SRTM 1Sec Grid as the Digital Elevation Model (DEM). Subsequently, the Sentinel-1 TOPS deburst and merge operations are performed to produce a single SLC image. The stack of images is then clipped within the in-situ measurement location to reduce the data volume and increase the computation speed.
These image subsets are multi-looked with a factor of 4 in range and 1 in azimuth directions to generate a square pixel. The final pixel resolution is ≈ 15 m . Finally, these multi-looked images are used to generate covariance elements images. A 5 × 5 boxcar filter is applied to further reduce the speckle information within the images. Next, the baseline information is deleted from the metadata, and the covariance images are exported in the PolSARpro format. Here, please note that the covariance elements have ensembled information that helps in applying the second-order statistics.
Then, we compute the target characterization parameter, θ xP using the PolSARtools software in the QGIS platform [48]. Subsequently, we geocode the images with UTM projected coordinate system. We further analyzed the in-situ measurement locations and extracted target characterization parameters from the geocoded products. From several Sentinel-1 images acquired during the campaign, four dual-polarization (VV and VH) C-band Sentinel-1 Single Look Complex (SLC) data were selected for use in this study (Table 1). The selection of Sentinel-1 datasets was based on acquisition dates that were near coincident with in-situ measurement periods.

2.3. Target Characterization Parameter

The scattering vector for a dual-polarized SAR data are represented as, k d = [ S X X S X Y ] T , where X and Y correspond to either horizontal (H) or vertical (V) polarization states. Using these scattering vectors, we define the 2 × 2 covariance matrix as
C 2 = k d · k d * T = | S X X | 2 S X X S X Y * S X Y S X X * | S X Y | 2
where · denotes ensemble average, and T denotes vector transpose. We define the elements of the 2 × 2 matrix as, C 11 = | S X X | 2 , C 22 = | S X Y | 2 and C 12 = S X X S X Y * .
Similar to the conventional degree of polarization, the 2D Barakat degree of polarization [49] given in Equation (2) ( 0 m xP 1 ) also characterizes the state of polarization (or purity) of an EM wave. For a completely polarized EM wave, m xP = 1 and for a completely unpolarized EM wave, m xP = 0 . In between these two extreme cases, the EM wave is said to be partially polarized, 0 < m xP < 1 :
m xP = 1 4 | C 2 | Tr ( C 2 ) 2 ,
where | · | is the determinant of a matrix and Tr is the trace of a matrix. Let us consider two auxiliary quantities defined as
tan η 1 = C 11 m xP Span and tan η 2 = C 22 m xP Span ,
where η 1 and η 2 are two free variables, and Span = C 11 + C 22 . The quantity C 11 m xP Span indicates the co-polarized scattering power with respect to the total polarized power, while C 22 m xP Span indicates the cross-polarized scattering power with respect to the total polarized power. Furthermore, using a simple relationship, we obtain:
tan θ xP = tan η 1 η 2 = m xP Span C 11 C 22 C 11 C 22 + m xP 2 Span 2 .
where θ xP [ 45 , 45 ] .
The eigen-decomposition of C 2 can be expressed as
C 2 = U 2 Σ U 2 1 ,
where Σ is a 2 × 2 diagonal matrix with non-negative real elements, λ 1 λ 2 0 , which are the eigenvalues of C 2 , and U 2 is 2 × 2 unitary matrix, where u i ’s are the unit orthogonal eigenvectors. We define the pseudo probabilities, p i , in terms of the eigenvalues as
p i = λ i k = 1 2 λ k ,
which we then use to define the scattering entropy as
H xP = k = 1 2 p k log 2 p k .
The degree of polarization m xP characterizes the degree of coherence of partially polarized waves. In contrast, the entropy H xP describes the degree of statistical disorder associated with partially polarized waves. The parameter θ xP describes the scattering characteristics of a target. Therefore, in this work, these parameters provide information about the growth stage or phenology of crops.

2.4. H xP / θ xP Bound

The feasible regions in the H xP / θ xP plot is represented by the bounding curves, Curve I ( C 2 ( I ) ) and Curve II ( C 2 ( I I ) ) as shown in Figure 2 and is given as
C 2 ( I ) = 1 1 + k 1 0 0 k , θ xP max = tan 1 k 2 2 k + 1 k 2 k + 1
C 2 ( I I ) = 1 k + 1 k 0 0 1 , θ xP min = tan 1 k 2 2 k + 1 k 2 k + 1
where k [ 0 , 1 ] denotes the scattering amplitude ratio, and therefore, θ xP [ 45 , 45 ] . One can observe from Equation (4) that when:
  • m xP = 0 (i.e., when there exists no polarization structure in the scattered EM wave), then θ xP = 0 characterize random scattering from targets.
  • m xP = 1 , and θ xP = 45 , characterize coherent scattering from deterministic targets (i.e., trihedral or dihedral).
  • m xP = 1 , and θ xP = 45 , characterize cross-polarized scattering from complex targets.
Therefore, θ xP [ 45 , 45 ] characterizes diverse scattering-type information in between these extreme cases, and hence it is suitable to characterize scattering-type information from various targets.

2.5. Example of Variation of α ^ and θ xP

In dual polarimetry, Cloude proposed a target characterization parameter, α ¯ [ 0 , 90 ] using the eigendecomposition approach [24]. In order to compare this with the proposed θ xP [ 45 , 45 ] , we have rescaled α ¯ to α ^ = 45 ° α ¯ , where both α ^ and θ xP vary from −45 ° to 45 ° .
To show the efficacy of θ xP over α ^ , we have considered three different land cover types, i.e., urban, oriented urban, and vegetation from the VV-VH data extracted from a full-pol C-band RADARSAT-2 data. The variations of θ xP and α ^ over these targets are shown in Figure 3. One can note from Figure 3 that the median value of α ^ over orthogonal urban area is approximately 35°. In contrast, the median value of θ xP is approximately 42°. The degree of polarization over the area is approximately 0.91 . This high degree of polarization essentially confirms the nearly coherent scattering characteristics of the target. In this regard, θ xP better characterizes the target to be coherent than α ^ .
However, over oriented urban area, we observe a large variation of both θ xP and α ^ . This large variation might be due to the diversity in the polarization randomness due to the orientation of buildings to the radar line of sight and orientation in the azimuth direction. Over vegetation surfaces, we observe a partially polarized scattered wave that consists of a major coherent part and a significant incoherent part. The coherent part might be due to the direct scattering from the nearly smooth soil surface, while the incoherent part is due to the scattering from the volume media. In this case, θ xP also correctly characterizes the coherent contribution of the target during scattering as compared to α ^ . A α ^ θ xP plot for RADARSAT-2 data over San-Francisco in VV-VH mode is shown in Appendix A.1. Therefore, from these analyses, we can observe the enhancement of target scattering characterization using θ xP . Thus, θ xP can be further utilized to characterize different phenological stages of diverse crop types.

2.6. Unsupervised Clustering Zones over Vegetative Surface

In this study, we are interested in observing the changes in the scattering mechanisms over different phenological stages of crops. One should note that, for natural targets, C 11 C 22 almost always, and hence, θ xP [ 0 , 45 ] . Let us now characterize m xP and θ xP for a few particular scattering scenarios:
  • For a pure diffused target, m xP = 0 , implies, θ xP = 0 .
  • For pure or point scatterer, m xP = 1 and η 2 = 0 , implies θ xP = 45 .
  • Infeasible scattering: C 11 = C 22 and m xP = 1 .
Therefore, in terms of vegetation development, for a highly random vegetative structure, m xP 0 and θ xP 0 . However, a slightly rough soil surface acts like a Bragg scatterer, and hence m xP 1 and θ xP 45 .
The variations of θ xP and m xP over bare field and vegetative field types are shown in Figure 4. It can be noted that θ xP and m xP are highly sensitive with the canopy structure. Interestingly, we observed that the sensitivity is much higher for the highly erectophile crop structures such as canola, corn, and wheat. In contrast, the sensitivity is low for the soybean crop. This high sensitivity of θ xP might be due to the volume decorrelation with an increased canopy structure.
In this study, we utilize the θ xP and H xP parameters to propose a new unsupervised clustering scheme to describe the target scattering behavior. The overall clustering scheme consists of 12 zones: Z1, Z2, Z3, Z4, Z5, Z6, Z7, Z8, Z9, Z10, Z11, and Z12. These divisions of the H xP / θ xP plane are based on certain scattering symmetry assumptions. In particular, the scattered Stokes vector, S s is a function of the target property, represented in terms of the Kennaugh matrix, K and the transmitted Stokes vector, S t . Moreover, the two Stokes vectors are related as S s = K S t . Furthermore, one can note that S s is also a function of the scattering order (n) [50]. Thus, the definition of scattering entropy, H [51,52] can be expressed in terms of n. Following this, we can observe that H increases with increasing n. For example, H = 0 for n = 0 ; H 0.3 for n = 1 ; H 0.5 for n = 2 , and, H 0.7 for n = 3 . However, for n > 3 , change in H is insignificant. Therefore, following this trend in the variation of the entropy, we divided the radial axis H xP into four sub-sections by considering the values: 0.3, 0.5, and 0.7. An illustration of the H xP / θ xP clustering plane with the zones is shown in Figure 5.

3. Results and Discussion

In this section, we analyze the temporal dynamics of crops using the proposed dualpolarimetric descriptors. Furthermore, we utilize the proposed clustering framework to assess the phenological stages of the two crops, i.e., canola and wheat, from the C-band Sentinel-1 dual-pol SLC SAR data.

3.1. Canola

The spatio-temporal changes of θ xP are shown in Figure 6. The temporal variations of θ xP and m xP are shown in Figure 7 and Figure 10. For temporal analysis, we considered 48 sampling points from three fields (Field No. 206, 208 and 224). We plot the variations of the proposed SAR descriptors with respect to the phenological changes of canola at different dates.
According to the in-situ information, canola was sown at the end of May. Hence, during June, the phenological stage of canola was bounded within the early to mid vegetative stage. In this regard, it is worth mentioning that canola is a broad-leaf crop with comparable leaf size to the wavelength of C-band SAR (≈5.6 cm). Therefore, due to the formation of dense rosette near the soil surface, the leaf structure considerably affects the SAR backscatter values. Thus, due to the high dynamic variations of crop morphology, canola is of particular interest in SAR research.
From Figure 7, we can observe that, on 13 June, the median value of m xP 0.78 . This high value of m xP typically suggests high polarized scattered components from canola fields. This trend might be due to the low vegetation density in the fields. From the in-situ information, we observe that the canola was in the stem elongation stage during this time. Consequently, low biomass and low Plant Area Index (PAI) are evident at this period. Moreover, the vertical crop structure slightly attenuated the transmitted vertical (V) polarization. Hence, we observe a median value of θ xP 40 ° . Apart from this, the low depolarization in the scattered wave might be due to the scattering from leaves and canopy structures combined.
Furthermore, Figure 8 provides information about the different morphological characteristics of canola at this particular period. We observe that the points from these canola fields are mainly distributed in two separate zones: Z4 and Z7 in Figure 8a. This information typically infers that, during 13 June, the crop density in some canola fields was low. As a result, low entropy pure scattering was evident from those fields (Z4). In contrast, other fields had comparatively dense foliage structures, which increased scattering randomness increased (Z7). From in-situ measurement, we observe that fields 206 and 208 had higher PAI and biomass values than field 224. In field 224, most canola crops were at the leaf development stage. In addition, we observe a point in the Z1 region and points in the Z10 zone, indicating more advanced phenological stages. These variations from Z1 to Z10 reflect the differences in sowing dates for canola, and, consequently, differences in crop emergence and development among fields early in the season.
On 7 July, canola crops reached their flowering stage. During this time, we can observe a sudden drop in the θ xP and m xP values. The median value of m xP 0.44 . This sudden decrease in the scattered polarization structure is due to the development of flowers, stems, and branches. During this time, the high attenuation of the V-polarized transmitted wave also lowered the θ xP value. Furthermore, the generation of high cross-polarization due to the complex crop canopy structure also aided in lowering θ xP values. During this time, the median value of θ xP is approximately 20°.
Similar changes in the clustering zones are evident from Figure 8b. We observe a shift of clustering zone from Z4 and Z7 to Z11. The formation of the dense cluster in Z11 indicates high randomness in the scattered wave due to the complex flower and branch structures. Moreover, during the flowering and early pod development stages, the leaf density drastically reduces. Canola crops drop their leaves during this period of pod and seed development. However, this dense cluster in Z11 also suggests that the number of flowers per plant was high, and the health of the crop was good. On the other hand, we observe some points in the comparatively low entropy region, Z10.
On 19 July, small changes in θ xP and m xP are observed as compared to 7 July. This effect might be due to the existence of a similar complex canopy structure due to which the depolarization is evident in the scattered wave. During this early to the mid-pod development stage, ramified stems and seeds appear. In the early stage, seeds remain watery and translucent, while seeds reach maximum dry weight at the advanced stage. Therefore, depending upon this particular phenomenon, changes in the scattering mechanisms might be observed. However, during the satellite acquisition, most fields reached the early to mid-stage of pod development. As a result, the variance in θ xP is low during this time. The median value of θ xP is ≈18°, while the median of m xP is approximately 0.4 .
We also observe the accumulation of H xP / θ xP cluster in Z11 region in Figure 8c. During this time, the complex vegetation canopy structure increased the randomness in the scattered wave. Consequently, we observe a marginal increase in H xP compared to 7 July. However, as stated earlier, most of the fields reached the early to mid-stage of pod development, and thus we observe the majority of the point cloud in the Z11 region.
Subsequently, on 24 August, we observe an increase in θ xP values. Similarly, the increase in m xP suggests that the polarized component in the scattered wave has increased during this period. This increase in the values of θ xP and m xP might be due to the changes in the canopy moisture content. However, during this period, high variance in the data is evident. During this period, the median value of θ xP is approximately 34°, while the median value of m xP is approximately 0.58 . In Figure 8d, we observe three distinct zones; Z7, Z10, and Z11. This phenomenon is due to the harvest of canola crops. The post-harvest residue might have generated marginal depolarization in the scattered wave. The percentages of pixels at each date are shown in Table 2.

3.2. Wheat

In this section, we analyse the temporal variation of θ xP and m xP over wheat fields. Wheat is an edible cereal grass with an erectophile canopy structure. Over the test site, sowing of wheat was performed during the beginning of May. This study has considered three wheat fields (Field No. 233, 220, 62) for the temporal evaluation. The temporal patterns of θ xP and m xP are shown in Figure 10. Alongside this, we represent the dynamic behaviour of unsupervised clusters in Figure 9.
On 13 June, wheat crops reached the tillering stage. This stage starts with the appearance of tillers and ends with mature leaf sheaths. During this time, the crop density remains low. As a result, the scattered wave remains significantly polarized. It can be observed from Figure 10 that the median value of m xP during this period is 0.78 . This essentially indicates that approximately 78% of the scattered wave is completely polarized, while approximately 22% is depolarized. This marginal depolarization effect might be due to the canopy interaction. Alongside this, the median value of θ xP is ≈39°, indicating that the target is a closely coherent scatterer. However, the high standard deviations in both θ xP and m xP represent the variations in the uneven phenological advancement within the fields due to variations in seeding dates.
The clustering zones on 13 June are shown in Figure 9a. Like canola, we can observe that the clusters are distributed in four distinct zones; Z4, Z7, Z10, and Z11. The existence of these distinct zones is due to the differences in the morphological characteristics of wheat crops. From in-situ measurements, it is observed that the plant density of Field No. 220 was low (≈ 100   m 2 ). In contrast, the plant density over Field Nos. 233 and 62 was ≈ 125   m 2 and ≈ 190   m 2 , respectively. Similarly, high PAI and biomass are evident over these two fields. As a result, a high attenuation of the transmitted wave is observed over Field No. 233 and 62. Hence, we observe clouds of clusters at Z10 and Z11. On the other hand, the attenuation of the scattered wave is comparatively lower over Field No. 220. Hence, cluster formation at Z4 and Z7 is observed in Figure 9a.
On 7 July, we observe a decrease in both θ xP and m xP values due to the appearance of more branches and leaves. The median value of θ xP is around 27°, while the median value of m xP is around 0.49 . At this stage, wheat has advanced to the early flowering stage. At this stage, flowers appear on the upper portion of the canopy layer. As a result, the dense wheat structure and the flowers increase the multiple interactions of the EM wave. Hence, a high proportion of depolarized components in the scattered wave is evident during this period. Similar changes in the clusters are evident from Figure 9b. At this period, the clusters shift to Z11 from Z4 and Z7 due to the increased scattering entropy due to the randomness in the scattered wave by the targets. However, a cloud of clusters is seen at Z10, which might be due to some wheat crops having late phenological development.
On 19 July, wheat crops advanced to the early dough stage. At this phenological period, wheat grains appear, and it remains milky. Therefore, the dielectric property of grains affects the scattering mechanism. Moreover, the denser wheat canopy increases multiple scattering. As a result, the values of m xP reduce further as compared to 7 July. Similarly, the coherent component in the scattered EM wave also reduces. The median value of m xP is 0.37 . The median value of θ xP is 20 ° . The high variation in the values of θ xP and m xP might be due to the randomly oriented wheat stems and heads. Similar scattering characteristics from the wheat during the dough stage are reported by Wu et al. [53]. In particular, during this period, the scattering phenomenon occurs primarily from two parts of the wheat canopy, i.e., thick upper canopy layer and relatively sparse lower canopy layer.
Further from Figure 9c, we can observe the formation of clusters in the Z11 region. However, we can also observe some points in Z12 due to additional complex scattering mechanisms from the wheat canopy subjected to partial lodging or related canopy anomalies. It is important to note that we observe much spread of points in the Z11 zone, unlike other phenological stages of canola and wheat. This suggests that, during this particular phenology period, the appearance of grains and their orientations might be responsible for breaking the reflection asymmetry condition during scattering from the wheat canopy. Moreover, as stated earlier, the complex canopy geometry, the spatial distribution of grains, and their random orientations increased the overall randomness in the scattered EM wave. As a result, for the majority of the points, we observe high entropy values.
During 24 August, the harvest of wheat started. At this stage, the interaction of EM waves mainly takes place with the soil surface and post-harvest residue. Consequently, a reverse trend is observed in θ xP and m xP values. From Figure 10, we can observe that the median value of m xP is approximately 0.58 , while the median value of θ xP is approximately 34. These values indicate the coherent component increase in the scattered EM wave from the wheat canopy. Additionally, Figure 9d shows three different clustering zones, Z7, Z10, and Z11, during this period. The percentages of temporal variations of clusters at different zones are shown in Table 3. In addition, the temporal variations of θ xP and α ^ are shown in Appendix A.2. One can observe that the temporal signature of these parameters differs for wheat and canola. Therefore, these descriptors can be efficiently utilized to distinguish wheat and canola fields using time series data. Furthermore, it might also be possible to distinguish these two crop fields from other land cover targets. The temporal variation of the clusters over the cultivation area is shown in Figure 11.

4. Conclusions

In this study, we propose a new polarimetric target characterization parameter, θ xP [ 45 , 45 ] for dual-pol Sentinel-1 Synthetic Aperture Radar (SAR) data. However, over natural surfaces, such as vegetation, θ xP [ 0 , 45 ] , where θ xP = 0 ° denotes totally incoherent targets and θ xP = 45 ° denotes totally coherent targets. Utilizing the scattering entropy parameter, H xP along with θ xP , we proposed a novel unsupervised classification scheme for dual-pol SAR data. We further utilized θ xP and the classification scheme to analyze the phenological development of canola and wheat over a Canadian test site. The results are promising and exhibit high sensitivity with the morphological changes of the crops at each phenological stage.
We observe a broad dynamic range of θ xP from the temporal analysis of both canola and wheat starting from their early vegetative stage to maturity and harvest stage. For both crops, the variation of θ xP within the desired phenology window is 40 ° to 15 ° . High values of θ xP during the early stages are due to the low crop density within the fields. In contrast, the low θ xP values during the flowering stage are due to the complex canopy structure and additional appearance of flowers and grains. Moreover, we observe a trend reversal for both canola and wheat on 24 August. This trend reversal is due to the harvest of both crops during this period.
Within the scope of this study, we have characterized diverse crop phenological stages in terms of the physical scattering of the electromagnetic wave. The unsupervised classification scheme comprises twelve distinct zones, which represent these different physical scattering mechanisms. Thus, these clustering zones are beneficial for identifying the morphological status of the crop when a priori information is not available. Alongside this, the clustering scheme can also capture the difference in the growth stages simultaneously. Moreover, during the initial growth period, most of the points lay in the low entropy zones, while those points shifted towards high entropy zones at advanced phenological stages. This phenomenon is due to the changes in the canopy complexity with the advancement of the crop phenological stages. We already know that the dual-pol configuration does not have either the VV or HH polarization information. Hence, by utilizing θ xP , we are unable to uniquely characterize scattering-type information, i.e., whether it is an even-bounce or odd-bounce mechanism. In addition, as stated earlier, if the SAR signal gets highly saturated due to the morphological characteristics of targets, θ xP might provide a similar saturation effect. This appearance is evident during the end of flowering to maturity stages of crops.
We can extend this study to several other crop types using different dual polarimetric SAR sensor configurations, such as HH-HV. The proposed parameter and the clustering framework will be valuable to analyze data from the upcoming dual-pol NASA-ISRO Synthetic Aperture Radar Mission (NISAR) and Copernicus Sentinel SAR missions for agricultural and environmental studies. The open-source code is available at: https://github.com/Subho07/dual_cross_pol_theta_clustering.

Author Contributions

Conceptualization, S.D., N.B., A.B. and S.H.; methodology, S.D. and A.B.; programming, S.D. and N.B.; validation, A.B., S.H. and H.M.; format analysis, all; writing, all; supervision, S.H., A.B. and H.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The Sentinel-1 SLC data were accessed from https://search.asf.alaska.edu/#/, accessed on 16 July 2021.
DateData ID
13 JuneS1A_IW_SLC__1SDV_20160613T001529_20160613T001556_011685_011E64_4083
7 JulyS1A_IW_SLC__1SDV_20160707T001530_20160707T001557_012035_01298D_9FD9
19 JulyS1A_IW_SLC__1SDV_20160719T001540_20160719T001604_012210_012F46_1DED
24 AugustS1A_IW_SLC__1SDV_20160824T001533_20160824T001600_012735_0140AB_4BC9

Acknowledgments

The authors would like to thank the European Space Agency (ESA) for providing Sentinel-1 images over the test site and the members of the SMAPVEX16 science team for providing ground truth information. The authors also acknowledge the GEO-AWS Earth Observation Cloud Credits Program, which supported the computation with Sentinel-1 on AWS cloud platform through the project: AWS4AgriSAR-Crop inventory mapping from SAR data on a cloud computing platform and formed the testbed for processing pipelines. N.B. and S.D. would like to acknowledge the support from the Ministry of Education (formerly the Ministry of Human Resource and Development-MHRD), Govt. of India, towards their doctoral research work.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Appendix A.1. Relationship with Cloude’s α i ’s

In this section, we derive the relationship of θ xP with Cloude α i ’s [54]. In this regard, please note that
C 11 Span = i = 1 2 p i cos 2 α i ,
where ( α i ) i = 1 , 2 are the individual scattering-type parameters obtained from Cloude target scattering vector [24]. Let us define
tan η 1 = C 11 m xP Span = i = 1 2 p i cos 2 α i 1 i = 1 2 2 p i ,
and
tan η 2 = C 22 m xP Span = 1 i = 1 2 p i cos 2 α i 1 i = 1 2 2 p i .
Therefore, using simple trigonometric formulation, we get
tan ( η 1 η 2 ) = 2 i = 1 2 p i cos 2 α i 1 1 i = 1 2 2 p i 1 + i = 1 2 p i cos 2 α i ( 1 i = 1 2 p i cos 2 α i ) 1 i = 1 2 2 p i ,
further simplifying the expression of θ xP given in terms of ( p i ) i = 1 , 2 and ( α i ) i = 1 , 2 as
tan θ xP = Num Den , in which Num = 2 i = 1 2 p i cos 2 α i 1 1 i = 1 2 2 p i , and Den = 1 i = 1 2 2 p i + i = 1 2 p i cos 2 α i 1 i = 1 2 p i cos 2 α i .
Therefore, one can note that θ xP is a function of i = 1 2 p i cos 2 α i , and the Barakat degree of polarization, 1 i = 1 2 2 p i . We have shown the α ^ θ xP scatter plot for VV-VH polarization modes in Figure A1. Here, α ^ = 45 ° α ¯
Figure A1. α ^ θ xP plot for RADARSAT-2 data over San Francisco in VV-VH mode.
Figure A1. α ^ θ xP plot for RADARSAT-2 data over San Francisco in VV-VH mode.
Remotesensing 13 04412 g0a1
In Figure A1, the deviation of data points from the 1:1 line in the plots can be observed in both VV-VH scatter plots. This deviation could be due to the presence of a partially polarized EM wave that is seized by θ xP through the utilization of m xP . On the contrary, α ^ misses this information that is obtained by eigen-decomposition of the covariance matrix. One should note that for completely polarized ( m xP = 1 ) EM wave from coherent targets, both α ^ and θ xP will either be equal to 45   or 45 .

Appendix A.2. Temporal Variations of θ xP and α ^

Table A1. Temporal variations of θ xP and α ^ over canola and wheat at different dates.
Table A1. Temporal variations of θ xP and α ^ over canola and wheat at different dates.
13 June7 July19 July24 August
Canola θ xP 40°20°18°34°
α ^ 37°18°16°32°
Wheat θ xP 39°27°20°34°
α ^ 37°26°15°28°

Appendix A.3. Temporal Variations of In-Situ Measurements of Crops

Table A2. Temporal variations of Plant Area Index (PAI), dry biomass (kg m 2 ), and Vegetation Water Content (VWC) (kg m 2 ) for canola and wheat at different dates.
Table A2. Temporal variations of Plant Area Index (PAI), dry biomass (kg m 2 ), and Vegetation Water Content (VWC) (kg m 2 ) for canola and wheat at different dates.
13 June7 July19 July24 August
CanolaPhenologyLeaf developmentFlowering stagePod developmentMaturity/ harvest
PAI1.82 ± 0.434.02 ± 0.626.32 ± 0.16N/A
dry biomass0.21 ± 0.080.43 ± 0.040.76 ±   0.06 N/A
VWC1.20 ± 0.135.82 ± 0.325.96 ±   0.20 N/A
WheatPhenologyTillering stageEarly flowering stageEarly dough stageMaturity/ harvest
PAI2.78 ± 0.315.92 ± 0.226.52 ± 0.11N/A
dry biomass0.23 ± 0.040.57 ± 0.020.98 ±   0.07 N/A
VWC2.21 ± 0.125.74 ± 0.246.11 ±   0.15 N/A

References

  1. Mandal, D.; Kumar, V.; Bhattacharya, A.; Rao, Y.S.; Siqueira, P.; Bera, S. Sen4Rice: A processing chain for differentiating early and late transplanted rice using time-series Sentinel-1 SAR data with Google Earth Engine. IEEE Geosci. Remote Sens. Lett. 2018, 15, 1947–1951. [Google Scholar] [CrossRef]
  2. Dey, S.; Mandal, D.; Robertson, L.D.; Banerjee, B.; Kumar, V.; McNairn, H.; Bhattacharya, A.; Rao, Y. In-season crop classification using elements of the Kennaugh matrix derived from polarimetric RADARSAT-2 SAR data. Int. J. Appl. Earth Obs. Geoinf. 2020, 88, 102059. [Google Scholar] [CrossRef]
  3. Mandal, D.; Kumar, V.; McNairn, H.; Bhattacharya, A.; Rao, Y. Joint estimation of Plant Area Index (PAI) and wet biomass in wheat and soybean from C-band polarimetric SAR data. Int. J. Appl. Earth Obs. Geoinf. 2019, 79, 24–34. [Google Scholar] [CrossRef]
  4. Mandal, D.; Kumar, V.; Lopez-Sanchez, J.M.; Bhattacharya, A.; McNairn, H.; Rao, Y. Crop biophysical parameter retrieval from Sentinel-1 SAR data with a multi-target inversion of Water Cloud Model. Int. J. Remote Sens. 2020, 41, 5503–5524. [Google Scholar] [CrossRef]
  5. Mandal, D.; Hosseini, M.; McNairn, H.; Kumar, V.; Bhattacharya, A.; Rao, Y.; Mitchell, S.; Robertson, L.D.; Davidson, A.; Dabrowska-Zielinska, K. An investigation of inversion methodologies to retrieve the leaf area index of corn from C-band SAR data. Int. J. Appl. Earth Obs. Geoinf. 2019, 82, 101893. [Google Scholar] [CrossRef]
  6. Mandal, D.; Ratha, D.; Bhattacharya, A.; Kumar, V.; McNairn, H.; Rao, Y.S.; Frery, A.C. A radar vegetation index for crop monitoring using compact polarimetric SAR data. IEEE Trans. Geosci. Remote Sens. 2020, 58, 6321–6335. [Google Scholar] [CrossRef]
  7. Mandal, D.; Kumar, V.; Ratha, D.; Lopez-Sanchez, J.M.; Bhattacharya, A.; McNairn, H.; Rao, Y.; Ramana, K. Assessment of rice growth conditions in a semi-arid region of India using the Generalized Radar Vegetation Index derived from RADARSAT-2 polarimetric SAR data. Remote Sens. Environ. 2020, 237, 111561. [Google Scholar] [CrossRef] [Green Version]
  8. Bhogapurapu, N.; Dey, S.; Bhattacharya, A.; Mandal, D.; Lopez-Sanchez, J.M.; McNairn, H.; López-Martínez, C.; Rao, Y.S. Dual-polarimetric descriptors from Sentinel-1 GRD SAR data for crop growth assessment. ISPRS J. Photogramm. Remote Sens. 2021, 178, 20–35. [Google Scholar] [CrossRef]
  9. Ferrazzoli, P.; Paloscia, S.; Pampaloni, P.; Schiavon, G.; Solimini, D.; Coppo, P. Sensitivity of microwave measurements to vegetation biomass and soil moisture content: A case study. IEEE Trans. Geosci. Remote Sens. 1992, 30, 750–756. [Google Scholar] [CrossRef]
  10. Davidson, M.W.; Le Toan, T.; Mattia, F.; Satalino, G.; Manninen, T.; Borgeaud, M. On the characterization of agricultural soil roughness for radar remote sensing studies. IEEE Trans. Geosci. Remote Sens. 2000, 38, 630–640. [Google Scholar] [CrossRef] [Green Version]
  11. Wiseman, G.; McNairn, H.; Homayouni, S.; Shang, J. RADARSAT-2 polarimetric SAR response to crop biomass for agricultural production monitoring. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2014, 7, 4461–4471. [Google Scholar] [CrossRef]
  12. 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. [Google Scholar]
  13. Kussul, N.; Lemoine, G.; Gallego, F.J.; Skakun, S.V.; Lavreniuk, M.; Shelestov, A.Y. Parcel-based crop classification in Ukraine using Landsat-8 data and Sentinel-1A data. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2016, 9, 2500–2508. [Google Scholar] [CrossRef]
  14. Nguyen, D.B.; Gruber, A.; Wagner, W. Mapping rice extent and cropping scheme in the Mekong Delta using Sentinel-1A data. Remote Sens. Lett. 2016, 7, 1209–1218. [Google Scholar] [CrossRef]
  15. 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. [Google Scholar] [CrossRef]
  16. Van Tricht, K.; Gobin, A.; Gilliams, S.; Piccard, I. Synergistic use of radar Sentinel-1 and optical Sentinel-2 imagery for crop mapping: A case study for Belgium. Remote Sens. 2018, 10, 1642. [Google Scholar] [CrossRef] [Green Version]
  17. Whelen, T.; Siqueira, P. Time-series classification of Sentinel-1 agricultural data over North Dakota. Remote Sens. Lett. 2018, 9, 411–420. [Google Scholar] [CrossRef]
  18. Minasny, B.; Shah, R.M.; Che Soh, N.; Arif, C.; Indra Setiawan, B. Automated Near-Real-Time Mapping and Monitoring of Rice Extent, Cropping Patterns, and Growth Stages in Southeast Asia Using Sentinel-1 Time Series on a Google Earth Engine Platform. Remote Sens. 2019, 11, 1666. [Google Scholar]
  19. Arias, M.; Campo-Bescós, M.Á.; Álvarez-Mozos, J. Crop Classification Based on Temporal Signatures of Sentinel-1 Observations over Navarre Province, Spain. Remote Sens. 2020, 12, 278. [Google Scholar] [CrossRef] [Green Version]
  20. Bousbih, S.; Zribi, M.; Lili-Chabaane, Z.; Baghdadi, N.; El Hajj, M.; Gao, Q.; Mougenot, B. Potential of Sentinel-1 radar data for the assessment of soil and cereal cover parameters. Sensors 2017, 17, 2617. [Google Scholar] [CrossRef] [Green Version]
  21. Kumar, P.; Prasad, R.; Gupta, D.; Mishra, V.; Vishwakarma, A.; Yadav, V.; Bala, R.; Choudhary, A.; Avtar, R. Estimation of winter wheat crop growth parameters using time series Sentinel-1A SAR data. Geocarto Int. 2018, 33, 942–956. [Google Scholar] [CrossRef]
  22. Dey, S.; Bhogapurapu, N.; Bhattacharya, A.; Mandal, D.; Lopez-Sanchez, J.M.; McNairn, H.; Frery, A.C. Rice phenology mapping using novel target characterization parameters from polarimetric SAR data. Int. J. Remote Sens. 2021, 42, 5519–5543. [Google Scholar] [CrossRef]
  23. Dey, S.; Bhattacharya, A.; Ratha, D.; Mandal, D.; McNairn, H.; Lopez-Sanchez, J.M.; Rao, Y. Novel clustering schemes for full and compact polarimetric SAR data: An application for rice phenology characterization. ISPRS J. Photogramm. Remote Sens. 2020, 169, 135–151. [Google Scholar] [CrossRef]
  24. Cloude, S. The dual polarisation entropy/alpha decomposition. In Proceedings of the 3rd International Workshop on Science and Applications of SAR Polarimetry and Polarimetric Interferometry, Frascati, Italy, 22–26 January 2007; pp. 22–26. [Google Scholar]
  25. Ainsworth, T.; Schuler, D.; Lee, J.S. Polarimetric SAR characterization of man-made structures in urban areas using normalized circular-pol correlation coefficients. Remote Sens. Environ. 2008, 112, 2876–2885. [Google Scholar] [CrossRef]
  26. Sugimoto, M.; Ouchi, K.; Nakamura, Y. On the similarity between dual-and quad-eigenvalue analysis in SAR polarimetry. Remote Sens. Lett. 2013, 4, 956–964. [Google Scholar] [CrossRef]
  27. Freeman, A.; Durden, S.L. A three-component scattering model for polarimetric SAR data. IEEE Trans. Geosci. Remote Sens. 1998, 36, 963–973. [Google Scholar] [CrossRef] [Green Version]
  28. Yamaguchi, Y.; Moriyama, T.; Ishido, M.; Yamada, H. Four-component scattering model for polarimetric SAR image decomposition. IEEE Trans. Geosci. Remote Sens. 2005, 43, 1699–1706. [Google Scholar] [CrossRef]
  29. Shan, Z.; Wang, C.; Zhang, H.; Chen, J. H-alpha decomposition and alternative parameters for dual Polarization SAR data. In Proceedings of the PIERS, Suzhou, China, 12–16 September 2011. [Google Scholar]
  30. Xie, L.; Zhang, H.; Wang, C.; Wu, F.; Zhang, B.; Tang, Y. Maritime application using Ha decomposition in compact and dual-pol SAR. In Proceedings of the 2013 Asia-Pacific Conference on Synthetic Aperture Radar (APSAR), Tsukuba, Japan, 23–27 September 2013; pp. 563–566. [Google Scholar]
  31. Sugimoto, M.; Ouchi, K.; Yang, C.S. On the eigenvalue analysis using HH-VV dual-polarization SAR data and its applications to monitoring of coastal oceans. In Ocean Sensing and Monitoring V; International Society for Optics and Photonics: Bellingham, WA, USA, 2013; Volume 8724, p. 87240G. [Google Scholar]
  32. Trudel, M.; Charbonneau, F.; Leconte, R. Using RADARSAT-2 polarimetric and ENVISAT-ASAR dual-polarization data for estimating soil moisture over agricultural fields. Can. J. Remote Sens. 2012, 38, 514–527. [Google Scholar]
  33. Periasamy, S. Significance of dual polarimetric synthetic aperture radar in biomass retrieval: An attempt on Sentinel-1. Remote Sens. Environ. 2018, 217, 537–549. [Google Scholar] [CrossRef]
  34. Mandal, D.; Kumar, V.; Ratha, D.; Dey, S.; Bhattacharya, A.; Lopez-Sanchez, J.M.; McNairn, H.; Rao, Y.S. Dual polarimetric radar vegetation index for crop growth monitoring using sentinel-1 SAR data. Remote Sens. Environ. 2020, 247, 111954. [Google Scholar] [CrossRef]
  35. Dey, S.; Chaudhuri, U.; Mandal, D.; Bhattacharya, A.; Banerjee, B.; McNairn, H. BiophyNet: A Regression Network for Joint Estimation of Plant Area Index and Wet Biomass From SAR Data. IEEE Geosci. Remote Sens. Lett. 2020, 18, 1701–1705. [Google Scholar] [CrossRef]
  36. Dey, S.; Chaudhuri, U.; Bhogapurapu, N.R.; Lopez-Sanchez, J.M.; Banerjee, B.; Bhattacharya, A.; Mandal, D.; Rao, Y. Synergistic Use of TanDEM-X and Landsat-8 Data for Crop-type Classification and Monitoring. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2021, 14, 8744–8760. [Google Scholar] [CrossRef]
  37. Nasrallah, A.; Baghdadi, N.; El Hajj, M.; Darwish, T.; Belhouchette, H.; Faour, G.; Darwich, S.; Mhawej, M. Sentinel-1 data for winter wheat phenology monitoring and mapping. Remote Sens. 2019, 11, 2228. [Google Scholar] [CrossRef] [Green Version]
  38. Wali, E.; Tasumi, M.; Moriyama, M. Combination of Linear Regression Lines to Understand the Response of Sentinel-1 Dual Polarization SAR Data with Crop Phenology—Case Study in Miyazaki, Japan. Remote Sens. 2020, 12, 189. [Google Scholar] [CrossRef] [Green Version]
  39. Schlund, M.; Erasmi, S. Sentinel-1 time series data for monitoring the phenology of winter wheat. Remote Sens. Environ. 2020, 246, 111814. [Google Scholar] [CrossRef]
  40. Lopez-Sanchez, J.M.; Cloude, S.R.; Ballester-Berman, J.D. Rice phenology monitoring by means of SAR polarimetry at X-band. IEEE Trans. Geosci. Remote Sens. 2011, 50, 2695–2709. [Google Scholar] [CrossRef]
  41. 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]
  42. McNairn, H.; Jiao, X.; Pacheco, A.; Sinha, A.; Tan, W.; Li, Y. Estimating canola phenology using synthetic aperture radar. Remote Sens. Environ. 2018, 219, 196–205. [Google Scholar] [CrossRef]
  43. Dey, S.; Bhattacharya, A.; Ratha, D.; Mandal, D.; Frery, A.C. Target Characterization and Scattering Power Decomposition for Full and Compact Polarimetric SAR Data. IEEE Trans. Geosci. Remote Sens. 2020, 59, 3981–3998. [Google Scholar] [CrossRef]
  44. Dey, S.; Bhattacharya, A.; Frery, A.C.; López-Martínez, C.; Rao, Y.S. A Model-Free Four Component Scattering Power Decomposition for Polarimetric SAR Data. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2021, 14, 3887–3902. [Google Scholar] [CrossRef]
  45. Barakat, R. n-fold polarization measures and associated thermodynamic entropy of N partially coherent pencils of radiation. Opt. Acta Int. J. Opt. 1983, 30, 1171–1182. [Google Scholar] [CrossRef]
  46. McNairn, H.; Tom, J.; Powers, J.; Bélair, J.; Berg, A.; Bullock, A.; Colliander, A.; Cosh, A.; Kim, M.; Ramata, S.; et al. Experimental Plan SMAP Validation Experiment 2016 in Manitoba, Canada (SMAPVEX16-MB). 2016. Available online: http://smapvex16-mb.espaceweb.usherbrooke.ca/ (accessed on 10 October 2021).
  47. Mandal, D.; Vaka, D.S.; Bhogapurapu, N.R.; Vanama, V.; Kumar, V.; Rao, Y.S.; Bhattacharya, A. Sentinel-1 SLC preprocessing workflow for polarimetric applications: A generic practice for generating dual-pol covariance matrix elements in SNAP S-1 toolbox. Preprints 2019. [Google Scholar] [CrossRef] [Green Version]
  48. Bhogapurapu, N.; Dey, S.; Mandal, D.; Bhattacharya, A.; Rao, Y.S. PolSAR tools: A QGIS plugin for generating SAR descriptors. J. Open Source Softw. 2021, 6, 2970. [Google Scholar] [CrossRef]
  49. Barakat, R. Degree of polarization and the principal idempotents of the coherency matrix. Opt. Commun. 1977, 23, 147–150. [Google Scholar] [CrossRef]
  50. Cloude, S. Polarisation: Applications in Remote Sensing; Oxford University Press: Oxford, UK, 2009. [Google Scholar]
  51. Brosseau, C. Polarization transfer and entropy transformation. Optik 1991, 88, 109–117. [Google Scholar]
  52. Bicout, D.; Brosseau, C. Multiply scattered waves through a spatially random medium: Entropy production and depolarization. J. Phys. I 1992, 2, 2047–2063. [Google Scholar] [CrossRef] [Green Version]
  53. Wu, L.k.; Moore, R.K.; Zoughi, R. Sources of scattering from vegetation canopies at 10 Ghz. IEEE Trans. Geosci. Remote Sens. 1985, 5, 737–745. [Google Scholar] [CrossRef]
  54. Cloude, S.R.; Pottier, E. An entropy based classification scheme for land applications of polarimetric SAR. IEEE Trans. Geosci. Remote Sens. 1997, 35, 68–78. [Google Scholar] [CrossRef]
Figure 1. Map of the study area showing the locations of the canola and wheat fields over the Manitoba province of Canada.
Figure 1. Map of the study area showing the locations of the canola and wheat fields over the Manitoba province of Canada.
Remotesensing 13 04412 g001
Figure 2. Feasible region of H xP / θ xP plane bounded by two curves (Curve I and Curve II). θ xP is represented in the angular direction and H xP is represented in the radial direction. The non-feasible scattering region is shaded in gray.
Figure 2. Feasible region of H xP / θ xP plane bounded by two curves (Curve I and Curve II). θ xP is represented in the angular direction and H xP is represented in the radial direction. The non-feasible scattering region is shaded in gray.
Remotesensing 13 04412 g002
Figure 3. Variations of θ xP and α ^ for different scattering target from RADARSAT-2 data over San Francisco, CA, USA.
Figure 3. Variations of θ xP and α ^ for different scattering target from RADARSAT-2 data over San Francisco, CA, USA.
Remotesensing 13 04412 g003
Figure 4. Variations of θ xP and H xP for bare field conditions and fully developed crops. These measurements are derived over in-situ fields.
Figure 4. Variations of θ xP and H xP for bare field conditions and fully developed crops. These measurements are derived over in-situ fields.
Remotesensing 13 04412 g004
Figure 5. Twelve clustering zones in the H xP / θ xP plane.
Figure 5. Twelve clustering zones in the H xP / θ xP plane.
Remotesensing 13 04412 g005
Figure 6. Temporal variation of θ xP over the study area. The dates are (a) 13 June, (b) 7 July, (c) 19 July and (d) 24 August.
Figure 6. Temporal variation of θ xP over the study area. The dates are (a) 13 June, (b) 7 July, (c) 19 July and (d) 24 August.
Remotesensing 13 04412 g006
Figure 7. Temporal variation of θ xP and m xP over the canola fields.
Figure 7. Temporal variation of θ xP and m xP over the canola fields.
Remotesensing 13 04412 g007
Figure 8. Temporal variation of clusters on H xP / θ xP plane over the canola fields. The dates are (a) 13 June, (b) 7 July, (c) 19 July, and (d) 24 August.
Figure 8. Temporal variation of clusters on H xP / θ xP plane over the canola fields. The dates are (a) 13 June, (b) 7 July, (c) 19 July, and (d) 24 August.
Remotesensing 13 04412 g008
Figure 9. Temporal variation of clusters on H xP / θ xP plane over the wheat fields. The dates are (a) 13 June; (b) 7 July; (c) 19 July and (d) 24 August.
Figure 9. Temporal variation of clusters on H xP / θ xP plane over the wheat fields. The dates are (a) 13 June; (b) 7 July; (c) 19 July and (d) 24 August.
Remotesensing 13 04412 g009
Figure 10. Temporal variation of θ xP and m xP over the wheat fields.
Figure 10. Temporal variation of θ xP and m xP over the wheat fields.
Remotesensing 13 04412 g010
Figure 11. Images of temporal variation of clusters on H xP / θ xP plane over the study area. The dates are (a) 13 June; (b) 7 July; (c) 19 July and (d) 24 August.
Figure 11. Images of temporal variation of clusters on H xP / θ xP plane over the study area. The dates are (a) 13 June; (b) 7 July; (c) 19 July and (d) 24 August.
Remotesensing 13 04412 g011
Table 1. Specification for Sentinel-1 data acquired for the Canadian test site.
Table 1. Specification for Sentinel-1 data acquired for the Canadian test site.
Acquisition
Date
Beam
Mode
Incidence Angle
Range (Deg.)
Orbit az ( m ) × rg ( m )
13 June 2016IW30.22–32.47Ascending 15 × 15
07 July 2016IW30.22–32.47Ascending 15 × 15
19 July 2016IW30.22–32.47Ascending 15 × 15
24 August 2016IW30.22–32.47Ascending 15 × 15
Table 2. Temporal variation in the percentage of data points in each zone for different phenology stages of canola. The zones with high percentages of points at a particular phenology stage are highlighted in bold.
Table 2. Temporal variation in the percentage of data points in each zone for different phenology stages of canola. The zones with high percentages of points at a particular phenology stage are highlighted in bold.
DatesZ1Z2Z3Z4Z5Z6Z7Z8Z9Z10Z11Z12
13 June2.00.00.058.00.00.035.40.00.04.60.00.0
7 July0.00.00.00.00.00.00.00.00.06.291.62.2
19 July0.00.00.00.00.00.00.00.00.06.293.80.0
24 August0.00.00.00.00.00.012.60.00.064.622.80.0
Table 3. Temporal variation in the percentage of data points in each zone for different phenology stages of wheat. The zones with high percentages of points at a particular phenology stage are highlighted in bold.
Table 3. Temporal variation in the percentage of data points in each zone for different phenology stages of wheat. The zones with high percentages of points at a particular phenology stage are highlighted in bold.
DatesZ1Z2Z3Z4Z5Z6Z7Z8Z9Z10Z11Z12
13 June0.00.00.027.00.00.047.90.00.010.414.70.0
7 July0.00.00.00.00.00.00.00.00.018.881.20.0
19 July0.00.00.00.00.00.00.00.00.00.083.316.7
24 August0.00.00.00.00.00.031.30.00.029.139.60.0
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Dey, S.; Bhogapurapu, N.; Homayouni, S.; Bhattacharya, A.; McNairn, H. Unsupervised Classification of Crop Growth Stages with Scattering Parameters from Dual-Pol Sentinel-1 SAR Data. Remote Sens. 2021, 13, 4412. https://0-doi-org.brum.beds.ac.uk/10.3390/rs13214412

AMA Style

Dey S, Bhogapurapu N, Homayouni S, Bhattacharya A, McNairn H. Unsupervised Classification of Crop Growth Stages with Scattering Parameters from Dual-Pol Sentinel-1 SAR Data. Remote Sensing. 2021; 13(21):4412. https://0-doi-org.brum.beds.ac.uk/10.3390/rs13214412

Chicago/Turabian Style

Dey, Subhadip, Narayanarao Bhogapurapu, Saeid Homayouni, Avik Bhattacharya, and Heather McNairn. 2021. "Unsupervised Classification of Crop Growth Stages with Scattering Parameters from Dual-Pol Sentinel-1 SAR Data" Remote Sensing 13, no. 21: 4412. https://0-doi-org.brum.beds.ac.uk/10.3390/rs13214412

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