imposm3/geom/ring.go

127 lines
3.6 KiB
Go

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
}