[OSM-dev] RuntimeError: need at least two points for linestring (way_id=619453148)

[Sarah Hoffmann]

Hi,

On Thu, Nov 01, 2018 at 12:13:40PM +0900, koji higuchi wrote:
> *I tried to extract data from .osm.pbf file and write to shapefile as
> follows:*
> 
> import os, osmium, fiona
> 
> fi = 'europe-latest.osm.pbf'
> fo = 'europe-latest.shp'
> 
> drv = 'ESRI Shapefile'
> 
> crs = {'no_defs': True, 'ellps': 'WGS84', 'datum': 'WGS84', 'proj':
> 'longlat'}
> 
>  schema = {'geometry': 'LineString',
>                    'properties': {'id': 'float', 'name' : 'str', 'kind' :
> 'str'}}
> 
> outfile = fiona.open(fo, 'w', driver=drv, crs=crs, schema=schema)
> 
> geomfab = osmium.geom.GeoJSONFactory()
> 
> class ShapeConverter(osmium.SimpleHandler):
>       def way(self, w):
>            if 'place' in w.tags:
>                rec = {'geometry' : eval(geomfab.create_linestring(w)),
>                       'properties' : {'id' : float(w.id),
>                                      'name' : w.tags.get('name'),
>                                        'kind' : w.tags['place']}}
>                outfile.write(rec)
> 
> ShapeConverter().apply_file(fi, locations=True)
> 
> I got the following error after extracting several contents:
> 
>  rec = {'geometry' : eval(geomfab.create_linestring(w)),
> RuntimeError: need at least two points for linestring (way_id=619453148)
> 
> How could I skip that erroneous id and extract data for other working ids?

You need to do this manually yourself in the handler. I recommend
simply catching the exception as there are some other error
conditions besides too few points:

def way(self, w):
	if 'place' in w.tags:
		try:
			geom = geomfab.create_linestring(w)
		except:
			print("Skipping way with bad geometry")
			return

	  rec = {'geometry' : eval(geom),
           'properties' : {'id' : float(w.id),
           'name' : w.tags.get('name'),
           'kind' : w.tags['place']}}
    outfile.write(rec)

Kind regards

Sarah