blob: c21a45af3841f1be457dcfe9ebe48727bc2d08f3 (
plain)
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
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
|
package de.lmu.ifi.dbs.elki.algorithm.outlier.lof.parallel;
/*
This file is part of ELKI:
Environment for Developing KDD-Applications Supported by Index-Structures
Copyright (C) 2015
Ludwig-Maximilians-Universität München
Lehr- und Forschungseinheit für Datenbanksysteme
ELKI Development Team
This program is free software: you can redistribute it and/or modify
it under the terms of the GNU Affero General Public License as published by
the Free Software Foundation, either version 3 of the License, or
(at your option) any later version.
This program is distributed in the hope that it will be useful,
but WITHOUT ANY WARRANTY; without even the implied warranty of
MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
GNU Affero General Public License for more details.
You should have received a copy of the GNU Affero General Public License
along with this program. If not, see <http://www.gnu.org/licenses/>.
*/
import de.lmu.ifi.dbs.elki.database.datastore.DataStore;
import de.lmu.ifi.dbs.elki.database.datastore.DoubleDataStore;
import de.lmu.ifi.dbs.elki.database.ids.DBIDRef;
import de.lmu.ifi.dbs.elki.database.ids.DBIDUtil;
import de.lmu.ifi.dbs.elki.database.ids.DoubleDBIDListIter;
import de.lmu.ifi.dbs.elki.database.ids.KNNList;
import de.lmu.ifi.dbs.elki.math.MathUtil;
import de.lmu.ifi.dbs.elki.parallel.Executor;
import de.lmu.ifi.dbs.elki.parallel.processor.AbstractDoubleProcessor;
import de.lmu.ifi.dbs.elki.parallel.variables.SharedDouble;
/**
* Processor for the "local reachability density" of LOF.
*
* @author Erich Schubert
* @since 0.7.0
*
* @apiviz.has Instance
*/
public class LRDProcessor extends AbstractDoubleProcessor {
/**
* KNN store
*/
private DataStore<? extends KNNList> knns;
/**
* k-distance store
*/
private DoubleDataStore kdists;
/**
* Constructor.
*
* @param knns k nearest neighbors
* @param kdists k distances
*/
public LRDProcessor(DataStore<? extends KNNList> knns, DoubleDataStore kdists) {
super();
this.knns = knns;
this.kdists = kdists;
}
@Override
public Instance instantiate(Executor master) {
return new Instance(master.getInstance(output));
}
/**
* Instance
*
* @author Erich Schubert
*/
private class Instance extends AbstractDoubleProcessor.Instance {
/**
* Constructor.
*
* @param output Output variable
*/
protected Instance(SharedDouble.Instance output) {
super(output);
}
@Override
public void map(DBIDRef id) {
KNNList knn = knns.get(id);
double lrd = 0.0;
int size = 0;
for(DoubleDBIDListIter n = knn.iter(); n.valid(); n.advance()) {
// Do not include the query object
if(DBIDUtil.equal(n, id)) {
continue;
}
lrd += MathUtil.max(kdists.doubleValue(n), n.doubleValue());
size += 1;
}
// Avoid division by 0:
output.set(lrd > 0 ? size / lrd : Double.POSITIVE_INFINITY);
}
}
}
|