胡忠超

摘? 要：傳統(tǒng)模糊C均值（Fuzzy C-Means， FCM）聚類(lèi)算法沒(méi)有充分考慮圖像中空間信息的約束作用，導(dǎo)致其對(duì)噪聲比較敏感。為了改善上述問(wèn)題，基于隱馬爾科夫（Hidden Markov Random Field， HMRF）模型的模糊聚類(lèi)算法被提出，該文將以遙感圖像為研究數(shù)據(jù)，對(duì)該算法的分類(lèi)結(jié)果進(jìn)行定性評(píng)價(jià)。

關(guān)鍵詞：圖像分類(lèi)? FCM? HMRF

中圖分類(lèi)號(hào)：TP7 ? ?文獻(xiàn)標(biāo)識(shí)碼：A 文章編號(hào)：1672-3791（2019）07（b）-0033-02

遙感圖像分類(lèi)作為遙感圖像處理中的重要環(huán)節(jié)之一被廣泛研究。傳統(tǒng)的圖像分類(lèi)方法，不能有效描述像素類(lèi)屬的不確定性。為此，將模糊集理論引入到聚類(lèi)算法中，但傳統(tǒng)的FCM算法對(duì)噪聲極其敏感。隱馬爾科夫隨機(jī)場(chǎng)理論（Hidden Markov Random Field， HMRF）[1]結(jié)合圖像構(gòu)建HMRF模型，由于HMRF缺乏明確的物理模型，為此，借助于Potts模型對(duì)其加以刻畫(huà)。此外，由于歐氏距離對(duì)噪聲極其敏感，假設(shè)圖像中像素服從高斯分布，并以其負(fù)對(duì)數(shù)函數(shù)定義像素與聚類(lèi)間非相似性測(cè)度。并且引入基于先驗(yàn)概率的K-L規(guī)則化項(xiàng)代替模糊因子控制算法模糊程度，賦予目標(biāo)函數(shù)以明確的物理意義。

1? 算法描述

假設(shè)圖像描述為z={zi，i=1，…，n}，其中zi為像素光譜測(cè)度值，取值范圍為0～255，i為像素索引，n為圖像中總像素?cái)?shù)。假設(shè)圖像中具有c個(gè)類(lèi)別，則隸屬度矩陣可表示為U=[uij]n×c，其中，j為聚類(lèi)索引。

基于HMRF模型模糊聚類(lèi)算法的目標(biāo)函數(shù)為：

JHRMF? ? ? ? ? ? ? ? ? ? ? ? ? ? ? （1）

其中，λ為模糊項(xiàng)系數(shù)，控制算法模糊程度，λ越大算法越模糊;dij為像素與聚類(lèi)非相似性測(cè)度。采用高斯分布的負(fù)對(duì)數(shù)定義：

其中，p（zi|μj，Σj）為多值高斯分布概率密度函數(shù)，w為數(shù)據(jù)維度，對(duì)于彩色遙感圖像一般取值為3，μ為均值，為協(xié)方差矩陣。

式（1）中的πij為先驗(yàn)概率，基于HMRF理論采用Potts模型定義。

其中，η為鄰域作用強(qiáng)度，值越大鄰域作用越強(qiáng);L={Li=l，i=1，…，n}為定義在圖像域中的標(biāo)號(hào)場(chǎng)，Li為第i個(gè)像素的標(biāo)號(hào)，l為其具體實(shí)現(xiàn)，且l∈{1，…，c}。Ni為像素i的8鄰域系統(tǒng)，i'為鄰域像素索引，且i′∈Ni。并通過(guò)求導(dǎo)得到均值、協(xié)方差及隸屬度矩陣的迭代公式。最后，通過(guò)反模糊化隸屬度矩陣得到分類(lèi)結(jié)果。

2? 實(shí)驗(yàn)結(jié)果與討論

選取像素為128×128大小的遙感，如圖1（a）所示，類(lèi)數(shù)為2。圖1（b）為基于HMRF模型的模糊聚類(lèi)算法分類(lèi)結(jié)果，其中λ=0.5，η=0.25。由圖1可知，引入HMRF模型的模糊聚類(lèi)方法能夠有效克服噪聲對(duì)分類(lèi)結(jié)果的影響，實(shí)驗(yàn)結(jié)果有效證明了該方法的有效性及魯棒性。

3? 結(jié)語(yǔ)

通過(guò)對(duì)遙感圖像進(jìn)行分類(lèi)實(shí)驗(yàn)，結(jié)果表明，基于HMRF的模糊聚類(lèi)方法具有較好的抗噪性能，并且能夠得到理想的分類(lèi)結(jié)果。

參考文獻(xiàn)

[1] Chatzis S P， Varvarigou T A. A fuzzy clustering approach toward hidden Markov random field models for enhanced spatially constrained image segmentation[J].IEEE Transactions on Fuzzy Systems，2008，16（5）：1351-1361.