DBSCAN算法

最新推荐文章于 2023-05-03 16:41:27 发布

小魏612

最新推荐文章于 2023-05-03 16:41:27 发布

阅读量839

点赞数 1

文章标签：算法聚类机器学习

本文链接：https://blog.csdn.net/qq_62945476/article/details/124509474

版权

本文深入解析了DBSCAN（基于密度的空间聚类算法），详细介绍了其算法流程、代码实现、效果演示及参数选择。DBSCAN能够处理噪声点并发现任意形状的簇，无需预设聚类数量，但对参数选择敏感，特别是邻域半径(Eps)和邻域点个数阈值(min_sample)。此外，文章展示了在不同数据集上的应用效果，突显了其在处理非均匀密度聚类时的挑战。

摘要由CSDN通过智能技术生成

定义

DBSCAN（Density-Based Spatial Clustering of Applications with Noise，具有噪声的基于密度的聚类方法）是一种基于密度的空间聚类算法。该算法将具有足够密度的区域划分为簇，并在具有噪声的空间数据库中发现任意形状的簇，它将簇定义为密度相连的点的最大集合。

算法流程

将所有点标为未标记点，依次遍历所有点，判断该点是否被标记，如果被标记则遍历下个点，反之计算该点邻域点，判断邻域点个数，如果小于min_sample(邻域半径内最小邻域点个数)则将该点标记为噪声点，反之标记为簇Ci(i为不同簇标号)，并将该点及其邻域点存入簇Ci，从Ci中选出一个点，判断该点是否为噪声点，若是将其标号改成簇Ci标号，反之判断该点邻域点个数是否大于min_sample，如果大于说明该点可继续扩充，将该点及其邻域点加入到Ci，继续遍历Ci中点，直到Ci中所有点被遍历完，则返回开头，遍历下一个点，如果该点邻域点个数小于min_sample，说明该点是簇Ci中的边界点，不可继续扩充，则标记该点为簇Ci点，继续遍历Ci中下一个点。直到原数据所有点被遍历完，算法结束，返回数据对应的簇标号。

代码

#寻找某点的邻域元素,返回对应的索引  便于后续判断该点的标签类别（即簇类别）
def find_nerighborhood(p, data, eps):
    N = []  #存放某个点邻域元素的索引
    for i in range(len(data)):
        distance = np.sqrt(np.sum((p - data[i])**2))
        if distance < eps:
            N.append(i)
    return N


#扩充簇元素
def expand_cluster(N, data, eps, min_sample, mask, clusters_label):
    #遍历簇中每个元素，找到符合扩充条件的元素 添加到簇中
    i = 0
    while i < len(N):
        #如果该点被标为噪声点，则将其改成簇标号 且不用计算该邻域元素因为以前被标为噪声点说明现在只是边缘点
        if mask[N[i]] == -1:
            mask[N[i]] = clusters_label
        #如果该点未被标记则标为簇标号
        elif mask[N[i]] == 0:
            mask[N[i]] = clusters_label
            #计算该点邻域点
            N_new = find_nerighborhood(data[N[i]], data, eps)
            if len(N_new) > min_sample:
                #将该点邻域存入簇
                N += N_new
        i += 1
    # return N



def DBSCAN_re(data, eps, min_sample):
    #所有数据标为未标记数据
    # mask = np.zeros_like(data[:, 0], dtype=bool)   #存放标记号
    mask = np.zeros(data.shape[0])
    # clusters = []   #存放某种簇所有元素
    clusters_label = 0  #簇类别号
    for p in range(len(data)):
        #如果该点被标记则跳过该点
        if mask[p] != 0:
            continue

        #如果该点未被处理则计算该点邻域点
        N = find_nerighborhood(data[p], data, eps)
        num_ner = len(N)
        #如果邻域元素过少则标为噪声点
        if num_ner < min_sample:
            mask[p] = -1
        else:
            #将该点标为第i类 i = 1, 2, ...
            clusters_label += 1
            mask[p] = clusters_label
            #扩充簇元素
            expand_cluster(N, data, eps, min_sample, mask, clusters_label)
    return mask  #返回每个点对应的类别标号