这么香的DBSCAN算法你确定不用？

lsekfe

机器学习中各式各样的聚类算法层出不穷，包括层次聚类、系统聚类、K中心聚类、DBSCAN聚类等等。DBSCAN(Density-Based Spatial Clustering of Applications with Noise)是一个比较有代表性的基于密度的聚类算法，它最大的优点是能够把具有足够高密度的区域划分为簇，并在具有噪声的空间中发现任意形状的聚类，因而被广泛应用。

　　01
　　传统DBSCAN聚类算法采用穷举搜索进行数据的划分，通过在聚类过程中多次遍历待划分数据得到划分结果。这种算法称为蛮力算法，时间复杂度达到了O(n2)，数据量较少时可以在较短时间内完成数据划分。

02
　　Kd-树是K-dimension tree的缩写，是对数据点在k维空间划分的一种数据结构，主要应用于多维空间关键数据的搜索。本质上说，Kd-树就是一种平衡二叉树。所以它的增加、删除、查询时间复杂度都是O(logn)。
　　首先必须搞清楚的是，kd-树是一种空间划分树，说白了，就是把整个空间划分为特定的几个部分，然后在特定空间的部分内进行相关搜索操作。想像一个三维空间，kd-树按照一定的划分规则把这个三维空间划分了多个空间，如下图所示：

03
　　为了解决传统DBSCAN聚类算法的低效率问题，我们将Kd-树数据结构与该算法相结合，在进行聚类之前，先对待聚类数据创建Kd-树，使用空间搜索算法更高效地进行数据划分。整个数据划分流程图如下图示：

04
　　通过流程图我们可以知道，Kd-树DBSCAN聚类算法会首先对待聚类数据创建kd-树，然后通过空间搜索算法检查数据中每点的r邻域来搜索簇。如果点p的r邻域包含的点多于MinPts，则创建一个以p为核心对象的新簇。通过迭代的聚集从这些核心对象直接密度可达的对象来达到数据划分的目的。当没有新的数据点可以添加到任何簇时，算法结束。

　　05
　　Kd-树创建伪代码：

1. kdtree (list of points pointList, int depth)
   
2. {
   
3.     // Select axis based on depth so that axis cycles through all valid values
   
4.     var int axis := depth mod k;
   
5. 
   
6.     // Sort point list and choose median as pivot element
   
7.     select median by axis from pointList;
   
8. 
   
9.     // Create node and construct subtrees
   
10.     var tree_node node;
    
11.     node.location := median;
    
12.     node.leftChild := kdtree(points in pointList before median, depth+1);
    
13.     node.rightChild := kdtree(points in pointList after median, depth+1);
    
14.     return node;
    
15. }

复制代码

DBSCAN聚类算法伪代码：

1. DBSCAN(D, eps, MinPts)
   
2.     C = 0                                          //类别标示
   
3.     for each unvisited point P in dataset D        //遍历
   
4.        mark P as visited                           //已经访问
   
5.        NeighborPts = regionQuery(P, eps)           //计算这个点的邻域     
   
6.        if sizeof(NeighborPts) < MinPts             //不能作为核心点
   
7.           mark P as NOISE                          //标记为噪音数据
   
8.        else                                        //作为核心点，根据该点创建一个类别
   
9.           C = next cluster
   
10.           expandCluster(P, NeighborPts, C, eps, MinPts)    //根据该核心店扩展类别
    
11.            
    
12. expandCluster(P, NeighborPts, C, eps, MinPts)
    
13.     add P to cluster C                                     //扩展类别，核心店先加入
    
14.     for each point P' in NeighborPts                       //然后针对核心店邻域内的点，如果该点没有被访问，
    
15.        if P' is not visited
    
16.           mark P' as visited                               //进行访问
    
17.           NeighborPts' = regionQuery(P', eps)              //如果该点为核心点，则扩充该类别
    
18.           if sizeof(NeighborPts') >= MinPts
    
19.              NeighborPts = NeighborPts joined with NeighborPts'
    
20.        if P' is not yet member of any cluster              //如果邻域内点不是核心点，并且无类别，比如噪音数据，则加入此类别
    
21.           add P' to cluster C
    
22.            
    
23. regionQuery(P, eps)                                       //计算邻域
    
24.     return all points within P's eps-neighborhood

复制代码

06
　　所以，各位机器学习、深度学习大佬们，不要再拘泥于传统DBSCAN聚类算法进行图像处理、数据分析啦，对于一个高分辨率图像，数以千百万级的像素点进行数据划分时，对于O(n2)的时间复杂度来说，传统算法带来的时间开销是相当大的。为何不将数据结构Kd-树与其加以融合，在完全相同的聚类效果面前，花费时间更少的算法它不香吗？哈哈~