locode: Use quadtree to find continent #830
11 changed files with 53 additions and 14 deletions
BIN
pkg/services/control/ir/service.pb.go
generated
BIN
pkg/services/control/ir/service.pb.go
generated
Binary file not shown.
BIN
pkg/services/control/ir/service_grpc.pb.go
generated
BIN
pkg/services/control/ir/service_grpc.pb.go
generated
Binary file not shown.
BIN
pkg/services/control/ir/types.pb.go
generated
BIN
pkg/services/control/ir/types.pb.go
generated
Binary file not shown.
BIN
pkg/services/control/service.pb.go
generated
BIN
pkg/services/control/service.pb.go
generated
Binary file not shown.
BIN
pkg/services/control/service_grpc.pb.go
generated
BIN
pkg/services/control/service_grpc.pb.go
generated
Binary file not shown.
BIN
pkg/services/control/types.pb.go
generated
BIN
pkg/services/control/types.pb.go
generated
Binary file not shown.
BIN
pkg/services/tree/service.pb.go
generated
BIN
pkg/services/tree/service.pb.go
generated
Binary file not shown.
BIN
pkg/services/tree/service_grpc.pb.go
generated
BIN
pkg/services/tree/service_grpc.pb.go
generated
Binary file not shown.
BIN
pkg/services/tree/types.pb.go
generated
BIN
pkg/services/tree/types.pb.go
generated
Binary file not shown.
|
@ -4,6 +4,8 @@ import (
|
||||||
"fmt"
|
"fmt"
|
||||||
"os"
|
"os"
|
||||||
|
|
||||||
|
"github.com/paulmach/orb/quadtree"
|
||||||
|
|
||||||
locodedb "git.frostfs.info/TrueCloudLab/frostfs-node/pkg/util/locode/db"
|
locodedb "git.frostfs.info/TrueCloudLab/frostfs-node/pkg/util/locode/db"
|
||||||
"github.com/paulmach/orb"
|
"github.com/paulmach/orb"
|
||||||
"github.com/paulmach/orb/geojson"
|
"github.com/paulmach/orb/geojson"
|
||||||
|
@ -36,22 +38,24 @@ func (db *DB) PointContinent(point *locodedb.Point) (*locodedb.Continent, error)
|
||||||
minDst float64
|
minDst float64
|
||||||
)
|
)
|
||||||
|
|
||||||
for _, feature := range db.features {
|
pointer := db.tree.Matching(planarPoint, func(p orb.Pointer) bool {
|
||||||
if multiPolygon, ok := feature.Geometry.(orb.MultiPolygon); ok {
|
return planar.PolygonContains(
|
||||||
if planar.MultiPolygonContains(multiPolygon, planarPoint) {
|
p.(*geojson.Feature).Geometry.(orb.Polygon),
|
||||||
|
planarPoint,
|
||||||
|
)
|
||||||
|
})
|
||||||
|
|
||||||
|
if pointer != nil {
|
||||||
|
continent = pointer.(*geojson.Feature).Properties.MustString(continentProperty)
|
||||||
|
}
|
||||||
|
|
||||||
|
if continent == "" {
|
||||||
|
for _, feature := range db.features {
|
||||||
|
distance := planar.DistanceFrom(feature.Geometry, planarPoint)
|
||||||
|
if minDst == 0 || minDst > distance {
|
||||||
|
minDst = distance
|
||||||
continent = feature.Properties.MustString(continentProperty)
|
continent = feature.Properties.MustString(continentProperty)
|
||||||
break
|
|
||||||
}
|
}
|
||||||
} else if polygon, ok := feature.Geometry.(orb.Polygon); ok {
|
|
||||||
if planar.PolygonContains(polygon, planarPoint) {
|
|
||||||
continent = feature.Properties.MustString(continentProperty)
|
|
||||||
break
|
|
||||||
}
|
|
||||||
}
|
|
||||||
distance := planar.DistanceFrom(feature.Geometry, planarPoint)
|
|
||||||
if minDst == 0 || minDst > distance {
|
|
||||||
minDst = distance
|
|
||||||
continent = feature.Properties.MustString(continentProperty)
|
|
||||||
}
|
}
|
||||||
}
|
}
|
||||||
|
|
||||||
|
@ -73,6 +77,38 @@ func (db *DB) init() error {
|
||||||
|
|
||||||
db.features = features.Features
|
db.features = features.Features
|
||||||
|
|
||||||
|
err = db.buildQuadtree()
|
||||||
|
if err != nil {
|
||||||
|
return fmt.Errorf("could not build quadtree: %w", err)
|
||||||
|
}
|
||||||
|
|
||||||
|
return nil
|
||||||
|
}
|
||||||
|
|
||||||
|
func (db *DB) buildQuadtree() error {
|
||||||
|
db.tree = quadtree.New(orb.Bound{
|
||||||
|
Min: orb.Point{-180, -180},
|
||||||
|
Max: orb.Point{180, 180},
|
||||||
|
})
|
||||||
|
|
||||||
|
for _, feature := range db.features {
|
||||||
|
if multiPolygon, ok := feature.Geometry.(orb.MultiPolygon); ok {
|
||||||
fyrchik marked this conversation as resolved
|
|||||||
|
for _, polygon := range multiPolygon {
|
||||||
|
newFeature := geojson.NewFeature(polygon)
|
||||||
|
newFeature.Properties = feature.Properties.Clone()
|
||||||
|
err := db.tree.Add(newFeature)
|
||||||
|
if err != nil {
|
||||||
|
return err
|
||||||
|
}
|
||||||
|
}
|
||||||
|
} else if _, ok := feature.Geometry.(orb.Polygon); ok {
|
||||||
|
err := db.tree.Add(feature)
|
||||||
|
if err != nil {
|
||||||
|
return err
|
||||||
|
}
|
||||||
|
}
|
||||||
|
}
|
||||||
|
|
||||||
return nil
|
return nil
|
||||||
}
|
}
|
||||||
|
|
||||||
|
|
|
@ -5,6 +5,7 @@ import (
|
||||||
"sync"
|
"sync"
|
||||||
|
|
||||||
"github.com/paulmach/orb/geojson"
|
"github.com/paulmach/orb/geojson"
|
||||||
|
"github.com/paulmach/orb/quadtree"
|
||||||
)
|
)
|
||||||
|
|
||||||
// Prm groups the required parameters of the DB's constructor.
|
// Prm groups the required parameters of the DB's constructor.
|
||||||
|
@ -31,6 +32,8 @@ type DB struct {
|
||||||
once sync.Once
|
once sync.Once
|
||||||
|
|
||||||
features []*geojson.Feature
|
features []*geojson.Feature
|
||||||
|
|
||||||
|
tree *quadtree.Quadtree
|
||||||
}
|
}
|
||||||
|
|
||||||
const invalidPrmValFmt = "invalid parameter %s (%T):%v"
|
const invalidPrmValFmt = "invalid parameter %s (%T):%v"
|
||||||
|
|
Loading…
Reference in a new issue
What about
var multiPolygon orb.MultiPolygon
?make(_, 0)
is usually never needed:nil
acts good enough as the default value.Fixed