esy-osm-shape 0.0.3

Convert OpenStreetMap objects to shapely objects.

esy-osm-shape

esy.osm.shape is a Python library to convert OpenStreetMap primitives to shapely objects.

Usage

The following examples operate on a historic dataset for Andorra from geofabrik. Let's download the dataset first:

>>> import urllib.request
>>> filename, headers = urllib.request.urlretrieve(
...     'https://download.geofabrik.de/europe/andorra-190101.osm.pbf',
...     filename='andorra.osm.pbf'
... )

Open the file and generate linestrings for each highway OpenStreetMap entry.

>>> import shapely, esy.osm.shape
>>> shape = esy.osm.shape.Shape('andorra.osm.pbf')
>>> highways = [
...     obj for obj in shape(lambda e: e.tags.get('highway') is not None)
...     if type(obj) is shapely.geometry.LineString
... ]

Using shapely objects it is also easy to compute geometric properties, like for example the length of the highways (note that the unit of this length is in longitude and latitude):

>>> sum(linestring.length for linestring in highways)
16.952160743015657

For slightly more details, jump to the documentation.