Next Article in Journal
The Impact of Different Value Types on Environmentally Responsible Behavior: An Empirical Study from Residents of National Park Communities in China
Next Article in Special Issue
Spatial and Temporal Variation Characteristics of Ecological Environment Quality in China from 2002 to 2019 and Influencing Factors
Previous Article in Journal
Unveiling Heavy Metal Links: Correlating Dust and Topsoil Contamination in Vilnius Schools
Previous Article in Special Issue
Investigating the Spatial Heterogeneity and Influencing Factors of Urban Multi-Dimensional Network Using Multi-Source Big Data in Hangzhou Metropolitan Circle, Eastern China
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Earth Observation Data and Geospatial Deep Learning AI to Assign Contributions to European Municipalities Sen4MUN: An Empirical Application in Aosta Valley (NW Italy)

by
Tommaso Orusa
1,2,*,
Annalisa Viani
3 and
Enrico Borgogno-Mondino
1
1
Department of Agricultural, Forest and Food Sciences (DISAFA), GEO4Agri DISAFA Lab, Università degli Studi di Torino, Largo Paolo Braccini 2, 10095 Grugliasco, Italy
2
IN.VA spa and Earth Observation Valle d’Aosta—eoVdA, Località L’Île-Blonde 5, 11020 Brissogne, Italy
3
Institute for Piedmont, Liguria, Aosta Valley (IZS PLV) S.C Valle d’Aosta—CeRMAS (National Reference Center for Wildlife Diseases), Località Amerique 7/G, 11020 Quart, Italy
*
Author to whom correspondence should be addressed.
Land 2024, 13(1), 80; https://doi.org/10.3390/land13010080
Submission received: 4 December 2023 / Revised: 8 January 2024 / Accepted: 9 January 2024 / Published: 10 January 2024
(This article belongs to the Special Issue Applying Earth Observation Data for Urban Land-Use Change Mapping)

Abstract

:
Nowadays, European program Copernicus’ Sentinel missions have allowed the development of several application services. In this regard, to strengthen the use of free satellite data in ordinary administrative workflows, this work aims to evaluate the feasibility and prototypal development of a possible service called Sen4MUN for the distribution of contributions yearly allocated to local municipalities and scalable to all European regions. The analysis was focused on the Aosta Valley region, North West Italy. A comparison between the Ordinary Workflow (OW) and the suggested Sen4MUN approach was performed. OW is based on statistical survey and municipality declaration, while Sen4MUN is based on geospatial deep learning techniques on aerial imagery (to extract roads and buildings to get real estate units) and yearly Land Cover map components according to European EAGLE guidelines. Both methods are based on land cover components which represent the input on which the financial coefficients for assigning contributions are applied. In both approaches, buffers are applied onto urban class (LCb). This buffer was performed according to the EEA-ISPRA soil consumption guidelines to avoid underestimating some areas that are difficult to map. In the case of Sen4MUN, this is applied to overcome Sentinel sensor limits and spectral mixing issues, while in the case of OW, this is due to limits in the survey method itself. Finally, a validation was performed assuming as truth the approach defined by law as the standard, i.e., OW, although it has limitations. MAEs involving LCb, road lengths and real estate units demonstrate the effectiveness of Sen4MUN. The developed approach suggests a contribution system based on Geomatics and Remote sensing to the public administration.

1. Introduction

The European Space program Copernicus, with its Sentinels missions, has allowed the creation of many research projects and prototypes aiming at developing several services, many of them already available based on earth observation data [1,2,3,4,5,6]. Unfortunately, there are still few available or prototype services based on an applied use of geomatics and remote sensing, despite an expansion of Space Economy. The distribution systems of regional contributions to municipalities are overwhelmingly based on territorial components. Therefore, the geospatial component and the technological transfer offered by the growing development of geomatics and remote sensing in this sector would be enormous, allowing for important impacts on administrative processes currently carried out with systems that are often onerous, as well as poorly efficient and inaccurate. To achieve real digitalization in public administration is necessary to develop services capable of responding to ordinary questions in a modern way, guaranteeing scalability and a certain degree of standardization. In fact, nowadays there is a growing need of standardization and validation procedures of geographical and earth observation data/product/services with special concerns about the expected roles they can have within monitoring/control actions from institutional subjects (e.g., CAP controls. natural hazards, etc.).
One of these, still under development, and one of the few that aims to use remote sensing in a highly economic key is the Sen4CAP system for the payment of contribution and control in agriculture, according to the Common Agricultural Policy (CAP) in light of EU regulation (N. 746/2018) [7,8,9]. The EU’s Common Agricultural Policy (CAP) seeks to raise agricultural productivity in Europe in a sustainable manner while maintaining a respectable level of living for EU farmers. The CAP, which has an annual budget of about EUR 59 billion, uses a variety of strategies such as direct payments, market reforms, and rural development, to improve the sustainability and competitiveness of agriculture in Europe [10]. The Integrated Administration and Control System (IACS) oversees and manages the majority of the CAP budget with the goal of protecting the program’s finances and assisting farmers in submitting their declarations. [11]. Earth observation (EO) satellites are regarded to play a bigger part in the CAP reform satellite’s enhancement and cost-effectiveness of the IACS. It is important to note that the CAP reform has established the role of EO data as being of utmost importance. Their adoption is expected to be required in 2024–2025, and may also extend to future planned missions like IRIDE, which the Italian Space Agency (ASI) and the European Space Agency (ESA) are working on together. Sentinel missions—Sentinel-1 and Sentinel-2 (hereinafter called Sentinels) in particular—play a major role. In the years after 2020, the initiative will give special consideration to demonstrating how data obtained from Sentinels may help modernize and streamline the CAP. Sen4CAP was established by the ESA in response to requests from European payment agencies like other EO services [12,13,14].
Also, other space agencies worldwide like NASA (US), JAXA (Japan), CNES (France), DLR (Germany), and more recently the ASI (Italy) with the ambitious program IRIDE, have developed or have planned to strengthen their services addressed to the private and public sector, not only for research purposes, but to promote technological transfer in many sectors [15], but at the present time fewer are related to contribution systems in different sectors which play a huge role. The space race has opened up new frontiers of investment for individuals and enterprises at different levels, not only the large economic giants of the ICT and geospatial services such as Planet, SatVu, Albedo, Maxar, e-Geos, etc. Despite several services being available from forest management monitoring to precision agriculture, to the management of migratory flows passing through the management of urban areas and planning of smart cities and much more [16], almost none are related to contributions system. The reason is related to the fact that, on one hand, many services and products offered are often too general and not specific, or worse, without solid scientific validation. On the other hand, public administration is too slow to embrace the new opportunities offered by applied sciences. It is often incapable of understanding them because they are not transferred in the correct way or suitable for use as a support and then a replacement for old methods without creating big changes (especially when it comes to money). Added to all this is bureaucracy and the need to legislate on new approaches so that they become operational and can complement and replace conventional ones. It is therefore not surprising that more bureaucratic countries are often slower to assimilate the news. Added to this is the need for foresight on the part of the political sector and of regional managers and technicians capable of understanding and getting involved, abandoning safe but often obsolete paths to support new ones that require hard work, experimentation, and continuous understanding [17].
Therefore, this work aims to create a tentative service based on European EO data and geospatially based processing about the distribution of contributions intended for municipalities called Sen4MUN, where Sen means Sentinels for municipalities. It is worth noting that in Italy, as in other EU countries, the municipalities receive income from different types according to different criteria [18]. Current incomes of taxes, contributions, and equalization nature are made up of four income items: taxes, duties, and similar income; tax sharing; equalization funds from central administrations; equalization funds from the region or autonomous province. Generally, equalization funds from central administrations and regions or autonomous provinces are the most important and the main core of the approach suggested [19,20]. Concerning contribution delivered by regions to municipalities, they are assigned following two possible approaches; the first one is based on the estimated expenditure, taking into account the previous year (the more you spend the more you receive money for local development within a budget defined at the regional and governmental level); the second is a little more rigorous, based on indicators of territorial development (in which, at least in Italy, the more the urbanized area develops, roads, housing units, factories and more in general areas for urban use, the higher the revenues) as a function of the resident population without going into the criticisms of the present approaches, which certainly should be reviewed from the point of view of environmental sustainability dealing with SDGs goals [21]. Since the second approach permits an exploitation of the possibility offered by EO data and GIS updated data; a system similar to the idea at the base of Sen4CAP has been suggested in this work. Urbanized components, like other territorial bio-physical surfaces, can potentially be mapped to the temporal resolution of Sentinels by translating the information produced by a medium-high resolution land cover into a datum that can be used by administrations at an economic level [22]. To date, the monitoring of land cover is carried out only from an environmental point of view (to map the territory or quantify land consumption or for other research purposes), but without, up to now, deriving any form of economic quantification. To try to bridge this gap, the Sen4MUN has been designed.
Sen4MUN aims to create a single and standardized approach for the item concerning government revenue or, as in the case of this study, regional revenues for local authorities considering also the environmental issues. Sen4MUN aims to replace the approach currently based on the estimation of territorial indicators obtained from statistical analyzes or very rough estimates used up to now by suggesting a more rigorous, efficient, and objective approach, based on the technology transfer offered by EO data trying to suggest new approaches by space economy and the actual and expected roles of geomatics within the next generation EU framework from science to public services.
In this regard, therefore, to strengthen the use of free European satellite data in Ordinary administrative workflows such as Sen4CAP, this work has been focused to assess the feasibility and prototypal development of a possible service called Sen4MUN for the distribution of contributions yearly allocated to local municipalities and scalable to all European regions.
The analysis was focused on Valle d’Aosta Region, North West Italy, considered complex in land cover mapping because of its geo-morphology and Sentinels geometrical resolution limits in alpine areas linked to SAR distortions [23] and spectral mixing, as well as phenology detection in pixels located in high slope degree in multispectral remote sensed data [24].
A comparison between the ordinary methodology based on the estimation of territorial indicators (OW), with the buffer zone retrieved with statistical surveys and the Sen4MUN approach, which is based on yearly land cover classification according to EAGLE guidelines, has been performed. Finally, due to the fact that some roads cannot be correctly mapped by Sentinels due to a GSD limiting factor, updated GIS geodatabases were included in the Sen4MUN prototypal approach.

2. Materials and Methods

2.1. Study Area

Sen4MUN was realized and tested within the Aosta Valley Region in the North West of Italy. Aosta Valley is a region in northwestern Italy which hosts the highest peaks in the whole alpine chain and in general in Europe. The region borders France in the Western part and Switzerland in the Northern. Figure 1 has reported Aosta Valley’s location in respect to Italy. Though it is the smallest in Italy, its geomorphology makes it one of the most complex. It is situated in the Western Alps [25].

2.2. Ordinary Workflow

This approach adopts municipality declaration and statistical data obtained from regional or national offices that starts from building practices or territorial surveys performed within a time range (for example every 5 years [26]. In the most likely cases, data comes from GIS databases (like the cadastre of buildings real estate units and municipalities streets lengths), and the land cover areas divided into three marco-classes derived from municipality declarations. These three land cover components are: urban areas which are the most important, semi-anthropic areas, and sterile areas, buffered 10–20 or 30 m (over-estimated) due to inhomogeneity data collection through the different times [22]. Nowadays, Ordinary Workflow (hereinafter called OW) does not use deep learning processing to update the GIS geodatabase and earth observation data. In fact, almost the entire OW data collection is based on municipality declaration and statistical surveys at different level with yearly temporal gaps. This is the case of the Aosta Valley region that follows the Cerutti’s approach since 1979, which is entirely based on statistical surveys and declarations [27,28,29]. This approach considered the following patterns: the municipal area, the urban and anthropic areas, the semi-anthropic areas (that includes all vegetated areas cultivated and not cultivated by humans), the sterile areas (that means all land unproductive surfaces, like water bodies and courses, snow and ice, rocks), roads lengths, and finally real estate units. All these components are properly weighted and normalized according to financial criteria defined at regional, national, and European levels. The financial criteria are based on the resources made available by the annual regional budget on the basis of a legal resolution, the expenses made in the current year, the actual resident population and the municipal area. These indicators, based on state and regional tax formulas, which vary over time because they are the result of political economy choices, define the coefficients to be applied to the territorial development components covered by this study and which characterize the contributions paid. In fact, the financial criteria are defined and regulated by national law and regional regulations, leaving greater freedom in their calculation to regions with special statutes such as the Valle d’Aosta Autonomous Region.
Therefore, OW is based on data collection campaigns for statistical purposes and declarations from municipalities. The municipal offices are required to communicate certain parameters such as the municipal area, resident population, roads belonging to the municipality and real estate units present in the municipality. The poor communication between offices at various levels often makes these data not always representative and homogeneous in qualitative terms for all municipalities.
It is worth noting that OW is based on a simple weighted sum of factors that sum to total land cover as reported in financial land components workflow in Figure 1 which is used also in Sen4MUN using a different collection data method. In this regard, the role of remote sensing and GIS is to properly map these areas instead of using periodic surveys as it happens nowadays and in the issue previously discussed. Equations (12) and (13) are the core of this approach, and each variable cab be obtained weighing each component (representing the input variable) properly mapped. Table 1 shows how each variable was obtained for the only purpose of complete clarity, although the procedure is as described previously, a simple weighing of weighed remote sensed variables.
The sub-conventional parameters, respectively obtained by Equations (7), (9) and (11) are used to obtain algorithms reported in Equations (12) and (13), respectively. These are adopted by the regional Aosta Valley Autonomous offices to assess the municipality contributions according to the financial budget foreseen annually.
  v = o + r + u  
where,
v = sum of all sub-conventional weighted parameters,
z = v × η
where,
η = a financial weight in this case 11.50%.

2.3. Sen4MUN

The Sen4MUN approach is based on remote sensing and, in particular, yearly land cover, GIS-updated geodatabases, and deep learning to assess the parameters necessary as input in the OW.

2.3.1. Earth Observation Data and Processing

In this regard, urban and anthropic areas, semi-anthropic areas, and sterile areas are computed from Aosta Valley yearly land cover classes aggregation available at SCT regional geoportal and eoVdA webpages. Road lengths are obtained by the GIS viability monthly update geodatabase, as well as real estate units also including deep learning (this part will be further discussed). The reference year 2020 was considered. Land cover (LC) complies with EAGLE guidelines and is realized according to the following methods [30] and it is based on Sentinels missions (S1–S2) [31,32,33,34] adopting Aosta Valley land cover [25,30]. It is worth to note that, to assist the semantic and technological foundation of a European harmonized information management capability for land monitoring, the EAGLE Group has been working on a solution and proof of concept since 2008. Mostly, but not only, in their capacities as Eionet members, land monitoring specialists from many European Environment Agency (EEA) member nations have formed the self-initiated and public EAGLE Group. As a result, the EAGLE Group uses a bottom-up methodology to compile information and experiences from current land cover (LC) and land use (LU) categorization methodologies and projects. The Copernicus Land Monitoring Service recognizes EAGLE as a key and necessary element to facilitate a broad change in emphasis from categorization to characterization. This has resulted in the enforcement of EAGLE compliance in new CLMS products. The EAGLE guidelines envisage dividing the land cover components according to a pyramidal hierarchical approach, starting from macro classes to reach detailed components.
Since LC classes are more detailed than those reported in the Cerutti’s approach (which are the same at EAGLE at macro-level), they were aggregated according to the Cerutti’s description to integrate the Sen4MUN approach into the OW as follows in Table 2.

2.3.2. Geospatial Deep Learning Data and Processing

Roads and real estate units were extracted both from cadastral maps and deep learning adopting open-source libraries and Python scripts integrated with ESRI ArcGIS Pro v.2.9 for object detection and classification [35,36,37,38]. Roads and building footprints [39,40,41,42,43] were extracted using Convolutional Neural Network (CNN) techniques [44,45,46,47] onto the AGEA (Agency for Disbursements in Agriculture) 2020 ortho-rectified imagery, yearly available at the national level. In particular, ArcGIS pretrained deep learning models were adopted to extract roads and real estate units (that has been assigned from cadastre and municipal declaration of habitability to building footprint extracted using these models). Concerning road extraction, the deep learning model to extract roads from high resolution satellite imagery named: Road Extraction—Global in ESRI ArcGIS Pro v.2.9 was adopted. The implementation is based on the Sat2Graph model by [48]. Sat2Graph relies on a novel encoding scheme. It is capable of creating a three-dimensional tensor from the road network graph. This capability is named as Graph Tensor Encoding (GTE). By combining the benefits of segmentation-based and graph-based techniques, this graph-tensor encoding scheme enables the training of a basic, non-recurrent neural network model to directly translate the input satellite/aerial images into the road network graph (i.e., edges and vertices). In the case of a road network graph G = {V, E} covering a region measuring W meters by H meters, GTE employs a W−λ × H−λ × (1 + 3 · Dmax) 3D-tensor (represented as T) to hold the graph’s encoding. In this case, Dmax is the maximum number of edges that can be encoded at each λ × λ grid, and λ is the spatial resolution, which limits the encoded graph so that no two vertices can be co-located within a λ × λ grid. The two spatial axes in the two-dimensional plane are represented by the first two dimensions of T. To encode the graph information, we utilize the vector at each spatial point, ux,y = [Tx,y,1, Tx,y,2,..., Tx,y,(1 + 3·Dmax)]T. Its first element (vertexness), pv ∈ [0, 1], specifies the likelihood that a vertex would exist at point (x, y). Dmax 3-element groups, which encode the data of a possible outgoing edge from location (x, y), come after the initial element. The first element pei ∈ [0, 1] (edgeness) of the i-th 3-element group encodes the likelihood of having an outgoing edge toward (dxi, dyi), that is, an edge pointing from (x, y) to (x + dxi, y + dyi). Since vertices with degrees greater than six are extremely uncommon in road network graphs, we have set Dmax to six in this instance. GTE only use the i-th 3-element group to encode edges pointing toward a 360 Dmax-degree sector from (i − 1) · 360 Dmax degrees to i · 360 Dmax degrees in order to minimize the number of possible distinct isomorphic encodings of the same input graph. It is simple to encode a road network graph into GTE. The encoding algorithm first interpolates the straight road segment in the road network graph for the purpose of road network extraction. To keep the distance between consecutive points under d meters, it chooses the bare minimum of equally spaced intermediate spots. By controlling the edge vector’s length in GTE, this interpolation technique stabilizes the training process.
Interpolation for stacked roads may result in vertices from two overlapped road segments at the same location. When this occurs, the pre-trained model in ArcGIS permits to move the endpoint vectors of the two edges using an iterative conflict-resolution process. The objective is to ensure that the separation between any two vertices (derived from the two edges that overlap). The GTE decoding technique returns a graph’s anticipated GTE, which is frequently noisy, to the standard graph format (G = {V,E}). There are two steps in the decoding algorithm: (1) vertex extraction and (2) edge connection. It has been considered vertices and edges with a probability larger than a threshold (referred to as pthr), since both the edgeness and verticeness predictions are real numbers between 0 and 1.
During the vertex extraction phase, the decoding algorithm locates the local maxima of the vertexness map in order to extract possible vertices. The algorithm takes into consideration just those local maxima whose vertexness exceeds pthr. The decoding algorithm joins the outgoing edges of each candidate vertex v ∈ V to other vertices in the edge connection stage. The following distance function (reported in Equation (14)) is used by the algorithm to calculate the distance of the i-th edge of vertex v ∈ V to all other neighboring vertices u.
d v , i , u = v x + d x i , v x + d y i u x , u y + + w · c o s d i s t d x i , d y i , u x v x , v y v y  
where w is the weight of the cosine distance in the distance function and cosdist(v1, v2) is the cosine distance of the two vectors. In this case, we set w to a high value, like 100, to prevent erroneous connections. Once this distance has been calculated, the decoding algorithm inserts an edge between v and u 0 and selects a vertex, u 0, that minimizes the distance function d(v, i, u). In order to prevent false edges from being added to the graph when there are no suitable candidate vertices nearby, we specified a maximum distance criterion of 15 m. Finally, Sat2Graph uses cross-entropy loss (denoted as LCE) and L2-loss. The vertexness channel (pv) and edgeness channels are subjected to the cross-entropy loss.
pei i ∈ {1, 2,..., Dmax}), and the edge vector channels ((dxi, dyi) i ∈ {1, 2,..., Dmax}) get the L2-loss. GTE varies across extended road sections.
In this instance, distinct ground truth labels in the GTE format can be mapped to the same road structure. We only calculate the losses for edgeness and edge vectors at position (x, y) when there is a vertex at that location in the ground truth due to this discrepancy. The overall loss function is displayed below in Equation (15) (ground truth is represented by T ^   p ^ v , p ^ e i , d x ^ i , and d y ^ i ).
L T , T ^ = x , y 1 . . W × 1 . . H ( L C E p v , p ^ v + + T ^ x , y , 1 × ( i = 1 D m a x ( L C E p e i , p ^ e i + L 2 ( d x i , d y i , d x ^ i , d y ^ i ) ) ) )    
Concerning building footprints, these were obtained using the deep learning model to extract building footprints from high-resolution aerial and satellite imagery named: Building Footprint Extraction—New Zealand. In this last case, real estate extraction is based on the Mask R-CNN model architecture implemented using ArcGIS API for Python [49].
Developed on top of Faster R-CNN, Mask R-CNN is a state-of-the-art model for instance segmentation. A region-based convolutional neural network called Faster R-CNN [50] provides bounding boxes together with a confidence score for each object’s class identification. Mask R-CNN works on two-stage mainly based on Faster R-CNN architecture: (a) (phase 1): there are two networks in the first stage; a region proposal network and a backbone network (ResNet, VGG, Inception, etc.) to provide a collection of region suggestions, these networks execute once for each image. The feature map’s regions that contain the object are called region proposals; (b) (phase 2): The network predicts object classes and bounding boxes for every suggested region that was acquired in phase one. While fully linked layers in the networks always need a constant size vector to produce predictions, each proposed region might have a variable size. The RoIAlign technique or RoI pool, which is quite similar to MaxPooling, are used to fix the size of these proposed regions. When segmenting scenes or objects with irregular borders, segmentation models may produce boundaries that are too smooth and may not be accurate. A point-based rendering neural network module named PointRend has been included as an upgrade into the pre-trained model in order to obtain a clear segmentation border. This module provides the segmentation problem with a rendering perspective, utilizing techniques from classical computer graphics. Labels are frequently predicted by image segmentation models on a regular grid with low resolution, such as 1/8th of the input. To upscale the forecasts to the original resolution in these models, interpolation is used. PointRend, on the other hand, upscales the predictions using an iterative subdivision technique by having a trained tiny neural network predict the labels of points at certain places. This technique efficiently produces output with great resolution [51].
In this work, in order to improve the quality of the results obtained and, in particular, the extraction of a building footprint to get real estate units, PointRend has been adopted into the pre-trained model. To enable PointRend within Mask R-CNN in the ESRI ArcGIS front-end or scripting console, the following parameters has been set as reported in Equation (16):
m o d e l = M a s k R C N N d a t a = d a t a , p o i n t r e n d = T r u e
An 8-bit, 4-band high resolution (0.20 m GSD) aerial AGEA imagery of 2020 of the whole Aosta Valley was adopted. Due to the high computation necessary to extract the features in the whole region, the processing was conducted in a workstation with the following main characteristics: 64 GB CPU memory; 2 TB SSD storage and with a graphic card Nvidia RTX A4000 16 GB GDDR6.
Then, from the roads extracted, they were checked with the viability geodatabase, and only municipality roads and areas were considered (as specifically reported in the regional laws) for the road length computation. Then, the land cover areas and roads were buffered 20 m for the following reasons: some sparse urban areas are difficult to map due to the GSD limits of Sentinels [52,53,54,55,56,57,58,59], then this threshold is adopted both in OW procedure and ISPRA Land Units in the case of soil consumption estimation [29]. Finally, the areas obtained were used as inputs in the OW equations reported in Table 1.

2.4. Validation

The validation of the Sen4MUN approach and in particular of the surface estimated adopting GIS [60,61,62,63,64] and earth observation data to retrieve the parameters adopted in the OW was realized by computing the Mean Absolute Error (MAE) according to Equation (17) here reported:
M A E = i = 1 n | p i o i | n
where pi is the prediction (Sen4MUN component area), oi is the OW component areas estimated without remote sensing methods, and n is the number of municipalities in the Aosta Valley autonomous area equal to 74).
It is worth noting that the EAGLE classification starts from the same macro-classes identified by Cerutti’s. The unique difference in the two approaches is represented by the names adopted (despite the description is perfectly the same). Therefore, the EAGLE approach starts from macro to arrive to detail with sub-classes while Cerutti’s stops to macro-level.
This is a key point because it permits comparability between the two methods even if different technical approaches are followed.
In order to sum up the suggested newer approach, a workflow of the Sen4MUN approach is provided in Figure 2.

3. Results

The Sen4MUN approach was validated by computing the MAE compared to the surfaces calculated with the traditional method from statistical surveys. It is worth noting that, for the reference year 2020, the main surfaces and inputs were computed both with OW and Sen4MUN for each municipality within the Aosta Valley. Furthermore, also deep learning pre-trained models concerning on roads and real estate units have been evaluated. The pre-trained ESRI ArcGIS model Building Footprint Extraction—New Zealand has produced an F1 score of 86.20 all over the region, with an overall precision of 0.885 and a recall of 0.861. Despite many building footprints being available in ArcGIS, this is one that works better in geomorphological complex areas, while the pre-trained ESRI ArcGIS model Road Extraction—Global showed a precision score of 0.804 in the area of study. However, including results without PointRend as reported in Equation (16), it had a score of 0.873. Excluding false positives due to rocks artifacts and after a manual refining due to photointerpretation, the final layer had an overall accuracy of 0.961. It is worth noting that a 4-bands aerial imagery was adopted to perform the extraction. In particular, to refine the accuracy of the models in ArcGIS Pro in order to reach a threshold score of 0.80, hyperparameter optimization has been performed, fine-tuning the hyperparameters of the model, such as learning rate, batch size, and optimizer settings, to improve the training process and enhance the model’s performance. Moreover, a repeat cycle approach has been followed by iteration through the training, evaluation, and refinement steps multiple times until the model’s performance reaches an acceptable level. Finally, MAEs were computed for land cover components, road length and real estate units.
Since urban areas play a major role in the present regional regulation in terms of township incomes, MAEs involving each municipality were computed also considering their road length and real estate units. The results obtained are reported in Table 3 and a general Table A1 in Appendix A and in the graphs concerning on the MAEs reported in Figure 3.
Table 3 reports the overall MAE computed considering the same component retrieved with OW and Sen4MUN.
Figure 3 shows the MAEs for the three territorial components used for the calculation of municipal contributions by the Valle d’Aosta Autonomous Region. In particular, the following are reported in the section: (A) the MAEs for urban and anthropic areas expressed in square km; (B) the MAEs for the length of municipal roads expressed in km; (C) the MAEs for real estate units expressed in units. In all three sections of Figure 3, it is interesting to note that the largest MAEs generally occur in municipalities with fragmented urban areas with scattered villages and houses and roads in steep areas due to the limitations of the sensor and processing techniques or errors inherent in the traditional method, such as to cause the imbalances as illustrated below.
As reported in Table 3, an overall MAE of 0.16 km2 was obtained involving urban and anthropic areas, while 0.81 km was obtained for road length and 11 units in the case of real estate units. An overall MAE of 0.82 km2 was computed involving all EAGLE Land cover classes, respecting those computed with Cerutti’s approach. All the values obtained are significant and seem to suggest a validity of Sen4MUN. The values obtained in Table 3 for the two approaches are not very distant as indicated by the MAEs. It should be highlighted that the standard for assigning contributions is currently OW because it is regulated at a legislative level. It is worth noting that the errors came from Cerutti’s approach which overestimate some areas due to limits in this approach itself. Furthermore, it is interesting to note that in all the components, for example urbanized areas, road lengths and real estate units, MAEs are associated with municipalities that underwent more changes in at least one of the components in the reference year and with sparse villages. Unlike what can be derived from the MAEs here computed, the Sen4MUN system is more effective in monitoring territorial changes because it is based on satellite and GIS data updated at a high temporal frequency than the ordinary system. The MAEs obtained allow the transferability of the new approach and possible replacement of the ordinary one even in an Alpine reality such as the one investigated. In fact, MAEs are generally higher in mountainous areas than those obtainable in lowland areas due to the limiting factors affecting some remote sensing applications (land cover mapping in particular).

4. Discussions

The results show the consistency of the methods. Furthermore, the procedure adopts free worldwide coverage data. The creation of land covers on an annual, biennial, or monthly basis or even on a time scale of Sentinels (potentially every 5 days) allows the continuous monitoring of territorial dynamics. Moreover, this activity is part of the actions that each European country carries out for the monitoring of soil consumption and that the regions or provinces carry out to offer detailed products according to the guidelines of the European Environmental Agency. Furthermore, the availability for the administration of a regional air flight with very high resolution imagery (like an AGEA flight), albeit for agricultural purposes, allows this type of data to be integrated into other service chains. Such as in the case of Sen4MUN, to extract the footprint of buildings and roads thus makes public spending more efficient through the expansion of derivable products and services. Certainly, it is complex to create universally valid models, and locally suitable models are always preferable. But it must also be said that to create standards such as the Sen4CAP adopted by the European Union, compromises must be reached. Sen4MUN was specifically developed in a complex area such as to make it scalable in a much simpler way to plain and hilly areas which are the majority at a European level. The fact that the approach returns interesting results in an alpine and geomorphologically complex context is certainly encouraged from a scalability perspective. It should be remembered that the largest MAEs occur in municipalities that have more scattered settlements with isolated houses mixed with woods. Areas by which their nature are complex to map at Sentinel scale are not so numerous. Finally, it should also be remembered that many of the differences are linked not to a lack of accuracy and precision of Sen4MUN, which from a purely photo-interpretative analysis proved to be more performing than the traditional OW approach, but precisely as a result of errors in the dataset from statistical investigations of OW. This dataset, despite the errors it suffers from, was taken as a reference base as current regulations require this to be the data collection tool for the purposes of assigning contributions. During the analyses in some municipalities with high MAEs, it was found that Sen4MUN, on the contrary, mapped components that had escaped OW.
Naturally, the critical points in the Sen4MUN approach result in a correct mapping and knowledge of geomatics and remote sensing and their limits, such as the most suitable approach (hierarchical, single-directed, data fusion, etc.), the classification algorithms, as well as the input data that best respond to the components to be mapped, etc. Furthermore, the definition of the optimal number of training areas validation sets according to the area to be mapped is crucial [65,66,67,68,69,70,71,72,73,74,75,76] in obtaining high accuracies and minimization of errors [77,78,79,80,81,82,83,84,85,86,87]. Despite this work, a ready land cover has been adopted with an overall accuracy upper than 0.94, the whole procedure considering all these issues has been tackled in previous scientific literature involving the Aosta Valley territory [88]. In fact, in the case of a stand-alone and not prototype of this service at an application level not only on a regional or national scale, but on a European scale given the characteristics of the continent, the optimal solution would be to implement a regulated procedure, but which gives a certain degree of freedom in the algorithms and training sets (providing only thresholds) so that chains capable of responding best according to the area to be mapped are developed. Sen4CAP has been excessively standardized also in the algorithms making it underperform in the mountain areas [89]. The need to standardize as much as possible a procedure that moves economic contributions is crucial for large-scale regulation [47], but this is performed with the support of researchers, technicians, and academic experts and after several experiments in UE, to better define the operational leans and freedoms useful for developing consistent products and services for each European reality without creating disparities and differences. As far as the authors are aware, there are no applications in scientific literature for the monetarization of land cover and usage or AI geospatial deep learning at a national, regional, or local contribution scale for services provided by the public administration, therefore a comparison of the results obtained to date is extremely difficult. At the same time, fewer are scientific works on monetary attribution to the components of the coverage for the estimation of ecosystem services and nature based solutions [90,91], but with the intent of estimates of mere research and not of actual transfer of technological application to the public sector. Furthermore, the works themselves are on a large and non-municipal scale and involve geographical areas different from those under study. Therefore, we hope that this work will stimulate the Italian, European and international scientific and technical community to deepen and explore the application potential of geomatics in the operational workflows of public administration by offering new services and comparing the results and approaches obtained for the purpose of continuous improvement.
The use of the prototypal Sen4MUN Aosta Valley has made it possible to make the procedure for assigning contributions to municipalities more objective through the use of Sentinel data and updated GIS geodatabases, favoring technology transfer and the implementation of a possible new service within the Copernicus program and other future programs, such as an IRIDE capable of offering increasingly high spatial and temporal resolution data [37] useful for direct applications to the public sector [92,93,94,95,96].

5. Conclusions

Sen4MUN can be used as a standard procedure for assigning contributions to municipalities. The procedure is consistent and in line with the ordinary territorial workflow based on statistical surveys without the use of earth observation data, deep learning, and geodatabase. Although the system has been tested and used on a prototype level in the Aosta Valley, it can be scaled up to other Italian and European regions. The hope is that this service will become operational and spread to all member countries of the European Union as well as other national realities. Sen4MUN could join other services; coming from various Earth observation programs such as Copernicus, favoring an ever more massive technology transfer to the public sector by rationalizing activities and processes with plural activities with a view to the ever-increasing importance of geomatics in management flows and implementation planners of local, national and European policies. Finally, Sen4MUN seems to be capable of providing useful data for contribution purposes to public administration.

Author Contributions

Conceptualization, T.O.; methodology, T.O. and A.V.; software, T.O.; validation, T.O.; formal analysis, T.O.; investigation, T.O.; resources, T.O.; data curation, T.O. and A.V.; writing—original draft preparation, T.O.; writing—review and editing, T.O. and A.V.; visualization, T.O.; supervision, T.O. and E.B.-M.; project administration, T.O. and E.B.-M.; funding acquisition, T.O. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

Data can be accessed through request at the Enti locali della Regione Autonoma Valle d’Aosta offices.

Acknowledgments

Thanks to the Regione Autonoma Valle d’Aosta and INVA spa and eoVdA for having made this work possible, and, in particular, the Regional Cartographic Office in particular Chantal Trèves, and to Enti Locali e Ufficio finanza e contabilità degli Enti locali della Regione Autonoma Valle d’Aosta and in particular, Tiziana Vallet, and also Stefania Fanizzi, Nicoletta Berno, Emanuela Oro, Alessandra Sibona and more in general, everyone working in the Aosta Valley offices. Finally, the Geo4Agri DISAFA Lab colleagues for the great support.

Conflicts of Interest

The authors declare no conflicts of interest.

Appendix A

The data used to realize the histogram available in Figure 3 are reported here below in Table A1.
Table A1. MAEs and areas computed with both OW and Sen4MUN (S4M).
Table A1. MAEs and areas computed with both OW and Sen4MUN (S4M).
Aosta Valley MunicipalitiesUrban Areas (km2)Roads Length (km)Real Estate Units
S4MOWMAES4MOWMAES4MOWMAE
Allein0.280.330.0516.3715.321.056436385
Antey-Saint-André1.041.000.0417.9815.672.313250320842
Aosta9.589.570.01106.26104.551.7165,07365,05023
Arnad1.681.410.2729.3426.732.612614260113
Arvier0.870.660.2122.7322.020.71197019711
Avise0.450.500.0512.5911.581.008998909
Ayas2.842.850.0121.2819.042.2410,02210,0211
Aymavilles1.371.390.0131.0530.031.02334933463
Bard0.260.200.063.042.560.493203155
Bionaz0.550.580.0315.2213.951.266266224
Brissogne2.110.791.3219.7818.461.321510149911
Brusson1.551.860.3128.9128.750.16522952227
Challand-Saint-Anselme0.850.930.0819.3015.493.81323132301
Challand-Saint-Victor0.630.700.0716.9315.661.27177617706
Chambave1.390.880.5224.3625.130.77180417995
Chamois0.170.330.163.153.150.005685626
Champdepraz0.940.650.2917.6015.781.81141114121
Champorcher0.691.000.3122.9721.721.25235523514
Charvensod1.551.420.139.959.900.05389738907
Châtillon2.972.280.6947.6343.324.31751575105
Cogne2.101.720.3825.2122.442.78552855271
Courmayeur3.372.900.4753.5252.900.6215,33415,30133
Donnas1.901.720.1826.8928.211.324069404821
Doues0.600.700.1024.4127.032.621329131712
Emarèse0.290.340.0613.5712.850.729849777
Etroubles0.700.570.1318.9014.894.011314130212
Fénis1.471.280.1827.1326.300.83340934054
Fontainemore0.661.470.8126.9126.740.161560153525
Gaby0.610.800.199.519.250.261394133757
Gignod1.371.230.1530.4430.140.302692267319
Gressan2.842.520.3227.4326.261.178330831317
Gressoney-La-Trinité0.820.880.064.183.800.381944191232
Gressoney-Saint-Jean1.861.970.1018.5616.322.25535053473
Hône1.160.890.2713.0014.511.52222622260
Introd0.560.530.0313.8713.330.54144414359
Issime0.771.020.2510.479.560.91126312603
Issogne1.501.080.4226.3530.834.482305229114
Jovençan0.720.430.2810.4513.322.8797495915
La Magdeleine0.280.360.085.465.620.161157110156
La Salle2.251.950.2942.8136.216.597445741926
La Thuile1.771.440.3327.2324.942.287223720518
Lillianes0.520.980.4619.4822.282.801160114713
Montjovet1.621.420.2139.4138.810.60327132629
Morgex2.331.660.6722.9719.413.57698269739
Nus2.441.880.5654.4051.672.73512551250
Ollomont0.420.550.138.807.681.129659605
Oyace0.270.260.013.501.931.574754732
Perloz0.460.910.4519.4521.952.511248123711
Pollein2.111.160.9510.2714.133.862217220611
Pontboset0.330.470.1410.3111.321.0180979118
Pontey1.170.500.676.804.722.07129612942
Pont-Saint-Martin2.001.610.3921.3921.490.10516551605
Pré-Saint-Didier1.100.910.1920.5819.461.13592559178
Quart3.902.561.3459.6958.830.86649564896
Rhêmes-Notre-Dame0.310.390.086.285.580.7071970712
Rhêmes-Saint-Georges0.310.390.086.195.570.6282080812
Roisan0.660.660.0113.4412.450.99145514478
Saint-Christophe2.981.931.0444.1844.320.15551955118
Saint-Denis0.490.530.049.986.923.06116311612
Saint-Marcel1.461.190.2734.9934.130.86230822999
Saint-Nicolas0.460.480.0219.1917.172.02118411759
Saint-Oyen0.350.220.147.685.721.965565551
Saint-Pierre2.221.820.4140.6737.223.45522552223
Saint-Rhémy-en-Bosses0.880.730.1519.4818.510.98127412668
Saint-Vincent2.322.130.2040.4839.291.198956894313
Sarre2.742.180.5641.5340.331.216616659125
Torgnon1.101.100.0021.1719.561.614332431319
Valgrisenche0.360.540.1714.0013.780.228598518
Valpelline0.600.620.0213.5814.250.671322130022
Valsavarenche0.480.660.188.469.290.831036102313
Valtournenche2.762.660.1037.5637.510.0513,87013,8619
Verrayes1.761.590.1726.2926.120.17321832108
Verrès2.141.270.8713.5010.742.77395639533
Villeneuve1.431.010.4228.3130.041.73219221875
TOTAL103.992.10.161686.71626.40.81303,049293,21411

References

  1. Matevosyan, H.; Lluch, I.; Poghosyan, A.; Golkar, A. A Value-Chain Analysis for the Copernicus Earth Observation Infrastructure Evolution: A Knowledgebase of Users, Needs, Services, and Products. IEEE Geosci. Remote Sens. Mag. 2017, 5, 19–35. [Google Scholar] [CrossRef]
  2. Žlebir, S. Copernicus Earth Observation Programme. In Proceedings of the 40th COSPAR Scientific Assembly, Moscow, Russia, 2–10 August 2014. [Google Scholar]
  3. Schroedter-Homscheidt, M.; Arola, A.; Killius, N.; Lefèvre, M.; Saboret, L.; Wandji, W.; Wald, L.; Wey, E. The Copernicus Atmosphere Monitoring Service (CAMS) Radiation Service in a Nutshell. In Proceedings of the 22nd SolarPACES Conference, Abu Dhabi, United Arab Emirates, 11–14 October 2016. [Google Scholar]
  4. Thépaut, J.-N.; Dee, D.; Engelen, R.; Pinty, B. The Copernicus Programme and Its Climate Change Service. In Proceedings of the IGARSS 2018—2018 IEEE International Geoscience and Remote Sensing Symposium, Valencia, Spain, 22–27 July 2018; pp. 1591–1593. [Google Scholar]
  5. Szantoi, Z.; Strobl, P. Copernicus Sentinel-2 Calibration and Validation; Taylor & Francis: Abingdon, UK, 2019; Volume 52. [Google Scholar]
  6. Peuch, V.-H.; Engelen, R.; Rixen, M.; Dee, D.; Flemming, J.; Suttie, M.; Ades, M.; Agustí-Panareda, A.; Ananasso, C.; Andersson, E.; et al. The Copernicus Atmosphere Monitoring Service: From Research to Operations. Bull. Am. Meteorol. Soc. 2022, 103, E2650–E2668. [Google Scholar] [CrossRef]
  7. Koetz, B.; Defourny, P.; Bontemps, S.; Bajec, K.; Cara, C.; de Vendictis, L.; Kucera, L.; Malcorps, P.; Milcinski, G.; Nicola, L.; et al. SEN4CAP Sentinels for CAP Monitoring Approach. In Proceedings of the 2019 JRC IACS Workshop, Valladolid, Spain, 10–11 April 2019. [Google Scholar]
  8. Sarvia, F.; Petris, S.D.; Orusa, T.; Borgogno-Mondino, E. MAIA S2 Versus Sentinel 2: Spectral Issues and Their Effects in the Precision Farming Context. In International Conference on Computational Science and Its Applications; Springer: Berlin/Heidelberg, Germany, 2021; pp. 63–77. [Google Scholar]
  9. Koontz, T.M. Money Talks? But to Whom? Financial versus Nonmonetary Motivations in Land Use Decisions. Soc. Nat. Resour. 2001, 14, 51–65. [Google Scholar]
  10. Sishodia, R.P.; Ray, R.L.; Singh, S.K. Applications of Remote Sensing in Precision Agriculture: A Review. Remote Sens. 2020, 12, 3136. [Google Scholar] [CrossRef]
  11. Lupia, F.; Antoniou, V. Copernicus Sentinels Missions and Crowdsourcing as Game Changers for Geospatial Information in Agriculture. GEOmedia 2018, 22, 32–35. [Google Scholar]
  12. Carella, E.; Orusa, T.; Viani, A.; Meloni, D.; Borgogno-Mondino, E.; Orusa, R. An Integrated, Tentative Remote-Sensing Approach Based on NDVI Entropy to Model Canine Distemper Virus in Wildlife and to Prompt Science-Based Management Policies. Animals 2022, 12, 1049. [Google Scholar] [CrossRef] [PubMed]
  13. Viani, A.; Orusa, T.; Borgogno-Mondino, E.; Orusa, R. Snow Metrics as Proxy to Assess Sarcoptic Mange in Wild Boar: Preliminary Results in Aosta Valley (Italy). Life 2023, 13, 987. [Google Scholar] [CrossRef]
  14. Ippoliti, C.; Candeloro, L.; Gilbert, M.; Goffredo, M.; Mancini, G.; Curci, G.; Falasca, S.; Tora, S.; Di Lorenzo, A.; Quaglia, M.; et al. Defining Ecological Regions in Italy Based on a Multivariate Clustering Approach: A First Step towards a Targeted Vector Borne Disease Surveillance. PLoS ONE 2019, 14, e0219072. [Google Scholar] [CrossRef]
  15. Amani, M.; Ghorbanian, A.; Ahmadi, S.A.; Kakooei, M.; Moghimi, A.; Mirmazloumi, S.M.; Moghaddam, S.H.A.; Mahdavi, S.; Ghahremanloo, M.; Parsian, S.; et al. Google Earth Engine Cloud Computing Platform for Remote Sensing Big Data Applications: A Comprehensive Review. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2020, 13, 5326–5350. [Google Scholar] [CrossRef]
  16. Orusa, T.; Mondino, E.B. Landsat 8 Thermal Data to Support Urban Management and Planning in the Climate Change Era: A Case Study in Torino Area, NW Italy. In Remote Sensing Technologies and Applications in Urban Environments IV; International Society for Optics and Photonics: Strasbourg, France, 2019; Volume 11157, p. 111570O. [Google Scholar]
  17. Gascon, F.; Cadau, E.; Colin, O.; Hoersch, B.; Isola, C.; Fernández, B.L.; Martimort, P. Copernicus Sentinel-2 Mission: Products, Algorithms and Cal/Val. In Earth Observing Systems XIX; International Society for Optics and Photonics: Bellingham, WA, USA, 2014; Volume 9218, p. 92181E. [Google Scholar]
  18. Louvin, R. Flessibilità Fiscale e Zone Franche. Profili Giuridici e Finanziari. In IUS Publicum Europaeum; Nuova Serie: Aosta, Italy, 2022; Volume 9, pp. 1–240. [Google Scholar]
  19. Louvin, R. Il Comitato Europeo Delle Regioni: Bilancio e Rilancio. In Il Comitato delle Regioni, Regioni e Regioni Alpine: Riflessioni ed Esperienze sul futuro dell’Unione Europea; Università degli Studi di Trento, Facoltà di Giurisprudenza: Trento, Italy, 2022; pp. 61–68. [Google Scholar]
  20. Louvin, R. L’évolution Des Compétences Communales En Italie: Analogies et Discordances Par Rapport Au Cadre Français. In Quelle (s) Commune (s) pour le XXIe Siècle? Approche de Droit Comparé; L’Harmattan: Paris, France, 2018; pp. 193–216. [Google Scholar]
  21. 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; Springer: Berlin/Heidelberg, Germany, 2018; pp. 191–219. [Google Scholar]
  22. Congedo, L.; Sallustio, L.; Munafò, M.; Ottaviano, M.; Tonti, D.; Marchetti, M. Copernicus High-Resolution Layers for Land Cover Classification in Italy. J. Maps 2016, 12, 1195–1205. [Google Scholar] [CrossRef]
  23. Samuele, D.P.; Filippo, S.; Orusa, T.; Enrico, B.-M. Mapping SAR Geometric Distortions and Their Stability along Time: A New Tool in Google Earth Engine Based on Sentinel-1 Image Time Series. Int. J. Remote Sens. 2021, 42, 9135–9154. [Google Scholar] [CrossRef]
  24. Orusa, T.; Viani, A.; Cammareri, D.; Borgogno Mondino, E. A Google Earth Engine Algorithm to Map Phenological Metrics in Mountain Areas Worldwide with Landsat Collection and Sentinel-2. Geomatics 2023, 3, 221–238. [Google Scholar] [CrossRef]
  25. Orusa, T.; Cammareri, D.; Borgogno Mondino, E. A Possible Land Cover EAGLE Approach to Overcome Remote Sensing Limitations in the Alps Based on Sentinel-1 and Sentinel-2: The Case of Aosta Valley (NW Italy). Remote Sens. 2022, 15, 178. [Google Scholar] [CrossRef]
  26. Ferrario, C.; Ferri, V. I Comuni Italiani e l’autonomia Finanziaria: Una Scomoda Opportunità? Sci. Reg. 2023, 23, 105–138. [Google Scholar]
  27. Cerutti, P. Uso Del Territorio e Forme Contributive, 1st ed.; UTET: Torino, Italy, 1979; Volume 1. [Google Scholar]
  28. Rosanò, A. La Riforma Della Legge Della Regione Autonoma Valle d’Aosta in Materia Di Attività Condotte Nell’ambito Delle Politiche Promosse Dall’Unione Europea. In Quaderni AISDUE 1/2023; Editoriale Scientifica: Naples, Italy, 2023; pp. 359–372. [Google Scholar]
  29. 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]
  30. Orusa, T.; Cammareri, D.; Borgogno Mondino, E. A Scalable Earth Observation Service to Map Land Cover in Geomorphological Complex Areas beyond the Dynamic World: An Application in Aosta Valley (NW Italy). Appl. Sci. 2022, 13, 390. [Google Scholar] [CrossRef]
  31. Berger, M.; Aschbacher, J. Preface: The Sentinel Missions—New Opportunities for Science. Remote Sens. Environ. 2012, 120, 1–2. [Google Scholar] [CrossRef]
  32. 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]
  33. Bereta, K.; Caumont, H.; Daniels, U.; Goor, E.; Koubarakis, M.; Pantazi, D.-A.; Stamoulis, G.; Ubels, S.; Venus, V.; Wahyudi, F. The Copernicus App Lab Project: Easy Access to Copernicus Data. In Proceedings of the 22nd International Conference on Extending Database Technology (EDBT), Lisbon, Portugal, 26–29 March 2019; pp. 501–511. [Google Scholar]
  34. 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]
  35. Zhang, X.; Li, X.; An, J.; Gao, L.; Hou, B.; Li, C. Natural Language Description of Remote Sensing Images Based on Deep Learning. In Proceedings of the 2017 IEEE International Geoscience and Remote Sensing Symposium (IGARSS), Fort Worth, TX, USA, 23–28 July 2017; pp. 4798–4801. [Google Scholar]
  36. Shakya, A.; Biswas, M.; Pal, M. Parametric Study of Convolutional Neural Network Based Remote Sensing Image Classification. Int. J. Remote Sens. 2021, 42, 2663–2685. [Google Scholar] [CrossRef]
  37. Bullock, E.L.; Healey, S.P.; Yang, Z.; Houborg, R.; Gorelick, N.; Tang, X.; Andrianirina, C. Timeliness in Forest Change Monitoring: A New Assessment Framework Demonstrated Using Sentinel-1 and a Continuous Change Detection Algorithm. Remote Sens. Environ. 2022, 276, 113043. [Google Scholar] [CrossRef]
  38. Cao, Y.; Niu, X.; Dou, Y. Region-Based Convolutional Neural Networks for Object Detection in Very High Resolution Remote Sensing Images. In Proceedings of the 2016 12th International Conference on Natural Computation, Fuzzy Systems and Knowledge Discovery (ICNC-FSKD), Changsha, China, 13–15 August 2016; pp. 548–554. [Google Scholar]
  39. Ji, S.; Zhang, C.; Xu, A.; Shi, Y.; Duan, Y. 3D Convolutional Neural Networks for Crop Classification with Multi-Temporal Remote Sensing Images. Remote Sens. 2018, 10, 75. [Google Scholar] [CrossRef]
  40. Baroud, S.; Chokri, S.; Belhaous, S.; Mestari, M. A Brief Review of Graph Convolutional Neural Network Based Learning for Classifying Remote Sensing Images. Procedia Comput. Sci. 2021, 191, 349–354. [Google Scholar] [CrossRef]
  41. Pan, X.; Zhao, J. A Central-Point-Enhanced Convolutional Neural Network for High-Resolution Remote-Sensing Image Classification. Int. J. Remote Sens. 2017, 38, 6554–6581. [Google Scholar] [CrossRef]
  42. Shirmard, H.; Farahbakhsh, E.; Heidari, E.; Beiranvand Pour, A.; Pradhan, B.; Müller, D.; Chandra, R. A Comparative Study of Convolutional Neural Networks and Conventional Machine Learning Models for Lithological Mapping Using Remote Sensing Data. Remote Sens. 2022, 14, 819. [Google Scholar] [CrossRef]
  43. Cresson, R. A Framework for Remote Sensing Images Processing Using Deep Learning Techniques. IEEE Geosci. Remote Sens. Lett. 2018, 16, 25–29. [Google Scholar] [CrossRef]
  44. Albawi, S.; Mohammed, T.A.; Al-Zawi, S. Understanding of a Convolutional Neural Network. In Proceedings of the 2017 International Conference on Engineering and Technology (ICET), Antalya, Turkey, 21–23 August 2017; pp. 1–6. [Google Scholar]
  45. Zhong, J.; Zhong, S.; Zhang, Q.; Maia, N.; Shen, Y.; Liu, S.; Yu, Y.; Peng, Z. Vision-Based System for Simultaneous Monitoring of Shaft Rotational Speed and Axial Vibration Using Non-Projection Composite Fringe Pattern. Mech. Syst. Signal Process. 2019, 120, 765–776. [Google Scholar] [CrossRef]
  46. Pickering, S.; Tanaka, S.; Yamada, K. The Impact of Municipal Mergers on Local Public Spending: Evidence from Remote-Sensing Data. J. East Asian Stud. 2020, 20, 243–266. [Google Scholar] [CrossRef]
  47. Basten, C.; von Ehrlich, M.; Lassmann, A. Income Taxes, Sorting and the Costs of Housing: Evidence from Municipal Boundaries in Switzerland. Econ. J. 2017, 127, 653–687. [Google Scholar] [CrossRef]
  48. He, S.; Bastani, F.; Jagwani, S.; Alizadeh, M.; Balakrishnan, H.; Chawla, S.; Elshrif, M.M.; Madden, S.; Sadeghi, M.A. Sat2graph: Road Graph Extraction through Graph-Tensor Encoding. In Proceedings of the Computer Vision–ECCV 2020: 16th European Conference, Glasgow, UK, 23–28 August 2020; Proceedings, Part XXIV 16. pp. 51–67. [Google Scholar]
  49. He, K.; Gkioxari, G.; Dollár, P.; Girshick, R. Mask R-Cnn. In Proceedings of the IEEE International Conference on Computer Vision, Venice, Italy, 22–29 October 2017; pp. 2961–2969. [Google Scholar]
  50. Ren, S.; He, K.; Girshick, R.; Sun, J. Faster R-Cnn: Towards Real-Time Object Detection with Region Proposal Networks. Adv. Neural Inf. Process. Syst. 2015, 28, 161–243. [Google Scholar] [CrossRef]
  51. Kirillov, A.; Wu, Y.; He, K.; Girshick, R. PointRend: Image Segmentation as Rendering. arXiv 2020, arXiv:1912.08193. [Google Scholar]
  52. Henriksson, G.; Weitoft, G.R.; Allebeck, P. Associations between Income Inequality at Municipality Level and Health Depend on Context–A Multilevel Analysis on Myocardial Infarction in Sweden. Soc. Sci. Med. 2010, 71, 1141–1149. [Google Scholar] [CrossRef]
  53. Esmailpour, N.; Dehghan Dehkordi, E.; Abdali, S.; Ilderemi, H. An Overview on Income Experiences and Procedures of Municipalities with an Emphasis on Their Sustainability. Eur. Online J. Nat. Soc. Sci. Proc. 2015, 4, 227. [Google Scholar]
  54. Viani, A.; Orusa, T.; Mandola, M.L.; Robetto, S.; Belvedere, M.; Renna, G.; Scala, S.; Borgogno-Mondino, E.; Orusa, R. Tick’s Suitability Habitat Maps and Tick-Host Relationships in Wildlife. A One Health Approach Based on Multitemporal Remote Sensed Data, Entropy and Meta® Population Dataset in Aosta Valley, NW Italy. In Proceedings of the GeoVet 2023 International Conference, Teramo, Italy, 19–21 September 2023. [Google Scholar]
  55. Büttner, G. CORINE Land Cover and Land Cover Change Products. In Land Use and Land Cover Mapping in Europe; Springer: Berlin/Heidelberg, Germany, 2014; pp. 55–74. [Google Scholar]
  56. Feranec, J.; Soukup, T.; Hazeu, G.; Jaffrain, G. European Landscape Dynamics: CORINE Land Cover Data; CRC Press: Boca Raton, FL, USA, 2016. [Google Scholar]
  57. ESCAP. Producing Land Cover Change Maps and Statistics: Step by Step Guide on the Use of QGIS and RStudio; 2020 ESA; ESCAP: Frascati, Italy, 2020. [Google Scholar]
  58. Comber, A.J.; Wadsworth, R.; Fisher, P. Using Semantics to Clarify the Conceptual Confusion between Land Cover and Land Use: The Example of ‘Forest’. J. Land Use Sci. 2008, 3, 185–198. [Google Scholar] [CrossRef]
  59. Orusa, T.; Orusa, R.; Viani, A.; Carella, E.; Borgogno Mondino, E. Geomatics and EO Data to Support Wildlife Diseases Assessment at Landscape Level: A Pilot Experience to Map Infectious Keratoconjunctivitis in Chamois and Phenological Trends in Aosta Valley (NW Italy). Remote Sens. 2020, 12, 3542. [Google Scholar] [CrossRef]
  60. Conrad, O.; Bechtel, B.; Bock, M.; Dietrich, H.; Fischer, E.; Gerlitz, L.; Wehberg, J.; Wichmann, V.; Böhner, J. System for Automated Geoscientific Analyses (SAGA) v. 2.1.4. Geosci. Model Dev. 2015, 8, 1991–2007. [Google Scholar] [CrossRef]
  61. QGIS Development Team. QGIS Geographic Information System. Open Source Geospatial Foundation Project. 2018. Available online: https://qgis.org/en/site/ (accessed on 7 January 2024).
  62. Grizonnet, M.; Michel, J.; Poughon, V.; Inglada, J.; Savinaud, M.; Cresson, R. Orfeo ToolBox: Open Source Processing of Remote Sensing Images. Open Geospat. Data Softw. Stand. 2017, 2, 15. [Google Scholar] [CrossRef]
  63. Inglada, J.; Christophe, E. The Orfeo Toolbox Remote Sensing Image Processing Software. In Proceedings of the 2009 IEEE International Geoscience and Remote Sensing Symposium, Cape Town, South Africa, 12–17 July 2009; Volume 4, p. IV-733. [Google Scholar]
  64. Racine, J.S. RStudio: A Platform-Independent IDE for R and Sweave; JSTOR: New York, NY, USA, 2012. [Google Scholar]
  65. Homewood, K.; Lambin, E.F.; Coast, E.; Kariuki, A.; Kikula, I.; Kivelia, J.; Said, M.; Serneels, S.; Thompson, M. Long-Term Changes in Serengeti-Mara Wildebeest and Land Cover: Pastoralism, Population, or Policies? Proc. Natl. Acad. Sci. USA 2001, 98, 12544–12549. [Google Scholar] [CrossRef] [PubMed]
  66. Seto, K.C.; Kaufmann, R.K. Modeling the Drivers of Urban Land Use Change in the Pearl River Delta, China: Integrating Remote Sensing with Socioeconomic Data. Land Econ. 2003, 79, 106–121. [Google Scholar] [CrossRef]
  67. Olioso, A.; Soria, G.; Sobrino, J.; Duchemin, B. Evidence of Low Land Surface Thermal Infrared Emissivity in the Presence of Dry Vegetation. IEEE Geosci. Remote Sens. Lett. 2007, 4, 112–116. [Google Scholar] [CrossRef]
  68. Parker, D.C.; Manson, S.M.; Janssen, M.A.; Hoffmann, M.J.; Deadman, P. Multi-Agent Systems for the Simulation of Land-Use and Land-Cover Change: A Review. Ann. Assoc. Am. Geogr. 2003, 93, 314–337. [Google Scholar] [CrossRef]
  69. Kennedy, R.E.; Townsend, P.A.; Gross, J.E.; Cohen, W.B.; Bolstad, P.; Wang, Y.; Adams, P. Remote Sensing Change Detection Tools for Natural Resource Managers: Understanding Concepts and Tradeoffs in the Design of Landscape Monitoring Projects. Remote Sens. Environ. 2009, 113, 1382–1396. [Google Scholar] [CrossRef]
  70. Kennedy, R.E.; Yang, Z.; Gorelick, N.; Braaten, J.; Cavalcante, L.; Cohen, W.B.; Healey, S. Implementation of the LandTrendr Algorithm on Google Earth Engine. Remote Sens. 2018, 10, 691. [Google Scholar] [CrossRef]
  71. 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]
  72. Rose, R.A.; Byler, D.; Eastman, J.R.; Fleishman, E.; Geller, G.; Goetz, S.; Guild, L.; Hamilton, H.; Hansen, M.; Headley, R.; et al. Ten Ways Remote Sensing Can Contribute to Conservation. Conserv. Biol. 2015, 29, 350–359. [Google Scholar] [CrossRef]
  73. Tamiminia, H.; Salehi, B.; Mahdianpari, M.; Quackenbush, L.; Adeli, S.; Brisco, B. Google Earth Engine for Geo-Big Data Applications: A Meta-Analysis and Systematic Review. ISPRS J. Photogramm. Remote Sens. 2020, 164, 152–170. [Google Scholar] [CrossRef]
  74. Brisco, B.; Brown, R.; Hirose, T.; McNairn, H.; Staenz, K. Precision Agriculture and the Role of Remote Sensing: A Review. Can. J. Remote Sens. 1998, 24, 315–327. [Google Scholar] [CrossRef]
  75. Green, D.G.; Leishman, T. Computing and Complexity—Networks, Nature and Virtual Worlds. In Philosophy of Complex Systems; Elsevier: Amsterdam, The Netherlands, 2011; pp. 137–161. [Google Scholar]
  76. Karra, K.; Kontgis, C.; Statman-Weil, Z.; Mazzariello, J.C.; Mathis, M.; Brumby, S.P. Global Land Use/Land Cover with Sentinel 2 and Deep Learning. In Proceedings of the 2021 IEEE International Geoscience and Remote Sensing Symposium IGARSS, Brussels, Belgium, 11–16 July 2021; pp. 4704–4707. [Google Scholar]
  77. Agarwal, C. A Review and Assessment of Land-Use Change Models: Dynamics of Space, Time, and Human Choice. Hum. Ecol. Risk Assess. Int. J. 2002, 346, 102–118. [Google Scholar]
  78. Long, H.; Li, Y.; Liu, Y.; Woods, M.; Zou, J. Accelerated Restructuring in Rural China Fueled by ‘Increasing vs. Decreasing Balance’ Land-Use Policy for Dealing with Hollowed Villages. Land Use Policy 2012, 29, 11–22. [Google Scholar] [CrossRef]
  79. Atzberger, C. Advances in Remote Sensing of Agriculture: Context Description, Existing Operational Monitoring Systems and Major Information Needs. Remote Sens. 2013, 5, 949–981. [Google Scholar] [CrossRef]
  80. Yeh, A.G.; Li, X. An Integrated Remote Sensing and GIS Approach in the Monitoring and Evaluation of Rapid Urban Growth for Sustainable Development in the Pearl River Delta, China. Int. Plan. Stud. 1997, 2, 193–210. [Google Scholar] [CrossRef]
  81. Li, X.; Yeh, A.G.-O. Analyzing Spatial Restructuring of Land Use Patterns in a Fast Growing Region Using Remote Sensing and GIS. Landsc. Urban Plan. 2004, 69, 335–354. [Google Scholar] [CrossRef]
  82. Tarasenko, M.V. Current Status of the Russian Space Programme. Space Policy 1996, 12, 19–28. [Google Scholar] [CrossRef]
  83. Xiao, J.; Shen, Y.; Ge, J.; Tateishi, R.; Tang, C.; Liang, Y.; Huang, Z. Evaluating Urban Expansion and Land Use Change in Shijiazhuang, China, by Using GIS and Remote Sensing. Landsc. Urban Plan. 2006, 75, 69–80. [Google Scholar] [CrossRef]
  84. Berger, M.M. Happy Birthday, Constitution: The Supreme Court Establishes New Ground Rules for Land-Use Planning. Urban Lawyer 1988, 20, 735–800. [Google Scholar]
  85. 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, S27–S31. [Google Scholar] [CrossRef]
  86. Heistermann, M.; Müller, C.; Ronneberger, K. Land in Sight?: Achievements, Deficits and Potentials of Continental to Global Scale Land-Use Modeling. Agric. Ecosyst. Environ. 2006, 114, 141–158. [Google Scholar] [CrossRef]
  87. Rindfuss, R.R.; Stern, P.C. Linking Remote Sensing and Social Science: The Need and the Challenges. In People and Pixels: Linking Remote Sensing and Social Science; National Academy of Science: Washington, DC, USA, 1998; pp. 1–27. [Google Scholar]
  88. Orusa, T.; Viani, A.; Moyo, B.; Cammareri, D.; Borgogno-Mondino, E. Risk Assessment of Rising Temperatures Using Landsat 4–9 LST Time Series and Meta® Population Dataset: An Application in Aosta Valley, NW Italy. Remote Sens. 2023, 15, 2348. [Google Scholar] [CrossRef]
  89. Kohler, T.; Giger, M.; Hurni, H.; Ott, C.; Wiesmann, U.; von Dach, S.W.; Maselli, D. Mountains and Climate Change: A Global Concern. Mt. Res. Dev. 2010, 30, 53–55. [Google Scholar] [CrossRef]
  90. Crossman, N.D.; Bryan, B.A.; de Groot, R.S.; Lin, Y.-P.; Minang, P.A. Land Science Contributions to Ecosystem Services. Curr. Opin. Environ. Sustain. 2013, 5, 509–514. [Google Scholar] [CrossRef]
  91. Hagedoorn, L.C.; Koetse, M.J.; van Beukering, P.J. Estimating Benefits of Nature-Based Solutions: Diverging Values from Choice Experiments with Time or Money Payments. Front. Environ. Sci. 2021, 9, 686077. [Google Scholar] [CrossRef]
  92. Napp, S.; Petrić, D.; Busquets, N. West Nile Virus and Other Mosquito-Borne Viruses Present in Eastern Europe. Pathog. Glob. Health 2018, 112, 233–248. [Google Scholar] [CrossRef] [PubMed]
  93. Caminade, C.; Morse, A.P.; Fevre, E.M.; Mor, S.; Baylis, M.; Kelly-Hope, L. Forecasting the Risk of Vector-Borne Diseases at Different Time Scales: An Overview of the CLIMate SEnsitive DISease (CLIMSEDIS) Forecasting Tool Project for the Horn of Africa. In Proceedings of the EGU23, 25th EGU General Assembly, Vienna, Austria, 23–28 April 2023. [Google Scholar]
  94. Viani, A.; Orusa, T.; Divari, S.; Lovisolo, S.; Zanet, S.; Borgogno-Mondino, E.; Orusa, R.; Bollo, E. Bartonella spp. Distribution Assessment in Red Foxes (Vulpes Vulpes) Coupling Geospatially-Based Techniques. 2023, 76° Convegno SISVET, Bari. Available online: https://hdl.handle.net/2318/1938950 (accessed on 7 January 2024).
  95. Bolton, D.K.; Gray, J.M.; Melaas, E.K.; Moon, M.; Eklundh, L.; Friedl, M.A. Continental-Scale Land Surface Phenology from Harmonized Landsat 8 and Sentinel-2 Imagery. Remote Sens. Environ. 2020, 240, 111685. [Google Scholar] [CrossRef]
  96. Bhaduri, B.; Bright, E.; Coleman, P.; Dobson, J. LandScan. Geoinformatics 2002, 5, 34–37. [Google Scholar]
Figure 1. Area of study, involving the Region Valle d’Aosta, NW Italy. False color Sentinel-2 imagery 2022 meteorological summer composite (NIR, Red, Green).
Figure 1. Area of study, involving the Region Valle d’Aosta, NW Italy. False color Sentinel-2 imagery 2022 meteorological summer composite (NIR, Red, Green).
Land 13 00080 g001
Figure 2. Sen4MUN workflow.
Figure 2. Sen4MUN workflow.
Land 13 00080 g002
Figure 3. MAEs computed per each municipality in Aosta Valley Autonomous Region comparing Sen4MUN with OW considering (A) Urban and anthropic areas; (B) Road length; (C) Real estate units.
Figure 3. MAEs computed per each municipality in Aosta Valley Autonomous Region comparing Sen4MUN with OW considering (A) Urban and anthropic areas; (B) Road length; (C) Real estate units.
Land 13 00080 g003
Table 1. Ordinary Workflow equations adopted to retrieve territorial inputs for the financial contribution counting. It is worth to note that these equations are adopted also in Sen4MUN approach.
Table 1. Ordinary Workflow equations adopted to retrieve territorial inputs for the financial contribution counting. It is worth to note that these equations are adopted also in Sen4MUN approach.
DescriptionAlgorithm
Land Cover Areas (LCA)d = e + g + i (1)
where,
e = urban and anthropic area
g = semi-anthropic areas
i = sterile areas
d = municipality’s administrative boundaries
LCA Weights f = e × α (2)
where,
e = urban and anthropic area
α = weight (in this case = 3)
f = urban weighted area
 
h = g × β (3)
where,
g = semi-anthropic area
β = weight (in this case = 1.5)
h = semi-anthropic weighted area
 
l = i × γ (4)
where,
i = sterile area
γ = weight (in this case = 0.5)
l = sterile weighted area
Weighted areasm = f + h + l (5)
where,
f = urban weighted area
h = semi-anthropic weighted area
l = sterile weighted area
 
n = m i = 1 n m × 100 (6)
where,
m = conventional municipality area
i = 1 n m = sum of all the municipalities in the regional areas
n = conventional weighted municipality area
 
o = n × δ (7)
where,
n = conventional weighted municipality area
δ = weight (in this case = 50%)
o = sub-conventional weighted municipality area
Roads length q = p i = 1 n p × 100 (8)
where,
p = roads length
i = 1 n p = sum of all the municipalities in the regional areas
q = roads weighted length
 
r = q × ε (9)
where,
q = roads weighted length
ε = weight (in this case = 30%)
r = sub-conventional weighted roads length
Real estate units t = s i = 1 n s × 100 (10)
where,
s = real estate units
i = 1 n s = sum of all the municipalities real estate units
t = real estate weighted units
 
u = t × ζ (11)
where,
t = real estate weighted units
ζ = weight (in this case = 20%)
u = sub-conventional weighted real estate units
Table 2. Comparison between Land Cover EAGLE and Cerutti’s classes.
Table 2. Comparison between Land Cover EAGLE and Cerutti’s classes.
Land Cover EAGLE ClassCerutti Class
Urban and anthropic areasurban and anthropic area
Shrubland and transitional woodssemi-anthropic areas
Woody cropssemi-anthropic areas
Water surfacessterile areas
Water coursessterile areas
Needle-leaved forestssemi-anthropic areas
Broad-leaved forestssemi-anthropic areas
Mixed forests and moorssemi-anthropic areas
Permanent snow and icesterile areas
Natural grasslands and alpine pasturessemi-anthropic areas
Lawn pasturessemi-anthropic areas
Bare rockssterile areas
Discontinuous herbaceous vegetation of medium-low altitudesemi-anthropic areas
Sparse herbaceous vegetation at high altitudessemi-anthropic areas
Alpine wetlandssterile areas
Table 3. Overall MAEs and areas computed with both OW and Sen4MUN.
Table 3. Overall MAEs and areas computed with both OW and Sen4MUN.
Urban & Anthropic Areas (km2)Road Length (km)Real Estate Units
Sen4MUNOWMAESen4MUNOWMAESen4MUNOWMAE
103.992.10.161686.71626.40.81303,049293,214133
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Orusa, T.; Viani, A.; Borgogno-Mondino, E. Earth Observation Data and Geospatial Deep Learning AI to Assign Contributions to European Municipalities Sen4MUN: An Empirical Application in Aosta Valley (NW Italy). Land 2024, 13, 80. https://doi.org/10.3390/land13010080

AMA Style

Orusa T, Viani A, Borgogno-Mondino E. Earth Observation Data and Geospatial Deep Learning AI to Assign Contributions to European Municipalities Sen4MUN: An Empirical Application in Aosta Valley (NW Italy). Land. 2024; 13(1):80. https://doi.org/10.3390/land13010080

Chicago/Turabian Style

Orusa, Tommaso, Annalisa Viani, and Enrico Borgogno-Mondino. 2024. "Earth Observation Data and Geospatial Deep Learning AI to Assign Contributions to European Municipalities Sen4MUN: An Empirical Application in Aosta Valley (NW Italy)" Land 13, no. 1: 80. https://doi.org/10.3390/land13010080

APA Style

Orusa, T., Viani, A., & Borgogno-Mondino, E. (2024). Earth Observation Data and Geospatial Deep Learning AI to Assign Contributions to European Municipalities Sen4MUN: An Empirical Application in Aosta Valley (NW Italy). Land, 13(1), 80. https://doi.org/10.3390/land13010080

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