-
Notifications
You must be signed in to change notification settings - Fork 0
/
DiffusionDistanceTest.java
227 lines (197 loc) · 7.09 KB
/
DiffusionDistanceTest.java
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
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
package apps;
import Jama.EigenvalueDecomposition;
import algebra.Matrix;
import algebra.Vec2;
import algebra.Vector;
import apps.utils.MyFrame;
import graph.Graph;
import graph.GraphLaplacian;
import graph.KnnGraph;
import numeric.MyMath;
import twoDimEngine.BoxEngine;
import twoDimEngine.elements.Line2D;
import twoDimEngine.elements.Point2D;
import twoDimEngine.shaders.FillShader;
import twoDimEngine.shaders.PaintMethod2D;
import java.awt.*;
import java.awt.event.KeyEvent;
import java.awt.event.MouseEvent;
import java.util.*;
public class DiffusionDistanceTest extends MyFrame {
private static final int numPoints = 500;
private static final double length = 5;
private static final double epsilon = 0.1;
private int knn = 5;
private KnnGraph<Vec2> knnGraph;
private BoxEngine engine;
private PaintMethod2D shader;
private ArrayList<Vec2> points;
// index starts at 1
private int vertexCoordinate = 1;
private Matrix eigenVectors;
private Vector eigenValues;
private double heatTime = -1;
/**
* mouse coordinates
*/
private int mx, my;
public DiffusionDistanceTest(String title, int width, int height) {
super(title, width, height);
engine = new BoxEngine(width, height);
shader = new FillShader(engine);
engine.setBackGroundColor(Color.white);
points = new ArrayList<>();
Random r = new Random();
for (int i = 0; i < numPoints; i++) {
points.add(new Vec2(2 * length * r.nextDouble() - length, 2 * length * r.nextDouble() - length));
}
drawData();
init();
}
public static void main(String[] args) {
new DiffusionDistanceTest("Teste", 800, 500);
}
private void addData2Engine() {
for (int i = 0; i < points.size(); i++) {
Point2D e = new Point2D(points.get(i));
e.setColor(Color.black);
e.setRadius(0.01);
engine.addToList(e, shader);
}
}
private void drawData() {
engine.removeAllElements();
if (this.knnGraph == null) {
this.knnGraph = new KnnGraph<>(points, knn, (x, y) -> Vec2.diff(x, y).norm());
}
for (int i = 0; i < points.size(); i++) {
knnGraph.putVertexProperty(i + 1, "pos", points.get(i));
}
drawKnnGraph(knnGraph);
GraphLaplacian graphLaplacian = new GraphLaplacian(knnGraph, x -> Math.exp(-(x * x) / 25));
Matrix laplacian = graphLaplacian.getL();
EigenvalueDecomposition eigenvalueDecomposition = new EigenvalueDecomposition(new Jama.Matrix(laplacian.getMatrix()));
eigenVectors = new Matrix(eigenvalueDecomposition.getV().getArray());
eigenValues = new Vector(eigenvalueDecomposition.getRealEigenvalues());
eigenValues.applyFunction((x) -> -x);
drawHeatKernel();
engine.buildBoundingBoxTree();
engine.setCameraAuto(1.25);
}
private void drawHeatKernel() {
Vector eVertex = new Vector(knnGraph.getNumVertex());
eVertex.setX(vertexCoordinate, 1.0);
Matrix transpose = Matrix.transpose(eigenVectors);
Vector v = new Vector(eigenValues);
if (heatTime < 0) {
// formula is t = -log(epsilon) / eigenvalue_k, for positive eigenvalues. but here all eigenvalues are negative(see drawData function)
heatTime = Math.log(0.1) / v.getX((int) (Math.floor(0.04 * v.size()) + 1));
}
v.applyFunction((x) -> Math.exp(x * heatTime));
Vector f = eigenVectors.prodVector(Matrix.diag(v).prodVector(transpose.prodVector(eVertex)));
double max = f.getMax().getX();
double min = f.getMin().getX();
for (int i = 0; i < points.size(); i++) {
Point2D e = new Point2D(points.get(i));
e.setRadius(0.1);
double intensity = (f.getX(i + 1) - min) / (max - min);
float[] heatColor = getHeatColor(intensity);
e.setColor(new Color(heatColor[0], heatColor[1], heatColor[2]));
engine.addToList(e, shader);
}
}
public float[] getHeatColor(double intensity) {
float[] ans = new float[3];
ans[0] = (float) MyMath.clamp((10.0 / 4.0) * intensity, 0, 1);
ans[1] = (float) MyMath.clamp((10.0 / 4.0) * (intensity - 0.4), 0, 1);
ans[2] = (float) MyMath.clamp(5 * (intensity - 0.8), 0, 1);
return ans;
}
@Override
public void reShape() {
engine.setImageSize(widthChanged, heightChanged);
}
@Override
public void updateDraw() {
engine.clearImageWithBackground();
drawData();
engine.drawElements();
engine.paintImage(this.getGraphics());
}
private void drawKnnGraph(Graph knnGraph) {
for (Integer u : knnGraph.getVertexSet()) {
for (Integer v : knnGraph.getAdjVertex(u)) {
Line2D line = new Line2D(knnGraph.getVertexProperty(u, "pos"), knnGraph.getVertexProperty(v, "pos"));
line.setColor(Color.black);
engine.addToList(line, shader);
}
}
}
@Override
public void keyPressed(KeyEvent arg0) {
Map<Integer, Runnable> switchMap = new HashMap<>();
switchMap.put(KeyEvent.VK_PLUS, () -> {
heatTime += epsilon;
System.out.println(heatTime);
});
switchMap.put(KeyEvent.VK_MINUS, () -> {
heatTime -= epsilon;
heatTime = Math.max(0, heatTime);
System.out.println(heatTime);
});
final int dx = arg0.getKeyCode() - KeyEvent.VK_1;
if (dx > 0 && arg0.getKeyCode() <= KeyEvent.VK_9) {
knn = dx + 1;
knnGraph = null;
} else {
Optional.ofNullable(switchMap.get(arg0.getKeyCode())).ifPresent(Runnable::run);
}
}
@Override
public void keyReleased(KeyEvent arg0) {
// TODO Auto-generated method stub
}
@Override
public void keyTyped(KeyEvent arg0) {
// TODO Auto-generated method stub
}
@Override
public void mouseClicked(MouseEvent arg0) {
// TODO Auto-generated method stub
}
@Override
public void mouseEntered(MouseEvent arg0) {
// TODO Auto-generated method stub
}
@Override
public void mouseExited(MouseEvent arg0) {
// TODO Auto-generated method stub
}
@Override
public void mousePressed(MouseEvent e) {
mx = e.getX();
my = e.getY();
Vec2 mouse = new Vec2(engine.inverseCoordX(mx), engine.inverseCoordY(my));
double min = Double.MAX_VALUE;
int minIndex = -1;
for (int i = 0; i < points.size(); i++) {
double norm = Vec2.diff(points.get(i), mouse).norm();
if (min > norm) {
min = norm;
minIndex = i;
}
}
vertexCoordinate = minIndex + 1;
}
@Override
public void mouseReleased(MouseEvent arg0) {
// TODO Auto-generated method stub
}
@Override
public void mouseDragged(MouseEvent e) {
}
@Override
public void mouseMoved(MouseEvent arg0) {
// TODO Auto-generated method stub
}
}