-
Notifications
You must be signed in to change notification settings - Fork 9
/
KMeansClusteringExample.cpp
68 lines (57 loc) · 1.83 KB
/
KMeansClusteringExample.cpp
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
/*=========================================================================
*
* Copyright David Doria 2012 daviddoria@gmail.com
*
* Licensed under the Apache License, Version 2.0 (the "License");
* you may not use this file except in compliance with the License.
* You may obtain a copy of the License at
*
* http://www.apache.org/licenses/LICENSE-2.0.txt
*
* Unless required by applicable law or agreed to in writing, software
* distributed under the License is distributed on an "AS IS" BASIS,
* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
* See the License for the specific language governing permissions and
* limitations under the License.
*
*=========================================================================*/
#include "KMeansClustering.h"
#include <iostream>
Eigen::MatrixXd GenerateData();
int main(int, char *[])
{
Eigen::MatrixXd points = GenerateData();
KMeansClustering kmeans;
kmeans.SetK(2);
kmeans.SetPoints(points);
kmeans.SetInitMethod(KMeansClustering::KMEANSPP);
kmeans.SetRandom(false); // for repeatable results
kmeans.SetRandom(true); // for real, random results
kmeans.Cluster();
std::vector<unsigned int> labels = kmeans.GetLabels();
std::cout << "Resulting cluster ids:" << std::endl;
for(unsigned int i = 0; i < labels.size(); ++i)
{
std::cout << labels[i] << std::endl;
}
return EXIT_SUCCESS;
}
Eigen::MatrixXd GenerateData()
{
// Create 6 2-D points
Eigen::MatrixXd points(2, 6);
Eigen::VectorXd p = Eigen::VectorXd::Zero(2);
p[0] = 10; p[1] = 10;
points.col(0) = p;
p[0] = 10.1; p[1] = 10.1;
points.col(1) = p;
p[0] = 10.2; p[1] = 10.2;
points.col(2) = p;
p[0] = 5; p[1] = 5;
points.col(3) = p;
p[0] = 5.1; p[1] = 5.1;
points.col(4) = p;
p[0] = 5.2; p[1] = 5.2;
points.col(5) = p;
return points;
}