package geom import ( "github.com/omniscale/imposm3/element" "github.com/omniscale/imposm3/geom/geos" ) type ring struct { ways []*element.Way refs []int64 nodes []element.Node geom *geos.Geom holes map[*ring]bool containedBy int area float64 outer bool inserted map[int64]bool } func (r *ring) isClosed() bool { return len(r.refs) >= 4 && r.refs[0] == r.refs[len(r.refs)-1] } func (r *ring) tryClose(maxRingGap float64) bool { return element.TryCloseWay(r.refs, r.nodes, maxRingGap) } func newRing(way *element.Way) *ring { r := ring{} r.ways = []*element.Way{way} r.refs = make([]int64, len(way.Refs)) r.nodes = make([]element.Node, len(way.Nodes)) r.containedBy = -1 r.holes = make(map[*ring]bool) copy(r.refs, way.Refs) copy(r.nodes, way.Nodes) return &r } func reverseRefs(refs []int64) { for i, j := 0, len(refs)-1; i < j; i, j = i+1, j-1 { refs[i], refs[j] = refs[j], refs[i] } } func reverseNodes(nodes []element.Node) { for i, j := 0, len(nodes)-1; i < j; i, j = i+1, j-1 { nodes[i], nodes[j] = nodes[j], nodes[i] } } func mergeRings(rings []*ring) []*ring { endpoints := make(map[int64]*ring) for _, ring := range rings { if len(ring.refs) < 2 { continue } left := ring.refs[0] right := ring.refs[len(ring.refs)-1] if origRing, ok := endpoints[left]; ok { // left node connects to.. delete(endpoints, left) if left == origRing.refs[len(origRing.refs)-1] { // .. right end origRing.refs = append(origRing.refs, ring.refs[1:]...) origRing.nodes = append(origRing.nodes, ring.nodes[1:]...) } else { // .. left end, reverse ring reverseRefs(origRing.refs) origRing.refs = append(origRing.refs, ring.refs[1:]...) reverseNodes(origRing.nodes) origRing.nodes = append(origRing.nodes, ring.nodes[1:]...) } origRing.ways = append(origRing.ways, ring.ways...) if rightRing, ok := endpoints[right]; ok && rightRing != origRing { // right node connects to another ring, close ring delete(endpoints, right) if right == rightRing.refs[0] { origRing.refs = append(origRing.refs, rightRing.refs[1:]...) origRing.nodes = append(origRing.nodes, rightRing.nodes[1:]...) } else { reverseRefs(rightRing.refs) origRing.refs = append(origRing.refs[:len(origRing.refs)-1], rightRing.refs...) reverseNodes(rightRing.nodes) origRing.nodes = append(origRing.nodes[:len(origRing.nodes)-1], rightRing.nodes...) } origRing.ways = append(origRing.ways, rightRing.ways...) right := origRing.refs[len(origRing.refs)-1] endpoints[right] = origRing } else { endpoints[right] = origRing } } else if origRing, ok := endpoints[right]; ok { // right node connects to.. delete(endpoints, right) if right == origRing.refs[0] { // .. left end origRing.refs = append(ring.refs[:len(ring.refs)-1], origRing.refs...) origRing.nodes = append(ring.nodes[:len(ring.nodes)-1], origRing.nodes...) } else { // .. right end, reverse ring reverseRefs(ring.refs) origRing.refs = append(origRing.refs[:len(origRing.refs)-1], ring.refs...) reverseNodes(ring.nodes) origRing.nodes = append(origRing.nodes[:len(origRing.nodes)-1], ring.nodes...) } origRing.ways = append(origRing.ways, ring.ways...) endpoints[left] = origRing } else { // ring is not connected (yet) endpoints[left] = ring endpoints[right] = ring } } uniqueRings := make(map[*ring]bool) for _, ring := range endpoints { uniqueRings[ring] = true } result := make([]*ring, 0, len(uniqueRings)) for ring, _ := range uniqueRings { result = append(result, ring) } return result }