tommy/libpostal
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

[polygons] Only fixing polygons in cases with inner rings

Browse Source
This commit is contained in:
Al
2015-11-01 12:36:35 -05:00
parent 477300c061
commit 90773294b9
1 changed files with 14 additions and 1 deletions
Download Patch File Download Diff File Expand all files Collapse all files

15
scripts/geodata/polygons/reverse_geocode.py
View File

@@ -23,6 +23,7 @@ import sys
import tempfile import tempfile
from functools import partial from functools import partial
from shapely.geos import TopologicalError
this_dir = os.path.realpath(os.path.dirname(__file__)) this_dir = os.path.realpath(os.path.dirname(__file__))
sys.path.append(os.path.realpath(os.path.join(os.pardir, os.pardir))) sys.path.append(os.path.realpath(os.path.join(os.pardir, os.pardir)))
@@ -687,7 +688,19 @@ class OSMReverseGeocoder(RTreePolygonIndex):
if poly is None or not poly.bounds or len(poly.bounds) != 4: if poly is None or not poly.bounds or len(poly.bounds) != 4:
continue continue
# Figure out which outer polygon contains each inner polygon # Figure out which outer polygon contains each inner polygon
interior = [p2 for p2 in inner if poly.contains(p2)] try:
interior = [p2 for p2 in inner if poly.contains(p2)]
except TopologicalError:
poly = cls.fix_polygon(poly)
fixed_inner = []
for p in inner:
p = cls.fix_polygon(p)
if p.type != 'MultiPolygon':
fixed_inner.append(p)
else:
fixed_inner.extend(p)
inner = fixed_inner
interior = [p2 for p2 in inner if poly.contains(p2)]
if interior: if interior:
# Polygon with holes constructor # Polygon with holes constructor