av在线精品一区二区,中文字幕人妻av,91精品国产综合婷婷香蕉

當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖

xiyue6911

新蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 421.2
帖子: 79
在線: 35小時
蟲號: 3927522

已領(lǐng)完

基于Kmeans的圖像聚類算法: 領(lǐng)取紅包 (小木蟲手機(jī)app專屬紅包)

掃一掃，下載小木蟲客戶端

使用方法：直接在matlab中運(yùn)行，但注意更改圖片路徑名；
主要分三步：
1，初始化聚類中心，我使用K等分點(diǎn)做的；
2，隨機(jī)性采樣，出于效率的考慮；
3，剩余點(diǎn)劃分，將未采集到的圖像剩余點(diǎn)劃分到簇類中。
clc;clear all;close all;
%% Kmeans聚類圖像
K = 4;
im = imread('penguin.jpg');
im = imresize(im,[321,321],'nearest');
[R c lmd] = size(im);
elenum = numel(im(:,:,1));
total = [1:elenum];%總體；
dt = 0.1;%聚類中心迭代閾值；
for i = 1:size(im,lmd) %Cij表示第i波段成像，第j個
min_c = min(min(im(:,:,i)));
max_c = max(max(im(:,:,i)));
for k = 1:K
      C(i,k) = min_c + k*((max_c - min_c)/(K + 1));
end
end
C = double(C);
C_orgi = C;%初始聚類中心；
%修改圖像的存儲模式
points(elenum) = struct('r',0,'g',0,'b',0,'cluster',0);
for i = 1:R
for j = 1:c
      points((j-1)*R + i).r = im(i,j,1);
      points((j-1)*R + i).g = im(i,j,2);
      points((j-1)*R + i).b = im(i,j,3);
      points((j-1)*R + i).cluster = 0;
end
end
smplrate = 1/10000;
smplseq = sort(randperm(elenum,ceil(elenum*smplrate)),'ascend');%隨機(jī)性采樣M*N個圖像樣本點(diǎn)
cluster(K,100) =  struct('r',1,'c',1,'data',points(1));%存儲簇類每行存儲該簇類的樣本點(diǎn)
counter(K) = 0;%記錄每個簇類的樣本點(diǎn)個數(shù)；
totalnum = 100;%聚類迭代最大次數(shù)；
for iteration = 1:totalnum
for i = smplseq
for j = 1:K
      ptvalue(1,1) = points(i).r;
      ptvalue(2,1) = points(i).g;
      ptvalue(3,1) = points(i).b;
      distclust(j) = norm(double(ptvalue) - double(C(:,j)));%存儲樣本點(diǎn)到聚類的距離；
end
clunum = find(distclust == min(distclust));%樣本點(diǎn)距離聚類中心最小的簇類編號；
points(i).cluster = clunum;%標(biāo)記該樣本點(diǎn)類型；
counter(clunum) = counter(clunum) + 1;%樣本計數(shù)增加一個
cluster(clunum,counter(clunum)).r = mod(i,R);
cluster(clunum,counter(clunum)).c = floor(i/R)+1;
cluster(clunum,counter(clunum)).data = points(i);
end
%更新每個聚類中心的
for i = 1:K
sum = repmat(double(0),[3,1]);
for j = 1:counter(i)
      ptvalue(1,1) = cluster(i,j).data.r;
      ptvalue(2,1) = cluster(i,j).data.g;
      ptvalue(3,1) = cluster(i,j).data.b;
      sum = sum + double(ptvalue);
end
if  i == 1
      C_new = sum/counter(i);
else
      C_new = cat(2,C_new,sum/counter(i));
end
end
if norm(abs(C-C_new)) < dt
      break;
end
C = C_new;
end
% 參數(shù)顯示；
%disp('隨機(jī)樣本個數(shù)'),disp(size(smplseq,2));disp('初始聚類中心：');disp(C_orgi);disp('Kmeans 聚類新中心');disp(C);disp('聚類迭代次數(shù)');disp(iteration);
smpldiff = setdiff(total,smplseq);
for i = smpldiff
for j = 1:K
      ptvalue(1,1) = points(i).r;
      ptvalue(2,1) = points(i).g;
      ptvalue(3,1) = points(i).b;
      distclust(j) = norm(double(ptvalue) - double(C_new(:,j)));%存儲樣本點(diǎn)到聚類的距離；
end
clunum = find(distclust == min(distclust));%樣本點(diǎn)距離聚類中心最小的簇類編號；
points(i).cluster = clunum;%標(biāo)記該樣本點(diǎn)類型；
counter(clunum) = counter(clunum) + 1;%樣本計數(shù)增加一個
cluster(clunum,counter(clunum)).r = mod(i,R);
cluster(clunum,counter(clunum)).c = floor(i/R)+1;
cluster(clunum,counter(clunum)).data = points(i);
end
clusmap = zeros(R,c);
count = 0;
for i = 1:numel(points)
r = mod((i-1),R)+1;
c = ceil(i/R);%一維索引轉(zhuǎn)換為二維索引，按列
if points(i).cluster~=0
      clusmap(c,r) = points(i).cluster;
      count = count + 1;
end
end
clusmap = clusmap';
stats = regionprops(clusmap);
imshow(clusmap,[]);

回復(fù)此樓

» 收錄本帖的淘帖專輯推薦

» 猜你喜歡

288求調(diào)劑一志愿哈工大材料與化工已經(jīng)有5人回復(fù)
289求調(diào)劑已經(jīng)有6人回復(fù)
調(diào)劑求院校招收已經(jīng)有7人回復(fù)
070300化學(xué)354求調(diào)劑已經(jīng)有10人回復(fù)
材料專業(yè)申博已經(jīng)有5人回復(fù)
287求調(diào)劑已經(jīng)有17人回復(fù)
求調(diào)劑：085600材料與化工，考材科基，總分319 已經(jīng)有13人回復(fù)
一志愿大連理工大學(xué)，機(jī)械工程學(xué)碩，341 已經(jīng)有3人回復(fù)
一志愿南昌大學(xué)324求調(diào)劑已經(jīng)有5人回復(fù)
266分，求材料相關(guān)專業(yè)調(diào)劑已經(jīng)有12人回復(fù)

» 搶金幣啦！回帖就可以得到:

查看全部散金貼

1樓 2015-10-25 22:05:02

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

青青呀呀

新蟲 (正式寫手)

應(yīng)助: 0 (幼兒園)
金幣: 275.7
帖子: 817
在線: 74.8小時
蟲號: 4154145

★
小木蟲: 金幣+0.5, 給個紅包，謝謝回帖

我當(dāng)初看著Origin自己起碼會用用了，就夸大其詞的說，聽說Matlab更加高大上，我也要學(xué)，看著樓主天書般的帖子，我真想說，還好我只是說說而已…

發(fā)自小木蟲IOS客戶端

贊一下

回復(fù)此樓

7樓2015-10-25 22:29:02

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

查看全部 11 個回答

wake_up

木蟲 (著名寫手)

應(yīng)助: 124 (高中生)
金幣: 4014.5
帖子: 1113
在線: 424.2小時
蟲號: 670119

★
小木蟲: 金幣+0.5, 給個紅包，謝謝回帖

Matlab自帶的Kmeans支持并行計算，快

贊一下(1人)

回復(fù)此樓

8樓2015-10-26 13:13:19

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

xiyue6911

新蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 421.2
帖子: 79
在線: 35小時
蟲號: 3927522

引用回帖:

8樓: Originally posted by wake_up at 2015-10-26 13:13:19
Matlab自帶的Kmeans支持并行計算，快

謝謝一會兒試一試

發(fā)自小木蟲Android客戶端

贊一下

回復(fù)此樓

9樓2015-10-26 14:47:27

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

xiyue6911

新蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 421.2
帖子: 79
在線: 35小時
蟲號: 3927522

引用回帖:

8樓: Originally posted by wake_up at 2015-10-26 13:13:19
Matlab自帶的Kmeans支持并行計算，快

%% Kmeans聚類圖像；
im = imread('penguin.jpg');
im = imresize(im,[256,256]);
im_new = size(im(:,:,1));
im_r = im(:,:,1);im_r = im_r(

;
im_g = im(:,:,2);im_g = im_g(

;
im_b = im(:,:,3);im_b = im_b(

;
datas = double([im_r im_g im_b]);
opts = statset('Display','final');
[idx,ctrs] = kmeans(datas,4,'Distance','sqEuclidean','Options',opts);
im_new = reshape(idx,[256 256]);

與我寫的滿足一樣的功能，使用matlab就是快，還是謝謝

贊一下

回復(fù)此樓

10樓2015-10-26 22:10:08

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

簡單回復(fù)

云湖莎莎3樓

2015-10-25 22:05 回復(fù)

xiyue6911(金幣+1): 謝謝參與

已獲得1個金幣發(fā)自小木蟲IOS客戶端

nqqtmm5樓

2015-10-25 22:05 回復(fù)

xiyue6911(金幣+1): 謝謝參與

已獲得1個金幣發(fā)自小木蟲IOS客戶端

wang47640842樓

2015-10-25 22:05 回復(fù)

xiyue6911(金幣+1): 謝謝參與

已獲得1個金幣發(fā)自小木蟲Android客戶端

土蛋要考研4樓

2015-10-25 22:05 回復(fù)

xiyue6911(金幣+1): 謝謝參與

已獲得1個金幣發(fā)自小木蟲Android客戶端

2188xinhiu6樓

2015-10-25 22:05 回復(fù)

xiyue6911(金幣+1): 謝謝參與

已獲得1個金幣發(fā)自小木蟲Android客戶端

最具人氣熱帖推薦 [查看全部]	作者	回/看	最后發(fā)表

[考研] 289求調(diào)劑 +6	BrightLL 2026-03-29	6/300	2026-03-31 11:22 by oooqiao
[考研] 320分，材料與化工專業(yè)，求調(diào)劑 +10	一定上岸aaa 2026-03-27	14/700	2026-03-31 10:47 by foria
[考研] 求調(diào)劑 +5	11ggg 2026-03-30	5/250	2026-03-31 10:46 by Jaylen.
[考研] 284求調(diào)劑 +5	小熊～～ 2026-03-31	5/250	2026-03-31 10:22 by sophie2180
[考研] 266求調(diào)劑 +3	哇呼哼呼哼 2026-03-29	3/150	2026-03-31 10:06 by cal0306
[考研] 269求調(diào)劑 +4	我想讀研11 2026-03-31	4/200	2026-03-31 10:04 by cal0306
[考研] 材料求調(diào)劑一志愿哈工大總分298分，前三科223分 +10	dongfang59 2026-03-27	10/500	2026-03-30 23:42 by 果果媽咪
[考研] 279求調(diào)劑 +12	j的立方 2026-03-29	12/600	2026-03-30 20:30 by dick_runner
[考研] 346求調(diào)劑一志愿070303有機(jī)化學(xué) +7	蘿卜燉青菜 2026-03-28	8/400	2026-03-30 19:59 by 無際的草原
[考研] 材料化工340求調(diào)劑 +3	jhx777 2026-03-30	3/150	2026-03-30 17:54 by JourneyLucky
[考研] 324求調(diào)劑 +9	hanamiko 2026-03-26	11/550	2026-03-30 14:27 by JourneyLucky
[考研] 071010 323 分求調(diào)劑 +3	Baekzhy 2026-03-27	3/150	2026-03-30 14:24 by andresqi
[考研] 求調(diào)劑 +7	青春裁為三截 2026-03-29	7/350	2026-03-30 13:14 by laoshidan
[考研] 317求調(diào)劑 +10	蛋黃咸肉粽 2026-03-26	10/500	2026-03-30 09:45 by longlotian
[考研] 327求調(diào)劑 +6	汲亦昊 2026-03-29	6/300	2026-03-29 13:40 by peike
[考研] 070300化學(xué)求調(diào)劑 +4	起個名咋這么難 2026-03-27	4/200	2026-03-27 21:39 by 83503孫老師
[考研] 一志愿211院校 344分東北農(nóng)業(yè)大學(xué)生物學(xué)學(xué)碩，求調(diào)劑 +5	丶風(fēng)雪夜歸人丶 2026-03-26	8/400	2026-03-27 19:22 by 丶風(fēng)雪夜歸人丶
[考研] 348求調(diào)劑 +4	小懶蟲不懶了 2026-03-27	5/250	2026-03-27 12:47 by 果果媽咪
[考研] 315調(diào)劑 +4	0860求調(diào)劑 2026-03-26	5/250	2026-03-27 11:23 by wangjy2002
[考研] 309求調(diào)劑 +4	gajsj 2026-03-25	5/250	2026-03-26 00:27 by Dyhoer