Next Article in Journal
Attribution Analysis of Climate and Anthropic Factors on Runoff and Vegetation Changes in the Source Area of the Yangtze River from 1982 to 2016
Next Article in Special Issue
An Alternative Method for the Generation of Consistent Mapping to Monitoring Land Cover Change: A Case Study of Guerrero State in Mexico
Previous Article in Journal
“No Entry into New South Wales”: COVID-19 and the Historic and Contemporary Trajectories of the Effects of Border Closures on an Australian Cross-Border Community
Previous Article in Special Issue
Spatiotemporal Change Analysis and Future Scenario of LULC Using the CA-ANN Approach: A Case Study of the Greater Bay Area, China
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Multispectral Sentinel-2 and SAR Sentinel-1 Integration for Automatic Land Cover Classification

1
Italian Institute for Environmental Protection and Research (ISPRA), Via Vitaliano Brancati 48, 00144 Rome, Italy
2
Department of Innovation in Biology, Agri-Food and Forest systems (DIBAF), University of Tuscia, Via San Camillo de Lellis SNC, 01100 Viterbo, Italy
3
Dipartimento di Scienze della Terra, Università degli Studi di Firenze, Via La Pira 4, 50121 Firenze, Italy
4
geoLAB-Laboratory of Forest Geomatics, Department of Agricultural, Food and Forestry Systems, University of Florence, Via San Bonaventura, 13, 50145 Firenze, Italy
5
Dipartimento di Bioscienze e Territorio, Università degli Studi del Molise, C/da Fonte Lappone, 86090 Pesche, Italy
*
Author to whom correspondence should be addressed.
Submission received: 29 April 2021 / Revised: 2 June 2021 / Accepted: 2 June 2021 / Published: 7 June 2021

Abstract

:
The study of land cover and land use dynamics are fundamental to understanding the radical changes that human activity is causing locally and globally and to analyse the continuous metamorphosis of landscape. In Europe, the Copernicus Program offers numerous territorial monitoring tools to users and decision makers, such as Sentinel data. This research aims at developing and implementing a land cover mapping and change detection methodology through the classification of Copernicus Sentinel-1 and Sentinel-2 satellite data. The goal is to create a versatile and economically sustainable algorithm capable of rapidly processing large amounts of data, allowing the creation of national-scale products with high spatial resolution and update frequency for operational purposes. Great attention was paid to compatibility with the main activities planned in the near future at the national and European level. In this sense, a land cover classification system consistent with the European specifications of the EAGLE group has been adopted. The methodology involves the definition of distinct sets of decision rules for each of the land cover macro-classes and for the land cover change classes. The classification refers to pixels’ spectral and backscatter characteristics, exploiting the main multi-temporal indices while proposing two new ones: the NDCI to distinguish between broad-leaved and needle-leaved trees, and the Burned Index (BI) to identify burned areas. This activity allowed for the production of a land cover map for 2018 and the change detection related to forest disturbances and land consumption for 2017–2018, reaching an overall accuracy of 83%.

1. Introduction

Land cover is defined as the biophysical cover of Earth’s surface, as artificial surfaces, agricultural areas, forests and semi-natural areas, water and wetlands by the Directive 2007/02 of the European Commission [1,2]. Different land cover types furnish various ecosystem services and support habitats; land cover is also considered a climate variable, an element for the description of the environment [3]. Human activities depend on land cover and its evolution [4]. Therefore, the comprehension of the processes that lead to the transformation of the territory are crucial for sustainable resource management, environmental protection, modelling of climate change scenarios, and for the assessment of the impacts on the earth system directly linked to human activities [5]. Remote sensing is important to the study of Earth phenomena, which can be monitored and assessed at different spatial and temporal resolutions [6,7].
Several user communities such as decision-makers, non-governmental organizations, European communities, scientists and researchers require different types of land cover information for their activities and policies [8]. For example, land cover information is used for evaluating progress towards UN Sustainable Development Goals (SDGs) targets [9], such as target 15.3, which is related to the reach of the Land Degradation Neutrality (LDN) by 2030 [10]. Despite the importance of land cover data in environmental monitoring and planning, the products currently available on a national scale are limited in number and their characteristics are not always adequate. The main national land cover and land use products are part of the Copernicus program (https://land.copernicus.eu, accessed on 28 February 2021). The CORINE Land Cover (CLC) guarantees a high thematic detail and a long historical series, but is limited from the point of view of spatial detail and updating frequency (https://land.copernicus.eu/pan-european/corine-land-cover, accessed on 28 February 2021). More recently, the High Resolution Layers (HRLs) have allowed the description with high spatial detail of the main land cover classes, maintaining a multi-year updating frequency (https://land.copernicus.eu/pan-european/high-resolution-layers, accessed on 28 February 2021). There are other Copernicus products on a national scale, with high thematic and spatial detail but related to specific areas, such as Urban Atlas, Riparian Zones, Coastal Zones and Natura 2000 (https://land.copernicus.eu/local, accessed on 28 February 2021). At the national level there is the National Land Consumption Map produced and updated by ISPRA and other land use and land cover maps on a regional scale. These regional products are often derived from CLC and are not up to date, with some exceptions such as in the Lombardy region. In this context, it is crucial to exploit the tools introduced by the Copernicus program, such as Sentinel data and DIAS computing resources. In 2014 the European Space Agency (ESA) launched Sentinel-1 and Sentinel-2 satellites, part of the Space component of the Copernicus program of Earth observation [11]. Sentinel-1 performs C-band synthetic aperture radar (SAR) imaging. It operates both day and night, acquiring images regardless of the weather and supporting the absence of optical images, with a revisit period of six days [12,13,14]. Sentinel-1 interacts with the structural properties of elements in different ways and through different signal polarizations (VV, VH) in relation to their roughness and moisture content [15,16]. The Sentinel-2 constellation (including Sentinel-2A and Sentinel-2B) has a 5-day revisit cycle and provides 13-band multispectral images from visible to Short Wave Infrared (SWIR), associated with three spatial resolutions from 10 to 60 m. Sentinel-2 data characteristics guarantee good results in the monitoring of land use and cover [17,18,19]. Much research has confirmed that the joint use of information provided by Sentinel-1 and Sentinel-2 lead to overcoming the limitations of using single data products in the land cover (LC) and land cover changes (LCG) classification [20], allowing the evaluation of different object characteristics, such as pigmentation, moisture, size and density [21,22,23]. The integration of these two data sources can improve land cover detection [24,25], for example for grassland [26], urban areas [27] and land cover changes [21] or hazard assessment [28].
Land cover and land use (LU) are interconnected and influence each other and their changes impact on Earth’s natural resources and balances [29]. The description of territory through its land use and land cover characteristics has led to the creation of many classification systems, conceived of different purposes and referring to various levels of detail and territorial areas. In many cases, a mixed approach has been adopted, which takes into consideration both the characteristics of the land cover and of land use, as in the case of CORINE Land Cover. The problem of comparing, interpreting and using data referring to different classification systems is well known [30]. However, the increase in information sources and classification techniques makes a standardization of classification systems more necessary. To satisfy the requirements proposed by European legislation for monitoring activities, both at the European and national level, a list of criteria for a new data model has been created, and the EIONET Action Group on Land monitoring in Europe (EAGLE) Concept [31] aims at being a tool for analytic class definitions and for linking recent and future nomenclatures. Moreover, it allows the harmonisation of European land monitoring systems and can be implemented as an object-oriented base for mapping and monitoring research. It is not a classification system, but a tool for unifying LC and LU information both in top-down and bottom-up approaches. The EAGLE concept is based on a clear distinction between land cover and land use and is summarized in a matrix composed by three separate descriptors: land cover components (LCC), land use attributes (LUA) and further characteristics (LCH). The descriptors can be combined in order to create specific classification systems for different needs, or to identify correspondences with existing classes, keeping the three descriptors independent [32].
There are many algorithms for investigating land cover and land use starting from the classification of satellite images. The choice of the most suitable method depends on aspects such as the type of data, distribution of classes, research purposes and interpretability of the classifier, and also the balance between objectives and the available resources. In general, the automatic classification methods are time consuming with the increase of the data dimension and volume, and sometimes data interpretation can become difficult [3,33]. Supervised classifications can be used to analyse a large amount of data; they are based on the selection of an adequate number of training samples with known values [34], which are used to predict unknown values of testing data [35]. For this reason, it is important to select examples which completely represent the variability of the characteristics of the analysed territory. In the past few years the supervised classification technique increased in popularity [36], thanks to the availability of ancillary data, which provide accurate information on ground truth [37]. Ground truth information is also important to define thresholds values for the classification; these values are important for isolating objects from the background, in order to improve classification precision [38]. Threshold methods are used for their simplicity and because they require few parameters [39,40]. The identification and use of thresholds related to the spectral characteristics and their variation over time are useful to assess both physical and biophysical characteristics of land cover [41,42], to classify cover change [43], to define crop phenology for determining start and end of growing season [40], and to classify tree species [44,45] or burned areas [46]. Some supervised classification methods, for example parametric supervised classifiers (maximum likelihood, minimum distance), are difficult to use for classifying large multi-temporal datasets because they have limited flexibility in decision boundaries. Since these classifiers assume a normal data distribution, they furnish excellent results when data distribution is unimodal but they can be difficult to apply for multi-modal dataset analysis [47]. Machine learning [48] approaches are largely used in land cover mapping, thanks to their capacity to model class signatures, to elaborate many input data and to produce higher accuracy compared to traditional parametric classifiers, especially for complex data with many predictor variables [49]. For example, decision trees and random forest are machine learning approaches which focus on decision rules on class boundaries which improve the classification of land cover types with unknown distribution and frequency [50,51,52]. On the other hand, machine learning is characterized by higher computational intensity, unknown decision rules (black box) and the need of input parameters [3,53]; moreover, it is not easy to apply to large areas, because it is characterized by difficult computing, such as the manual choice of a large number of training areas [54] and it is negatively influenced by noisy observations and over-fitting [3].
This research aims to develop a methodology that can provide land cover and land cover change products for operational purposes which can achieve the following project requirements:
-
allow the production of a 10 m land cover map on a national scale, exploiting the potential offered by free Sentinel-1 and Sentinel-2 data;
-
allow at least annual updates for the entire Italian territory;
-
ensure accuracy in line with Copernicus products;
-
adopt a structure that allows the updating and expansion of the methodology;
-
ensure compatibility with the main European and national initiatives in the field of land monitoring, such as Copernicus program and the Mirror Copernicus initiative of the Italian Ministry of Economic Development.
The main challenge of this research regards the creation of a tool capable of meeting all the indicated requirements, and this is also demonstrated by the absence of national and European products with these characteristics. In this sense, proven classification techniques were used, obtaining a methodology based on training areas and decision rules [55,56] for the production of pixel-based land cover classification products [57].
The study is in continuity with other previous analysis conducted by this working group on land consumption [58,59] and forest areas [60]. For the identification of the land cover classes mapped by the algorithm, an EAGLE-compliant classification system has been developed. In detail, eight land cover classes have been identified and characterized:
  • unvegetated abiotic surfaces (subdivided into artificial surfaces and natural surfaces);
  • vegetated biotic surfaces (distinguishing between woody vegetation, subdivided into broad-leaved and needle-leaved, and permanent and periodically herbaceous vegetation);
  • water surfaces (water bodies and permanent snow and ice).
Change detection classes regarded variations in artificial surfaces (land consumption and restoration) and forest disturbances, distinct into fires and other types of disturbances.
In order to test the methodology, this was applied to the Sentinel-1 and Sentinel-2 data series relating to 2018 for the classification of the land cover, and to the period 2017–2018 for the detection of land cover changes.
The methodology allowed the identification of land cover classes by exploiting only Sentinel-1 and Sentinel-2 data. The only ancillary data used is the National Land Consumption Map of ISPRA for the distinction of artificial abiotic surfaces.

2. Materials and Methods

2.1. Overview

The methodology presented in this study uses Sentinel-1 GRD and Sentinel-2 images to map land cover and land cover changes in Italy, on a yearly basis (Figure 1).
Decision rules are defined to identify land cover classes at pixel size of 10 m, setting fixed threshold values on composites (e.g., median, maximum) of multitemporal data. The rules are based on experience, on deep knowledge of literature and on previous analysis of temporal trend of spectral signatures, indexes and backscatter coefficients of each class. In this sense, the study is in continuity with other analysis published in previous scientific papers related with land consumption [58,59] and forest area monitoring [18,60]. Furthermore, the thresholds for vegetation classes were also defined using a collection of training areas. The output classifications were validated through the photo-interpretation of very high-resolution images.
The methodology was developed with Google Earth Engine (Javascript language program) but was also adapted for QGIS Semi-Automatic Classification Plugin and on a DIAS environment (ONDA https://www.onda-dias.eu, accessed on 4 March 2021). This DIAS experimentation refers to the agreement between the Italian Space Agency (ASI) and ISPRA for habitat mapping.

2.2. Study Area

The study area is the Italian territory (Figure 2), characterized by a total surface of 301,338 km2. Geographically it is divided into three parts: The Alpine arc in the north, the peninsula in the south central part and the insular area that includes Sardinia and Sicily. The climate in the Italian territory is alpine in the north, continental in the peninsula and Mediterranean along the coast, with drier summers and milder winters. The average annual rainfall is higher in the mountain areas with a maximum of 2500–3500 mm and reaches a minimum of 500 mm in Apulia and Molise.
The territory is mainly composed of mountains (35.0%) or hills (41.6%) and plains (23.2%). The mountain belts are the northern area, forming the Alps with peaks reaching over 4000 m in the western part. The Appennines chain extends over the whole peninsula, reaching its highest peak with Gran Sasso (2912 m), located in Abruzzo. Due to these characteristics, the northern sector and the central chain are often covered by clouds or affected by perturbations that make it difficult to find cloud-free optical images for various periods of the year.
Land cover is characterized by forest, grassland, built-up areas and water (both liquid water and ice). Herbaceous vegetation is dominated by cropland and pastures and is mainly concentrated on coastal areas and, in Padana, plains. The orography conditions the distribution of the settlements, which are concentrated in the plains and on the coastal areas.

2.3. Classification System

The classification system defined for this research is consistent with the EAGLE group specifications [32,61] and it is composed of eight land cover classes, organized into three levels (Table 1), and four land cover change classes (Table 2).
The first classification level distinguishes three main land cover classes. It consists of land cover classes at the first two levels, while the “biotic-vegetated” class is then further articulated through a third classification level based on biophysical characteristics (broad-leaved, needle-laved) and temporal parameters (permanent herbaceous, temporary herbaceous).
1.
Abiotic-non vegetated: The class includes any unvegetated surfaces, either covered with man-made artificial structures or geologically natural material surfaces (with or without anthropogenic influence or impact). It includes consolidated or unconsolidated unvegetated natural surfaces, such as rocky regions and sands.
At the second classification level, the class is subdivided into artificial and natural surfaces.
-
Artificial abiotic surfaces include impervious surfaces but also those of reversible land consumption [62]. Impervious and sealed surfaces are mainly covered with features with a certain height above ground (artificial buildings and constructions) or features without a specific height above ground (flat waterproof surfaces) such as artificial surface pavements (e.g., asphalt, concrete, tarmacadam). Reversible land consumption includes all areas where natural surface material has been replaced by man-made material or where natural material has been removed, forming a non-impermeable and undeveloped surface. This includes soil compaction, excavation or temporary impervious cover. For example: unpaved roads, construction sites or courtyards or sports fields, permanent deposits of material, photovoltaic fields, and quarries not yet restored. The EAGLE classification system includes quarries and extraction sites in the 1.2.1.1 bare rock class, but in this research they are classified as artificial abiotic surfaces since quarries and extraction sites are strongly modified due to compaction and morphology changes caused by material removal. The physical characteristics of these areas are more similar to those of the EAGLE 1.1.2 “Non-Sealed Artificial Surfaces” class, which includes any open areas where natural surface material has been replaced by artificial material or where natural material has been removed from its place of origin as result of human activity forming an unsealed and non-built-up surface.
-
Natural abiotic surfaces are any kind of surface material that remains in its natural consistence or form, either with or without anthropogenic influence. Consolidated and unconsolidated surfaces include: unvegetated rocky mountainous regions, sand, and bare soil.
2.
Biotic-vegetated: The class includes any terrestrial surface with spontaneous, semi-natural or artificial vegetation (e.g., crops and urban parks), with or without anthropogenic influence.
At the second classification level, biotic-vegetated surfaces are divided into woody vegetation and herbaceous vegetation.
-
Woody vegetation includes perennial woody plants with a single, self-supporting main stem or trunk, containing woody tissue and branching into smaller branches and shoots.
A third classification level has been defined applying to the land cover class “trees” a characteristic related to the “leaf form”, allowing for distinguishing needle-leaved and broad-leaved trees.
-
Herbaceous vegetation includes annual, biennial or perennial plants such as most ferns and grasses that do not have a persistent woody stem above ground. On the opposite to woody plants, the buds of herbaceous plants die at the end of the growing season, to regenerate from the tissues left above or below the ground (e.g., bulbs, rhizomes, tubers, seeds).
A third classification level has been defined applying to the land cover class “herbaceous vegetation”, a characteristic related to the “temporal parameter” of “duration”. In this way the class is divided on the basis of the persistence of the coverage throughout the year.
Permanent herbaceous areas are characterized by a continuous vegetation cover throughout a year. No bare soil occurs within a year. These areas are either unmanaged or extensively managed natural grasslands or permanently managed grasslands, or arable areas with a permanent vegetation cover or even set-aside land in agriculture.
Periodically herbaceous areas are characterized by at least one land cover change between bare soil and herbaceous vegetation within one year. Depending on the management intensities these areas can change between these two land cover components within a year. Normally these areas are arable areas.
3.
Water surfaces class includes water in its liquid or solid state of aggregation (snow and ice), both of natural formation (water basins, rivers, streams, stagnant waters, glaciers) and of artificial origin.
-
Water bodies are water in a liquid state of aggregation regardless of location, shape, salinity and origin (natural or artificial). This includes all types of inland water surfaces without direct interference or interchange with open sea water, regardless of salinity and origin (natural or artificial). Watercourses (flowing water surfaces) come from: rivers, streams, canals, waterways. Stagnant waters (water surfaces of non-current waters, mainly lakes and ponds, or meanders of cut rivers) come from: natural lakes (both fresh and salt water), ponds, reservoirs, oxbow lakes, bottom pools, irrigation ponds, ponds for the production of artificial snow, river banks for the production of hydroelectric energy, and ponds for firefighting.
-
Permanent snow and ice includes the snow cover that persists all year round, above the climatic snow limit and the persistent ice cover formed by the accumulation of snow.
The possible land cover change classes can be identified through a transformation matrix obtained starting from the considered land cover classes (Table 3). This paper focused on the changes that can occur with a good probability within a year, between 2017 and 2018 (code 2 in Table 3). These changes refer primarily to variations of artificial abiotic surfaces in terms of expansion (i.e., land consumption), reduction (i.e., restoration), and the loss of woody vegetation associated with forest disturbances (transformation from a woody vegetation class to natural abiotic surfaces or herbaceous vegetation). Code 1 of Table 3 indicates changes that can be identified with good approximation by exploiting a series of multi-year images. The changes associated with these land cover classes have not been analysed in this paper, as the considered period is not wide enough to allow monitoring of transformations such as increases in woody vegetation increase or to distinguish seasonal variations from real changes in herbaceous classes. Finally, code 0 indicates transitions that are not within the objectives of the project.
A definition of the four land cover change classes analysed in this paper is provided below.
1.
Land consumption is the replacement of a non-artificial land cover to an artificial abiotic surface, both permanent and reversible. Artificial abiotic surfaces that have been changed by, or are under the influence of human activities resulting in a land consumption process can be sealed or unsealed.
2.
Restoration is the replacement of an artificial abiotic surface with a semi-natural land cover where permeable land is back, such as herbaceous (permanent and periodic), woody, or urban green.
3.
Forest disturbances are drastic decreases or disappearances of vegetation in areas classified as woody vegetation, not associated with land consumption. Forest disturbances are further distinguished into two classes, through the introduction of two elements relating to the status segment of the landscape characteristics of the EAGLE matrix.
-
Burnt areas: The class includes natural woody vegetation affected by recent fires. This class includes recently burnt areas of forests, moors and heathlands, sclerophyllous vegetation, transitional forest-shrub formations, and areas with sparse vegetation.
-
Other disturbances: The class “other disturbances” identifies the removal of all or most of the trees, large and small, on a surface, following a disturbance event other than a fire. Virtually all woody vegetation is removed from the site in preparation for the settlement of new wood or another class of biotic vegetation.

2.4. Pre-Processing of Images and Collection of Training Areas

Before the classification process, preliminary steps are required to collect training areas and preprocess the input images, in order to obtain the backscatter values for Sentinel-1 GRDs and reflectance values for Sentinel-2 images (Figure 1).
Sentinel-2 L1C images were used in this study. The atmospheric disturbance affecting L1C pixels can increase the uncertainty of spectral signatures, therefore decreasing classification accuracy; however, Sentinel-2 L2A historical archive was not available at the time of research. The input images were cloud masked and pixels affected by clouds excluded from calculations. The cloud mask was made through a simple algorithm based on the quality assessment band as described in Google Earth Engine website (https://developers.google.com/earth-engine/datasets/catalog/COPERNICUS_S2, accessed on 25 March 2021). It is worth highlighting that cloud mask products provided with Sentinel-2 images generally underestimate clouds [63].
Sentinel-1 GRD images provided in Google Earth Engine are already converted to backscatter values, following the steps described in the website (https://developers.google.com/earth-engine/sentinel1, accessed on 25 March 2021), which are:
-
application of orbit file;
-
removal of GRD border noise (low intensity and invalid data);
-
thermal noise removal to reduce discontinuities between sub-swaths;
-
backscatter intensity calculation using radiometric calibration;
-
terrain correction (orthorectification using the SRTM 30 m DEM);
-
conversion of backscatter coefficient to dB; due to the wide dynamic range of SAR data, decibels were used to amplify the small variations and make the data more readable [64].
The same preprocessing can be applied in the DIAS environment.
Sentinel-1 GRD and Sentinel-2 grids are not aligned and have different spatial resolution: in detail, the Sentinel-1 GRD have ground range geometry, and the pre-processing involves the geocoding to map coordinates. Sentinel-2 images are provided in WGS 84 UTM coordinates and cover the Italian territory with about 80 granules. In order to use both images in the same workflow, the spatial resampling to a common coordinate reference system and spatial resolution is required. In this study, Sentinel-1 data were aligned to the Sentinel-2 grid, in the WGS 84 UTM reference system.
Another preliminary operation is the photo-interpretation of 905 training areas. These are used for setting the thresholds in classifying woody vegetation and for the distinction between broad-leaved and needle leaved vegetation.

2.5. Classification Process

The production of the land cover map involves the calculation of multitemporal indices and decision rules that allows for the separation of land cover classes. The classification is pixel based: a pixel is assigned to a specific class if it has spectral and geometric characteristics that satisfy all the rules and thresholds defined for that class. The classification primarily regards the assignment to each pixel of one of the three classes at the first classification level (Table 1). In detail, the water surfaces pixels were first identified. In the areas not classified as water surfaces, the presence of abiotic-non vegetated areas was evaluated, and finally the remaining pixels were attributed to the class of biotic-vegetated surfaces. This operation provides three masks, which are then further classified with reference to the second or third classification level, according to the methodology described below. For further details, please refer to Appendix A.
Table 4 shows the spectral indices used for classification and calculations on the base of the Sentinel-2 images.

2.5.1. Water Surfaces

Water classes are required input for the classification of abiotic land cover classes; therefore, they must be classified before the others (Figure A1).

Permanent Snow and Ice

The classification of permanent snow and ice is based on the definition of rules starting from the NDVI index and the “snow raster”. The “snow raster” derives from the methodology developed by ESA [69] to detect snow surfaces, which are identified by pixels that have band reflectance lower than threshold values reported below:
N D S I 0.2
and
B a n d   8 0.15
and
B a n d   2 0.28
and
B a n d   2 B a n d   4 > 0.85
The algorithm is based on the NDSI index, which exploits the property of snow surfaces to be very bright in the visible and dark in shortwave infrared bands (cloud and snow have similar reflectance in band 3 but band 11 reflectance for clouds is higher than snow). This index allows for the detection of snow and differentiating it from clouds, with values varying from −1 to 1. As a general rule, positive values define areas probably with a snow cover. The NDSI is flanked by three “filters” relating to band 8 (NIR), band 2 (blue) and the ratio between band 2 and 4. The thresholds are those indicated by ESA analysis. The snow raster represents the percentage of times that all the aforementioned conditions are satisfied during a specific reporting period.
The following three conditions are defined for the identification of permanent snow and ice.
1.
Snow raster ≥ 10%
The first condition to indicate the presence of permanent snow and ice in the pixel refers to the Sentinel-2 for the whole year. It defines that snow raster conditions must be verified in the pixel for at least 10% of valid observations. The threshold of 10% allows the elimination of outlier pixels.
2.
Snow raster summer > 1%
The second condition on the snow raster refers to Sentinel-2 images for the summer period (1 July–31 September) and a less conservative threshold of 1% has been adopted for the elimination of outliers.
3.
Maximum NDVI < 0.4
For each pixel, the maximum NDVI value is calculated with respect to the annual series of valid Sentinel-2 acquisitions. A threshold of 0.4 was defined to exclude areas with the presence of vegetation.

Water Bodies

Five conditions were considered in order to identify water bodies; each of them contributes to reduce possible mistakes between the spectral signature of water and other classes. Three rules are based on indices derived from Sentinel-2 images and two used Sentinel-1 GRD data sets.
1.
High NDWI raster ≥ 5%
It was observed by [67] that water class has NDWI > 0 while non-water classes have NDWI < 0; a threshold value of 0.3 of NDWI was considered to better isolate water from the built up area and other surfaces [70]. The first condition to indicate the presence of water bodies is based on NDWI values referring to Sentinel-2 for the whole year. It defines that NDWI ≥ 0.3 conditions must be verified in the pixel for at least 5% of valid observations. The threshold of 5% allows the elimination of outlier pixels.
2.
Median NDVI < 0.3
From the NDVI multi-temporal series for the whole year, the median of the NDVI for each pixel was calculated. A threshold of 0.3 was then defined in order to eliminate areas with vegetation.
3.
Snow raster < 20%
As for the class of perennial ice and snow, a “snow raster” was elaborated, relating to the annual Sentinel-2 series. The potential presence of water bodies was referred to pixels where snow raster conditions are verified in less than 20% of the annual valid observations.
4.
Median VH polarization < −20 dB
The backscatter values used in this classification are derived from tests carried out on sample areas and from values reported in the literature [24,71,72]. The VH band discriminates quite well between the different classes presenting backscatter values around −20 dB, which remains constant throughout the year. Variations can be due to the wind, which makes the water surface rougher and modifies its backscatter. The potential presence of water bodies was referred to pixels with an annual median of the VH polarization < −20 dB.
5.
Low VV backscatter < −15%
The VV values are generally higher than VH ones. In the case of VV backscatter the value −10 dB was determined as a good threshold for identifying water [24,72]. A binary raster was produced starting from the monthly median (i.e., divided 12), attributing value 0 for values below −10 dB and 1 elsewhere. It was considered water class when the percentage of 1 throughout the year is less than 15%.

2.5.2. Abiotic Classes

The workflow for abiotic classes is based on the definition of two rules and requires only Sentinel-2 input images for the whole year. The abiotic surfaces are then distinguished between artificial and natural surfaces (Figure A2).
1.
Maximum NDVI ≤ 0.35
A threshold of 0.35 is defined on the maximum annual NDVI value, in order to exclude areas that may present vegetation during the year. This calculation is in common with the processing of water classes; therefore, the same data were calculated only once.
2.
Not Water
Both the abiotic surface class and the water surfaces have low maximum NDVI values. To avoid commission errors, this second rule is defined, which excludes pixels previously classified as water surfaces.

Artificial Abiotic Surfaces

The distinction of natural and artificial abiotic classes is based on the National Land Consumption Map elaborated by ISPRA. The map was produced through supervied classification of RapidEye images using maximum likelihood algorithms [73] and the methodology was then revised to use of Sentinel data. The map is updated every year through photointerpretation, with the help of masks of potential changes. The masks allow to reduce the area inspected by the photo interpreters, limiting the time and costs of the updating activity [62]. This paper describes the results of the updating methodology, based on the use of Sentinel-1 and Sentinel-2 data; for details on the methodology for the production of the map, see Strollo et al. [59] and ISPRA [74].

Natural Abiotic Surfaces

Artificial surfaces and natural abiotic surfaces are very similar from a spectral point of view, so it is difficult to distinguish them only through the use of satellite images. The natural abiotic surfaces class has been identified considering the pixels classified at the first level as abiotic surfaces that have not already been classified as artificial abiotic surfaces.

2.5.3. Biotic Classes

Biotic vegetation pixels are those that are not included in the classes of water and abiotic surfaces.
Woody vegetation was first identified and then divided into needle-leaved and broad-leaved trees. The land cover of the remaining areas was divided into periodically or permanent herbaceous, according to the EAGLE definition [61].
The classification of biotic surfaces required the Sentinel-1 GRD dataset for the whole year and Sentinel-2 images for the whole year, for the summer period (1 June–31 August) and for the winter period (1 January–30 April and 1–31 December). Furthermore, 905 training areas were collected throughout the national territory, for the identification of the woody vegetation and for the distinction of needle-leaved and broad-leaved trees.

Woody Vegetation

Woody vegetation was identified through the definition of three rules. The first two rules are related to Sentinel-2 optical images, while the third is based on Sentinel-1 GRD data (Figure A3).
1.
High NDVI summer ≥ 70%
The first condition to indicate the presence of woody vegetation is based on NDVI values referring to Sentinel-2 for the summer period (between 1 June to 31 August). It defines that NDVI ≥ 0.5 conditions must be verified in the pixel for at least 70% of the valid observations.
2.
Maximum NDVI summer ≥ Threshold on minimum value of the Maximum NDVI summer
The maximum NDVI value for the summer period (between 1 June to 31 August) was calculated for each pixel. This value was compared with a threshold obtained from the training area. In detail, the maximum summer NDVI was calculated in each pixel of the training areas and the lowest of them was taken as the reference threshold. Pixels with a maximum NDVI above the training area threshold are associated with the presence of woody vegetation.
3.
High vegetation
The backscatter of VH polarization for the whole year was used, as it can partially improve the detection between high vegetation (trees) and low vegetation (e.g., crops) [75]. Since the acquisition orbit influences the angle of view, the ascending and descending orbits are considered separately. A conservative threshold of −20 dB has been defined for distinguishing low vegetation and fallow land (that generally have backscatter values lower than −11 dB [76]), from trees (that generally have higher backscatter values). Pixels with backscatter values > −20 dB in at least 2% of valid observations, both for VH ascending and VH descending, were considered as woody vegetation. The 2% threshold is used to exclude possible outliers from time series of backscatter values.
Broad-Leaved and Needle-Leaved
Pixels classified as woody vegetation can be further distinguished in broad-leaved and needle-leaved using Sentinel-2 summer and winter images and training areas.
1.
Woody vegetation raster = 1
The conditions for distinguishing broad-leaved and needle-leaved are applied only to pixels classified as woody vegetation.
2.
Mean SWIR summer ≤ Threshold Mean SWIR summer (needle-leaved)
Mean SWIR summer > Threshold Mean SWIR summer (broad-leaved)
The spectral characteristics of needle-leaved and broad-leaved are very different in their wavelengths of short wave infrared (SWIR, band 11 of Sentinel-2 images). Starting from the training area of the needle-leaved pixels, the mean of the SWIR was calculated on the summer period (1 June–31 August) and the maximum value was set as a threshold. Needle-leaved pixels will have mean SWIR values for summer period lower than or equal to the threshold, while broad-leaved pixels will have higher.
3.
High NDCI winter > Threshold high NDCI winter (needle-leaved)
The third condition exploits a new empirical index developed in this research and called the “Normalized Difference Coniferous Index” (NDCI). The index is calculated on the winter subset of Sentinel-2 images (1 January–30 April and 1–31 December) and is the product of two normalized difference indexes: the first is based on band 6 (red edge, 740 nm) and band 12 (short wave infrared band, 2190 nm), and the second is based on band 8 (infrared, 842 nm) and band 11 (short wave infrared band, 1610 nm).
B 6 B 12 B 6 + B 12 B 8 B 11 B 8 + B 11
These bands show a more pronounced difference between needle-leaved and broad-leaved areas in winter than in summer. This is probably due to the drastic fall in the spectral response of broad-leaved trees compared to needle-leaved which do not drop their leaves in a senescence period [44].
Starting from the needle-leaved training areas, the NDCI was calculated and a threshold of 0.3 was set. The percentage of times that NDCI > 0.3 in the winter period was then calculated for each pixel of the needle-leaved training area. The lowest of these percentages was set as a reference threshold for the distinction of needle-leaved and broad-leaved trees. Needle-leaved pixels will have frequencies of exceeding the threshold higher than or equal to the training area one, while broad-leaved will have lower.

Herbaceous Vegetation

Herbaceous vegetation relates to pixels not previously assigned to other land cover classes (Figure A4). These areas meet the following conditions:
1.
Maximum NDVI > 0.35
As seen for the identification of abiotic surfaces, the condition distinguishes the vegetated areas from the unvegetated ones
2.
Not woody vegetation
Areas with vegetation cover but which do not respect the rules for falling into woody vegetation are considered herbaceous.
Permanent and Periodically Herbaceous
The pixels of herbaceous vegetation can be divided into permanent and periodically herbaceous, by defining conditions on the NDVI values calculated on the Sentinel-2 annual series:
1.
Low NDVI ≤ 5% (permanent herbaceous)
Low NDVI > 5% (periodically herbaceous)
Starting from the Sentinel-2 annual image series, the frequency with which NDVI < 0.35 (low NDVI) is calculated for each pixel. Pixels with NDVI values < 0.35 in less than 5% of valid observations have been classified as “permanent herbaceous”, since this condition implies a very constant presence of vegetation. Otherwise the pixel is classified as “periodically herbaceous”, since the vegetated cover was replaced by non-vegetated cover for several periods of the year. It is worth noticing that non-vegetated cover could also be snow cover.

2.6. Land Cover Changes

The algorithm allows the user to directly identify changes, through specific procedures, for the different change classes. This approach limits the propagation of errors that would occur if changes were sought by comparing two complete land cover classifications [77]. The identification of changes of artificial abiotic surfaces involved the creation of masks of potential changes, used for the annual update of the National Land Consumption Map. In detail, Sentinel-1 and Sentinel-2 images for 2017 and 2018 were used. For Sentinel-1 data, thresholds on ascending and descending backscatter variations were defined, while for Sentinel-2 data thresholds regard the maximum NDVI value in the second year and the NDVI difference between the two years. Areas with high backscatter variations and/or strong reductions in NDVI between the two years were included in a raster of potential changes. The raster was used as a basis for the photo-interpretation activities, carried out for the updating of the National Land Consumption Map. The details of this study are described in another paper [58].
A specific methodology for mapping forest disturbance was developed, which also allows distinguishing burned areas from other kinds of changes such as forest harvesting (Figure A5). The methodology is applied to areas classified as woody vegetation in the land cover map. This class acts as a mask for the delimitation of the forest area of year 1, within which the search for changes takes place. The input data used are the woody vegetation class for year 1 and the Sentinel-2 summer images (from 1 June to 31 August) for the two reference years. In this research, the changes that took place between 2017 and 2018 were analysed. Two spectral indices were calculated for the two reference years: NDVI and NBR [78] (Table 4).
NBR is useful for identifying areas that have been burned in the reference period. Notably, healthy vegetation has very high near-infrared reflectance and low reflectance in the shortwave infrared portion of the spectrum. Recently burned areas have a low reflectance in near infrared and a high reflectance in the short-wave infrared band.
NDVI is high in vegetated areas, while it decreases if there are reductions in vegetation cover.
To identify a forest disturbance, at least one of the following two conditions must be verified.
1.
NDVI Difference ≤ −0.2
The median of the summer NDVI is calculated for the two reference years and then the difference between the two values is considered:
N D V I   d i f f e r e n c e = M e d i a n   N D V I   y e a r   1     M e d i a n   N D V I   y e a r   2
If the difference is greater than 0.2, the change in NDVI indicates the presence of a change.
2.
NBR Difference ≤ −0.2
The difference between the median of the summer NBR of the two reference years is considered:
N B R   d i f f e r e n c e = M e d i a n   N B R   y e a r   1     M e d i a n   N B R   y e a r   2
The threshold value 0.2 was identified empirically as optimal value for detecting tree loss and avoiding errors caused by phenological variations between years.

Distinction between “Burnt Areas” and “Other Distubances”

Starting from the pixels classified as forest disturbance, a distinction was made between “burned trees” and “other forest changes”, through an innovative burned index (BI) (Table 4). The BI is higher for dark or brown pixels, as is in case of burned areas.
The median of BI of year 2 was calculated, obtaining a raster of “Median Burned Index Year 2”. Pixels already classified as forest disturbance can be further detailed in fires if Median BI Year 2 ≥ 0.45
The threshold value 0.45 was determined empirically as the optimal value for distinguishing burned areas and avoiding errors due to shadow areas where pixels tend to be dark.

2.7. Accuracy Assessment

The validation involved the 2018 land cover map with an eight-class classification system, integrated with the four land cover change classes.
First, a qualitative verification was carried out through a systematic visual search for macroscopic errors. An accuracy assessment on the base of Olofsson methodology was then performed [79]; this approach is widely used in the literature [79,80,81,82].
The validation was carried out by identifying, ad hoc, a random sample of points which was then photo-interpreted using very high resolution images. The points were compared with the values of the land cover map at the same locations. The sample size (n) was calculated using the equation reported by Olofsson [79] and referred to 5.25 in [83].
W i S i 2 [ S O ^ ] 2
where:
Wi—is area proportion of each classes derived from the map classification
Si—standard deviation of stratum i, Si = √(Ui(1−Ui)) ([83], Eq. 5.55)
Ui—User accuracy of class i
S(Ô)—is the standard target standard error
Si influences the sample size, and its value is related to the Ui, which is unknown. A conservative scenario was assumed, considering Ui equal to 0.6 for all twelve classes.
The target standard error for overall accuracy was assumed to be 0.01 as suggested by Olofsson [79], which corresponds to a confidence interval of 1%.
The application of (2) gives a sample of size 2400 (Table 5).
Sample size calculation according to Olofsson [79]. The required inputs are the areas of each and the expected user’s accuracy for the sample. The points were divided among the land cover classes through a stratified sampling. The number of points attributed to each class was assumed to be equal to the average between equal and area-proportional distribution (see https://fromgistors.blogspot.com/2019/09/Accuracy-Assessment-of-Land-Cover-Classification.html accessed on 4 April 2021) (Table 6).
The points were photo-interpreted with very high resolution images from Google Earth, considering 2018 as the reference year for the land cover classes and 2017 and 2018 for the land cover change classes. The photo-interpretation followed a pixel based approach, assigning to each point the class on which it falls exactly.

3. Results

The land cover map shown in Figure 3 was obtained from the application of the methodology. The map shows the eight land cover classes for 2018 and the four land cover change classes for 2017–2018 described above.
Accuracy assessment was done on the map of Figure 3 and provides the results shown in Figure 4 and in Table 7.
The map has 83% overall accuracy (Table 7). This value offers a general estimation of the accuracy, but it is influenced by the distribution of the error between the classes, especially those of land cover change. This circumstance also affects the spatial distribution of the classification errors detected by the sample points of Figure 4, since the land cover change classes do not have a specific distribution within the national territory.
Water bodies (98%), abiotic artificial surfaces (92%) and permanent herbaceous (92%) have the best user accuracy, while the land cover change classes show the lowest values, except for “land consumption” which presents a user accuracy of 81%.
As regarding the producer’s accuracy, four classes have values of 100%, and 9 out of 12 greater than 81%. The lowest values regard the land cover change class of “other disturbances” (71%) and the land cover class of permanent herbaceous (62%). Analysing the land cover data for 2018 and land cover change that occurred between 2017 and 2018 (Table 8), 10.2% of the Italian territory is represented by abiotic surfaces, of which about 70% are artificial surfaces. The rest of the territory is divided between woody vegetation (45.2%, with a prevalence of broad-leaved trees) and herbaceous vegetation (42.8%, especially periodic herbaceous). The remaining 1.49% is represented by water bodies (1.1%) and perennial ice and snow (0.4%).
The four classes of change occupy 0.3% of the national territory and about 90% are forest disturbances.
Table 9 shows that between 2017 and 2018, land consumption (4231 hectares, i.e., 64.1% of total land consumption) and restoration (636 hectares, equal to 70.2% of restorations) occurred mainly in the periodically herbaceous class. The artificial abiotic surface is the class with more changes (both land consumption and restoration) compared to the class extension.
The decrease in vegetation caused by forest disturbances affects 89,596 hectares (Table 10). 12.4% (11,106 ha) of forest disturbances are burnt areas, while 87.6% (78,490 ha) had other causes.
The “other disturbances” class was further analysed, showing that more than half of the disturbances are attributable to forest harvesting.

4. Discussion

The high overall accuracy of the map supports the rule-based classification method chosen, which allows the classification of the entire national territory with a relatively small number of training samples and frequent and fast updates of the map. Maxwell et al. and Thanh Noi et al. [84,85] assessed that the machine learning classification method achieves a high accuracy level; however, these methods are usually applied to limited areas, at a regional or local scale [86,87], or to a few land cover classes [45,88]. Moreover, as illustrated by Zhang et al. [89], machine learning approaches such as random forest need to be trained with many samples, which must be balanced in order to avoid misclassification errors [90], making the map update process slower.
The developed methodology allowed for mapping the land cover and the changes in Italy, although with different results depending on the land cover class. The artificial abiotic surfaces and the changes related to land consumption were mapped using the National Land Consumption Map produced and updated by ISPRA. The use of ancillary data allowed for overcoming the issue related to spectral similarity between the natural abiotic surfaces and the artificial ones using Sentinel-2 images (Figure 5). The data have high accuracy, due to periodic updates and improvements through photo-interpretation (Figure 6). Hyperspectral images (such as PRISMA) could be used to improve the distinction of abiotic surfaces exploiting the bands in the SWIR region [91].
The class of permanent water bodies has an accuracy of over 90%, as the characteristics of the class can be well discriminated through the optical and radar instruments used in this research (Figure 7). The few omission errors regard boundary areas between water bodies and other land cover classes, such as river and lake banks, and some wetlands.
The class of permanent snow and ice has high accuracy (Figure 8), although it is slightly lower than permanent water bodies (86%). The accuracy of the class is affected by commission errors related to the seasonal variation of snow cover.
The classes of woody vegetation correctly distinguished between needle-leaved and broad-leaved trees in areas with homogeneous coverage, while in mixed forest pixels there are classification errors due to the spatial resolution of the data, causing mixed spectral signatures, and to the spectral similarity between the different types of trees (Figure 9) [6,7,92,93].
The results obtained for the herbaceous cover are related with the definition of the classes of periodic (Figure 10) and permanent (Figure 11) herbaceous, which consider the persistence of the herbaceous cover in the pixel throughout the year. In detail, pixels showing herbaceous cover in at least 95% of valid observations were classified as permanent herbaceous. This threshold is conservative and aligned with the definitions of CLC+ products [61], but it could be adapted to different temporal definitions.
Forest disturbances (Figure 12 and Figure 13) are mainly affected by commission errors, due to the adoption of wide thresholds for the identification of the classes. Commissions are mainly related with permanent crops, since they have a spectral behavior similar to forest disturbances during the year.
The mapping of burned areas shows results aligned with the European Forest Fire Information System (EFFIS) (https://effis.jrc.ec.europa.eu/, accessed on 25 February 2021).
The user accuracy of the class “other disturbances” is also low because it has been validated assimilating the class to “forest harvesting” only. The choice was made to evaluate the algorithm’s effectiveness in identifying this type of disturbance. Recently, several studies refer to the use of Sentinel data for monitoring forest disturbances, even if for limited areas [94]. The size of the study area conditioned the identification of changes, since cuttings are conducted with different techniques which vary according to the geographical area; therefore, in some cases it is not possible to clearly associate the spectral variations to a certain type of disturbance. The detection of forest harvesting is more effective in Central Italy, while there are classification errors where selective logging method is practiced, such as on the Alps.
The observation of the areas of potential disturbance over a longer period can allow the improvement of the data [60,95].
The main objective of this research regards the formulation of a methodology that would allow the mapping and monitoring of land cover for the development of operational services, in compliance with the project requirements in terms of resolution, study area extension and update frequency. This paper presents a first product for 2018 that allowed the achievement of the project objectives and that actually lends itself to future refinements and improvements.
The EAGLE-based classification system allows for the introduction of new classes, such as shrubland, and can be modified for specific purposes, for example with land use attributes to identify agricultural activities in the herbaceous vegetation classes. The identification of shrubland will make possible a better description of Mediterranean macchia, permanent crops such as vineyards, transition areas and woodland boundaries. In this sense, it is necessary to conduct an in-depth analysis of the spectral behavior of shrubs, which are similar to scattered trees or to a mixed herbaceous and arboreal cover.
The use of Sentinel data constitutes an added value of the method, as they offer excellent results in the study of land cover at high spatial and temporal resolution, achieving high accuracy level with various classification algorithms [57]. Furthermore, the free access policy also allows its use in contexts where financial resources for the acquisition of remotely sensed data are limited. Of course, additional classes would require additional methodologies, in the implementation of which the integration of hyperspectral data will be fundamental. The new HyperSpectral Precursor of the Application Mission (PRISMA) sensor, developed by the Italian Space Agency, has shown encouraging results in the distinction of forest types [96]. This technology will be particularly useful with the launch of the “Hyperspectral Imaging Mission for the Environment” (CHIME) [97], whose data will be freely accessible in the coming years.
A characteristic of the proposed method is the definition of thresholds and a set of rules which can be modified and adapted to different conditions. The use of a decision rules classification method has made it possible to develop ad-hoc rules for each land cover class, while the definition of fixed thresholds for the entire national territory makes the method easily controllable and implementable. On the other hand, the Italian territory presents a strong heterogeneity from the morphological and climatic points of view that cannot be neglected. In this sense, an important development of the methodology will regard the adaptation of the thresholds according to the characteristics of the territory, considering in the first place the influence of hydrometeorological variables and seasonal climatic variations on the indices. This topic is addressed by the working group in the publication by Spadoni et al. [18], which involved a study area consisting of three Sentinel-2 granules and will be studied in depth. In this way it is possible to better analyse the phenological evolution of vegetation in order to improve the distinction between needle-leaved and broad-leaved trees and the identification of the different types of herbaceous classes and forest disturbances.
A further improvement regards the use of elevation models. At the moment this data has been used in the monitoring of changes associated with land consumption described by Luti et al. [58] and will be able to guarantee important improvements in the land cover classification of mountain areas, especially in the shaded slopes of Alpine areas. In these areas it will be useful also to analyse the spatial correlation between pixels. This deepening can address the interpolation of missing or incomplete data [98], enhancing some critical areas and allowing users to overcome classification errors due to input data characteristics (such as image shadows, or cloud coverage).
The versatility of the algorithm and of the classification system make the methodology suitable for future integrations and for adaptation for specific purposes, while remaining consistent with the requirements proposed at the national and European level.

5. Conclusions

The research has allowed the creation of a land cover classification methodology capable of providing nationwide products with high spatial resolution and accuracy. The methodology shows satisfactory results, allowing the achievement of all design requirements. It was conceived with the aim of obtaining products that can be used for operational purposes, with a high update frequency. Indeed, they fit well into the Mirror Copernicus initiative of the Space Economy Strategic Plan of the Italian Ministry of Economic Development, which provides for the creation of frequently updated products at a national scale, for the monitoring of the territory.
The approach based on decision rules guarantees the control of the choices adopted for the classification, and a high level of modularity. In this sense, the rules can be refined by inserting variable thresholds according to the characteristics of the area, such as latitude and phytoclimatic zone, hydrometeorological variables or elevation, but also they can be accompanied by other rules for the identification of additional classes. In this sense, the classification system based on EAGLE has an easily expandable structure, remaining consistent with the European indications in terms of classification systems. The methodology can be implemented quickly and easily and has a low computational cost, allowing the production of land cover maps that can be updated annually and with high spatial resolution.
The maps produced by the application of the methodology can also be used as ancillary data for the development of algorithms for monitoring and characterizing specific coverage classes, such as the detailed distinction of forest types, the characterization of agricultural areas or a more detailed analysis of forest disturbances. The methodology is therefore promising, thanks to the simplicity of the approach on which it is based and the ability to be applied to large areas while maintaining high spatial resolution and accuracy levels. In this sense, it is unique in the panorama of European and national land cover monitoring products.

Author Contributions

Conceptualization, P.D.F., L.C., T.L. and M.M. (Michele Munafò); methodology, P.D.F., L.C. and T.L.; software, P.D.F., L.C., P.D. and T.L.; validation, P.D.F., T.L., P.D., C.G. and A.C.; formal analysis, P.D.F. and L.C.; investigation, P.D.F., L.C., T.L. and A.C.; resources, M.M. (Michele Munafò), M.M. (Marco Marchetti) and G.C.; data curation, P.D.F. and L.C.; writing—original draft preparation, P.D.F. and T.L.; visualization, P.D.F., A.C., P.D. and C.G.; supervision, M.M. (Michele Munafò), M.M. (Marco Marchetti) and G.C.; funding acquisition, M.M. (Michele Munafò); All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Italian Institute for Environmental Protection and Research (ISPRA) and the Italian Space Agency (ASI) agreement for the Habitat Mapping project (Agr. number F82F17000010005).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data presented in this study are available on request from the corresponding author. The data are not publicly available because they are part of ongoing research.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

This appendix describes the workflows of the classification methodology for each land cover class, organized in macro-class of water surfaces, abiotic surfaces and biotic surfaces (distinguished between woody vegetation and herbaceous vegetation).
Each workflow is structured in the following sections:
-
“DATA” shows the input data used for the classification.
-
“PRE-PROCESSING” illustrates the calculations of spectral indices and selected backscatter polarizations.
-
“BINARY RASTER” reports the conditions and thresholds defined for the various inputs in order to create binary rasters (whose values are 0 or 1). A binary raster is calculated by attributing to each pixel the value 1 if the conditions defined by the threshold are verified and 0 elsewhere. A binary raster is calculated for each valid image acquisition (i.e., the number of valid acquisitions can be lower or equal to the total number of satellite acquisitions, depending on cloud cover).
-
“VARIABLES” include the raster calculations derived from preprocessed data and binary rasters, which summarize the multi-temporal properties in a unique raster using statistics such as the median, or the percentage of occurrence calculated as:
P C l a s s =   B R N c l o u d   100
where B R is the summation of binary rasters and N c l o u d is the total number of valid acquisitions.
-
The conditions for the classification of every land cover class are described in the lower boxes, which are based on the above “VARIABLES” and thresholds that are fixed or calculated from training input data.
Figure A1. Workflow for water surfaces classification.
Figure A1. Workflow for water surfaces classification.
Land 10 00611 g0a1
Figure A2. Workflow for abiotic surfaces classification.
Figure A2. Workflow for abiotic surfaces classification.
Land 10 00611 g0a2
Figure A3. Workflow for woody vegetation classification.
Figure A3. Workflow for woody vegetation classification.
Land 10 00611 g0a3
Figure A4. Workflow for herbaceous vegetation classification.
Figure A4. Workflow for herbaceous vegetation classification.
Land 10 00611 g0a4
Figure A5. Workflow for forest disturbance detection.
Figure A5. Workflow for forest disturbance detection.
Land 10 00611 g0a5

References

  1. Sallustio, L.; Munafò, M.; Riitano, N.; Lasserre, B.; Fattorini, L.; Marchetti, M. Integration of land use and land cover inventories for landscape management and planning in Italy. Environ. Monit. Assess. 2016, 188, 1–20. [Google Scholar] [CrossRef] [PubMed]
  2. European Parliament, Council of the E.U. Directive 2007/2/EC of the European Parliament and of the Council of 14 March 2007 establishing an Infrastructure for Spatial Information in the European Community (INSPIRE). Off. J. Eur. Union 2007, 1–14.
  3. Gómez, C.; White, J.C.; Wulder, M.A. Optical remotely sensed time series data for land cover classification: A review. ISPRS J. Photogramm. Remote Sens. 2016, 116, 55–72. [Google Scholar] [CrossRef] [Green Version]
  4. Steinhausen, M.J.; Wagner, P.D.; Narasimhan, B.; Waske, B. Combining Sentinel-1 and Sentinel-2 data for improved land use and land cover mapping of monsoon regions. Int. J. Appl. Earth Obs. Geoinf. 2018, 73, 595–604. [Google Scholar] [CrossRef]
  5. EEA. The European Environment—State and Outlook 2020: Knowledge for Transition to a Sustainable Europe; European Environment Agency: Copenhagen, Denmark, 2019; pp. 115–131. [Google Scholar]
  6. Del Río-Mena, T.; Willemen, L.; Vrieling, A.; Nelson, A. Understanding intra-annual dynamics of ecosystem services using satellite image time series. Remote Sens. 2020, 12, 710. [Google Scholar] [CrossRef] [Green Version]
  7. Kuenzer, C.; Ottinger, M.; Wegmann, M.; Guo, H.; Wang, C.; Zhang, J.; Dech, S.; Wikelski, M. Earth observation satellite sensors for biodiversity monitoring: Potentials and bottlenecks. Int. J. Remote Sens. 2014, 35, 6599–6647. [Google Scholar] [CrossRef] [Green Version]
  8. Reba, M.; Seto, K.C. A systematic review and assessment of algorithms to detect, characterize, and monitor urban land change. Remote Sens. Environ. 2020, 242, 111739. [Google Scholar] [CrossRef]
  9. UNCCD. The Global Land Outlook, 1st ed.; United Nations Convention to Combat Desertification: Bonn, Germany, 2017; pp. 1–336. ISBN 9789295110489. [Google Scholar]
  10. Wunder, S.; Kaphengst, T.; Frelih-Larsen, A. Implementing Land Degradation Neutrality (SDG 15.3) at National Level: General Approach, Indicator Selection and Experiences from Germany. In International Yearbook of Soil Law and Policy 2017; Ginzky, H., Dooley, E., Heuser, I.L., Kasimbazi, E., Markus, T., Qin, T., Eds.; Springer International Publishing: Cham, Germany, 2018; pp. 191–219. ISBN 978-3-319-68885-5. [Google Scholar]
  11. Berger, M.; Aschbacher, J. Preface: The Sentinel missions-new opportunities for science. Remote Sens. Environ. 2012, 120, 1–2. [Google Scholar] [CrossRef]
  12. Torres, R.; Snoeij, P.; Geudtner, D.; Bibby, D.; Davidson, M.; Attema, E.; Potin, P.; Rommen, B.Ö.; Floury, N.; Brown, M.; et al. GMES Sentinel-1 mission. Remote Sens. Environ. 2012, 120, 9–24. [Google Scholar] [CrossRef]
  13. Chatziantoniou, A.; Petropoulos, G.P.; Psomiadis, E. Co-Orbital Sentinel 1 and 2 for LULC mapping with emphasis on wetlands in a mediterranean setting based on machine learning. Remote Sens. 2017, 9, 1259. [Google Scholar] [CrossRef] [Green Version]
  14. Bhattarai, R.; Rahimzadeh-Bajgiran, P.; Weiskittel, A.; Meneghini, A.; MacLean, D.A. Spruce budworm tree host species distribution and abundance mapping using multi-temporal Sentinel-1 and Sentinel-2 satellite imagery. ISPRS J. Photogramm. Remote Sens. 2021, 172, 28–40. [Google Scholar] [CrossRef]
  15. Nezry, E. Adaptive Speckle Filtering in Radar Imagery. In Land Applications of Radar Remote Sensing; IntechOpen: London, UK, 2014; pp. 1–55. [Google Scholar]
  16. Muro, J.; Varea, A.; Strauch, A.; Guelmami, A.; Fitoka, E.; Thonfeld, F.; Diekkrüger, B.; Waske, B. Multitemporal optical and radar metrics for wetland mapping at national level in Albania. Heliyon 2020, 6, e04496. [Google Scholar] [CrossRef]
  17. Drusch, M.; Del Bello, U.; Carlier, S.; Colin, O.; Fernandez, V.; Gascon, F.; Hoersch, B.; Isola, C.; Laberinti, P.; Martimort, P.; et al. Sentinel-2: ESA’s Optical High-Resolution Mission for GMES Operational Services. Remote Sens. Environ. 2012, 120, 25–36. [Google Scholar] [CrossRef]
  18. Spadoni, G.L.; Cavalli, A.; Congedo, L.; Munafò, M. Analysis of Normalized Difference Vegetation Index (NDVI) multi-temporal series for the production of forest cartography. Remote Sens. Appl. Soc. Environ. 2020, 20, 100419. [Google Scholar] [CrossRef]
  19. Teodoro, A.; Amaral, A. A statistical and spatial analysis of portuguese forest fires in summer 2016 considering landsat 8 and sentinel 2A data. Environments 2019, 6, 36. [Google Scholar] [CrossRef] [Green Version]
  20. Malenovský, Z.; Rott, H.; Cihlar, J.; Schaepman, M.E.; García-Santos, G.; Fernandes, R.; Berger, M. Sentinels for science: Potential of Sentinel-1, -2, and -3 missions for scientific observations of ocean, cryosphere, and land. Remote Sens. Environ. 2012, 120, 91–101. [Google Scholar] [CrossRef]
  21. Ienco, D.; Interdonato, R.; Gaetano, R.; Ho Tong Minh, D. Combining Sentinel-1 and Sentinel-2 Satellite Image Time Series for land cover mapping via a multi-source deep learning architecture. ISPRS J. Photogramm. Remote Sens. 2019, 158, 11–22. [Google Scholar] [CrossRef]
  22. Erinjery, J.J.; Singh, M.; Kent, R. Mapping and assessment of vegetation types in the tropical rainforests of the Western Ghats using multispectral Sentinel-2 and SAR Sentinel-1 satellite imagery. Remote Sens. Environ. 2018, 216, 345–354. [Google Scholar] [CrossRef]
  23. Joshi, N.; Baumann, M.; Ehammer, A.; Fensholt, R.; Grogan, K.; Hostert, P.; Jepsen, M.R.; Kuemmerle, T.; Meyfroidt, P.; Mitchard, E.T.A.; et al. A review of the application of optical and radar remote sensing data fusion to land use mapping and monitoring. Remote Sens. 2016, 8, 70. [Google Scholar] [CrossRef] [Green Version]
  24. Mansaray, L.R.; Huang, W.; Zhang, D.; Huang, J.; Li, J. Mapping rice fields in urban Shanghai, southeast China, using Sentinel-1A and Landsat 8 datasets. Remote Sens. 2017, 9, 257. [Google Scholar] [CrossRef] [Green Version]
  25. Clerici, N.; Valbuena Calderón, C.A.; Posada, J.M. Fusion of sentinel-1a and sentinel-2A data for land cover mapping: A case study in the lower Magdalena region, Colombia. J. Maps 2017, 13, 718–726. [Google Scholar] [CrossRef] [Green Version]
  26. Dusseux, P.; Corpetti, T.; Hubert-Moy, L.; Corgne, S. Combined use of multi-temporal optical and Radar satellite images for grassland monitoring. Remote Sens. 2014, 6, 6163–6182. [Google Scholar] [CrossRef] [Green Version]
  27. Iannelli, G.C.; Gamba, P. Jointly exploiting Sentinel-1 and Sentinel-2 for urban mapping. In Proceedings of the IGARSS 2018—2018 IEEE International Geoscience and Remote Sensing Symposium, Valencia, Spain, 22–27 July 2018; pp. 8209–8212. [Google Scholar]
  28. Colson, D.; Petropoulos, G.P.; Ferentinos, K.P. Exploring the Potential of Sentinels-1 & 2 of the Copernicus Mission in Support of Rapid and Cost-effective Wildfire Assessment. Int. J. Appl. Earth Obs. Geoinf. 2018, 73, 262–276. [Google Scholar]
  29. Duong, P.C.; Trung, T.H.; Nasahara, K.N.; Tadono, T. JAXA high-resolution land use/land cover map for Central Vietnam in 2007 and 2017. Remote Sens. 2018, 10, 1406. [Google Scholar] [CrossRef] [Green Version]
  30. Anderson, J.R. Land use and land cover changes-A framework for monitoring. J. Res. U. S. Geol. Surv. 1976, 5, 142–152. [Google Scholar]
  31. Arnold, S.; Kosztra, B.; Banko, G.; Smith, G.; Hazeu, G.; Bock, M. The EAGLE concept—A vision of a future European Land Monitoring Framework. In Proceedings of the 33rd EARSeL Symposium “Towards Horizon 2020”, Matera, Italy, 3–6 June 2013; pp. 551–568. [Google Scholar]
  32. Arnold, S.; Kosztra, B.; Banko, G.; Milenov, P.; Smith, G.; Hazeu, G. Explanatory Content Documentation of the EAGLE Concept 2021, Version 3.1. Available online: https://land.copernicus.eu/eagle/content-documentation-of-the-eagle-concept/manual/eagle-explanatory-documentation-v3-1-version-2021 (accessed on 4 March 2021).
  33. Chen, Y.; Gong, P. Clustering based on eigenspace transformation—CBEST for efficient classification. ISPRS J. Photogramm. Remote Sens. 2013, 83, 64–80. [Google Scholar] [CrossRef]
  34. Zhu, Z.; Gallant, A.L.; Woodcock, C.E.; Pengra, B.; Olofsson, P.; Loveland, T.R.; Jin, S.; Dahal, D.; Yang, L.; Auch, R.F. Optimizing selection of training and auxiliary data for operational land cover classification for the LCMAP initiative. ISPRS J. Photogramm. Remote Sens. 2016, 122, 206–221. [Google Scholar] [CrossRef] [Green Version]
  35. Holloway, J.; Helmstedt, K.J.; Mengersen, K.; Schmidt, M. A decision tree approach for spatially interpolating missing land cover data and classifying satellite images. Remote Sens. 2019, 11, 1796. [Google Scholar] [CrossRef] [Green Version]
  36. Thanh Noi, P.; Kappas, M. Comparison of random forest, k-nearest neighbor, and support vector machine classifiers for land cover classification using Sentinel-2 imagery. Sensors 2018, 18, 18. [Google Scholar] [CrossRef] [Green Version]
  37. Colditz, R.R. An evaluation of different training sample allocation schemes for discrete and continuous land cover classification using decision tree-based algorithms. Remote Sens. 2015, 7, 9655–9681. [Google Scholar] [CrossRef] [Green Version]
  38. Liu, K.; Shi, W.; Zhang, H. A fuzzy topology-based maximum likelihood classification. ISPRS J. Photogramm. Remote Sens. 2011, 66, 103–114. [Google Scholar] [CrossRef]
  39. Xin, Q.; Li, J.; Li, Z.; Li, Y.; Zhou, X. Evaluations and comparisons of rule-based and machine-learning-based methods to retrieve satellite-based vegetation phenology using MODIS and USA National Phenology Network data. Int. J. Appl. Earth Obs. Geoinf. 2020, 93, 102189. [Google Scholar] [CrossRef]
  40. Huang, X.; Liu, J.; Zhu, W.; Atzberger, C.; Liu, Q. The optimal threshold and vegetation index time series for retrieving crop phenology based on a modified dynamic threshold method. Remote Sens. 2019, 11, 2725. [Google Scholar] [CrossRef] [Green Version]
  41. Otukei, J.R.; Blaschke, T. Land cover change assessment using decision trees, support vector machines and maximum likelihood classification algorithms. Int. J. Appl. Earth Obs. Geoinf. 2010, 12, 27–31. [Google Scholar] [CrossRef]
  42. Kobayashi, N.; Tani, H.; Wang, X.; Sonobe, R. Crop classification using spectral indices derived from Sentinel-2A imagery. J. Inf. Telecommun. 2020, 4, 67–90. [Google Scholar] [CrossRef]
  43. Sinha, P.; Kumar, L. Independent two-step thresholding of binary images in inter-annual land cover change/no-change identification. ISPRS J. Photogramm. Remote Sens. 2013, 81, 31–43. [Google Scholar] [CrossRef]
  44. Persson, M.; Lindberg, E.; Reese, H. Tree species classification with multi-temporal Sentinel-2 data. Remote Sens. 2018, 10, 1794. [Google Scholar] [CrossRef] [Green Version]
  45. Agrillo, E.; Filipponi, F.; Pezzarossa, A.; Casella, L.; Smiraglia, D.; Orasi, A.; Attorre, F.; Taramelli, A. Earth observation and biodiversity big data for forest habitat types classification and mapping. Remote Sens. 2021, 13, 1231. [Google Scholar] [CrossRef]
  46. Ngadze, F.; Mpakairi, K.S.; Kavhu, B.; Ndaimani, H.; Maremba, M.S. Exploring the utility of Sentinel-2 MSI and Landsat 8 OLI in burned area mapping for a heterogenous savannah landscape. PLoS ONE 2020, 15, e0232962. [Google Scholar] [CrossRef]
  47. 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. [Google Scholar] [CrossRef]
  48. Lary, D.J.; Alavi, A.H.; Gandomi, A.H.; Walker, A.L. Machine learning in geosciences and remote sensing. Geosci. Front. 2016, 7, 3–10. [Google Scholar] [CrossRef] [Green Version]
  49. Maxwell, A.E.; Warner, T.A.; Vanderbilt, B.C.; Ramezan, C.A. Land cover classification and feature extraction from national agriculture imagery program (NAIP) Orthoimagery: A Review. Photogramm. Eng. Remote Sens. 2017, 83, 737–747. [Google Scholar] [CrossRef]
  50. Foody, G.M.; Mathur, A. The use of small training sets containing mixed pixels for accurate hard image classification: Training on mixed spectral responses for classification by a SVM. Remote Sens. Environ. 2006, 103, 179–189. [Google Scholar] [CrossRef]
  51. Kulkarni, A.D.; Lowe, B. Random Forest Algorithm for Land Cover Classification. Pattern Recognit. Lett. 2016, 27, 294–300. [Google Scholar]
  52. Liu, Y.; Gong, W.; Hu, X.; Gong, J. Forest type identification with random forest using Sentinel-1A, Sentinel-2A, multi-temporal Landsat-8 and DEM data. Remote Sens. 2018, 10, 946. [Google Scholar] [CrossRef] [Green Version]
  53. Rodriguez-Galiano, V.F.; Chica-Olmo, M.; Abarca-Hernandez, F.; Atkinson, P.M.; Jeganathan, C. Random Forest classification of Mediterranean land cover using multi-seasonal imagery and multi-seasonal texture. Remote Sens. Environ. 2012, 121, 93–107. [Google Scholar] [CrossRef]
  54. Sun, Z.; Xu, R.; Du, W.; Wang, L.; Lu, D. High-resolution urban land mapping in China from sentinel 1A/2 imagery based on Google Earth Engine. Remote Sens. 2019, 11, 752. [Google Scholar] [CrossRef] [Green Version]
  55. Houhoulis, P.F.; Michener, W.K. Detecting wetland change: A rule-based approach using NWI and SPOT-XS data. Photogramm. Eng. Remote Sens. 2000, 66, 205–211. [Google Scholar]
  56. Berhane, T.M.; Lane, C.R.; Wu, Q.; Autrey, B.C.; Anenkhonov, O.A.; Chepinoga, V.V.; Liu, H. Decision-tree, rule-based, and random forest classification of high-resolution multispectral imagery for wetland mapping and inventory. Remote Sens. 2018, 10, 580. [Google Scholar] [CrossRef] [Green Version]
  57. Phiri, D.; Simwanda, M.; Salekin, S.; Nyirenda, V.R.; Murayama, Y.; Ranagalage, M. Sentinel-2 Data for Land Cover/Use Mapping: A Review. Remote Sens. 2020, 12, 2291. [Google Scholar] [CrossRef]
  58. Luti, T.; De Fioravante, P.; Marinosci, I.; Strollo, A.; Riitano, N.; Falanga, V.; Mariani, L.; Congedo, L.; Munafò, M. Land Consumption Monitoring with SAR Data and Multispectral Indices. Remote Sens. 2021, 13, 1586. [Google Scholar] [CrossRef]
  59. Strollo, A.; Smiraglia, D.; Bruno, R.; Assennato, F.; Congedo, L.; De Fioravante, P.; Giuliani, C.; Marinosci, I.; Riitano, N.; Munafò, M. Land consumption in Italy. J. Maps 2020, 16, 113–123. [Google Scholar] [CrossRef]
  60. Giannetti, F.; Pegna, R.; Francini, S.; McRoberts, R.E.; Travaglini, D.; Marchetti, M.; Scarascia Mugnozza, G.; Chirici, G. A New Method for Automated Clearcut Disturbance Detection in Mediterranean Coppice Forests Using Landsat Time Series. Remote Sens. 2020, 12, 3720. [Google Scholar] [CrossRef]
  61. Kleeschulte, S.; Banko, G.; Smith, G.; Arnold, S.; Scholz, J.; Kosztra, B.; Maucha, G. Technical Specifications for Implementation of a New Land-Monitoring Concept Based on EAGLE, D5: Design Concept and CLC+ Backbone, Technical Specifications, CLC+ Core and CLC+ Instances Draft Specifications, Including Requirements Review; Version 5.4; European Environment Agency: Copenhagen, Denmark, 2020; pp. 1–78. [Google Scholar]
  62. Munafò, M. Consumo di Suolo, Dinamiche Territoriali e Servizi Ecosistemici; Edizione 2020; SNPA: Rome, Italy, 2020; p. 291. ISBN 9788844810139. [Google Scholar]
  63. Coluzzi, R.; Imbrenda, V.; Lanfredi, M.; Simoniello, T. A first assessment of the Sentinel-2 Level 1-C cloud mask product to support informed surface analyses. Remote Sens. Environ. 2018, 217, 426–443. [Google Scholar] [CrossRef]
  64. Flores, A.; Herndon, K.; Thapa, R.; Cherrington, E. (Eds.) The SAR Handbook: Comprehensive Methodologies for Forest Monitoring and Biomass Estimation; NASA: Washington, DC, USA, 2019. [Google Scholar]
  65. Rouse, J.W., Jr.; Deering, D.W.; Schell, J.A.; Harlan, J.C. Monitoring the Vernal Advancement and Retrogradation (Green Wave Effect) of Natural Vegetation; NASA: Greenbelt, MD, USA, 1974. [Google Scholar]
  66. Key, C.H.; Benson, N.C. The Normalized Burn Ratio (NBR): A Landsat TM Radiometric Measure of Burn Severity; USA Geological Survey Northern Rocky Mountain Science Center: Bozeman, MT, USA, 1999.
  67. McFeeters, S.K. The use of the Normalized Difference Water Index (NDWI) in the delineation of open water features. Int. J. Remote Sens. 1996, 17, 1425–1432. [Google Scholar] [CrossRef]
  68. Dozier, J. Spectral signature of alpine snow cover from the Landsat Thematic Mapper. Remote Sens. Environ. 1989, 28, 9–22. [Google Scholar] [CrossRef]
  69. ESA Level-2A Algorithm Overview. Available online: https://sentinels.copernicus.eu/web/sentinel/technical-guides/sentinel-2-msi/level-2a/algorithm (accessed on 4 March 2021).
  70. McFeeters, S.K. Using the normalized difference water index (ndwi) within a geographic information system to detect swimming pools for mosquito abatement: A practical approach. Remote Sens. 2013, 5, 3544–3561. [Google Scholar] [CrossRef] [Green Version]
  71. Gulácsi, A.; Kovács, F. Sentinel-1-imagery-based high-resolutionwater cover detection on wetlands, aided by google earth engine. Remote Sens. 2020, 12, 1–20. [Google Scholar] [CrossRef]
  72. Zhou, T.; Pan, J.; Zhang, P.; Wei, S.; Han, T. Mapping winter wheat with multi-temporal SAR and optical images in an urban agricultural region. Sensors 2017, 17, 1210. [Google Scholar] [CrossRef]
  73. Luti, T.; Segoni, S.; Catani, F.; Munafò, M.; Casagli, N. Integration of remotely sensed soil sealing data in landslide susceptibility mapping. Remote Sens. 2020, 12, 1486. [Google Scholar] [CrossRef]
  74. ISPRA Il consumo di suolo in Italia—Edizione 2015. Available online: https://www.isprambiente.gov.it/it/pubblicazioni/rapporti/il-consumo-di-suolo-in-italia-edizione-2015 (accessed on 18 March 2021).
  75. Holtgrave, A.K.; Röder, N.; Ackermann, A.; Erasmi, S.; Kleinschmit, B. Comparing Sentinel-1 and -2 data and indices for agricultural land use monitoring. Remote Sens. 2020, 12, 2919. [Google Scholar] [CrossRef]
  76. Formaggio, A.R.; Epiphanio, J.C.N.; dos Santos Simões, M. Radarsat backscattering from an agricultural scene. Pesqui. Agropecuária Bras. 2001, 36, 823–830. [Google Scholar] [CrossRef] [Green Version]
  77. Lefebvre, A.; Sannier, C.; Corpetti, T. Monitoring urban areas with Sentinel-2A data: Application to the update of the Copernicus High Resolution Layer Imperviousness Degree. Remote Sens. 2016, 8, 606. [Google Scholar] [CrossRef] [Green Version]
  78. Chirici, G.; Giannetti, F.; Mazza, E.; Francini, S.; Travaglini, D.; Pegna, R.; White, J.C. Monitoring clearcutting and subsequent rapid recovery in Mediterranean coppice forests with Landsat time series. Ann. For. Sci. 2020, 77, 1–14. [Google Scholar] [CrossRef]
  79. Olofsson, P.; Foody, G.M.; Herold, M.; Stehman, S.V.; Woodcock, C.E.; Wulder, M.A. Good practices for estimating area and assessing accuracy of land change. Remote Sens. Environ. 2014, 148, 42–57. [Google Scholar] [CrossRef]
  80. FAO. Map Accuracy Assessment and Area Estimation: A Practical Guide; National Forest Monitoring Assessment Working Paper; Food and Agriculture Organization of the United Nations: Rome, Italy, 2016. [Google Scholar]
  81. Stehman, S.V.; Foody, G.M. Key issues in rigorous accuracy assessment of land cover products. Remote Sens. Environ. 2019, 231, 111199. [Google Scholar] [CrossRef]
  82. Stehman, S.V.; Czaplewski, R.L. Design and Analysis for Thematic Map Accuracy Assessment—An application of satellite imagery. Remote Sens. Environ. 1998, 64, 331–344. [Google Scholar] [CrossRef]
  83. Cochran, W.G.; William, G. Sampling Techniques; John Wiley& Sons: New York, NY, USA, 1977. [Google Scholar]
  84. Maxwell, A.E.; Warner, T.A.; Fang, F. Implementation of machine-learning classification in remote sensing: An applied review. Int. J. Remote Sens. 2018, 39, 2784–2817. [Google Scholar] [CrossRef] [Green Version]
  85. Vuolo, F.; Neuwirth, M.; Immitzer, M.; Atzberger, C.; Ng, W.T. How much does multi-temporal Sentinel-2 data improve crop type classification? Int. J. Appl. Earth Obs. Geoinf. 2018, 72, 122–130. [Google Scholar] [CrossRef]
  86. Pesaresi, S.; Mancini, A.; Quattrini, G.; Casavecchia, S. Mapping mediterranean forest plant associations and habitats with functional principal component analysis using Landsat 8 NDVI time series. Remote Sens. 2020, 12, 1132. [Google Scholar] [CrossRef] [Green Version]
  87. Rüetschi, M.; Small, D.; Waser, L.T. Rapid detection of windthrows using sentinel-1 c-band sar data. Remote Sens. 2019, 11, 115. [Google Scholar] [CrossRef] [Green Version]
  88. Zimmermann, N.E.; Edwards, T.C., Jr.; Moisen, G.G.; Frescino, T.S.; Blackard, J.A. Remote sensing-based predictors improve distribution models of rare, early successional and broadleaf tree species in Utah. J. Appl. Ecol. 2007, 44, 1057–1067. [Google Scholar] [CrossRef] [Green Version]
  89. Zhang, F.; Yang, X. Improving land cover classification in an urbanized coastal area by random forests: The role of variable selection. Remote Sens. Environ. 2020, 251, 112105. [Google Scholar] [CrossRef]
  90. Millard, K.; Richardson, M. On the importance of training data sample selection in Random Forest image classification: A case study in peatland ecosystem mapping. Remote Sens. 2015, 7, 8489–8515. [Google Scholar] [CrossRef] [Green Version]
  91. Pepe, M.; Pompilio, L.; Gioli, B.; Busetto, L.; Boschetti, M. Detection and Classification of Non-Photosynthetic Vegetation from PRISMA Hyperspectral Data in Croplands. Remote Sens. 2020, 12, 3903. [Google Scholar] [CrossRef]
  92. Fassnacht, F.E.; Latifi, H.; Stereńczak, K.; Modzelewska, A.; Lefsky, M.; Waser, L.T.; Straub, C.; Ghosh, A. Review of studies on tree species classification from remotely sensed data. Remote Sens. Environ. 2016, 186, 64–87. [Google Scholar] [CrossRef]
  93. Modzelewska, A.; Fassnacht, F.E.; Stereńczak, K. Tree species identification within an extensive forest area with diverse management regimes using airborne hyperspectral data. Int. J. Appl. Earth Obs. Geoinf. 2020, 84, 101960. [Google Scholar] [CrossRef]
  94. Puletti, N.; Bascietto, M. Towards a tool for early detection and estimation of forest cuttings by remotely sensed data. Land 2019, 8, 58. [Google Scholar] [CrossRef] [Green Version]
  95. Banskota, A.; Kayastha, N.; Falkowski, M.J.; Wulder, M.A.; Froese, R.E.; White, J.C. Forest monitoring using Landsat time series data: A review. Can. J. Remote Sens. 2014, 40, 362–384. [Google Scholar] [CrossRef]
  96. Vangi, E.; D’Amico, G.; Francini, S.; Giannetti, F.; Lasserre, B.; Marchetti, M.; Chirici, G. The New Hyperspectral Satellite PRISMA: Imagery for Forest Types Discrimination. Sensors 2021, 21, 1182. [Google Scholar] [CrossRef]
  97. Rast, M.; Ananasso, C.; Bach, H.; Ben-Dor, E.; Chabrillat, S.; Colombo, R.; Del Bello, U.; Feret, J.B.; Giardino, C.; Green, R.O. Copernicus Hyperspectral Imaging Mission for the Environment: Mission Requirements Document; European Space Agency: Paris, France, 2019. [Google Scholar]
  98. Zakeri, F.; Mariethoz, G. A review of geostatistical simulation models applied to satellite remote sensing: Methods and applications. Remote Sens. Environ. 2021, 259, 112381. [Google Scholar] [CrossRef]
Figure 1. Workflow of the methodology.
Figure 1. Workflow of the methodology.
Land 10 00611 g001
Figure 2. Study area.
Figure 2. Study area.
Land 10 00611 g002
Figure 3. Land cover 2018 and land cover change 2017–2018 in Italy.
Figure 3. Land cover 2018 and land cover change 2017–2018 in Italy.
Land 10 00611 g003
Figure 4. Result of the accuracy assessment with reference to the sample of 2400 photo-interpreted points.
Figure 4. Result of the accuracy assessment with reference to the sample of 2400 photo-interpreted points.
Land 10 00611 g004
Figure 5. Artificial abiotic surfaces, in black (on the left), compared with a Sentinel-2 image (on the right).
Figure 5. Artificial abiotic surfaces, in black (on the left), compared with a Sentinel-2 image (on the right).
Land 10 00611 g005
Figure 6. Land consumption, highlighted in red on the third image.
Figure 6. Land consumption, highlighted in red on the third image.
Land 10 00611 g006
Figure 7. Water bodies.
Figure 7. Water bodies.
Land 10 00611 g007
Figure 8. Permanent snow and ice.
Figure 8. Permanent snow and ice.
Land 10 00611 g008
Figure 9. Woody vegetation, the distinction between broad-leaved (brighter) and needle-leaved (darker).
Figure 9. Woody vegetation, the distinction between broad-leaved (brighter) and needle-leaved (darker).
Land 10 00611 g009
Figure 10. Periodic herbaceous.
Figure 10. Periodic herbaceous.
Land 10 00611 g010
Figure 11. Permanent herbaceous.
Figure 11. Permanent herbaceous.
Land 10 00611 g011
Figure 12. Other disturbances—Forest harvesting, highlighted in red on the right image.
Figure 12. Other disturbances—Forest harvesting, highlighted in red on the right image.
Land 10 00611 g012
Figure 13. Burnt area, highlighted in grey on the right image.
Figure 13. Burnt area, highlighted in grey on the right image.
Land 10 00611 g013
Table 1. Land cover classification system.
Table 1. Land cover classification system.
Land Cover
I LevelII LevelIII Level
Abiotic-non vegetatedArtificial abiotic surfaces
Natural abiotic surfaces
Biotic-vegetatedWoody vegetationBroad-leaved
Needle-leaved
Herbaceous vegetationPermanent
Periodically
Water surfacesWater bodies
Permanent snow and ice
Table 2. Land cover change classification system.
Table 2. Land cover change classification system.
Land Cover Change
Land consumption
Restoration
Forest disturbancesBurnt areas
Other disturbances
Table 3. Land cover transformation matrix on yearly basis. The following transformations are considered: 3—Likely (Stability), 2—Probable (change can happen any time), 1—Possible (change can happen in specific cases), 0—Not considered.
Table 3. Land cover transformation matrix on yearly basis. The following transformations are considered: 3—Likely (Stability), 2—Probable (change can happen any time), 1—Possible (change can happen in specific cases), 0—Not considered.
2018
Artificial Abiotic Surfaces Natural Abiotic Surfaces Broad-Leaved Needle-Leaved Permanent Herbaceous Periodic Herbaceous Water Bodies Permanent Snow and Ice
2017 Artificial abiotic surfaces 32222200
Natural abiotic surfaces 23111100
Broad-leaved 22302200
Needle-leaved 22032200
Permanent herbaceous 21113100
Periodic herbaceous 21111300
Water bodies 00000030
Permanent snow and ice 00000003
Table 4. Spectral indices used for land cover classification and change detection.
Table 4. Spectral indices used for land cover classification and change detection.
SatelliteIndex NameBrief Index FormulaReference
Sentinel-2Normalized Difference Vegetation Index N D V I = B 8 B 4 B 8 + B 4 [65]
Normalized Burn Ratio N B R = B 8 B 12 B 8 + B 12 [66]
Normalized Difference Water Index N D W I = B 2 B 8 B 2 + B 8 [67]
Normalized Difference Snow Index N D S I = B 3 B 11 B 3 + B 11 [68]
Normalized Difference Coniferous Index N D C I = B 6 B 12 B 6 + B 12 B 8 B 11 B 8 + B 11
Burned Index B I = 1 B 3 + B 4 + B 8 1 + B 3 + B 4 + B 8
Table 5. Calculation of sample size.
Table 5. Calculation of sample size.
Land Cover ClassArea (ha)WiUiSiWi*Si
Artificial Abiotic 2,127,9990.0710.60.4900.035
Natural Abiotic 953,5740.0320.60.4900.015
Broad-leaved 12,006,5000.3980.60.4900.195
Needle-leaved1,607,5820.0530.60.4900.026
Permanent herb.12,657,9000.4200.60.4900.206
Periodic herb.239,8340.0080.60.4900.004
Water bodies337,7370.0110.60.4900.005
Permanent snow and ice111,7730.0040.60.4900.002
Other disturbance78,4900.0030.60.4900.001
Burned areas11,1060.0000.60.4900.000
Land consumption66000.0000.60.4900.000
Restoration9060.0000.60.4900.000
Total30,139,974 1
S(Ô) overall accuracy0.01
Total number of samples2400
Table 6. Allocation of the sample between the considered classes.
Table 6. Allocation of the sample between the considered classes.
ClassesAllocation
EqualProportionalFinal
Artificial Abiotic 200169185
Natural Abiotic 20076138
Broad-leaved 200956578
Needle-leaved200128164
Permanent herb.2001008604
Periodic herb.20019110
Water bodies20027113
Permanent snow and ice2009104
Other disturbance 2006103
Burned areas 2001100
Land consumption 2001100
Restoration 2000100
Total240024002400
Table 7. Confusion matrix obtained from the accuracy assessment.
Table 7. Confusion matrix obtained from the accuracy assessment.
Overall Accuracy
0.83
Land Cover
Class nameUser’s
accuracy
Producer’s
accuracy
Artificial abiotic surfaces 0.920.88
Natural abiotic surfaces 0.830.78
Water bodies 0.980.90
Permanent snow and ice 0.861.00
Broad-leaved 0.870.81
Needle-laved 0.900.88
Permanent herbaceous 0.920.62
Periodic herbaceous 0.860.81
Land Cover Change
Class nameUser’s
accuracy
Producer’s
accuracy
Other disturbance0.350.71
Burned areas 0.671.00
Land consumption 0.811.00
Restoration 0.501.00
Table 8. Distribution of land cover classes for 2018 and land cover change for 2017–2018 on the Italian territory.
Table 8. Distribution of land cover classes for 2018 and land cover change for 2017–2018 on the Italian territory.
ha% Total% Class
Abiotic Surfaces3,081,57310.22-
Artificial surfaces2,127,9997.0669.06
Natural surfaces953,5743.1630.94
Woody Vegetation13,614,06445.17-
Broad-leaved12,006,48239.8488.19
Needle-leaved1,607,5825.3311.81
Herbaceous Vegetation12,897,72442.79-
Periodic12,657,89042.0098.14
Permanent239,8340.801.86
Water and Ice449,5101.49-
Water bodies337,7371.1275.13
Permanent snow and ice 111,7730.3724.87
Land Cover Change97,1020.32
Italy30,139,974100.00-
Table 9. Distribution of land consumption between 2017 and 2018 among land cover classes.
Table 9. Distribution of land consumption between 2017 and 2018 among land cover classes.
Land ConsumptionRestoration
ha%m2chg/haha%%m2chg/ha
Abiotic Surfaces131119.94.29110.00.3
Artificial surfaces------
Natural surfaces131119.913.69110.00.9
Woody Vegetation104915.90.816918.60.
Broad-leaved101915.40.816618.30.1
Needle-leaved300.50.230.30.0
Herbaceous Vegetation424064.23.364771.40.5
Periodic423164.13.363670.20.5
Permanent90.10.3111.20.3
Water and Ice00.00.000.00.0
Water bodies00.00.000.00.0
Permanent snow and ice 00.00.000.00.0
Italy6600100.02.2906100.00.3
Table 10. Distribution of forest disturbances occurring between 2017 and 2018 among woody vegetation classes.
Table 10. Distribution of forest disturbances occurring between 2017 and 2018 among woody vegetation classes.
2017–18ha%
Burnt areasTotal11,10612.4
Broad-leaved661059.5
Needle-leaved449640.5
Other disturbancesTotal78,49087.6
Broad-leaved71,09690.6
Needle-leaved73949.4
Total 89,596100.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

De Fioravante, P.; Luti, T.; Cavalli, A.; Giuliani, C.; Dichicco, P.; Marchetti, M.; Chirici, G.; Congedo, L.; Munafò, M. Multispectral Sentinel-2 and SAR Sentinel-1 Integration for Automatic Land Cover Classification. Land 2021, 10, 611. https://0-doi-org.brum.beds.ac.uk/10.3390/land10060611

AMA Style

De Fioravante P, Luti T, Cavalli A, Giuliani C, Dichicco P, Marchetti M, Chirici G, Congedo L, Munafò M. Multispectral Sentinel-2 and SAR Sentinel-1 Integration for Automatic Land Cover Classification. Land. 2021; 10(6):611. https://0-doi-org.brum.beds.ac.uk/10.3390/land10060611

Chicago/Turabian Style

De Fioravante, Paolo, Tania Luti, Alice Cavalli, Chiara Giuliani, Pasquale Dichicco, Marco Marchetti, Gherardo Chirici, Luca Congedo, and Michele Munafò. 2021. "Multispectral Sentinel-2 and SAR Sentinel-1 Integration for Automatic Land Cover Classification" Land 10, no. 6: 611. https://0-doi-org.brum.beds.ac.uk/10.3390/land10060611

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