| 5 | 1/1 | 返回列表 |
| 查看: 676 | 回復(fù): 3 | |||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | |||
jgj030106銀蟲(chóng) (正式寫(xiě)手)
|
[求助]
matlab中實(shí)現(xiàn)c-means聚合分類算法 已有2人參與
|
||
|
如何在matlab中實(shí)現(xiàn)c-means聚合分類算法,最好給個(gè)源代碼,必有重謝 發(fā)自小木蟲(chóng)Android客戶端 |
程序 |

新蟲(chóng) (初入文壇)
|
%%%%%%%%%%%%%%% FCM算法分割圖像 %%%%%%%%%%%%%% %matlab function clusterResult = FCM(imagePath, C, V, M, iter, epsm) % 模糊C均值(FCM)聚類算法分割圖像 % clusterResult = FCM(imagePath, C, V, M, iter, epsm) % Example: clusterResult = FCM('E:\Image\lena.bmp') % clusterResult = FCM('E:\Image\lena.bmp',3,[0 127 255]) % Input: % imagePath: 圖像路徑 % C: 類別數(shù),缺省值為2 % V: 初始化聚類中心,缺省值為[0 255] % M: 加權(quán)指數(shù),缺省值為2 % iter: 迭代次數(shù),缺省值為100 % epsm: 迭代停止閾值,缺省值為1.0e-2 % Output: % clusterResult: 聚類中心結(jié)果 % Note: % C的值要與V的初始化聚類中心個(gè)數(shù)相同 % 設(shè)定缺省值 I= imread('1024.bmp'); if nargin < 6 epsm = 1.0e-2; end if nargin < 5 iter = 100; end if nargin < 4 M = 2; end if nargin < 3 [x,y,z]=impixel(I); V(1)=z(1); V(2)=z(2); V(3)=z(3); end if nargin < 2 C = 3; end % 讀入圖像及其信息 figure, imshow(I); title('原圖像'); [row col] = size(I); grayHist = imhist(I); % figure, imhist(I); % title('直方圖'); histProb = grayHist / (row * col); len = length(histProb); % tic % FCM迭代過(guò)程 cnt = 0; while(cnt < iter) % 計(jì)算隸屬度函數(shù)(注意要特殊考慮某個(gè)像素點(diǎn)和聚類中心一樣的情況) for i = 1 : len flag = 0; for j = 1 : C if i == V(j) U(j, i) = 1.0; if j == 1 U(j + 1 : C, i) = 0.0; elseif j == 3 U(1 : C - 1, i) = 0.0; elseif j == 2 U(1 : j - 1, i) = 0.0; U(j + 1 : C, i) = 0.0; end flag = 1; break; end end if flag == 0 u = (1.0 ./ ((i - V) .^ 2)) .^ (1.0 / (M - 1)); uSum = sum(u); U(1 : C, i) = u' / uSum; end end % 計(jì)算更新各類聚類中心 for j = 1 : C i = linspace(1, len, len); v = sum(histProb' .* i .* (U(j, .^ M));vSum = sum(histProb' .* (U(j, .^ M));if vSum == 0 clusterResult(j) = 0; else clusterResult(j) = v / vSum; end end % 計(jì)算誤差并判斷算法迭代是否停止 diff = sum((clusterResult - V) .^ 2); if diff <= epsm break; else V = clusterResult; end cnt = cnt + 1; end % toc % 分割圖像 for i = 1 : row for j = 1 : col temp = (double(I(i, j)) - clusterResult) .^ 2; [fmin pos] = min(temp); I(i, j) = pos * 255 / C; end end figure, imshow(uint8(I)); title('分類后的圖像'); imshow(uint8(I),'Colormap',hsv(255)); %disp('迭代次數(shù):iterTimes = '); disp(cnt); % end of FCM.m |
銀蟲(chóng) (小有名氣)
| CodeForge的這個(gè)代碼或許可以解答你的問(wèn)題:http://www.codeforge.cn/article/289338 |
新蟲(chóng) (初入文壇)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿哈爾濱工業(yè)大學(xué)材料與化工方向336分 +9 | 辰沐5211314 2026-03-26 | 9/450 |
|
|---|---|---|---|---|
|
[考研] 調(diào)劑310 +7 | 溫柔的晚安 2026-03-25 | 8/400 |
|
|
[考研] 322求調(diào)劑 +7 | 宋明欣 2026-03-27 | 7/350 |
|
|
[考研] 299求調(diào)劑 +7 | 嗯嗯嗯嗯2 2026-03-27 | 7/350 |
|
|
[考研] 【求調(diào)劑】085601材料工程專碩 | 總分272 | +6 | 腳滑的守法公民 2026-03-27 | 6/300 |
|
|
[考研] 0856,材料與化工321分求調(diào)劑 +12 | 大饞小子 2026-03-27 | 13/650 |
|
|
[考研] 一志愿211院校 344分 東北農(nóng)業(yè)大學(xué)生物學(xué)學(xué)碩,求調(diào)劑 +5 | 丶風(fēng)雪夜歸人丶 2026-03-26 | 8/400 |
|
|
[考研] 0703化學(xué)338求調(diào)劑! +6 | Zuhui0306 2026-03-26 | 7/350 |
|
|
[考研] 求調(diào)劑 +6 | 林之夕 2026-03-24 | 6/300 |
|
|
[考研] 081200-11408-276學(xué)碩求調(diào)劑 +3 | 崔wj 2026-03-26 | 3/150 |
|
|
[考研] 340求調(diào)劑 +3 | Amber00 2026-03-26 | 3/150 |
|
|
[考研]
|
WWW西西弗斯 2026-03-24 | 8/400 |
|
|
[考研] 一志愿 南京郵電大學(xué) 288分 材料考研 求調(diào)劑 +3 | jl0720 2026-03-26 | 3/150 |
|
|
[考研] 材料專碩 335 分求調(diào)劑 +4 | 拒絕冷暴力 2026-03-25 | 4/200 |
|
|
[考研] 考研一志愿蘇州大學(xué)初始315(英一)求調(diào)劑 +3 | sbdksD 2026-03-24 | 4/200 |
|
|
[考研] 340求調(diào)劑 +5 | 話梅糖111 2026-03-24 | 5/250 |
|
|
[考研] 生物學(xué)學(xué)碩求調(diào)劑 +7 | 小羊睡著了? 2026-03-23 | 10/500 |
|
|
[考研] 284求調(diào)劑 +3 | yanzhixue111 2026-03-23 | 6/300 |
|
|
[考研] 269求調(diào)劑 +4 | 我想讀研11 2026-03-23 | 4/200 |
|
|
[考研] 336化工調(diào)劑 +4 | 王大坦1 2026-03-23 | 5/250 |
|