Complementarity of Sentinel-1 and Sentinel-2 Data for Mapping Agricultural Areas in Senegal ()
1. Introduction
Senegal is a country where agriculture is mainly rainfed (with only 5% of the land being irrigated [1], making agricultural production subject to climatic hazards. In recent years, several public and private actions have been initiated to improve agricultural productivity, but the agricultural system is especially vulnerable to climate change. The lack of reliable and recent data is a barrier to good management of the sector’s challenges. In this context, remote sensing can provide considerable benefits to improve the accuracy of agricultural data and statistics.
In recent years, many studies have shown the value of using remote sensing for agricultural monitoring in divers situations [2] - [9]. Satellite remote sensing techniques for monitoring and mapping crops are mainly based on the approach of spectral treatment using thematic indices such as NDVI (Standardized Difference Vegetation Index) and LAI (Leaf Area Index) [10] [11] [12]. The main limitation of this approach has the low discrimination in the context of very heterogeneous micro-plots of culture [13]. In Senegal, high crop diversity and small parcel size are important challenges to crop mapping at the plot scale using Earth observation data.
In the recent years, the availability of time series of optical and radar images time series with high spatial and temporal resolutions has brought particular interest for agricultural mapping at the plot scale. This is especially true since the different crops to be discriminated do not have the same phenological cycles [14] [15] [16].
With the high sensitivity of certain wavelengths to chlorophyl activity and also to plant characteristics (density, biological cycle, etc.), optics is an effective means of monitoring vegetation [17] [18] [19] [20] [21]. Sentinel-2 data, with 13 spectral bands provides a good solution for distinguishing between different crops. However, the sensitivity of optical sensors to atmospheric conditions is often a limitation for monitoring vegetation during the rainy season due to cloud cover.
Conversely, radar data are not disturbed by atmospheric conditions but their sensitivity to different surface parameters (soil roughness, soil moisture and vegetation) limit their efficiency in the monitoring and mapping of herbaceous vegetation [22]. However, analysing the temporal profile of these data under different polarizations can provide important information to discriminate types of crops, especially because of differences in phenological cycles plant shapes and sizes during the various stages of the growing season [23].
This work aims to test the complementarity of optical and radar images, on the one hand, and the temporal evolution of these two types of signals during the agricultural season on the other hand. By combining time series for Sentinel-2 (optical) and Sentinel-1 (C-band radar), we expect to be able to discriminate between different crops at the plot scale, enabling better monitoring of agricultural areas in the Senegalese context.
2. Materials and Methods
2.1. Study Area
The study area, about 400 km2, is in the department of Nioro du Rip (Figure 1). It is located between latitudes 13.64˚ and 13.82˚ north, and longitudes −15.74˚ and −15.92˚ west and belongs to the groundnut production basin, which is the most important agricultural region of Senegal.
The climate is of the Sahelian type, with two seasons: a dry season (mid-October to mid-June) and a rainy one (mid-June to mid-October) marked by a strong interannual variability. Each year, between 900 - 1100 mm is recorded during the three months of the rainy season, causing strong hydric erosion often causing extensive gullying.
Demographic pressure has strongly modified the agricultural system, notably with deforestation and the abandonment of fallow cropping cycles [24].
The agricultural system is shared between cash crops dominated by peanuts and cereal crops mainly composed of millet and maize grown on small parcels of about 2 hectares.
2.2. Field Data
The field data consisted of 168 control plots (59 of peanut, 55 of maize and 54 of millet), serving as a training and validation plots (Figure 2). Field missions were
Figure 2. Spatial repartition of control plots.
carried out during the crop growth period to delimit control plots contours (digitization of contours) with GPS (Global Position System). The choice of plots was based on several criteria including the spatial representativeness of each crop, the level of growth and density. This to have a good spatial representativeness of each type of culture. In order to avoid edge effects and eliminate mixed pixels between contiguous plots of different crops, plots borders were adjusted from Google Earth.
2.3. Satellite Image Data
Two datasets from multiple dates during the 2016 rainy season (June to October) were used for Sentinel-2A (Optical) and Sentinel-1A (Radar), as shown in Table 1.
Sentinel-2A multispectral images are available on the ESA (European Space Agency) website (https://sentinel.esa.int/web/sentinel). These are level-1C products, with resolutions from 10 meters (visible + near-infrared bands) to 60 meters, organized in 100 km × 100 km ortho-rectified tiles in UTM WGS84 projections. Data from the T28PDA tile, which fully covers the study area were used. Although the satellite imaging periodicity is 10 days, only three images had less than 10% cloud cover rate: 0! And 18 of August and 07 of October, with 0.22%, 5.8% and 0% of cloud cover, respectively.
Sentinel-1 data, are available on the https://scihub.copernicus.eu/dhub website and on the Sentinel product exploitation platform (PEPS) http://scihub//peps.cnes.fr. Sentinel-1 is a C-band Synthetic Aperture Radar (SAR) sensor that provides images at 10 meters of spatial resolution in two VV and VH polarizations mode. As radar data are not affected by cloud cover, a total of 14 images were analysed.
Table 1. Satellite data acquisition date.
3. Methods
The work was done in four steps (Figure 3).
The first step in the work was to correct the satellite data. Downloaded optical data (Level 1C) were corrected using ESA’s Sen2Cor algorithm to suppress atmospheric effects and convert numerical values to reflectance (Level-2A). Radar data were corrected for geometric distortions (geocoding, georeferencing), and filtered to reduce speckle. The data were then normalized (radiometric calibration) in order to be able to perform the multi-temporal analyses. This required transforming the backscattered signal (Digital Number) into a dimensionless quantity called sigma backscattering coefficient (σ˚ in dB) proportional to the ratio of power received and power transmitted at the antenna [25].
Secondly pixel values were extracted for control plots for all the images (optical and radar), with plot boundaries being adjusted to avoid contiguous pixels.
In the third step, data was analysed to explore their ability to differentiate between the three types of crops. For optical data, the spectral signature of the three cultures for the three dates (08 and 18 of August, 07 of October) was carried out to determine the bands and the periods for optimal discrimination. Pair-wise comparison of individual parcels across the three dates where used to determine which date would provide optimal discrimination for subsequent NDVI analysis. The dispersion of radar data was evaluated by comparing the mean and median values of the backscattering coefficients of each crop to see the influence of the extreme values. This is usually caused by the sensitivity of the radar signal to humidity and soil roughness, especially at the beginning of the growing season when the plant cover is low [22] [26]. Subsequently, the temporal profiles of the signals of the three cultures following under two polarizations (VV and VH) were analysed and compared in order to look for differences in backscattering between the cultures during their phenological cycle. These analyses showed that the relationship between the two polarizations can improve the differentiation of cultures. An index called IG (IG = 2 * Exp(VH/VV) was calculated for each date and its temporal profile analysed. the purpose of this index is to highlight the spectral differences of crops that the polarizations taken individually do not able to show. Finally, the differences between the signals of the cultures taken two by two were analysed to identify dates exhibiting significant differences (more than 5 dB). This allowed the selection of the most discriminating images.
Finally, the most discriminating optical and radar data were combined to map different crops. Two thirds of the control plots were randomly selected to provide a training base. Segmentation was then used to delineate the contours of the
homogeneous sets grouping pixels which had similar criteria and form them into subsets, homogeneous regions or segments [27]. This was undertaken with Idrisi Silva software which enables several images with different characteristics (spatial resolution, sensors, etc) to be combined. The combination of segmentation and classification results was used to produce yielded a crop map which accuracy was then tested using the one third of the control plots.
4. Results and Discussions
4.1. Optical Data
The analysis of the spectral signature of the three crops for the three dates showed a strong confusion in the visible, especially during the growth phase (August). However, towards the end of the cropping season (October images), differences between peanut and cereal crops (millet and maize), were noted especially in the near infrared (PRI, 730 nm) as shown in Figure 4. There were also small differences in the red band (665 nm). This shows that the NDVI could differentiate the three crops during this period (October). For this reason, investigations were then limited to NDVI in the rest of the optical data analysis.
NDVI values for August were strongly correlated (Figure 5(a)) but overall maize values were low (less than 0.18 on August 8th and 0.3 on August 18th). In October, the NDVI of the peanut remained generally higher than that of millet (Figure 5(b) and Figure 5(c)). It can be deduced that August was not suitable to discriminate crops. The NDVIs values of October, showed a good discrimination of the three crops, in particular between peanut and millet. Further analysis was based solely on the image of October 7th.
Figure 4. Spectral signature of the three crops in: (a) August 08; (b) August 18 and (c) October 07, 2016.
Figure 5. Comparison of the NDVIs of crops of the three dates.
4.2. Radar Data
The radar signals were sensitive to seasonal variability for all three crops. This is consistent with other research obtained for C-band radar over the entire Sahelian band [28] [29]. On the other hand, the VV polarization signal remained higher than that in VH. This is because the penetration ability of a wave through a plant layer depends on the geometry, the orientation, the size, the density as well as and the dielectric constant of the elementary diffusers constituting them [30].
The signals for each of the three crops were similar during the growth phase (Figure 6(a)). This situation is explained by the fact that during this period, the contribution of soil (roughness and humidity) was often greater than that of vegetation because of the low vegetation cover [22]. It was therefore difficult to discriminate crops during the first part of the rainy season. However, the IG index showed a difference between millet and the two other crops at the 192nd and 204th day of year corresponding to July 11 and July 23. This difference could be the result of cultural practices that differ from one crop to the other. Millet is sown before the onset of rains unlike maize and peanut where the crop requires deep plowing. As a result, the roughness of the soil is very different between millet fields and those of maize and peanut. This explains the fact that the VH/VV ratio amplified by the exponential function differentiated millet from the other two crops.
During the period of maturity and harvest, significant differences were observed between peanut and cereal crops (Figure 6(a)). This situation is explained by the fact that harvesting the peanut crop leaves a bare soil, unlike the harvest of the two cereal crops (millet and maize) which only removes the ears and leaves the stalks in place to slowly dry and die. Thus, peanut signal drops sharply after the harvest, when millet and maize are slowly declining.
Pairwise comparison of the signals for the different crops showed that the difference between the peanut signal and the two cereals can reach 2 dB in VV and up to 5 dB in VH towards the end of the rainy season (326th day of year corresponding to November 20th). This shows that it is possible to discriminate peanut from cereal crops in that period. However, the difference between the millet and maize signal was very small and rarely exceeded 1 dB (Figure 7).
In conclusion, the images of July 11 (193rd Julian days) and November 20 (324th Julian days) were the most discriminating. In fact, the IG index of July 11 differentiated the millet from the two other crops and the image of November 20th discriminated peanut from millet and maize in VV and VH polarization.
Figure 6. Temporal evolution of three crops signal (a) in VV and VH polarisation (b) IG index (2 * Exp(VH/VV)).
Figure 7. Temporal evolution of the differences between crops signal taken in pairs in (a) VV polarisation and (b) VH polarization.
4.3. Combining of Optical and Radar Data
The false color compositions between the different images (NDVI of October 7th, VV and VH polarization of July 11 and November 20) make it possible to highlight differences in land cover (Figure 8):
NDVI combined with VH polarization (July 11th and November 20th Figure 8(a)) or with VV (July 11th and November 20th Figure 8(b)) highlighted habitations (yellow) and water (black);
The combination of the NDVI with the VV and VH of July 11 (Figure 8(c)) brought out the natural vegetation zones (green);
The NDVI with VV and VH of November 20 (Figure 8(d)) showed differences between crops (black for peanut, green for millet and yellow maize);
The image of July 11th (VV and VH polarization) with the VH of November 20th (Figure 8(e)) or with VV of November 20th (Figure 8(f)) highlighted the heterogeneity of agricultural areas, although there is significant noise was observed;
The image of November 20th (VV and VH) combined with VH of July 11th (Figure 8(g)) or with the VV of July 11th (Figure 8(h)) showed differences between crops.
Combined optical and radar data combination reduced the noise. This shows once again the importance of associating these two types of data for crop mapping. The radar showed a sensitivity of the different crops at different dates but with a background noise on the images while the optical image (NDVI October 07) is able to reduce noise and facilitated the detection of contours.
The segmentation performed from the five combined images (NDVI of October 7th, VV and VH of July 11th and VV and VH of November 20th) allowed the distribution the pixels in homogeneous segments. The segmentation results were merged with the result of the supervised classification to smooth and refine the result (Figure 9). Finally, to assess the validity of the result, a confusion matrix was calculated using 1/3 of the remaining data (Table 2). The overall error of the classification was 14% with a Kappa coefficient of 0.83. This error is mainly due to the strong confusion between the two cereals where 32% of the maize fields were classified as millet and 30% of millet fields were considered as maize. Thus, maize presented a global classification error of 38% and millet has a 34% error, while the error in classifying peanut only 7%. Non-agricultural areas were well identified with much lower errors for water (0.5%), natural vegetation (4.6%), habitation (2.4%) and bare soil (6.2%).
Optical signals for the three crops were poorly distinguished during the growing season (August). This is partly because during this period the coverage rates are almost identical for all crops [23] and so, all fields reflected in almost the same way making the differentiation of crops difficult. Only the infrared wavelength detected a slight difference for maize or the reflectance remained lower than that of millet and peanut. In October, the peanut NDVI stood out clearly from that of millet with a sort of barrier around 0.55. This is explained by
Figure 8. Color composite images created by combining optical and radar data.
Figure 9. Land cover of the area obtained.
Kappa = 0.8251.
the late harvest of peanut which began in the last 10-day period of October when the millet and maize are already harvested. Thus, during the first 2/3 of October, peanut was still in a phase of intense chlorophyllous activity.
The radar data confirmed the differences noted with the optical data but the time series provides greater precision (14 radar images against only 3 optics). Until day 253 of the year (09 September) radar signals of the three crops remained confused under both VV and VH polarizations. However, the IG index (amplified ratio VH/VV) showed a slight difference between millet and the two other crops on day 192 and day 204 of year (July 11 and 23). This difference could result from cropping practices where peanut and maize plowing is often much more marked than for millet, as the two polarizations do not have the same sensitivities to roughness [31]. Between the day 253 and day 300 of the year (September 9th and October 27th), the signals of millet and maize were lower than that of peanut under both polarizations. This, because this period corresponded to the harvest of cereal crops (millet and maize), while peanut has not yet completed its phenological cycle. However, there was a reversal of signals beyond October 27 (day 300 of the year) with a sharp drop of peanut signal corresponding to the beginning of its harvest during which all the plant is cut to leave a completely bare soil. During this period, the difference between the peanut and cereal (millet and maize) backscatter coefficient reached 2 dB in VV and 5 dB in VH polarization, indicating the possibility of discriminating this crop at the end of the rainy season.
Combining optical and radar images show strong potential for distinguishing between the three crops with more than 85% of precision. However, the differentiation of millet and maize remained difficult especially with the radar where the difference in signals rarely exceeds 1 dB throughout the rainy season. It’s only in December that there is a real difference (up to 2 dB). Only the IG index of 11 July and the NDVI of October 07 discriminated millet and maize but with small differences.
5. Conclusions
The paper examined the capacity of satellite data to map crops (peanut, millet and maize) in the Sahelian context where parcels are small (less than 2 hectares).
The purpose of this study is to achieve a multi-temporal and multi-sensor analysis method that allows mapping small crop parcels by combining optical (Sentinel-2) and radar (Sentinel-1) data.
Optical satellite data gives a better result towards the end of rainy season, period during which we note differences between the NDVI of millet and that of peanut
Radar data show, by combining the two bands (VV and VH), differences between crops at the beginning of season. This may relate to differences in agricultural practices. Indeed, peanut and maize are usually sown by plowing, thereby affecting soil roughness which influences the radar signal in VV and VH.
Significant differences in the two polarizations (VV and VH) were noted towards the end of the rainy season (October) due to the fact that the three crops were not harvested at the same period and in the same way.
Combination of optic and radar data made it possible to considerably improve the discrimination of the three crops with a global accuracy of the 82%. However, the capacity of satellite data to discriminate the three crops depends on the period during the growing season.
In order to improve the results once again, mathematical and classification algorithms are being developed.
Acknowledgements
This work is carried out within the framework of the research project funded by Centre d’excellence Africain en Mathématique, Informatique et Technologie de l’Information et de la Communication CEA-MITIC of Gaston Berger University of Saint-Louis, Senegal. We warmly thank CEA-MITIC for accepting this work.