BE_OSM_ADMIN: BE_OSM_ADMIN: Administrative boundaries of Belgium based on...

Description References Examples

Description

BE_OSM_ADMIN: SpatialPolygonsDataFrame with administrative boundaries of Belgium. Extracted on 2015/11/19. It contains polygons of all administrative boundaries (levels 2, 4, 6, 7, 8, 9). Polygons are stored in degrees latitude/longitude (EPSG:4326 WGS 84). The data slot in the SpatialPolygonsDataFrame contains the following fields:

The OpenStreetMap data of Europe (europe-latest.osm.pbf as of 2015-11-18T22:23:02Z) was clipped with the bounding box 51.6, 49.4, 2.3, 6.5 using osmosis and the resulting file was imported with osm2pgsql in PostGIS and converted to a shapefile with pgsql2shp to obtain shape files of polygons of administrative boundaries of administrative levels 2, 4, 6, 7, 8, 9.

SQL: select osm_id, boundary, admin_level, name, way_area, way, id, way_off, rel_off, parts, members, rels.tags from (select * from planet_osm_polygon where boundary = 'administrative' and admin_level IN ('2', '4', '6', '7', '8', '9')) polygon left outer join (select * from planet_osm_rels) as rels on -polygon.osm_id = rels.id;

The data was further joined based on the INS code with aggregated TF_SOC_POP_STRUCT_2015 data from the BelgiumStatistics package and spatially grouped at different administrative levels which were retained from the TF_SOC_POP_STRUCT_2015 data namely: municipality / district / province / region / country as well as clipped to the Belgium boundary. All the administrative boundaries are made available in the BE_OSM_ADMIN dataset which can be found inside the package.

Mark that Brussels is not considered a a province and that in this dataset at the province level the region level information is taken from the TF_SOC_POP_STRUCT_2015 dataset.

References

http://wiki.openstreetmap.org/wiki/Tag:boundary=administrative, http://download.geofabrik.de/europe.html, https://en.wikipedia.org/wiki/World_Geodetic_System

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
## Not run: 
data(BE_OSM_ADMIN)
str(as.data.frame(BE_OSM_ADMIN))

library(sp)
plot(BE_OSM_ADMIN)
plot(subset(BE_OSM_ADMIN, admin.level %in% "2"))
plot(subset(BE_OSM_ADMIN, admin.level %in% "4"))
plot(subset(BE_OSM_ADMIN, admin.level %in% "6"))
plot(subset(BE_OSM_ADMIN, admin.level %in% "7"))
plot(subset(BE_OSM_ADMIN, admin.level %in% "8"))
plot(subset(BE_OSM_ADMIN, admin.level %in% "9"))

## End(Not run)

weRbelgium/BelgiumMaps.Admin documentation built on May 4, 2019, 2:04 a.m.