DBSCAN密度聚類

2022-06-11 17:18:14 字數 3715 閱讀 7671

**

dbscan是一種基於密度的聚類演算法,這類密度聚類演算法一般假定類別可以通過樣本分佈的緊密程度決定。同一類別的樣本,他們之間的緊密相連的,也就是說,在該類別任意樣本周圍不遠處一定有同類別的樣本存在。

通過將緊密相連的樣本劃為一類,這樣就得到了乙個聚類類別。通過將所有各組緊密相連的樣本劃為各個不同的類別,則我們就得到了最終的所有聚類類別結果。

在上一節我們定性描述了密度聚類的基本思想,本節我們就看看dbscan是如何描述密度聚類的。dbscan是基於一組鄰域來描述樣本集的緊密程度的,引數(ϵϵ, minpts)用來描述鄰域的樣本分佈緊密程度。其中,ϵϵ描述了某一樣本的鄰域距離閾值,minpts描述了某一樣本的距離為ϵϵ的鄰域中樣本個數的閾值。

假設我的樣本集是d=(x1,x2,...,xm)(x1,x2,...,xm),則dbscan具體的密度描述定義如下:

1) ϵϵ-鄰域:對於xj∈dxj∈d,其ϵϵ-鄰域包含樣本集d中與xjxj的距離不大於ϵϵ的子樣本集,即nϵ(xj)=nϵ(xj)=, 這個子樣本集的個數記為|nϵ(xj)||nϵ(xj)| 

2) 核心物件:對於任一樣本xj∈dxj∈d,如果其ϵϵ-鄰域對應的nϵ(xj)nϵ(xj)至少包含minpts個樣本,即如果|nϵ(xj)|≥minpts|nϵ(xj)|≥minpts,則xjxj是核心物件。 

3)密度直達:如果xixi位於xjxj的ϵϵ-鄰域中,且xjxj是核心物件,則稱xixi由xjxj密度直達。注意反之不一定成立,即此時不能說xjxj由xixi密度直達, 除非且xixi也是核心物件。

4)密度可達:對於xixi和xjxj,如果存在樣本樣本序列p1,p2,...,ptp1,p2,...,pt,滿足p1=xi,pt=xjp1=xi,pt=xj, 且pt+1pt+1由ptpt密度直達,則稱xjxj由xixi密度可達。也就是說,密度可達滿足傳遞性。此時序列中的傳遞樣本p1,p2,...,pt−1p1,p2,...,pt−1均為核心物件,因為只有核心物件才能使其他樣本密度直達。注意密度可達也不滿足對稱性,這個可以由密度直達的不對稱性得出。

5)密度相連:對於xixi和xjxj,如果存在核心物件樣本xkxk,使xixi和xjxj均由xkxk密度可達,則稱xixi和xjxj密度相連。注意密度相連關係是滿足對稱性的。

從下圖可以很容易看出理解上述定義,圖中minpts=5,紅色的點都是核心物件,因為其ϵϵ-鄰域至少有5個樣本。黑色的樣本是非核心物件。所有核心物件密度直達的樣本在以紅色核心物件為中心的超球體內,如果不在超球體內,則不能密度直達。圖中用綠色箭頭連起來的核心物件組成了密度可達的樣本序列。在這些密度可達的樣本序列的ϵϵ-鄰域內所有的樣本相互都是密度相連的。

有了上述定義,dbscan的聚類定義就簡單了。

dbscan的聚類定義很簡單:由密度可達關係匯出的最大密度相連的樣本集合,即為我們最終聚類的乙個類別,或者說乙個簇。

這個dbscan的簇裡面可以有乙個或者多個核心物件。如果只有乙個核心物件,則簇裡其他的非核心物件樣本都在這個核心物件的ϵϵ-鄰域裡;如果有多個核心物件,則簇裡的任意乙個核心物件的ϵϵ-鄰域中一定有乙個其他的核心物件,否則這兩個核心物件無法密度可達。這些核心物件的ϵϵ-鄰域裡所有的樣本的集合組成的乙個dbscan聚類簇。

那麼怎麼才能找到這樣的簇樣本集合呢?dbscan使用的方法很簡單,它任意選擇乙個沒有類別的核心物件作為種子,然後找到所有這個核心物件能夠密度可達的樣本集合,即為乙個聚類簇。接著繼續選擇另乙個沒有類別的核心物件去尋找密度可達的樣本集合,這樣就得到另乙個聚類簇。一直執行到所有核心物件都有類別為止。

基本上這就是dbscan演算法的主要內容了,是不是很簡單?但是我們還是有三個問題沒有考慮。

第乙個是一些異常樣本點或者說少量游離於簇外的樣本點,這些點不在任何乙個核心物件在周圍,在dbscan中,我們一般將這些樣本點標記為噪音點。

第二個是距離的度量問題,即如何計算某樣本和核心物件樣本的距離。在dbscan中,一般採用最近鄰思想,採用某一種距離度量來衡量樣本距離,比如歐式距離。這和knn分類演算法的最近鄰思想完全相同。對應少量的樣本,尋找最近鄰可以直接去計算所有樣本的距離,如果樣本量較大,則一般採用kd樹或者球樹來快速的搜尋最近鄰。如果大家對於最近鄰的思想,距離度量,kd樹和球樹不熟悉,建議參考之前寫的另一篇文章k近鄰法(knn)原理小結。

第三種問題比較特殊,某些樣本可能到兩個核心物件的距離都小於ϵϵ,但是這兩個核心物件由於不是密度直達,又不屬於同乙個聚類簇,那麼如果界定這個樣本的類別呢?一般來說,此時dbscan採用先來後到,先進行聚類的類別簇會標記這個樣本為它的類別。也就是說dbscan的演算法不是完全穩定的演算法。

下面我們對dbscan聚類演算法的流程做乙個總結。

輸入:樣本集d=(x1,x2,...,xm)(x1,x2,...,xm),鄰域引數(ϵ,minpts)(ϵ,minpts), 樣本距離度量方式

輸出: 簇劃分c. 

1)初始化核心物件集合ω=∅ω=∅, 初始化聚類簇數k=0,初始化未訪問樣本集合γγ = d,  簇劃分c = ∅∅

2) 對於j=1,2,...m, 按下面的步驟找出所有的核心物件:

a) 通過距離度量方式,找到樣本xjxj的ϵϵ-鄰域子樣本集nϵ(xj)nϵ(xj)

b) 如果子樣本集樣本個數滿足|nϵ(xj)|≥minpts|nϵ(xj)|≥minpts, 將樣本xjxj加入核心物件樣本集合:ω=ω∪ω=ω∪

3)如果核心物件集合ω=∅ω=∅,則演算法結束,否則轉入步驟4.

4)在核心物件集合ωω中,隨機選擇乙個核心物件oo,初始化當前簇核心物件佇列ωcur=ωcur=, 初始化類別序號k=k+1,初始化當前簇樣本集合ck=ck=, 更新未訪問樣本集合γ=γ−γ=γ−

5)如果當前簇核心物件佇列ωcur=∅ωcur=∅,則當前聚類簇ckck生成完畢, 更新簇劃分c=, 更新核心物件集合ω=ω−ckω=ω−ck, 轉入步驟3。

6)在當前簇核心物件佇列ωcurωcur中取出乙個核心物件o′o′,通過鄰域距離閾值ϵϵ找出所有的ϵϵ-鄰域子樣本集nϵ(o′)nϵ(o′),令δ=nϵ(o′)∩γδ=nϵ(o′)∩γ, 更新當前簇樣本集合ck=ck∪δck=ck∪δ, 更新未訪問樣本集合γ=γ−δγ=γ−δ,  更新ωcur=ωcur∪(δ∩ω)−o′ωcur=ωcur∪(δ∩ω)−o′,轉入步驟5.

輸出結果為: 簇劃分c=

和傳統的k-means演算法相比,dbscan最大的不同就是不需要輸入類別數k,當然它最大的優勢是可以發現任意形狀的聚類簇,而不是像k-means,一般僅僅使用於凸的樣本集聚類。同時它在聚類的同時還可以找出異常點,這點和birch演算法類似。

那麼我們什麼時候需要用dbscan來聚類呢?一般來說,如果資料集是稠密的,並且資料集不是凸的,那麼用dbscan會比k-means聚類效果好很多。如果資料集不是稠密的,則不推薦用dbscan來聚類。

下面對dbscan演算法的優缺點做乙個總結。

dbscan的主要優點有:

1) 可以對任意形狀的稠密資料集進行聚類,相對的,k-means之類的聚類演算法一般只適用於凸資料集。

2) 可以在聚類的同時發現異常點,對資料集中的異常點不敏感。

3) 聚類結果沒有偏倚,相對的,k-means之類的聚類演算法初始值對聚類結果有很大影響。

dbscan的主要缺點有:

1)如果樣本集的密度不均勻、聚類間距差相差很大時,聚類質量較差,這時用dbscan聚類一般不適合。

2) 如果樣本集較大時,聚類收斂時間較長,此時可以對搜尋最近鄰時建立的kd樹或者球樹進行規模限制來改進。

3) 調參相對於傳統的k-means之類的聚類演算法稍複雜,主要需要對距離閾值ϵϵ,鄰域樣本數閾值minpts聯合調參,不同的引數組合對最後的聚類效果有較大影響。

DBSCAN 密度聚類

dbscan 密度聚類 dbscan演算法是一種基於密度的聚類演算法 聚類的時候不需要預先指定簇的個數 最終的簇的個數不定 dbscan演算法將資料點分為三類 核心點 在半徑eps內含有超過minpts數目的點 邊界點 在半徑eps內點的數量小於minpts,但是落在核心點的鄰域內 噪音點 既不是核...

密度聚類演算法DBScan

利用密度聚類dbscan對樣本進行分類。基本思路 1 讀取兩組不同資料訓練 2 分別對兩組資料進行聚類 3 將聚類結果視作乙個多維空間的點,計算其到原點的歐氏距離 4 根據ans1和ans2的歐氏距離找到合適的閾值 幾個必要概念 鄰域 對於樣本集中的xj,它的 鄰域為樣本集中與它距離小於 的樣本所構...

(十八)密度聚類DBSCAN

dbscan是一種基於密度的聚類演算法,這類密度聚類演算法一般假定類別可以通過樣本分佈的緊密程度決定。同一類別的樣本,他們之間的緊密相連的,也就是說,在該類別任意樣本周圍不遠處一定有同類別的樣本存在。通過將緊密相連的樣本劃為一類,這樣就得到了乙個聚類類別。通過將所有各組緊密相連的樣本劃為各個不同的類...