-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathknn.go
67 lines (57 loc) · 1.21 KB
/
knn.go
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
package hdb
import (
"github.com/intdxdt/heap"
"github.com/intdxdt/mbr"
"github.com/TopoSimplify/node"
)
func predicate(_ *KObj) (bool, bool) {
return true, false
}
func (tree *Hdb) Knn(
query mbr.MBR, limit int, score func(*mbr.MBR, *KObj) float64,
predicates ...func(*KObj) (bool, bool)) []*node.Node {
var predFn = predicate
if len(predicates) > 0 {
predFn = predicates[0]
}
var nd = &tree.data
var result []*node.Node
var child *dbNode
var stop, pred bool
var queue = heap.NewHeap(kobjCmp, heap.NewHeapType().AsMin())
for !stop && (nd != nil) {
for i := range nd.children {
child = &nd.children[i]
var o = &KObj{
dbNode: child,
MBR: &child.bbox,
IsItem: len(child.children) == 0,
Distance: -1,
}
o.Distance = score(&query, o)
queue.Push(o)
}
for !queue.IsEmpty() && queue.Peek().(*KObj).IsItem {
var candidate = queue.Pop().(*KObj)
pred, stop = predFn(candidate)
if pred {
result = append(result, candidate.GetNode())
}
if stop {
break
}
if limit != 0 && len(result) == limit {
return result
}
}
if !stop {
var q = queue.Pop()
if q == nil {
nd = nil
} else {
nd = q.(*KObj).dbNode
}
}
}
return result
}