OpenStreetMap+ Land Use / Land Cover classes and administrative regions of Europe

This dataset contains 23 30m resolution raster data of continental Europe land use / land cover classes extracted from OpenStreetMap, as well as administrative areas, and a harmonized building dataset based on OpenStreetMap and Copernicus HRL Imperviousness data.

The land use / land cover classes are:

1. buildings.commercial
2. buildings.industrial
3. buildings.residential
4. cemetery
5. construction.site
6. dump.site (landfill)
7. farmland
8. farmyard
9. forest
10. grass
11. greenhouse
12. harbour
13. meadow
14. military
15. orchard
16. quarry
17. railway
18. reservoir
19. road
20. salt
21. vineyard

The land use / land cover data was generated by extracting OSM vector layers from https://download.geofabrik.de/). These were then transformed into a 30 m density raster for each feature type. This was done by first creating a 10 m raster where each pixel intersecting a vector feature was assigned the value 100. These pixels were then aggregated to 10 m resolution by calculating the average of every 9 adjacent pixels. This resulted in a 0—100 density layer for the three feature types. Although the digitized building data from OSM offers the highest level of detail, its coverage across Europe is inconsistent. To supplement the building density raster in regions where crowd-sourced OSM building data was unavailable, we combined it with Copernicus High Resolution Layers (HRL) (obtained from https://land.copernicus.eu/pan-european/ high-resolution-layers), filling the non-mapped areas in OSM with the Impervious Built-up 2018 pixel values, which was averaged to 30 m. The probability values produced by the averaged aggregation were integrated in such a way that values between 0—100 refer to OSM (lowest and highest probabilities equal to 0 and 100 respectively), and the values between 101—200 refer to Copernicus HRL (lowest and highest probability equal to 200 and 101 respectively). This resulted in a raster layer where values closer to 100 are more likely to be buildings than values closer to 0 and 200. Structuring the data in this way allows us to select the higher probability building pixels in both products by the single boolean expression: Pixel > 50 AND pixel <150.

This dataset is part of the OpenStreetMap+ was used to pre-process the LUCAS/CORINE land use / land cover samples (https://doi.org/10.5281/zenodo.4740691) used to train machine learning models in Witjes et al., 2022 (https://doi.org/10.21203/rs.3.rs-561383/v4)

Each layer can be viewed interactively on the Open Data Science Europe data viewer at maps.opendatascience.eu.