DBSCAN

Posted 各种控恩恩恩

tags:

篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了DBSCAN相关的知识,希望对你有一定的参考价值。

[+]
  1. 一基于密度的聚类算法的概述
  2. 二DBSCAN算法的原理
    1. 基本概念
    2. 算法流程
  3. 三实验仿真
  4. 参考文献

一、基于密度的聚类算法的概述

    最近在Science上的一篇基于密度的聚类算法《Clustering by fast search and find of density peaks》引起了大家的关注(在我的博文“论文中的机器学习算法——基于密度峰值的聚类算法”中也进行了中文的描述)。于是我就想了解下基于密度的聚类算法,熟悉下基于密度的聚类算法与基于距离的聚类算法,如K-Means算法之间的区别。

    基于密度的聚类算法主要的目标是寻找被低密度区域分离的高密度区域。与基于距离的聚类算法不同的是,基于距离的聚类算法的聚类结果是球状的簇,而基于密度的聚类算法可以发现任意形状的聚类,这对于带有噪音点的数据起着重要的作用。


二、DBSCAN算法的原理


1、基本概念

    DBSCAN(Density-Based Spatial Clustering of Application with Noise)是一种典型的基于密度的聚类算法,在DBSCAN算法中将数据点分为一下三类:

  • 核心点。在半径Eps内含有超过MinPts数目的点
  • 边界点。在半径Eps内点的数量小于MinPts,但是落在核心点的邻域内
  • 噪音点。既不是核心点也不是边界点的点

在这里有两个量,一个是半径Eps,另一个是指定的数目MinPts。

    一些其他的概念

  1. Eps邻域。简单来讲就是与点的距离小于等于Eps的所有的点的集合,可以表示为
  2. 直接密度可达。如果在核心对象的Eps邻域内,则称对象从对象出发是直接密度可达的。
  3. 密度可达。对于对象链:是从关于Eps和MinPts直接密度可达的,则对象是从对象关于Eps和MinPts密度可达的。


2、算法流程


(流程)


三、实验仿真

    在实验中使用了两个测试数据集,数据集的原始图像如下:


(数据集1)

(数据集2)

数据集1相对比较简单。显然我们可以发现数据集1共有两个类,数据集2有四个类,下面我们通过DBSCAN算法实现数据点的聚类:


MATLAB代码

主程序



[plain] view plain copy


print?

  1. %% DBSCAN  
  2. clear all;  
  3. clc;  
  4.   
  5. %% 导入数据集  
  6. % data = load('testData.txt');  
  7. data = load('testData_2.txt');  
  8.   
  9. % 定义参数Eps和MinPts  
  10. MinPts = 5;  
  11. Eps = epsilon(data, MinPts);  
  12.   
  13. [m,n] = size(data);%得到数据的大小  
  14.   
  15. x = [(1:m)' data];  
  16. [m,n] = size(x);%重新计算数据集的大小  
  17. types = zeros(1,m);%用于区分核心点1,边界点0和噪音点-1  
  18. dealed = zeros(m,1);%用于判断该点是否处理过,0表示未处理过  
  19. dis = calDistance(x(:,2:n));  
  20. number = 1;%用于标记类  
  21.   
  22. %% 对每一个点进行处理  
  23. for i = 1:m  
  24.     %找到未处理的点  
  25.     if dealed(i) == 0  
  26.         xTemp = x(i,:);  
  27.         D = dis(i,:);%取得第i个点到其他所有点的距离  
  28.         ind = find(D<=Eps);%找到半径Eps内的所有点  
  29.           
  30.         %% 区分点的类型  
  31.           
  32.         %边界点  
  33.         if length(ind) > 1 && length(ind) < MinPts+1  
  34.             types(i) = 0;  
  35.             class(i) = 0;  
  36.         end  
  37.         %噪音点  
  38.         if length(ind) == 1  
  39.             types(i) = -1;  
  40.             class(i) = -1;  
  41.             dealed(i) = 1;  
  42.         end  
  43.         %核心点(此处是关键步骤)  
  44.         if length(ind) >= MinPts+1  
  45.             types(xTemp(1,1)) = 1;  
  46.             class(ind) = number;  
  47.               
  48.             % 判断核心点是否密度可达  
  49.             while ~isempty(ind)  
  50.                 yTemp = x(ind(1),:);  
  51.                 dealed(ind(1)) = 1;  
  52.                 ind(1) = [];  
  53.                 D = dis(yTemp(1,1),:);%找到与ind(1)之间的距离  
  54.                 ind_1 = find(D<=Eps);  
  55.                   
  56.                 if length(ind_1)>1%处理非噪音点  
  57.                     class(ind_1) = number;  
  58.                     if length(ind_1) >= MinPts+1  
  59.                         types(yTemp(1,1)) = 1;  
  60.                     else  
  61.                         types(yTemp(1,1)) = 0;  
  62.                     end  
  63.                       
  64.                     for j=1:length(ind_1)  
  65.                        if dealed(ind_1(j)) == 0  
  66.                           dealed(ind_1(j)) = 1;  
  67.                           ind=[ind ind_1(j)];     
  68.                           class(ind_1(j))=number;  
  69.                        end                      
  70.                    end  
  71.                 end  
  72.             end  
  73.             number = number + 1;  
  74.         end  
  75.     end  
  76. end  
  77.   
  78. % 最后处理所有未分类的点为噪音点  
  79. ind_2 = find(class==0);  
  80. class(ind_2) = -1;  
  81. types(ind_2) = -1;  
  82.   
  83. %% 画出最终的聚类图  
  84. hold on  
  85. for i = 1:m  
  86.     if class(i) == -1  
  87.         plot(data(i,1),data(i,2),'.r');  
  88.     elseif class(i) == 1  
  89.         if types(i) == 1  
  90.             plot(data(i,1),data(i,2),'+b');  
  91.         else  
  92.             plot(data(i,1),data(i,2),'.b');  
  93.         end  
  94.     elseif class(i) == 2  
  95.         if types(i) == 1  
  96.             plot(data(i,1),data(i,2),'+g');  
  97.         else  
  98.             plot(data(i,1),data(i,2),'.g');  
  99.         end  
  100.     elseif class(i) == 3  
  101.         if types(i) == 1  
  102.             plot(data(i,1),data(i,2),'+c');  
  103.         else  
  104.             plot(data(i,1),data(i,2),'.c');  
  105.         end  
  106.     else  
  107.         if types(i) == 1  
  108.             plot(data(i,1),data(i,2),'+k');  
  109.         else  
  110.             plot(data(i,1),data(i,2),'.k');  
  111.         end  
  112.     end  
  113. end  
  114. hold off  

%% DBSCANclear all;clc;%% 导入数据集% data = load('testData.txt');data = load('testData_2.txt');% 定义参数Eps和MinPtsMinPts = 5;Eps = epsilon(data, MinPts);[m,n] = size(data);%得到数据的大小x = [(1:m)' data];[m,n] = size(x);%重新计算数据集的大小types = zeros(1,m);%用于区分核心点1,边界点0和噪音点-1dealed = zeros(m,1);%用于判断该点是否处理过,0表示未处理过dis = calDistance(x(:,2:n));number = 1;%用于标记类%% 对每一个点进行处理for i = 1:m %找到未处理的点 if dealed(i) == 0 xTemp = x(i,:); D = dis(i,:);%取得第i个点到其他所有点的距离 ind = find(D<=Eps);%找到半径Eps内的所有点 %% 区分点的类型 %边界点 if length(ind) > 1 && length(ind) < MinPts+1 types(i) = 0; class(i) = 0; end %噪音点 if length(ind) == 1 types(i) = -1; class(i) = -1; dealed(i) = 1; end %核心点(此处是关键步骤) if length(ind) >= MinPts+1 types(xTemp(1,1)) = 1; class(ind) = number; % 判断核心点是否密度可达 while ~isempty(ind) yTemp = x(ind(1),:); dealed(ind(1)) = 1; ind(1) = []; D = dis(yTemp(1,1),:);%找到与ind(1)之间的距离 ind_1 = find(D<=Eps); if length(ind_1)>1%处理非噪音点 class(ind_1) = number; if length(ind_1) >= MinPts+1 types(yTemp(1,1)) = 1; else types(yTemp(1,1)) = 0; end for j=1:length(ind_1) if dealed(ind_1(j)) == 0 dealed(ind_1(j)) = 1; ind=[ind ind_1(j)]; class(ind_1(j))=number; end end end end number = number + 1; end endend% 最后处理所有未分类的点为噪音点ind_2 = find(class==0);class(ind_2) = -1;types(ind_2) = -1;%% 画出最终的聚类图hold onfor i = 1:m if class(i) == -1 plot(data(i,1),data(i,2),'.r'); elseif class(i) == 1 if types(i) == 1 plot(data(i,1),data(i,2),'+b'); else plot(data(i,1),data(i,2),'.b'); end elseif class(i) == 2 if types(i) == 1 plot(data(i,1),data(i,2),'+g'); else plot(data(i,1),data(i,2),'.g'); end elseif class(i) == 3 if types(i) == 1 plot(data(i,1),data(i,2),'+c'); else plot(data(i,1),data(i,2),'.c'); end else if types(i) == 1 plot(data(i,1),data(i,2),'+k'); else plot(data(i,1),data(i,2),'.k'); end endendhold off


距离计算函数

[plain] view plain copy


print?

  1. %% 计算矩阵中点与点之间的距离  
  2. function [ dis ] = calDistance( x )  
  3.     [m,n] = size(x);  
  4.     dis = zeros(m,m);  
  5.       
  6.     for i = 1:m  
  7.         for j = i:m  
  8.             %计算点i和点j之间的欧式距离  
  9.             tmp =0;  
  10.             for k = 1:n  
  11.                 tmp = tmp+(x(i,k)-x(j,k)).^2;  
  12.             end  
  13.             dis(i,j) = sqrt(tmp);  
  14.             dis(j,i) = dis(i,j);  
  15.         end  
  16.     end  
  17. end  

%% 计算矩阵中点与点之间的距离function [ dis ] = calDistance( x ) [m,n] = size(x); dis = zeros(m,m); for i = 1:m for j = i:m %计算点i和点j之间的欧式距离 tmp =0; for k = 1:n tmp = tmp+(x(i,k)-x(j,k)).^2; end dis(i,j) = sqrt(tmp); dis(j,i) = dis(i,j); end endend


epsilon函数

[plain] view plain copy


print?

  1. function [Eps]=epsilon(x,k)  
  2.   
  3. % Function: [Eps]=epsilon(x,k)  
  4. %  
  5. % Aim:   
  6. % Analytical way of estimating neighborhood radius for DBSCAN  
  7. %  
  8. % Input:   
  9. % x - data matrix (m,n); m-objects, n-variables  
  10. % k - number of objects in a neighborhood of an object  
  11. % (minimal number of objects considered as a cluster)  
  12.   
  13.   
  14.   
  15. [m,n]=size(x);  
  16.   
  17. Eps=((prod(max(x)-min(x))*k*gamma(.5*n+1))/(m*sqrt(pi.^n))).^(1/n);  

function [Eps]=epsilon(x,k)% Function: [Eps]=epsilon(x,k)%% Aim: % Analytical way of estimating neighborhood radius for DBSCAN%% Input: % x - data matrix (m,n); m-objects, n-variables% k - number of objects in a neighborhood of an object% (minimal number of objects considered as a cluster)[m,n]=size(x);Eps=((prod(max(x)-min(x))*k*gamma(.5*n+1))/(m*sqrt(pi.^n))).^(1/n);





最终的结果



(数据集1的聚类结果)


(数据集2的聚类结果)

在上面的结果中,红色的点代表的是噪音点,点代表的是边界点,十字代表的是核心点。不同的颜色代表着不同的类。



参考文献

[1] M. Ester, H. Kriegel, J. Sander, X. Xu, A density-based algorithm for discovering clusters in large spatial databases with noise,  www.dbs.informatik.uni-muenchen.de/cgi-bin/papers?query=--CO
[2] M. Daszykowski, B. Walczak, D. L. Massart, Looking for Natural Patterns in Data. Part 1: Density Based Approach

来源:http://blog.csdn.net/google19890102/article/details/37656733

以上是关于DBSCAN的主要内容,如果未能解决你的问题,请参考以下文章

初探DBSCAN聚类算法

python实现DBSCAN算法

DBSCAN密度聚类

DBSCAN算法

机器学习--聚类系列--DBSCAN算法

Python:3 维空间中的 DBSCAN