Next Article in Journal
Approach to Facilitating Geospatial Data and Metadata Publication Using a Standard Geoservice
Next Article in Special Issue
Highlighting Current Trends in Volunteered Geographic Information
Previous Article in Journal
Forecasting Urban Vacancy Dynamics in a Shrinking City: A Land Transformation Model
Previous Article in Special Issue
User-Generated Geographic Information for Visitor Monitoring in a National Park: A Comparison of Social Media Data and Visitor Survey
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Generating Up-to-Date and Detailed Land Use and Land Cover Maps Using OpenStreetMap and GlobeLand30

1
Department of Mathematics, University of Coimbra, Largo D. Dinis, 3001-501 Coimbra, Portugal
2
INESC Coimbra, Rua Sílvio Lima, Pólo II, 3030-290 Coimbra, Portugal
3
Department of Civil and Environmental Engineering, Politecnico di Milano, Piazza Leonardo da Vinci 32, 20133 Milan, Italy
4
Hellenic Military Academy, Leof. Varis—Koropiou, 16673 Vari, Greece
5
Hellenic Military Geographical Service, 4, Evelpidon Str., 11362 Athens, Greece
6
International Institute for Applied Systems Analysis (IIASA), Schlossplatz 1, A2361 Laxenburg, Austria
7
School of Rural and Surveying Engineering, National Technical University of Athens, 9 H. Polytechniou, 15780 Zografou, Greece
*
Author to whom correspondence should be addressed.
ISPRS Int. J. Geo-Inf. 2017, 6(4), 125; https://0-doi-org.brum.beds.ac.uk/10.3390/ijgi6040125
Submission received: 4 March 2017 / Revised: 17 April 2017 / Accepted: 17 April 2017 / Published: 22 April 2017
(This article belongs to the Special Issue Volunteered Geographic Information)

Abstract

:
With the opening up of the Landsat archive, global high resolution land cover maps have begun to appear. However, they often have only a small number of high level land cover classes and they are static products, corresponding to a particular period of time, e.g., the GlobeLand30 (GL30) map for 2010. The OpenStreetMap (OSM), in contrast, consists of a very detailed, dynamically updated, spatial database of mapped features from around the world, but it suffers from incomplete coverage, and layers of overlapping features that are tagged in a variety of ways. However, it clearly has potential for land use and land cover (LULC) mapping. Thus the aim of this paper is to demonstrate how the OSM can be converted into a LULC map and how this OSM-derived LULC map can then be used to first update the GL30 with more recent information and secondly, enhance the information content of the classes. The technique is demonstrated on two study areas where there is availability of OSM data but in locations where authoritative data are lacking, i.e., Kathmandu, Nepal and Dar es Salaam, Tanzania. The GL30 and its updated and enhanced versions are independently validated using a stratified random sample so that the three maps can be compared. The results show that the updated version of GL30 improves in terms of overall accuracy since certain classes were not captured well in the original GL30 (e.g., water in Kathmandu and water/wetlands in Dar es Salaam). In contrast, the enhanced GL30, which contains more detailed urban classes, results in a drop in the overall accuracy, possibly due to the increased number of classes, but the advantages include the appearance of more detailed features, such as the road network, that becomes clearly visible.

1. Introduction

Land cover is an essential climate variable [1] as it strongly influences current and future climate, particularly with rapid changes in the landscape due to human activities [2]. Baseline information on land use and land cover (LULC) is also a key input to many different types of models (e.g., economic land use models; species distribution modelling) and environmental assessments; the latter includes international conventions such as the Convention on Biodiversity, monitoring of the Sustainable Development Goals (SDGs), monitoring of environmental directives at the EU level and for evidence-based policy-making that is focused on issues such as land take.
Some countries and political organizations have their own detailed high resolution land cover products, e.g., the EU, the USA and Australia, but this is not the case for many developing countries. The Africover project from the Food and Agriculture Organization of the United Nations (FAO) produced LULC maps for a number of African countries [3] but these are now out of date and at a coarse resolution of 100 m. In other cases, products do exist but they are not openly shared. Land cover products are also available globally but these have traditionally been produced at coarse resolutions, i.e., between 300 m to 1 km, due to the satellite sensors used and because of their original intended purpose, i.e., to serve the needs of the climate modelling community. More recently, with the opening up of the Landsat archive, global land cover products have been created at a 30-m resolution. The FROM-GLC (Finer Resolution Observation and Monitoring of Global Land Cover) product has been created using fully automated approaches [4] while the GlobeLand30 (GL30) product has been developed using a combination of automated and manual methods [5]. However, these high resolution land cover products still require further validation at the international level to determine if they are useful for different applications.
Recently new data sources derived from citizen science, geo-crowdsourcing and applications of volunteered geographic information (VGI) [6] have shown great potential for deriving LULC information (see e.g., [7,8]). Among these, the OpenStreetMap (OSM) project, which consists of a freely-licensed, global geospatial database, has attracted considerable research interest from the academic community [9]. OSM represents a potentially valuable source of up-to-date LULC information due to the richness and variety of its map contents [10]. A recent development at the University of Heidelberg is OSMLandUse.org, which represents OSM data as a level 2 CORINE Land Cover (CLC) map [11]. Many OSM-based humanitarian mapping projects have produced large and detailed geospatial databases in the world’s less developed countries and where no other sources of accurate spatial data are available. The availability of detailed OSM data in these areas represents a source of up-to-date information that could be used to generate LULC maps that are more accurate or detailed than currently available global products. However, a major issue with OSM data is incomplete spatial coverage, so the key questions are: can we convert OSM data to LULC and then merge it with a high resolution global product such as GL30 to provide an updated LULC product for any location? Can we further enhance this updated product by adding more detail from OSM on urban areas? The aim of this paper is to present a methodology that can answer these questions, i.e., demonstrate how to create two LULC maps from OSM using the nomenclatures from the GL30, produced by the National Geomatic Center of China, and the Urban Atlas (UA), produced by the European Environment Agency. The former is used to produce an updated version of GL30 with more up-to-date information. The latter, after being merged with this updated version, is used to produce an enhanced version of GL30, which contains a more detailed characterization of urban areas. The procedure is tested on two study areas in places where OSM coverage is good but where detailed authoritative LULC maps are not readily available.

2. Data Sources

In this section, the two main datasets that are used in this paper are described, i.e., the OSM and the GL30. In addition, a brief description of the UA is provided since the nomenclature is used for deriving the LULC product from OSM.

2.1. OpenStreetMap (OSM)

The OSM project was initiated in 2004 based on Steve Coast’s vision that a global map of the world could be crowdsourced using the extensive local knowledge of people living and working in the areas to be mapped. After 13 years since the start of the project, OSM has indeed become the largest, most diverse, most complete and most up-to-date open access geospatial database in the world. Any object that is physically located in a position on the Earth’s surface (from forests, lakes and buildings to detailed features such as benches, drinking fountains and manholes) can be added to the OSM database. Rather than producing map products or cartographic outputs, the focus of the OSM project is to maintain a living geospatial database of the world. OSM data are available under the Open Database License (ODbL) [12], an inclusive license allowing anyone to freely copy, distribute, share and adapt the database provided that attribution is made to the OSM project and its contributors and that derived datasets are released under the same license.
The OSM database is comprised of vector data (i.e., the geometry and the associated attributes for each feature). The data model is made up of three kinds of geometric primitives: nodes (which are encoded points), ways (sets of nodes used for encoding linear and polygonal features) and relations (logical collections of two or more nodes, ways or other relations). Each of these primitives is associated with one or more attributes or tags. A tag is comprised of a key and a value. Although OSM contributors are free to use their own tags, there is a quasi-official collection of tags that has been established and agreed upon over the years by the global OSM community. The starting point of this collection is the Map Features Wiki page [13] that links to numerous other detailed pages on the OSM Wiki and explains the likely usage and use-case scenarios of each tag. For example, to identify a river, the tag “waterway = river” should be used, where “waterway” is the key and “river” is the value.
Anyone can contribute data to OSM provided that a user account is first created. At the time of writing (February 2017), more than 3 million users had registered in the OSM [14], although some studies have shown that most of the mapping is performed by a very small fraction of volunteers [15]. Users can contribute data to the project in three ways: (1) by digitizing features from satellite imagery (armchair mapping); (2) by inserting/uploading elements that have been physically surveyed in the field, for instance using a GPS (Global Positioning System) receiver (field mapping); and (3) by uploading other datasets that have been released under an open and compatible license. Access to OSM data via application programming interfaces (APIs) has also enabled the development of a wealth of other OSM-based software, websites, services and applications that make direct use of OSM data for a variety of purposes (see [16] for an extensive and up-to-date list).

Humanitarian Applications of OSM

The digital empowerment of both affected communities and volunteer networks are reshaping humanitarian responses in the twenty-first century [17]. Over the last few years, an increase in mapping and exploitation of OSM data has been associated with humanitarian applications in the less developed countries of the world, which were typically absent or underrepresented in OSM [18]. This is a result of the efforts of the Humanitarian OpenStreetMap Team (HOT), an international, nonprofit organization leading the collaborative mapping in OSM when a disaster strikes anywhere in the world. HOT was formed after the tragic earthquake in Haiti in 2010, when the remote mapping by volunteers from around the world became crucial in assisting the humanitarian aid work in the field [19]. Another initiative, the Missing Maps (MM) project [20], is a consortium of organizations (including HOT) that supports the creation of maps in the most vulnerable areas of the world that are subject, e.g., to epidemics, political crises, natural hazards and other risks, and where maps do not currently exist. Maps are created remotely from volunteers and then enriched by local organizations working in the field.
Remote volunteers often produce OSM data through the organization of so-called mapathons (literally “map marathons”), i.e., social events where experienced and novice contributors meet together in a room (e.g., at a university) and focus their efforts on a specific area that needs armchair mapping [21]. In contrast to mapathons, mapping parties are field mapping events where volunteers aim at improving the local OSM map of an area [22,23,24]. Both types of mapping actions may have humanitarian purposes. Latif et al. [25] describe the challenges of stimulating voluntary mapping efforts in the local community of a disaster-prone country like Bangladesh. Feinmann [26] presents an ambitious project from MM that is being undertaken in collaboration with Médecins Sans Frontières (MSF), HOT, and the British and American Red Cross. The aim is to add 200 M addresses to the OSM world map in two years via numerous mapathons held around the world. The OpenStreetMap Haiti project, which focused on the remote mapping of Haiti after the earthquake, is presented by Shemak [27], while Moeller and Furhmann [18] describe the mapathons organized by the American Red Cross and HOT (both partners of MM) to fight the 2014 Ebola outbreak in Western Africa. Ebrahim et al. [21] describe an experience of humanitarian mapathon for Swaziland co-organized with HOT and MM and performed by 200 10-year old children.
Humanitarian mapping in OSM has often proven to provide the only available source of map data and thus becomes a crucial source of information in emergency situations as well as for other preventive or planning needs such as census taking, risk mapping, control of outbreaks, etc. As an example, Figure 1 provides a meaningful comparison between Google Maps and OSM for the Kibera slum in Nairobi (Kenya) obtained using the Map Compare tool from Geofabrik [28]. While this slum is totally absent in the former, the latter shows the impressively high level of detail that has resulted from the work of the local volunteers in the frame of the Map Kibera project [29]. Situations like this one clearly show the potential of using OSM data to derive accurate map products like LULC maps, as discussed in Section 4.

2.2. GlobeLand30

GL30 has been produced by the National Geomatic Center of China (NGCC) as one of the first global 30-m land cover products produced from freely available Landsat imagery [5]. Baseline maps for the years 2000 and 2010 were created using more than 10,000 Landsat images to produce maps with 10 high-level land cover classes as follows: water bodies, wetland, artificial surfaces; cultivated land; permanent snow and ice; forest; shrubland; grassland; bareland and tundra. Pixel- and object-based classification approaches were combined to extract the classes in a hierarchical manner starting with water bodies and ending with tundra. These classes were then combined with ancillary layers such as other global, regional and national land cover products, OpenStreetMap and very high resolution satellite imagery from Google Earth, among others, via a customized web-based information service in order to label the objects and verify the results with knowledge related to specific geographical areas.
The GL30 product was validated globally for 2010 and was found to have an overall accuracy greater than 80% [5]. For artificial surfaces, it was also compared to other global land cover products for eight test areas around the world and accuracies ranged from 79% to 97%, outperforming both CLC and the FROM-GLC product of Gong et al. [4]. There have been a few other studies that have independently compared GL30 to authoritative products in Italy [30], Germany [31] and water bodies in Scandinavian countries [32] and more recently, land cover in Iran [33], with high agreements, i.e., greater than 78%. Other than the comparison for Iran, there has been little validation of this product outside Europe, and particularly in less developed countries, which may benefit from such a product if no other land cover maps are available. Table 1 shows the GL30 nomenclature along with the MMU, which varies by class.

2.3. Urban Atlas

The Global Monitoring for Environment and Security Urban Atlas, hereafter referred to as the Urban Atlas (UA) is a pan-European product that provides a detailed land cover and land use classification of European cities with a population greater than 100 K in the EU [34] although more cities were added in 2012 with populations greater than 50 K. The UA is intended as an input to evidence-based policy making, and allows for comparison between, and benchmarking of, major European cities. For example, studies related to the amount of greenspace or accessibility studies can be undertaken with such a product. More than 300 cities were mapped for the baseline year of 2006 using high resolution satellite imagery including other reference data such as commercial navigation products; cadastral and zoning data; local city maps; and aerial photographs, where available. Google Earth has also been used for verification of the locational and thematic accuracy of the classes.
The UA has a minimum mapping unit (MMU) of 0.25 ha (0.0025 km2) for artificial classes and 1 ha (0.01 km2) for other land cover classes. The positional accuracy is ±5 m with a minimum thematic accuracy of 85% for artificial surface classes and 80% for other land cover types [34]. The reason for introducing the UA in this paper is because of the detailed urban nomenclature that is available. This nomenclature is used to convert OSM data to a LULC map and then to produce an updated and enhanced version of the GL30 LULC map, as described in Section 4.
The UA is organized into four levels of hierarchy with increasing levels of detail. Table 2 lists the classes of the first three levels.

3. Study Areas

3.1. Kathmandu, Nepal

Between April and May 2015, Nepal was hit by a severe earthquake followed by many strong aftershocks, which caused widespread damage in the capital city of Kathmandu and the surrounding regions. Driven by HOT and locally coordinated by the nonprofit Kathmandu Living Labs (KLL) [35], the humanitarian mapping efforts have seen the participation of thousands of volunteers from around the world. The produced maps, which in many regions were completely missing before the event, were then used to support relief, rescue and rebuilding operations. Thanks to the coordination and training efforts of KLL, impressive field mapping was also performed locally, which has currently made the OSM Kathmandu map highly detailed and up-to-date. The 18 km × 18 km study area, which includes Kathmandu city, is shown in Figure 2.

3.2. Dar es Salaam, Tanzania

With a population of over 4 million people, Dar es Salaam in Tanzania is the largest city in East Africa [36]. The city is particularly prone to flooding, which every year during the rainy season causes many deaths and millions of dollars of damages, which could be prevented with adequate planning. However, since much of the city is comprised of unplanned and informal settlements, after the severe flood event of May 2015, the Dar Ramani Huria community-based mapping project [37] was started to raise awareness of flood resilience. They have trained university students and local citizens to create highly accurate OSM maps of the most flood-prone areas of the city, which were previously unmapped [38]. Maps are also delivered in printed form to the local governing bodies. The project, still ongoing, is managed by the Tanzanian Commission for Science and Technology, with partners including the City Council of Dar es Salaam, Ardhi University, University of Dar es Salaam, Buni Innovation Hub and HOT, with support from the Red Cross, the Global Facility for Disaster Reduction and Recovery and the World Bank. As a result of this project, which has also involved remote volunteers through the HOT channels (see e.g., [37]), Dar es Salaam is now one of the most densely mapped areas in the whole OSM database. Figure 3 shows the 18 km by 18 km study area in Dar es Salaam.

4. Methodology

4.1. Conversion from OSM to LULC

There is a large diversity of tags in OSM, and, as mentioned in Section 2.1, volunteers may even create their own tags. One of the keys available in OSM is “landuse”. Some of the values used for this key may have a direct conversion to LULC classes, such as “forest”, “vineyard” or “residential”. However, a large percentage of other proposed keys may also provide information on LULC, such as the keys “building”, “highway” or “amenity”. Hence, a methodology was developed to convert OSM features into LULC classes, which is shown schematically in Figure 4. This methodology was implemented for the set of tags listed in the OSM Map Features Wiki page that were considered to be the most relevant for extracting LULC information. However, as contributors are free to create new tags, whenever the conversion is made, an analysis of the available tags needs to be done for that specific study area and time of download, as additional important tags may exist. In a previous study by Fonte et al. [39], this conversion was already made for the UA and CLC nomenclatures.
The implemented methodology requires not only the identification of the nomenclature and the mapping of the OSM features into it, but also the use of rules and parameters that enable, among other aspects: (1) the conversion of linear features, such as roads and waterways (which are represented as polylines in OSM) into polygons; (2) the identification of the most appropriate class to be assigned to some tags which may be associated to more than one LULC class, such as grass (which may represent urban gardens, agricultural fields, pastures, or natural vegetation) or water (which may correspond to an urban lake inside a park, a river or the ocean); and (3) the solution of OSM data inconsistencies, which assign different non-compatible classes to the same location. Finally, as it may be convenient to create LULC maps with a particular MMU, the maps obtained through this process can be generalized, which requires the definition of additional rules. However, generalization was not considered in this paper.
From the technical perspective, the procedure described above was implemented in a Free and Open Source Software for Geospatial (FOSS4G) environment. The technologies used include GRASS GIS (Geographic Resources Analysis Support System Geographic Information System) [40], GDAL/OGR (Geospatial Data Abstraction Library) [41], and PostgreSQL [42] with its spatial extension PostGIS [43] as well as some specific tools to manage OSM data such as osm2pgsql [44] for PostgreSQL/PostGIS, and Osmosis [45]. Python was the main language used to integrate and manipulate the spatial processing of the data. The details about how these FOSS4G (Free and Open Source Software for Geospatial) technologies were used to accomplish the steps shown in Figure 4 are described in detail by Fonte et al. [39].
To be able to compare the LULC maps extracted from OSM with GL30, the tags associated with OSM features were also mapped to the classes of the GL30. Table 3 shows the mapping between the OSM tags considered for the study areas to both the UA and the GL30 nomenclature.

4.2. Comparison of OSM-Derived LULC Map with GlobeLand30

The comparison between the LULC map derived from OSM using the GL30 nomenclature and GL30 was made through: (1) computation of the percentage of unmapped regions in the OSM-derived map relative to the study area; and (2) a direct comparison of the maps for those regions with data in both maps (OSM-derived and GL30). In order to undertake this comparison, the vector LULC map derived from OSM data was first rasterized using a pixel size of 30 m to match the resolution of GL30. A confusion matrix was then calculated. The overall agreement of the overlapping regions was computed, as well as the marginal proportions of agreement [46]. In addition, a binary disagreement map was computed, which shows the spatial distribution of the classification agreement/disagreement.

4.3. Update of GlobeLand30 through the OSM-Derived LULC Map

The OSM-derived LULC map (with the GL30 nomenclature) was used to update the original GL30 product, which corresponds to the year 2010. To do this, the OSM-derived LULC map was rasterized to match the same 30 m grid as the GL30 for Kathmandu and Dar es Salaam. To minimize the loss of detail in the OSM data, the rasterization was performed using the prevalence method, i.e., the final 30 m pixel is assigned the class that has the largest vector presence in the pixel. The updated GL30 is then produced using the following rule, which is applied to each pixel:
  • If a pixel contains OSM-derived LULC data, update the GL30 with this information;
  • If a pixel contains no OSM-derived LULC data (i.e., a null value), then retain the original value of the GL30.
The original and the updated versions of GL30 were then compared by analyzing the difference in the relative proportions of the different LULC classes in order to understand what the additional contribution of OSM data are in the new updated map.

4.4. Enhancing GlobeLand30 through the Use of More Detailed OSM-Derived Data

In this section we go one step further and consider how to enhance the information content of the GL30 by adding additional detailed information that is available in OSM. The GL30 has 10 high-level classes as shown previously in Table 1, with only one class for urban areas (or artificial surfaces) while the UA nomenclature contains a number of detailed classes on urban areas (Table 2). Hence it is possible to enhance the artificial surfaces class of the GL30 with more detailed classes from the OSM-derived LULC map. The updated GL30 map (generated with the approach described in Section 4.3) is used as the starting layer. The non-urban classes, i.e., cultivated land (10), forest (20), grassland (30), shrubland (40), wetland (50) and water bodies (60), remain unchanged. Four new classes, which correspond to the UA level 2 nomenclature as listed in Table 2, are added to the GL30, replacing the single artificial surfaces class with these four new ones (Table 4). In the situation where a pixel contains no OSM-derived LULC data for the urban classes 80.1 to 80.4 but was classified as artificial surfaces (80) in the updated GL30 map, then the enhanced GL30 map pixel is assigned a value of 80.1 corresponding to the urban fabric class.

4.5. Validation of the Derived Maps

To validate the updated and enhanced GL30 maps and to compare them with the original GL30 LULC map, a random stratified sampling design was implemented to create a reference database of 50 points per class, considering the classes of the enhanced GL30 map (described in the previous section) as strata. The validation was then performed using the LACO-Wiki tool, which is a free, online portal that offers standardized land cover validation [47]. Interpretation in LACO-Wiki was implemented using what is referred to as blind validation, which means that the user is not shown the map class at each location but instead, interprets the very high resolution satellite imagery from Google Maps and Bing using the enhanced GL30 classification. For the validation of the GL30 and updated GL30, the four artificial surface classes (80.1 to 80.4) were grouped together. Once the validation was completed in LACO-Wiki, the reference database was downloaded for the calculation of confusion matrices and accuracy indices, which were computed using the methodology proposed by Card [48], where the area occupied by each class is considered in the computation.

5. Results

5.1. Conversion of OSM into LULC Classes

Figure 5 and Figure 6 show the original GL30 LULC maps, the OSM-derived LULC map using the GL30 nomenclature and the updated GL30 map for the two study areas.
Figure 7 shows binary disagreement maps computed for the pixels having non-null values in both the LULC maps obtained for the two study areas. For Kathmandu the regions with more agreement correspond to the main urban area and the forested areas, while most agricultural areas of the original GL30 are missing in OSM so the urban areas may have expanded since 2010. For Dar es Salaam there is a large agreement in the urban area, where most of the disagreement occurs in the areas of wetland, which are to a large extent absent from GL30, as well as in the classification of a natural region shown on the upper left hand side.
Table 5 and Table 6 show the confusion matrices between the original GL30 and the OSM-derived LULC map, the overall agreement (OA) and the row and column marginal proportions of agreement (RMPA and CMPA, respectively) for the study areas of Kathmandu and Dar es Salaam, respectively.
The results are considerably different for the two study areas. The overall agreement between the original GL30 and the OSM-derived LULC map is 56% for Kathmandu and 81% for Dar es Salaam. In the former case, the class having the highest RMPA and CMPA is 20 (Forest); class 10 (Cultivated land) shows a low RMPA but a very high CMPA, while the opposite happens for class 80 (Artificial surfaces). The remaining classes are almost absent in both LULC maps. In the case of Dar es Salaam, the highest agreement is found for class 80 (artificial surfaces). Class 20 (forest) and class 50 (wetlands) have a high CMPA and low RMPA (most of them are mapped as artificial surfaces in the original GL30) while class 30 (grassland), which is the one occupying the largest area in the original GL30, is almost absent in the OSM-derived map.
In general, classification disagreement between the pixels of the OSM-derived map and the GL30 could be due to: (1) classification errors in the original GL30 map; (2) LULC changes that occurred between 2010 (date of the original GL30 map) and the date of the OSM derived data, which were captured in the OSM data; (3) tagging errors in the original OSM data from which the LULC map is obtained; and (4) classification errors introduced by the conversion procedure.

5.2. Evaluation of the Updated GlobeLand30 Map

Table 7 and Table 8 show the relative distribution of classes and the percentage of non-null pixels for the same series of LULC maps, again for Kathmandu and Dar es Salaam, respectively.
As shown in Table 8, in the Dar es Salaam study area, there are pixels with null values in the original GL30, and as a consequence also in the updated GL30. This is due to the presence of the sea, which is not mapped in the GL30 LULC classification (see Figure 6). Conversely, for the Kathmandu study area, pixels with null values are only found in the OSM-derived LULC map. These pixels correspond to regions where no OSM data are available or no OSM data were used to generate the LULC map (because the tags were not considered appropriate for LULC purposes). However, it is worth noting that the percentage of non-null pixels in the OSM-derived LULC map is similar for both Kathmandu and Dar es Salaam (71.71% and 62.28%, respectively). In both study areas, the class with the largest percentage in the OSM-derived product is 80 (artificial surfaces), which reflects the fact that the highest mapping efforts in OSM are focused on built-up structures, mainly roads and buildings.
The effects of updating GL30 with the OSM-derived LULC map are also different for the two case studies. For Kathmandu, there is a significant increase in artificial surfaces from 23.42% to 48.64% (see Table 7), which results in an extension of the city area at the expense of cultivated land, as is clearly visible from Figure 5. For Dar es Salaam, the proportion among the classes is instead quite unchanged, except for a slight increase in class 50 (wetlands) and a decrease in class 30 (grassland). In particular, a visible difference between the original GL30 and the OSM-derived LULC map is the presence in the latter of the wetlands around the Msimbazi River and its tributaries. This area, which is the one typically suffering from flooding, is not mapped in GL30 due either to its MMU of 0.0009 km2 for rivers and 0.0729 km2 for wetlands or simply due to classification errors. Summarizing, the update of GL30 with the more detailed and small-scale LULC information derived from OSM mainly results in an extension of the urban area of Kathmandu and an increase in wetlands for Dar es Salaam.

5.3. Evaluation of the Enhanced GlobeLand30 Map

Figure 8 shows the enhanced GL30 maps for the two study areas produced after the application of the procedure described in Section 4.4. As is clearly visible from the figure, the main advantage of these maps compared to the updated GL30 (see the lower left part of Figure 5 and Figure 6) is the more detailed characterization of the urban areas according to the UA level 2 nomenclature. In particular, the details of the road network (included in the new class 80.2) are clearly visible. Supplementing Table 7 and Table 8, which list the proportions of the single urban class for Kathmandu and Dar es Salaam study areas, Table 9 shows the relative proportions for the 4 urban classes of the enhanced GL30. In both study areas, classes 80.3 and 80.4 are almost absent. The two remaining classes, 80.1 and 80.2, are roughly equally distributed in Kathmandu, with proportions of 53.56% and 45.89%, respectively. Conversely, class 80.1 is by far the most commonly occurring class in Dar es Salaam (proportion of 62.68% against 36.72% for class 80.2).
Finally, the binary maps in Figure 9 show the provenance of each pixel of the enhanced GL30 maps, i.e., they describe whether the pixel information comes from the original GL30 or the OSM-derived LULC map. It is easy to see that most of the pixel values (71.71% for Kathmandu and 74.87% for Dar es Salaam) originate from the OSM database. As already mentioned, these pixels correspond, in particular, to the urban and most populated areas, thus reflecting the primary focus of the OSM mapping effort.

5.4. Validation

The validation of the three maps was made using the reference database created as described in Section 4.5. Table 10 shows the user’s (UAc), producer’s (PAc) and overall accuracy (OAc) obtained for the original GL30 and updated GL30 for the two study areas.
In Section 2.2, the accuracy of GL30 was reported to vary between 79% and 97%. The results found here are much lower. This may be caused by the methodology used to create the reference database (i.e., the rules used to assign a class to the sample pixels), the fact that the reference database used in this study was created using only photo interpretation, and therefore in some cases the differentiation between, for example, agriculture and natural vegetation may be difficult, and also because the aerial and satellite images for the photo interpretation are more recent than GL30. For example, the Bing imagery used for Dar Es Salaam is dated 2013 (or occasionally 2010) while it is dated 2016 for Dar Es Salaam. Google Earth imagery is 2016 for both cities. The unavailability of imagery with the same date as the date of creation of GL30 prevents the estimation of the real amount of change as well as the identification of the effects on the accuracy assessment of the actual change and the influence of the accuracy assessment procedure used.
When comparing the results obtained for both study areas, the OAc increased from the original GL30 to the updated version, from 54% to 65% for Kathmandu (+11%) and from 61% to 69% for Dar es Salaam (+8%). Regarding the UAc and PAc, for Kathmandu, most differences are found for: class 60 (water bodies), since these are practically non-existent in the original GL30 and were mapped in the updated version; class 30 (grassland), as the UAc decreased from 33% to 4%, and class 40 (shrubland), where the opposite occurred, with an increase from 0% in the original GL30 to 38% in the updated version. This may result from actual changes in the vegetation or interpretation difficulties in the creation of the reference dataset. For Dar es Salaam, classes 50 (wetlands) and 60 (water bodies) are the most affected, as large regions around waterlines are classified in OSM as wetlands. These, however, do not correspond to regions with water in the satellite images, but likely to regions frequently flooded, as they also include some urban areas. This explains the relatively low values of UAc and large values of PAc for the updated GL30 map.
Table 11 shows the user’s (UAc), producer’s (PAc) and overall accuracy (OAc) obtained for the enhanced LULC map for the two study areas. For Kathmandu, the OAc of the enhanced LULC map is lower than the OAc of the updated map, but it is higher than the OAc of the original GL30. However, for Dar es Salaam the OAc is lower than the OAc of both the original GL30 and its updated version using OSM data. The decrease in accuracy is partly explained by an increase in the number of classes in the enhanced GL30 map since there are now four urban classes instead of one so the potential for error has increased. What has been gained is a more detailed urban characterization that clearly shows the road patterns in these cities.

6. Conclusions

In this paper a methodology is presented to demonstrate the use of OSM data for updating and enhancing the GL30 LULC product. GL30 exists at the global scale but it only has a small number of high level classes, e.g., only one urban class. In addition, the product is static from the year 2010 although a more recent version is currently being developed. OSM, on the other hand, is a dynamic data source created by citizens that is continually being updated. The results showed that, for the chosen study areas of Kathmandu, Nepal and Dar es Salaam, Tanzania, the conversion of OSM data into the GL30 classes, and the integration of the resultant maps with the original GL30, produced more up-to-date maps, with a higher OAc. A visual analysis of the updated products shows the expansion of urban areas and a more detailed map of water bodies. The richness of OSM data, which enables the use of a more detailed nomenclature, especially for the urban environment, was exploited by replacing the single urban class of GL30 with four more detailed urban classes using the UA level 2 nomenclature to create an enhanced map with a higher number of classes and thus greater detail. More than 70% of the content of the enhanced map for the two study areas is derived from the OSM database. Although the OAc of the resultant map was lower than the OAc of the original GL30 in the case of Dar es Salaam, a visual analysis shows clear identification of the road networks in both study areas.
The creation of hybrid maps, i.e., merging GL30 and OSM, has the advantage of obtaining a more detailed and updated product, especially related to urban development. However, the OSM data are not validated (even though there is a continuous validation made by the crowd), which may lead to some potential errors in the database. For example, there may be different tagging practices for some classes, in particular the natural classes, such as grassland (GL), shrubland (SL) or forest (F), which may create inconsistencies in the resulting LULC classes.
We have considered two “extreme” cases of mapping where many people have contributed, i.e., mainly remote mappers after the earthquake in Nepal or through an internationally funded project to address flooding in Dar es Salaam. However, this does not mean that areas in developing countries “need” to be hit by a disaster or funded by a project to be very well mapped in OSM. On the contrary, there are many areas worldwide that are well mapped simply thanks to local mapping groups or volunteers (i.e., without any contribution from remote mappers); see e.g., Abidjan city in the Ivory Coast (http://osm.org/go/ardJDOX). Thus, the same procedure can be applied anywhere that OSM data are available.
In the future we plan to create a Web service to apply the procedure outlined here in real time to any user defined area (without having to install software and run scripts). This service will be based on the Web processing service (WPS) specifications issued by the Open Geospatial Consortium (OGC). Such a service would be linked to an interface that would allow users to choose the optimal values of the parameters (e.g., the buffer values) so that users familiar with the area to be classified can use their local knowledge; otherwise default values would be used. This would bring an advantage compared to the current OSMLanduse.org service developed at the University of Heidelberg [11], which generates the final OSM-derived LULC map without allowing users to tune the parameter values used by the algorithm. We also plan to improve the algorithm with the use of other VGI, such as geotagged photographs from publicly available repositories such as Flickr and Panoramio (or its replacement) to complement LULC classification in areas where there is no OSM, or where the OSM data are not clearly mappable into a specific LULC target class.

Acknowledgments

The authors would like to acknowledge the support and contribution of EU COST Action TD1202 ‘Mapping and the Citizen Sensor’ (http://www.citizensensor-cost.eu) and EU COST Action IC1203 ‘European Network Exploring Research into Geospatial Information Crowdsourcing: software and methodologies for harnessing geographic information from the crowd (ENERGIC)’ (http://vgibox.eu/), and FCT project grant UID/MULTI/00308/2013 do INESC Coimbra. The authors would also like to acknowledge partial support from the EU-funded ERC CrowdLand project (No. 617754) and Horizon2020 LandSense project (No. 689812).

Author Contributions

C.C.F. and J.P. developed the automated methodology and J.P. created the OSM-derived LULC maps. C.C.F., M.M., V.A. and A.S. compared the OSM LULC map with GlobeLand30, created the reference database and calculated the evaluation measures for all of the different products outlined in the paper. L.S. coordinated the writing of the paper with the others, who all contributed sections and helped to edit the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
CLCCORINE Land Cover
CMPAColumn Marginal Proportions of Agreement
GL30GlobeLand30
HOTHumanitarian OpenStreetMap Team
KLLKathmandu Living Labs
LULCLand use land cover
MMUMinimum mapping unit
OGCOpen Geospatial Consortium
OSMOpenStreetMap
OAcOverall accuracy
PAcProducer’s accuracy
RMPARow marginal proportions of agreement
SDGsSustainable Development Goals
UAUrban Atlas
UAcUser’s accuracy
VGIVolunteered Geographic Information
WPSWeb processing service

References

  1. Hollmann, R.; Merchant, C.J.; Saunders, R.; Downy, C.; Buchwitz, M.; Cazenave, A.; Chuvieco, E.; Defourny, P.; de Leeuw, G.; Forsberg, R.; et al. The ESA Climate Change Initiative: Satellite Data Records for Essential Climate Variables. Bull. Am. Meteorol. Soc. 2013, 94, 1541–1552. [Google Scholar] [CrossRef]
  2. Feddema, J.J.; Oleson, K.W.; Bonan, G.B.; Mearns, L.O.; Buja, L.E.; Meehl, G.A.; Washington, W.M. The Importance of Land-Cover Change in Simulating Future Climates. Science 2005, 310, 1674–1678. [Google Scholar] [CrossRef] [PubMed]
  3. Di Gregorio, A.; Latham, J. AFRICOVER land cover classification and mapping project. In Land Use, Land Cover and Soil Sciences—Vol. I.; Verheye, W.H., Ed.; UNESCO-Encylopedia of Life Support Systems (EOLSS): Paris, France, 2009; pp. 236–251. [Google Scholar]
  4. Gong, P.; Wang, J.; Yu, L.; Zhao, Y.; Zhao, Y.; Liang, L.; Niu, Z.; Huang, X.; Fu, H.; Liu, S.; et al. Finer resolution observation and monitoring of global land cover: first mapping results with Landsat TM and ETM+ data. Int. J. Remote Sens. 2013, 34, 2607–2654. [Google Scholar] [CrossRef]
  5. Chen, J.; Chen, J.; Liao, A.; Cao, X.; Chen, L.; Chen, X.; He, C.; Han, G.; Peng, S.; Lu, M.; et al. Global land cover mapping at 30m resolution: A POK-based operational approach. ISPRS J. Photogramm. Remote Sens. 2015, 103, 7–27. [Google Scholar] [CrossRef]
  6. See, L.; Mooney, P.; Foody, G.; Bastin, L.; Comber, A.; Estima, J.; Fritz, S.; Kerle, N.; Jiang, B.; Laakso, M.; et al. Crowdsourcing, citizen science or Volunteered Geographic Information? The current state of crowdsourced geographic information. ISPRS Int. J. Geo-Inf. 2016, 5, 55. [Google Scholar] [CrossRef]
  7. Jokar Arsanjani, J.; Helbich, M.; Bakillah, M.; Hagenauer, J.; Zipf, A. Toward mapping land-use patterns from volunteered geographic information. Int. J. Geogr. Inf. Sci. 2013, 27, 2264–2278. [Google Scholar] [CrossRef]
  8. Antoniou, V.; Fonte, C.; See, L.; Estima, J.; Arsanjani, J.; Lupia, F.; Minghini, M.; Foody, G.; Fritz, S. Investigating the feasibility of geo-tagged photographs as sources of land cover input data. ISPRS Int. J. Geo-Inf. 2016, 5, 64. [Google Scholar] [CrossRef]
  9. Jokar Arsanjani, J.; Zipf, A.; Mooney, P.; Helbich, M. An introduction to OpenStreetMap in Geographic Information Science: Experiences, research, and applications. In OpenStreetMap in GIScience; Springer International Publishing: Cham, Switzerland, 2015; pp. 1–18. [Google Scholar]
  10. Fonte, C.; Minghini, M.; Antoniou, V.; See, L.; Brovelli, M.A.; Patriarca, J.; Milčinski, G. An automated methodology for converting OSM data into a land use/cover map. In Proceedings of the 6th International Conference on Cartography & GIS, Albena, Bulgaria, 13–17 June 2016. [Google Scholar]
  11. Zipf, A. Introducing OSMLanduse.org: OpenStreetMap Landuse Landcover (LULC) WebMap is Online. Available online: http://k1z.blog.uni-heidelberg.de/2016/11/28/introducing-osmlanduseorg-openstreetmap-landuse-landcover-lulc-webmap-is-online/ (accessed on 15 February 2017).
  12. OpenStreetMap Copyright and License. Available online: https://www.openstreetmap.org/copyright (accessed on 30 September 2016).
  13. OpenStreetMap Wiki Map Features. Available online: https://www.openstreetmap.org/copyright (accessed on 30 September 2016).
  14. OpenStreetMap Wiki Stats. Available online: http://wiki.openstreetmap.org/wiki/Stats (accessed on 7 February 2017).
  15. Neis, P.; Zielstra, D. Recent developments and future trends in volunteered geographic information research: The case of OpenStreetMap. Future Internet 2014, 6, 76–106. [Google Scholar] [CrossRef]
  16. Mooney, P.; Minghini, M. A review of OpenStreetMap data. In Mapping and the Citizen Sensor; Foody, G.M., See, L., Fritz, S., Fonte, C.C., Mooney, P., Olteanu-Raimond, A.-M., Antoniou, V., Eds.; Ubiquity Press: London, UK, 2017. [Google Scholar]
  17. Meier, P. New information technologies and their impact on the humanitarian sector. Int. Rev. Red Cross 2011, 93, 1239–1263. [Google Scholar] [CrossRef]
  18. Moeller, M.S.; Furhmann, S. Mapping the world: A new approach for Volunteered Geographic Information in the cloud. ISPRS Int. Arch. Photogramm. Remote Sens. Spat. Inf. Sci. 2015, XL-6/W1, 9–10. [Google Scholar] [CrossRef]
  19. Soden, R.; Palen, L. From crowdsourced mapping to community mapping: The post-earthquake work of OpenStreetMap Haiti. In Proceedings of the 11th International Conference on the Design of Cooperative Systems, Cham, Switzerland, 27–30 May 2014. [Google Scholar]
  20. OpenStreetMap Foundation Missing Maps. Available online: http://www.missingmaps.org/ (accessed on 31 October 2016).
  21. Ebrahim, M.; Minghini, M.; Molinari, M.E.; Torrebruno, A. MiniMapathon: Mapping the world at 10 years old. In Proceedings of the 8th Annual International Conference on Education and New Learning Technologies (EDULEARN 2016), Barcelona, Spain; 2016; pp. 4200–4208. [Google Scholar]
  22. Perkins, C.; Dodge, M. The potential of user-generated cartography: A case study of the OpenStreetMap project and Manchester mapping party. North West Geogr. 2008, 8, 19–32. [Google Scholar]
  23. Hristova, D.; Quattrone, G.; Mashhadi, A.; Capra, L. The life of the party: Impact of social mapping in OpenStreetMap. In Proceedings of the Seventh International AAAI Conference on Weblogs and Social Media, Boston, MA, USA, 8–11 July 2013. [Google Scholar]
  24. Mooney, P.; Minghini, M.; Stanley-Jones, F. Observations on an OpenStreetMap mapping party organised as a social event during an open source GIS conference. Int. J. Spat. Data Infrastruct. Res. 2015, 10, 138–150. [Google Scholar]
  25. Latif, S.; Islam, K.M.R.; Khan, M.M.I.; Ahmed, S.I. OpenStreetMap for the disaster management in Bangladesh. In Proceedings of the 2011 IEEE Conference on Open Systems, Langkawi, Malaysia, 25–28 September 2011. [Google Scholar]
  26. Feinmann, J. How MSF is mapping the world’s medical emergency zones. BMJ 2014, 349, 7540. [Google Scholar] [CrossRef] [PubMed]
  27. Shemak, A. The cartographic dimensions of humanitarianism: Mapping refugee spaces in post-earthquake Haiti. Cult. Dyn. 2014, 26, 251–275. [Google Scholar] [CrossRef]
  28. GeoFabrik Map Compare. Available online: http://tools.geofabrik.de/mc (accessed on 31 October 2016).
  29. Map Kibera Map Kibera. Available online: http://mapkibera.org/ (accessed on 30 September 2016).
  30. Brovelli, M.A.; Molinari, M.E.; Hussein, E.; Chen, J.; Li, R. The first comprehensive accuracy assessment of GlobeLand30 at a national level: Methodology and results. Remote Sens. 2015, 7, 4191–4212. [Google Scholar] [CrossRef]
  31. Jokar Arsanjani, J.; See, L.; Tayyebi, A. Assessing the suitability of GlobeLand30 for mapping land cover in Germany. Int. J. Digit. Earth 2016, 9, 873–891. [Google Scholar] [CrossRef]
  32. Ban, Y.; See, L.; Haas, J.; Jacob, A.; Fritz, S. Validation of the Water Layers of Global Land Cover Products Using Geo-Wiki and National Land Cover Maps. In Proceedings of the 35th EARSel Symposium, Stockholm, Sweden, 15–19 June 2015. [Google Scholar]
  33. Jokar Arsanjani, J.; Tayyebi, A.; Vaz, E. GlobeLand30 as an alternative fine-scale global land cover map: Challenges, possibilities, and implications for developing countries. Habitat Int. 2016, 55, 25–31. [Google Scholar] [CrossRef]
  34. Urban Atlas; European Environment Agency: Copenhagen, Denmark, 2010; DAT-107-en.
  35. Wiseman, A. Kathmandu Living Labs is working to respond in Nepal: help them out! Available online: https://hotosm.org/updates/2015-05-01_kathmandu_living_labs_is_working_to_respond_in_nepal_help_them_out (accessed on 31 October 2016).
  36. OpenStreetMap Wiki Dar es Salaam. Available online: http://wiki.openstreetmap.org/wiki/Dar_es_Salaam (accessed on 22 April 2017).
  37. Dar Ramani Huria. Available online: http://ramanihuria.org (accessed on 31 October 2016).
  38. Felix, B. Building Infrastructure with Mapillary in Dar es Salaam, Tanzania. Available online: http://blog.mapillary.com/update/2015/12/17/worldbank.html (accessed on 3 October 2016).
  39. Fonte, C.; Patriarca, J.; Minghini, M.; Antoniou, V.; See, L.; Brovelli, M. Using OpenStreetMap to Create Land Use and Land Cover Maps: Development of an Application. In Volunteered Geographic Information and the Future of Geospatial Data; Calazans Campelo, C.E., Bertolotto, M., Corcoran, P., Eds.; Advances in Geospatial Technologies, IGI Global: Hershey, PA, USA, 2017. [Google Scholar]
  40. GRASS Development Team Geographic Resources Analysis Support System (GRASS) Software, Version 7.2. Open Source Geospatial Foundation. Available online: http://grass.osgeo.org (accessed on 31 October 2016).
  41. OS Geo GDAL—Geospatial Data Abstraction Library. Available online: http://www.gdal.org/ (accessed on 31 October 2016).
  42. The PostgreSQL Global Development Group PostgreSQL. Available online: https://www.postgresql.org/ (accessed on 31 October 2016).
  43. PostGIS Project Steering Committee PostGIS. Available online: http://postgis.net/ (accessed on 31 October 2016).
  44. OpenStreetMap Wiki Osm2pgsql. Available online: http://wiki.openstreetmap.org/wiki/Osm2pgsql (accessed on 22 April 2017).
  45. OpenStreetMap Wiki Osmosis. Available online: http://wiki.openstreetmap.org/wiki/Osmosis (accessed on 31 October 2016).
  46. Stehman, S.V. Sampling designs for accuracy assessment of land cover. Int. J. Remote Sens. 2009, 30, 5243–5272. [Google Scholar] [CrossRef]
  47. See, L.; Perger, C.; Hofer, M.; Weichselbaum, J.; Dresel, C.; Fritz, S. LACO-WIKI: An open access online portal for land cover validation. ISPRS Ann. Photogramm. Remote Sens. Spat. Inf. Sci. 2015, II-3/W5, 167–171. [Google Scholar] [CrossRef]
  48. Card, D.H. Using known map category marginal frequencies to improve estimates of thematic map accuracy. Photogramm. Eng. Remote Sens. 1982, 48, 431–439. [Google Scholar]
Figure 1. Comparison between OpenStreetMap (OSM) (top) and Google Maps (bottom) showing the Kibera slum in Nairobi. Copyright: © OpenStreetMap contributors.
Figure 1. Comparison between OpenStreetMap (OSM) (top) and Google Maps (bottom) showing the Kibera slum in Nairobi. Copyright: © OpenStreetMap contributors.
Ijgi 06 00125 g001
Figure 2. Location of the study area in Kathmandu, Nepal. Copyright: © OpenStreetMap contributors.
Figure 2. Location of the study area in Kathmandu, Nepal. Copyright: © OpenStreetMap contributors.
Ijgi 06 00125 g002
Figure 3. Location of the study are in Dar es Salaam, Tanzania. Copyright: © OpenStreetMap contributors.
Figure 3. Location of the study are in Dar es Salaam, Tanzania. Copyright: © OpenStreetMap contributors.
Ijgi 06 00125 g003
Figure 4. Scheme of the conversion of OSM data into a land use and land cover (LULC) map.
Figure 4. Scheme of the conversion of OSM data into a land use and land cover (LULC) map.
Ijgi 06 00125 g004
Figure 5. Original GL30 (upper left), OSM-derived LULC map (upper right) and updated GL30 (lower left) for the Kathmandu study area.
Figure 5. Original GL30 (upper left), OSM-derived LULC map (upper right) and updated GL30 (lower left) for the Kathmandu study area.
Ijgi 06 00125 g005
Figure 6. Original GL30 (upper left), OSM-derived LULC map (upper right) and updated GL30 (lower left) for the Dar es Salaam study area.
Figure 6. Original GL30 (upper left), OSM-derived LULC map (upper right) and updated GL30 (lower left) for the Dar es Salaam study area.
Ijgi 06 00125 g006
Figure 7. Disagreement maps between the original GL30 and the OSM-derived LULC map for Kathmandu (left) and Dar es Salaam (right).
Figure 7. Disagreement maps between the original GL30 and the OSM-derived LULC map for Kathmandu (left) and Dar es Salaam (right).
Ijgi 06 00125 g007
Figure 8. Enhanced GL30 for Kathmandu (left) and Dar es Salaam (right) study areas.
Figure 8. Enhanced GL30 for Kathmandu (left) and Dar es Salaam (right) study areas.
Ijgi 06 00125 g008
Figure 9. The source of the data of the enhanced GL30 maps for Kathmandu (left) and Dar es Salaam (right) study areas.
Figure 9. The source of the data of the enhanced GL30 maps for Kathmandu (left) and Dar es Salaam (right) study areas.
Ijgi 06 00125 g009
Table 1. GlobeLand30 (GL30) nomenclature. MMU: minimum mapping unit.
Table 1. GlobeLand30 (GL30) nomenclature. MMU: minimum mapping unit.
Class CodeClass NameClass DescriptionMMU (km2)
10Cultivated land (CL)Arable land (cropland): dry land, paddy field, land for greenhouses, vegetable fields, artificial tame pastures, economic cropland in which shrub crops or herbaceous crops are planted, and land abandoned with the reclamation of arable land0.0324
20Forest (F)Broadleaved deciduous forest, evergreen broad-leaf forest, deciduous coniferous forest, evergreen coniferous forest, mixed broadleaf-conifer forest0.0576
30Grassland (GL)Typical grassland, meadow grassland, alpine grassland, desert grassland, grass0.09
40Shrubland (SL)Desert scrub, mountain scrub, deciduous and evergreen shrubs0.09
50Wetland (WL)Lake swamp, river flooding wetlands, seamarsh, shrub/forest wetlands, mangrove forest, tidal flats/salt marshes0.0729
60Water bodies (WB)Open water: lakes, reservoirs/fishponds, rivers0.0009 (Rivers)
0.0081 (Lakes)
70Tundra (T)Brush tundra, poaceae tundra, wet tundra, bare tundra, mixed tundraNot provided
80Artificial surfaces (AS)Settlement place, industrial and mining area, traffic facilities0.0144
90Bareland (BL)Saline-alkali land, sand, gravel, rock, microbiotic crust0.0324
100Permanent snow/ice (SI)Permanent snow, ice sheet and glacier0.0081
Table 2. Nomenclature of the Urban Atlas (UA).
Table 2. Nomenclature of the Urban Atlas (UA).
Level 1Level 2Level 3
1 Artificial surfaces (AS)1.1 Urban fabric (UF)1.1.1 Continuous urban fabric
1.1.2 Discontinuous urban fabric
1.1.3 Isolated structures
1.2 Industrial, commercial, public, military, private and transport units (ICP)1.2.1 Industrial, commercial, public, military and private units
1.2.2 Road and rail network and associated land
1.2.3 Port areas
1.2.4 Airports
1.3 Mine, dump and construction sites (MDC)1.3.1 Mineral extraction and dump sites
1.3.3 Construction sites
1.3.4 Land without current use
1.4 Artificial non-agricultural vegetated areas (AV)1.4.1 Green urban areas
1.4.2 Sports and leisure facilities
2 Agricultural areas, semi-natural areas and wetlands (A)
3 Forests (F)
5 Water (W)
Table 3. Combined nomenclature of UA and GL30.
Table 3. Combined nomenclature of UA and GL30.
OSMUA ClassGL30 Class
KeyKey Values
Amenity (polygon)animal_shelter, arts_centre, bank, bar, brothel, cafe, car_rental, car_wash, casino, cinema, clinic, college, community_centre, courthouse, crematorium, crypt, dentist, dive_centre, driving_school, embassy, fast_food, ferry_terminal, fire_station, food_court, fuel, grave_yard, gym, hospital, internet_cafe,1.2 (ICP)80 (AS)
Amenity (point)bar, hospital, pub, restaurant, school, university1.2 (ICP)80 (AS)
building (polygon)apartments, house, garage, garages, residential, yes, “null”1.1 (UF)80 (AS)
cathedral, chapel, church, civic, commercial, hangar, hospital, hotel, industrial, kiosk, mosque, office, public, retail, school, shrine, stadium, synagogue, temple, train_station, transportation, warehouse1.2 (ICP)80 (AS)
building (point)hotel, office, commercial, hospital, industrial, retail, warehouse, cathedral, chapel, church, mosque, temple, synagogue, school, stadium, train_station, transportation, public, shrine, civic, hangar, kiosk1.2 (ICP)80 (AS)
highway (line)living_street, motorway, motorway_link, pedestrian, primary, primary_link, residential, road, secondary, secondary_link, service, tertiary, tertiary_link, trunk, trunk_link, unclassified1.2 (ICP)80 (AS)
landuse (polygon)residential1.1 (UF)80 (AS)
cemetery, commercial, industrial, military, railway, retail1.2 (ICP)80 (AS)
brownfield, construction, depot, landfill, quarry, 1.3 (MDC)80 (AS)
scrub1.3 (MDC)40 (SL)
beer-garden1.4 (AV)80 (AS)
flowers, grass, greenfield, meadow, plants, recreation_ground, recreational_area, recreational, village_green1.4 (AV)30 (GL)
allotments, farm, farmland, farmyard, grass, greenhouse-horticulture, meadow, orchard, vineyard, allotments, farm, farmland, farmyard, meadow2.0 (A)10 (CL)
grass, meadow2.0 (A)30 (GL)
forest3.0 (F)20 (F)
reservoir, pond5.0 (W)60 (WB)
leisure (polygon)adult_gaming_centre, amusement_arcade, dance, hackerspace, ice_rink, sports_centre, stadium, swimming_pool1.2 (ICP)80 (AS)
dog_park, garden, park, golf_course, miniature_golf, pitch, playground, summer_camp, track, track, water_park1.4 (AV)30 (GL)
marina1.4 (AV)80 (AS)
nature_reserve2.0 (A) or 3.0 (F)20 (F), 30 (GL), 40 (SL), 50 (WL) or 90 (BL)
beach_resort2.0 (A)90 (BL)
swimming_area5.0 (W)60 (WB)
natural (polygon)forest (if area < 10,000 m2), grassland, park1.4 (AV)20 (F)
grassland2.0 (A)30 (GL)
fell, grassland, heath, scrub2.0 (A)90 (BL), 30 (GL) or 40 (SL)
beach, sand, shingle, bare_rock, scree, glacier2.0 (A)90 (BL)
mud, wetland2.0 (A)50 (WL)
forest, wood3.0 (F)20 (F)
bay, water, riverbank5.0 (W)60 (WB)
railway (line)rail1.2 (ICP)80 (AS)
waterway (line)river, riverbank, stream5.0 (W)60 (WB)
dock1.2 (ICP)80 (AS)
Table 4. Nomenclature used in the enhanced GL30 map.
Table 4. Nomenclature used in the enhanced GL30 map.
Class CodeNameMapping Rules Applied to Each Pixel
80.1Urban fabric (UF)If a pixel contains OSM-derived LULC data of class 80.1 to 80.4, then change updated GL30 to this class
ELSE
If a pixel was assigned to class 80 and contains no OSM-derived LULC data (i.e., a value of null), then assign a value of 80.1
80.2Industrial, commercial, public, military, private and transport units (ICP)
80.3Mine, dump and construction sites (MDC)
80.4Artificial non-agricultural vegetated areas (AV)
10Cultivated land (CL)If a pixel was not assigned to any of the classes 80.1 to 80.4 and belongs to one of these classes, then the class value remains unchanged
20Forest (F)
30Grassland (GL)
40Shrubland (SL)
50Wetland (WL)
60Water bodies (WB)
Table 5. Confusion matrix and indices of agreement for the Kathmandu study area: row marginal proportion of agreement (RMPA), column marginal proportion of agreement (CMPA) and overall agreement (OA).
Table 5. Confusion matrix and indices of agreement for the Kathmandu study area: row marginal proportion of agreement (RMPA), column marginal proportion of agreement (CMPA) and overall agreement (OA).
GL30 LULC Map
10 (CL)20 (F)30 (GL)40 (SL)50 (WL)60 (WB)80 (AS)RMPAi
OSM-derived LULC map10 (CL)14,048 2906077960.94
20 (F)777648,8231261498005680.83
30 (GL)499101078310
40 (SL)224420300140.01
50 (WL)40000000
60 (WB)236211000316770.01
80 (AS)75,68923225134405058,1110.43
CMPAi0.140.9500-0.330.95OA = 0.56
Table 6. Confusion matrix and indices of agreement for the Dar es Salaam study area.
Table 6. Confusion matrix and indices of agreement for the Dar es Salaam study area.
GL30 LULC Map
10 (CL)20 (F)30 (GL)40 (SL)50 (WL)60 (WB)80 (AS)RMPAi
OSM-derived LULC map10 (CL)4147022251640023770.01
20 (F)1721783850662044690.21
30 (GL)0111741410013200.11
40 (SL)00514008240.02
50 (WL)58641081130410061910,6420.25
60 (WB)073262361213720.00
80 (AS)109859666329402160162,7800.94
CMPAi 0.030.640.010.010.9500.89OA = 0.81
Table 7. Distribution of LULC classes and percentage of non-null pixels in the original GL30, the OSM-derived LULC map and the updated GL30 for Kathmandu study area.
Table 7. Distribution of LULC classes and percentage of non-null pixels in the original GL30, the OSM-derived LULC map and the updated GL30 for Kathmandu study area.
ClassOriginal GL30OSM-Derived LULC MapUpdated GL30
10 (CL)55.57%6.93%26.97%
20 (F)19.98%27.35%22.50%
30 (GL)0.11%0.62%0.50%
40 (SL)0.88%0.13%0.35%
50 (WL)000
60 (WB)0.04%1.43%1.04%
80 (AS)23.42%63.53%48.64%
non-null pixels100%71.71%100%
Table 8. Distribution of LULC classes and percentage of non-null pixels in the original GL30, the OSM-derived LULC map and the updated GL30 for Dar es Salaam study area.
Table 8. Distribution of LULC classes and percentage of non-null pixels in the original GL30, the OSM-derived LULC map and the updated GL30 for Dar es Salaam study area.
ClassOriginal GL30OSM-Derived LULC MapUpdated GL30
10 (CL)0.82%2.52%2.27%
20 (F)1.82%5.06%4.38%
30 (GL)9.59%0.79%5.09%
40 (SL)0.76%0.40%0.55%
50 (WL)1.68%7.99%6.16%
60 (WB)0.25%0.84%0.65%
80 (AS)85.08%82.40%80.90%
non-null pixels83.31%62.28%83.44%
Table 9. Distribution of urban LULC classes in the enhanced GL30 for Kathmandu and Dar es Salaam study areas.
Table 9. Distribution of urban LULC classes in the enhanced GL30 for Kathmandu and Dar es Salaam study areas.
Class CodeName% of Class 80 Kathmandu% of Class 80 Dar es Salaam
80.1Urban fabric (UF)53.56%62.68%
80.2Industrial, commercial, public, military, private and transport units (ICP)45.89%36.71%
80.3Mine, dump and construction sites (MDC)0.10%0.10%
80.4Artificial non-agricultural vegetated areas (AV)0.45%0.51%
Table 10. User’s accuracy (UAc), producer’s accuracy (PAc) and overall accuracy (OAc) obtained for the GL30 and updated GL30 for both study areas.
Table 10. User’s accuracy (UAc), producer’s accuracy (PAc) and overall accuracy (OAc) obtained for the GL30 and updated GL30 for both study areas.
Classes
10 (CL)20 (F)30 (GL)40 (SL)50 (WL)60 (WB)80 (AS)
Kathmandu GL30UAc28873300090
PAc9690150047
OAc54
Updated GL30UAc586043804973
PAc63921204673
OAc65
Dar es SalaamGL30UAc056344698065
PAc12352612294578
OAc61
Updated GL30UAc12352612294578
PAc2071141821089
OAc69
Table 11. User’s accuracy (UAc), producer’s accuracy (PAc) and overall accuracy (OAc) obtained for the enhanced LULC map for both study areas.
Table 11. User’s accuracy (UAc), producer’s accuracy (PAc) and overall accuracy (OAc) obtained for the enhanced LULC map for both study areas.
Classes
80.1 (UF)80.2 (ICP)80.3 (MDC)80.4 (AV)10 (CL)20 (F)30 (GL)40 (SL)50 (WL)60 (WB)
Kathmandu UAc80220445886438050
PAc547403719414051
OAc61
Dar es SalaamUAc6632840144226123046
PAc6358221587414108314
OAc49

Share and Cite

MDPI and ACS Style

Fonte, C.C.; Minghini, M.; Patriarca, J.; Antoniou, V.; See, L.; Skopeliti, A. Generating Up-to-Date and Detailed Land Use and Land Cover Maps Using OpenStreetMap and GlobeLand30. ISPRS Int. J. Geo-Inf. 2017, 6, 125. https://0-doi-org.brum.beds.ac.uk/10.3390/ijgi6040125

AMA Style

Fonte CC, Minghini M, Patriarca J, Antoniou V, See L, Skopeliti A. Generating Up-to-Date and Detailed Land Use and Land Cover Maps Using OpenStreetMap and GlobeLand30. ISPRS International Journal of Geo-Information. 2017; 6(4):125. https://0-doi-org.brum.beds.ac.uk/10.3390/ijgi6040125

Chicago/Turabian Style

Fonte, Cidália Costa, Marco Minghini, Joaquim Patriarca, Vyron Antoniou, Linda See, and Andriani Skopeliti. 2017. "Generating Up-to-Date and Detailed Land Use and Land Cover Maps Using OpenStreetMap and GlobeLand30" ISPRS International Journal of Geo-Information 6, no. 4: 125. https://0-doi-org.brum.beds.ac.uk/10.3390/ijgi6040125

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