基于改進(jìn)譜聚類算法在圖像分割中的應(yīng)用

王貝貝 楊明 燕慧超

摘要：為解決譜聚類算法應(yīng)用于圖像分割時(shí)，相似矩陣內(nèi)存占用較大甚至滿溢以及后續(xù)計(jì)算量大的問題，利用Nystrom方法隨機(jī)獲取一部分樣本點(diǎn)，根據(jù)樣本點(diǎn)和樣本點(diǎn)、樣本點(diǎn)和非樣本點(diǎn)2種相似關(guān)系近似表征所有像素點(diǎn)的相似性，得到原圖像的近似相似矩陣。在構(gòu)建上述所需2種相似關(guān)系的相似矩陣時(shí)，距離度量采用余弦函數(shù)。結(jié)果表明，采用近鄰傳播聚類算法代替kmeans算法對得到的低維向量子空間聚類，克服了聚類過程對初始值的敏感性，得到的分割結(jié)果較穩(wěn)定，4幅真實(shí)圖片也驗(yàn)證了研究算法的優(yōu)越性。改進(jìn)的譜聚類算法為圖像分割的穩(wěn)定性研究提供了依據(jù)。

關(guān)鍵詞：圖像處理；圖像分割；譜聚類；Nystrom方法；余弦距離；AP算法

中圖分類號：TP391文獻(xiàn)標(biāo)志碼：Adoi： 10.7535/hbgykj.2018yx01010

Application of improved spectral clustering algorithm in image segmentation

WANG Beibei， YANG Ming， YAN Huichao

（School of Science， North University of China， Taiyuan， Shanxi 030051， China）

Abstract：To solve the problem of large usage of the similar matrix memory， even overflowing，and the large amount of calculation when spectral clustering algorithm is applied in image segmentation， a part of sample points are obtained by using the Nystrom method， according to the two similarity relations between sample points and sample points， sample points and nonsample points to get the similarity relation of all pixel points， then the approximate similarity matrix of the original image is got. In building the two block similar matrixs required in Nystrom， cosine function is used to overcome the absoluteness of Euclidean distance in the distance measurement . In this paper， the nearest neighbor propagation clustering algorithm is used to replace the kmeans algorithm to overcome the sensitivity to the initial value in the clustering process， thus obtaining a more stable effect. The experimental results of the 4 real images show that the method is far superior， and the improved spectral clustering algorithm provides reference for the stable study of image segmentation.

Keywords：image processing； image segmentation； spectral clustering； Nystrom method； cosine distance； AP algorithm

圖像分割是模式識別中的一個(gè)重要研究領(lǐng)域，它的目標(biāo)是根據(jù)對圖像的像素屬性（如灰度、顏色和紋理等）進(jìn)行劃分從而實(shí)現(xiàn)對圖像的分割。在實(shí)際應(yīng)用中，對同一幅圖像，不同的研究人員可能需要提取不同的信息，人們需要也即感興趣的區(qū)域稱為前景，其余區(qū)域稱為背景。圖像分割就是把圖像分成各具特性的區(qū)域并提取出感興趣目標(biāo)的技術(shù)和過程[12]。

聚類算法在圖像分割中應(yīng)用得較為廣泛，其中譜聚類算法是依托譜圖劃分理論演化而來的，優(yōu)于kmeans算法[3]、FCM算法[4]等傳統(tǒng)算法，這是因?yàn)槠洳皇軜颖究臻g形狀的限制，能收斂于全局最優(yōu)解的特性。在進(jìn)行圖像分割時(shí)[5]，譜聚類算法將圖像的每個(gè)像素點(diǎn)視為對應(yīng)無向圖的一個(gè)節(jié)點(diǎn)，用像素矩陣來表達(dá)圖像，這樣就把圖像的分割問題轉(zhuǎn)換為對數(shù)據(jù)的聚類問題，此時(shí)，再根據(jù)某種距離準(zhǔn)則構(gòu)造新數(shù)據(jù)集的相似矩陣，新矩陣中的元素表示任意兩像素點(diǎn)間的相似度，若需要分割的圖像像素點(diǎn)個(gè)數(shù)為m×n，則它對應(yīng)的相似矩陣的維數(shù)為（m×n） ×（m×n），計(jì)算量大，占用內(nèi)存也較大。為解決該問題，文獻(xiàn)＼[6＼]基于Nystrom方法提出了一種近似方法估計(jì)相似矩陣，首先解決一個(gè)小的隨機(jī)像素子集的分組問題，然后將該解外推到圖像或圖像序列中的全部像素組中，這樣能簡單、高效地估計(jì)出相似矩陣的主特征向量，分割的結(jié)果也較穩(wěn)定。JIAO等[7]將譜聚類算法與網(wǎng)格結(jié)合，通過將原始網(wǎng)格模型嵌入到光譜空間中，來解決分割問題。首先根據(jù)最小規(guī)則確定網(wǎng)格凹區(qū)，然后通過考慮網(wǎng)格顯著性和曲率信息定義拉普拉斯矩陣，最后用特征分解法計(jì)算拉普拉斯矩陣的前k個(gè)特征向量，將原始網(wǎng)格嵌入到一個(gè)k維的光譜空間中，利用混合高斯方法實(shí)現(xiàn)最終的分割。李楊等[8]研究了在不同的顏色空間中，不同的距離測度構(gòu)造的相似矩陣對分割結(jié)果的影響，通過分析比較6種不同相似度構(gòu)造方法的分割效果，得出無論在RGB還是HSV顏色空間下，采用余弦距離公式得到的分割效果都是最優(yōu)的。李小斌等[9]設(shè)計(jì)了一個(gè)遞歸的多尺度分層分割算法，解決了譜聚類應(yīng)用于圖像分割時(shí)權(quán)值矩陣的特征值和特征向量難以計(jì)算的問題，首先定義像素點(diǎn)與類之間的距離，在研究對圖像的像素進(jìn)行抽樣的基礎(chǔ)上得到一個(gè)抽樣數(shù)定理，對少量的采樣點(diǎn)采用SCAWM聚類，然后再按距離最近原則對大量的非采樣點(diǎn)聚類，引入尺度因子使圖像從“粗”到“細(xì)”完成分割。

1經(jīng)典譜聚類算法研究

在數(shù)據(jù)聚類和圖像分割中，結(jié)合譜圖理論的譜聚類算法一直表現(xiàn)很好，并已成功用于解決數(shù)據(jù)聚類和圖像分割問題，而對圖像的分割實(shí)質(zhì)也是對數(shù)據(jù)的聚類，因?yàn)閷D像的分割是對代表圖像信息的像素點(diǎn)集X=（x1，x2，…，xn）的聚類。譜聚類算法在數(shù)據(jù)聚類中的高性能和較易實(shí)現(xiàn)，近年來吸引了越來越多人的關(guān)注。其本質(zhì)是將聚類問題轉(zhuǎn)化為圖的最優(yōu)劃分問題，將圖像映射為帶權(quán)無向圖，圖的每個(gè)節(jié)點(diǎn)對應(yīng)一個(gè)像素，每對像素之間的相似性作為無向圖的邊，通過邊上的權(quán)值構(gòu)造相似性矩陣。大致過程可歸納為3步。

1）根據(jù)某種相似性度量構(gòu)建待聚類數(shù)據(jù)集的相似矩陣W=[wij]；

2）計(jì)算規(guī)范化拉普拉斯矩陣L=D-12WD-12的前k個(gè)特征值及其對應(yīng)的特征向量，其中D=diag（d1，d2，…，di，…，dn），di=∑nj=1wij，構(gòu)建特征空間的向量矩陣V=[v1，v2，…，vk]∈Rn×k；

3）對上述得到的向量矩陣用經(jīng)典聚類算法進(jìn)行聚類。

然而，在傳統(tǒng)的譜聚類算法中依然存在如下問題：

1）常用的基于相似性度量的高斯核函數(shù)：

W（xi，xj）=exp（-‖xi-xj‖2/2σ2）（1）

不能充分反映數(shù)據(jù)集的復(fù)雜空間分布，當(dāng)集群是復(fù)雜的流形結(jié)構(gòu)時(shí)不適用[10]；

2）當(dāng)數(shù)據(jù)集的規(guī)模N比較大時(shí)，總的時(shí)間復(fù)雜度和空間復(fù)雜度分別達(dá)到了O（N3）和O（N2），儲存和分解一個(gè)大的相似矩陣就比較困難，尤其是對一個(gè)圖像。

2改進(jìn)譜聚類算法

2.1基于Nystrom方法的譜聚類

譜聚類算法的瓶頸之一是相似矩陣的構(gòu)建及其計(jì)算，首先解決計(jì)算量大的問題。

因?yàn)橄嗨凭仃囍性氐膫€(gè)數(shù)是樣本數(shù)的平方，例如：一幅由512×512個(gè)像素組成的圖像，需要計(jì)算262 144×262 144維的相似矩陣，計(jì)算量較大，計(jì)算時(shí)間長，很可能會引起內(nèi)存滿溢的情況。有研究表明，基于Nystrom方法的譜聚類能在很大程度上提升譜聚類算法的性能。

Nystrom方法是一種基于積分方程，計(jì)算矩陣特征值和特征向量的近似方法，適用于對稱正定矩陣。其基本思想是采用隨機(jī)采樣的方式獲得一定的樣本點(diǎn)，用樣本點(diǎn)與非樣本點(diǎn)之間的相似性來近似估計(jì)非樣本點(diǎn)之間的相似性[11]。若原圖像的像素總數(shù)為N，從中隨機(jī)采集m個(gè)樣本像素點(diǎn)（mN），則非樣本點(diǎn)的個(gè)數(shù)為n=N-m，相似性矩陣W可以表示為

式（2）中A∈Rm×m，B∈Rn×m，C∈Rn×n。矩陣A表示隨機(jī)采集的m個(gè)樣本像素點(diǎn)之間的相似矩陣，矩陣B表示隨機(jī)采集的m個(gè)樣本像素點(diǎn)與n個(gè)非樣本像素點(diǎn)之間的相似矩陣，矩陣C表示n個(gè)非樣本像素點(diǎn)之間的相似矩陣。因?yàn)閙N，所以矩陣C計(jì)算量如同未采樣前的矩陣的計(jì)算量一樣巨大，故通過Nystrom采樣方法用A和B估計(jì)C[12]。

先將矩陣A對角化為A=UΛUT，其中Λ是一個(gè)對角矩陣，非零元素對應(yīng)矩陣A的特征值，再將特征向量U擴(kuò)展為相似矩陣W的近似特征向量：

由式（3）可知Nystrom采樣方法是用BTA-1B估計(jì)C，即C=BTA-1B。由此，便可以用較少的樣本像素點(diǎn)逼近全部點(diǎn)的相似矩陣W。

筆者得到了相似矩陣和其主特征向量的近似估計(jì)和。但是應(yīng)用到譜聚類中，需要的是特征系統(tǒng)D-1/2WD-1/2的主特征向量。所以接下來，計(jì)算的行和，將其對應(yīng)為對角矩陣D中對角線上的元素，估計(jì)矩陣的行和為

若矩陣A正定，接下來便解決估計(jì)特征向量的正交化問題。令A(yù)1/2表示A的對稱正定平方根，并定義Q=A+A-1/2BBTA-1/2，將Q對角化為Q=UQΛQUTQ，則W可正交化為=VΛVT，其中：

經(jīng)過上述變換，則可取矩陣V的前k個(gè)向量構(gòu)成矩陣E，再對矩陣E的行向量聚類，由此便可減少在譜聚類算法中的計(jì)算量。

2.2基于余弦函數(shù)的相似矩陣

在上述過程中，筆者已經(jīng)用Nystrom方法解決了相似矩陣計(jì)算量大的問題，該節(jié)解決的則是相似矩陣的構(gòu)建問題，也即采取合適的高斯度量來表示相似矩陣。圖像的相似度矩陣由該圖像中各像素之間的相似度構(gòu)成，而高斯核函數(shù)中的歐氏距離度量只考慮絕對距離，在考慮相對距離時(shí)，余弦距離比其他距離表現(xiàn)出更好的實(shí)際意義。

假設(shè)待分割的圖像包含N個(gè)像素點(diǎn)，每個(gè)像素有p個(gè)特征數(shù)，則該圖像可表示為A∈RN×p，這里A=[a1，a2，…，aN]T，其中ai∈Rp，1≤i≤N，用余弦函數(shù)得到的圖的相似度矩陣可表示為W=[wij]N×N，其中：

如果對數(shù)據(jù)集A進(jìn)行歸一化，即有‖ai‖=‖aj‖=1，則像素點(diǎn)xi到像素點(diǎn)xj的余弦相似度可簡化為wij=aiaj，相應(yīng)地相似矩陣W簡化為W=[wij]N×N=AAT。由此，圖像的相似度矩陣可由圖像矩陣A和其轉(zhuǎn)置相乘得到。

2.3近鄰傳播聚類算法

在2.1節(jié)最后所得的矩陣E的每一行都對應(yīng)一個(gè)像素點(diǎn)，這節(jié)需要做的則是將新得到的矩陣相應(yīng)的數(shù)據(jù)點(diǎn)進(jìn)行聚類處理，傳統(tǒng)譜聚類算法的最后一步在采用kmeans聚類算法時(shí)，對初始值較敏感，需要人為初始化類的個(gè)數(shù)，分割結(jié)果不穩(wěn)定，即取不同的初始值得到的分割效果也不同。為解決該問題，本研究引入近鄰傳播算法（AP算法）[13]，該算法是2007年刊登在《科學(xué)》期刊上的的一篇文章提出的一種新型的聚類算法，事先不需要知道類的數(shù)目，通過對矩陣進(jìn)行迭代收斂，最終更新到穩(wěn)定的類中心，然后計(jì)算每一個(gè)數(shù)據(jù)點(diǎn)與所有類中心的相似度的大小，將其劃分到相似度最大的類中心，每個(gè)類中心及其劃分到該類中心的其他數(shù)據(jù)點(diǎn)即組成一個(gè)類別。

近鄰傳播算法需要首先構(gòu)建新的數(shù)據(jù)集中點(diǎn)與點(diǎn)之間的相似度矩陣，相似度用歐氏距離度量即可。Sij=-‖yi-yj‖2，i≠j，然后把相似度矩陣S（i，j）=-‖yi-yj‖作為輸入，交換數(shù)據(jù)點(diǎn)間的信息，交換方式主要依靠2種度量：吸引度和歸屬度。吸引度r（i，k）是數(shù)據(jù)點(diǎn)yi發(fā)送到候選類中心yk的信息，表示點(diǎn)yi對候選類中心點(diǎn)yk的支持程度，該值越大，則表明點(diǎn)yk作為點(diǎn)yi的類中心的可能性越大；歸屬度a（i，k）是候選類中心yk發(fā)送給數(shù)據(jù)點(diǎn)yi的信息，表示候選類中心yk作為點(diǎn)yi的類中心的合適程度，該值越大，則表明點(diǎn)yi選擇點(diǎn)yk作為其類中心的可能性越大。在算法迭代過程中，每個(gè)數(shù)據(jù)點(diǎn)都會進(jìn)行2種信息的交替更新，更新公式為

當(dāng)每個(gè)數(shù)據(jù)點(diǎn)都進(jìn)行了上述2次的信息更新后，算法也迭代了1次。在運(yùn)算過程中，引入阻尼系數(shù)λ∈[0，1），來避免AP算法的大幅度震蕩，迭代次數(shù)和每次迭代的結(jié)果都會隨著阻尼系數(shù)的改變而改變：λ取值較小時(shí)，迭代的次數(shù)會較少，但是運(yùn)算過程中的數(shù)值波動(dòng)性卻較大，此時(shí)得到的較大圖像的收斂結(jié)果不理想；λ取值較大時(shí)，迭代的次數(shù)雖會較多，但圖像的收斂效果較理想。本文取λ=0.9，對應(yīng)的吸引度矩陣和歸屬度矩陣更新公式如下：

其中t表示迭代次數(shù)。數(shù)據(jù)點(diǎn)yi的聚類中心即是令r（i，k）+a（i，k）取最大值的點(diǎn)yk。

2.4改進(jìn)譜聚類算法

譜聚類算法在進(jìn)行圖像分割，初始化相似矩陣時(shí)會出現(xiàn)內(nèi)存滿溢的情況，無法進(jìn)行后續(xù)工作。為了解決該問題，本研究提出一種改進(jìn)的譜聚類算法——INJW算法，對預(yù)處理的圖像，研究利用Nystrom方法隨機(jī)獲取一定的樣本點(diǎn)。先構(gòu)建樣本點(diǎn)與樣本點(diǎn)間的相似矩陣、樣本點(diǎn)與非樣本點(diǎn)間的相似矩陣，再根據(jù)式（4），用這2個(gè)分塊矩陣來近似表征所有像素點(diǎn)間的相似矩陣。在構(gòu)造上述所需的分塊矩陣時(shí)，距離度量采用余弦函數(shù)，克服歐氏距離的絕對性，由此便完整的得到了譜聚類算法第1步所需的相似矩陣；第2步所需規(guī)范化拉普拉斯矩陣的主特征向量為式（8）的前k個(gè)向量；最后對這個(gè)低維向量子空間聚類時(shí)采用近鄰傳播聚類算法代替kmeans算法，克服了聚類過程對初始值的敏感性。

1）預(yù)處理輸入圖像G，隨機(jī)抽取1%的樣本像素點(diǎn)，根據(jù)式（9）的余弦函數(shù)構(gòu)造像素點(diǎn)的分塊相似性矩陣A和B，得到待處理圖像的相似矩陣為

2）計(jì)算拉普拉斯矩陣L=D-1/2D-1/2的前k個(gè)特征向量，即根據(jù)式（8）計(jì)算出矩陣V，并取V的前k個(gè)最大特征值對應(yīng)的特征向量，得到新的特征矩陣為E。

3）將矩陣E的行向量規(guī)范化為單位矩陣Zij，

4）矩陣Zij的每一行對應(yīng)原圖像的一個(gè)像素點(diǎn)，用AP算法對其進(jìn)行聚類。

5）根據(jù)步驟4）得到的聚類結(jié)果確定所有像素點(diǎn)的類別：當(dāng)矩陣Zij的第i行對應(yīng)第j類時(shí)，像素點(diǎn)Ei對應(yīng)第j類。

3結(jié)果與分析

為了驗(yàn)證本研究算法的性能，在操作平臺為64位Windows 7系統(tǒng)、CPU為Intel（R） Core（TM）i52450M（2.50 GHz）、4 GB內(nèi)存的計(jì)算機(jī)和Matlab R2014a中進(jìn)行實(shí)驗(yàn)。實(shí)驗(yàn)圖片來自美國南加州大學(xué)的USISIPI圖像數(shù)據(jù)庫和美國加州大學(xué)伯克利分校的BSDB標(biāo)準(zhǔn)測試圖像數(shù)據(jù)庫中的4幅圖像，對比算法為k均值聚類算法（kmeans）、基于規(guī)范化拉普拉斯矩陣的譜聚類算法（NJW），實(shí)驗(yàn)對比圖見圖1—圖4。

圖1是數(shù)據(jù)庫中4幅圖的原始圖像：圖1 a）和圖1 b）分別是人物圖lena和cameraman，圖1 c）是動(dòng)物圖ducks，圖1 d）是實(shí)物圖airplane。圖2是kmeans算法的分割結(jié)果，由于該算法對初始值較敏感，且取不同的分割數(shù)時(shí)結(jié)果也不相同，故在編程時(shí)設(shè)置k分別取2，3，4，5，6，且在程序運(yùn)行20次后選擇最優(yōu)直觀分割圖，其中l(wèi)ena圖的最優(yōu)分割為k=3，cameraman圖的最優(yōu)分割為k=4，ducks圖的最優(yōu)分割為k=6，airplane圖的最優(yōu)分割為k=4。 圖3是基于規(guī)范化拉普拉斯矩陣的譜聚類算法（NJW）的分割結(jié)果圖。圖4是本文算法INJW——改進(jìn)的譜聚類算法的分割結(jié)果圖。

從總體上來看，kmeans算法分割出來的圖較模糊，前景和背景區(qū)域都含有很多噪聲點(diǎn)，NJW算法的分割結(jié)果圖較純凈，本文算法的分割效果圖最清晰，不管是人物還是動(dòng)物，前景和背景區(qū)分都較明顯，區(qū)域?qū)Ρ榷?、區(qū)域內(nèi)部均勻性都較前2種算法明顯提高。

1）人物lena圖kmeans算法的分割圖中人物的帽子、頭發(fā)和五官都已模糊化，且前景、背景混為一體；NJW算法的分割圖雖然前背景明顯，但人物的五官不是很清楚，鼻梁和臉頰并沒有分割開；本研究算法運(yùn)行出來的分割圖中帽子和頭發(fā)都更加細(xì)膩，五官清晰。

2）人物cameraman圖kmeans算法的分割圖中右側(cè)有少量噪聲點(diǎn)，前景人物和背景海水混為一體；NJW算法的分割圖則把前景人物分割出來，但是人臉五官不明顯，背景的輪廓沒有分割出來；本研究算法較好地分割開前景和背景，海上的建筑物高樓雖然并沒有完全分割出來，但總體上還是達(dá)到較理想的效果。

3）動(dòng)物ducks圖3種算法都分割出前背景，樹干上的涂白劑都和未噴涂部分分割開來，鴨子和海水也分割開。但是很明顯，kmeans算法分割出來的鴨子輪廓不明顯，且遠(yuǎn)處的3只已經(jīng)模糊化；NJW算法則把靠岸的3只鴨子和岸邊混為一體；本研究算法細(xì)節(jié)分割的更加明顯，鴨子的嘴巴也能顯露出來。

4）實(shí)物airplane圖該幅圖不僅要把前景的飛機(jī)和背景的天空、雪山分割開，還要準(zhǔn)確地分割出飛機(jī)上的圖案和字母。kmeans和NJW算法都沒能達(dá)到理想的效果，本研究算法能分割出機(jī)尾的字母和數(shù)字F16，以及飛機(jī)上的五角星標(biāo)志。

4結(jié)論

每一幅圖像都能用對應(yīng)的像素矩陣表示，所以對圖像的分割本質(zhì)上是對代表其像素點(diǎn)的數(shù)據(jù)的分割。譜聚類算法用于分割圖像時(shí)，最大的困難是相似矩陣的構(gòu)造，經(jīng)常會出現(xiàn)內(nèi)存不足甚至滿溢的情況，即使能計(jì)算出，時(shí)間成本也較高。本研究結(jié)合Nystrom采樣方法，很大程度上減少了計(jì)算的復(fù)雜度，且對相似矩陣中距離度量做了改進(jìn)，用余弦距離代替?zhèn)鹘y(tǒng)歐氏距離。最后階段，以近鄰傳播算法取代kmeans算法對k個(gè)特征主向量進(jìn)行聚類。實(shí)驗(yàn)結(jié)果驗(yàn)證了本文算法的穩(wěn)定性和有效性，分割出來的圖像質(zhì)量有了很大提升，接下來的工作重點(diǎn)將放在對前背景較復(fù)雜圖像的分割。

參考文獻(xiàn)/References：

[1]周莉莉， 姜楓. 圖像分割方法綜述研究[J]. 計(jì)算機(jī)應(yīng)用研究， 2017， 34（7）：19211928.

ZHOU Lili， JIANG Feng. Survey on image segmentation methods [J]. Application Research of Computers， 2017， 34（7）：19211928.

[2]章國紅，辛斌杰.圖像處理技術(shù)在紗線毛羽檢測方面的應(yīng)用[J].河北科技大學(xué)學(xué)報(bào)，2016，37（1）：7682.

ZHANG Guohong， XIN Binjie. Application of image processing technology in yarn hairiness detection[J]. Journal of Hebei University of Science and Technology， 2016，37（1）：7682.

[3]SHI B， PANG Z F， XU J. Image segmentation based on the hybrid total variation model and the kmeans clustering strategy[J]. Inverse Problems & Imaging， 2016， 10（3）：807828.

[4]DUBEY Y K， MUSHRIF M M. FCM clustering algorithms for segmentation of brain MR images[J]. Advances in Fuzzy Systems， 2016（3）：114.

[5]BAILLARD C， HELLIER P， BARILLOT C. Segmentation of brain 3D MR images using level sets and dense registration[J]. Medical Image Analysis， 2001， 5（3）：185194.

[6]FOWLKES C， BELONGIE S， CHUNG F， et al. Spectral grouping using the Nystrm method [J]. Pattern Analysis and Machine Intelligence IEEE Transactions， 2004， 26（2）：214225.

[7]JIAO X， WU T， QIN X. Mesh segmentation by combining mesh saliency with spectral clustering[J]. Journal of Computational & Applied Mathematics， 2018，329：134146.

[8]李揚(yáng)， 陸璐， 崔紅霞. 譜聚類圖像分割中相似度矩陣構(gòu)造研究[J]. 計(jì)算機(jī)技術(shù)與發(fā)展， 2016， 26（7）：5558.

LI Yang， LU Lu， CUI Hongxia. Research on similarity matrix structure in spectral clustering image segmentation [J]. Computer Technology and Development， 2016， 26（7）：5558.

[9]李小斌， 田錚. 基于譜聚類的圖像多尺度隨機(jī)樹分割[J]. 中國科學(xué)， 2007， 37（8）：10731085.

[10]YANG Y， WANG Y， XUE X. A novel spectral clustering method with superpixels for image segmentation[J]. Optik International Journal for Light and Electron Optics， 2016， 127（1）：161167.

[11]管濤， 李玉玲. 大規(guī)模矩陣降維的隨機(jī)逼近方法[J]. 數(shù)學(xué)的實(shí)踐與認(rèn)識， 2016，46（24）：184193.

GUAN Tao， LI Yuling. Stochastic approximation method for largescale matrix reduction[J]. Mathematics in Practice and Theory， 2016，46（24）：184193.

[12]BERGH F V D， ENGELBRECHT A P. Training product unit networks using cooperative particle swarm optimisers[J]. International Joint Conference on Neural Networks， 2001，1：126131.

[13]FREY B J， DUECK D. Clustering by passing messages between data points[J]. Science， 2007，315（5814）：972976.