/// <summary>
/// Performs the DBSCAN clustering algorithm.
/// </summary>
/// <param name="allPoints">Dataset</param>
/// <param name="epsilon">Desired region ball radius</param>
/// <param name="minPts">Minimum number of points to be in a region</param>
/// <param name="clusters">returns sets of clusters, renew the parameter</param>
public void ComputeClusterDbscan(MyCustomDatasetItem[] allPoints, double epsilon, int minPts, out HashSet <MyCustomDatasetItem[]> clusters)
{
DbscanPoint <MyCustomDatasetItem>[] allPointsDbscan = allPoints.Select(x => new DbscanPoint <MyCustomDatasetItem>(x)).ToArray();
int clusterId = 0;
for (int i = 0; i < allPointsDbscan.Length; i++)
{
DbscanPoint <MyCustomDatasetItem> p = allPointsDbscan[i];
if (p.IsVisited)
{
continue;
}
p.IsVisited = true;
DbscanPoint <MyCustomDatasetItem>[] neighborPts = null;
MN_RegionQuery(allPointsDbscan, p, epsilon, out neighborPts);
if (neighborPts.Length < minPts)
{
p.ClusterId = (int)ClusterIds.Noise;
}
else
{
clusterId++;
ExpandCluster(allPointsDbscan, p, neighborPts, clusterId, epsilon, minPts);
}
}
clusters = new HashSet <MyCustomDatasetItem[]>(
allPointsDbscan
.Where(x => x.ClusterId > 0)
.GroupBy(x => x.ClusterId)
.Select(x => x.Select(y => y.ClusterPoint).ToArray())
);
}
private void MN_RegionQuery(DbscanPoint <MyCustomDatasetItem>[] allPoints, DbscanPoint <MyCustomDatasetItem> point, double epsilon, out DbscanPoint <MyCustomDatasetItem>[] neighborPts)
{
//neighborPts = allPoints.Where(x => _metricFunc(point, x.ClusterPoint) >= epsilon).ToArray();
List <DbscanPoint <MyCustomDatasetItem> > tmp = new List <DbscanPoint <MyCustomDatasetItem> >();
foreach (var item in point.ClusterPoint.Columns)
{
if (item.Value >= epsilon)
{
tmp.Add(allPoints[contains_MDS(allPoints, item.Key)]);
}
}
neighborPts = tmp.ToArray();
}
/// <summary>
///
/// </summary>
/// <param name="allPoints">Dataset</param>
/// <param name="point">point to be in a cluster</param>
/// <param name="neighborPts">other points in same region with point parameter</param>
/// <param name="clusterId">given clusterId</param>
/// <param name="epsilon">Desired region ball range</param>
/// <param name="minPts">Minimum number of points to be in a region</param>
private void ExpandCluster(DbscanPoint <MyCustomDatasetItem>[] allPoints, DbscanPoint <MyCustomDatasetItem> point, DbscanPoint <MyCustomDatasetItem>[] neighborPts, int clusterId, double epsilon, int minPts)
{
point.ClusterId = clusterId;
for (int i = 0; i < neighborPts.Length; i++)
{
DbscanPoint <MyCustomDatasetItem> pn = neighborPts[i];
if (!pn.IsVisited)
{
pn.IsVisited = true;
DbscanPoint <MyCustomDatasetItem>[] neighborPts2 = null;
MN_RegionQuery(allPoints, pn, epsilon, out neighborPts2);
if (neighborPts2.Length >= minPts)
{
neighborPts = neighborPts.Union(neighborPts2).ToArray();
}
}
if (pn.ClusterId == (int)ClusterIds.Unclassified)
{
pn.ClusterId = clusterId;
}
}
}