def build(self):
try:
time_start = time.time()
ways = self.fetch_ways()
time_ways = time.time() - time_start
if not ways:
raise IncompletePolygonError('no ways found')
time_start = time.time()
rings = self.build_rings(ways)
time_rings = time.time() - time_start
if (imposm.config.imposm_multipolygon_max_ring
and len(rings) > imposm.config.imposm_multipolygon_max_ring):
log.warn('skipping relation %d with %d ways (%.1fms) and %d rings (%.1fms): too many rings',
self.relation.osm_id, len(ways), time_ways*1000, len(rings), time_rings*1000)
raise IncompletePolygonError('skipping too large multipolygon')
time_start = time.time()
self.build_relation_geometry(rings)
time_relations = time.time() - time_start
if time_ways + time_rings + time_relations > imposm.config.imposm_multipolygon_report:
log.warn('building relation %d with %d ways (%.1fms) and %d rings (%.1fms) took %.1fms',
self.relation.osm_id, len(ways), time_ways*1000, len(rings), time_rings*1000, time_relations*1000)
except InvalidGeometryError, ex:
log.debug(ex)
raise IncompletePolygonError(ex)
def fetch_ways(self):
ways = []
for member in self.relation.members:
# skip label nodes, relations of relations, etc
if member[1] != 'way': continue
way = self.ways_cache.get(member[0])
if way is None:
log.debug('way not found %s:%s', self.relation.osm_id,
member[0])
if imposm.config.import_partial_relations:
continue
else:
raise IncompletePolygonError(
'way not found %s:%s' %
(self.relation.osm_id, member[0]))
if way.partial_refs:
log.warn('multiple linestrings in way %s (relation %s)',
member[0], self.relation.osm_id)
raise IncompletePolygonError()
way.coords = self.fetch_way_coords(way)
if way.coords is None:
if not imposm.config.import_partial_relations:
raise IncompletePolygonError()
else:
ways.append(way)
return ways
def build_rings(self, ways):
rings = []
incomplete_rings = []
for ring in (Ring(w) for w in ways):
if ring.is_closed():
ring.geom = self.polygon_builder.build_checked_geom(ring, validate=self.validate_rings)
rings.append(ring)
else:
incomplete_rings.append(ring)
merged_rings = self.build_ring_from_incomplete(incomplete_rings)
if len(rings) + len(merged_rings) == 0:
raise IncompletePolygonError('linestrings from relation %s have no rings' % (self.relation.osm_id, ))
return rings + merged_rings
time_start = time.time()
self.build_relation_geometry(rings)
time_relations = time.time() - time_start
if time_ways + time_rings + time_relations > imposm.config.imposm_multipolygon_report:
log.warn('building relation %d with %d ways (%.1fms) and %d rings (%.1fms) took %.1fms',
self.relation.osm_id, len(ways), time_ways*1000, len(rings), time_rings*1000, time_relations*1000)
except InvalidGeometryError, ex:
log.debug(ex)
raise IncompletePolygonError(ex)
except IncompletePolygonError:
raise
except Exception, ex:
log.warn('error while building multipolygon:')
log.exception(ex)
raise IncompletePolygonError(ex)
class UnionRelationBuilder(RelationBuilderBase):
def build_relation_geometry(self, rings):
"""
Build relation geometry from rings.
"""
rings.sort(key=lambda x: x.geom.area, reverse=True)
# add/subtract all rings from largest
polygon = rings[0]
rel_tags = relation_tags(self.relation.tags, polygon.tags)
polygon.mark_as_inserted(rel_tags)
geom = polygon.geom