i need a rough k means code anyone help me plzz ...the following code has some errors.

Question

saradindu rana 2019-2-19

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/445804-i-need-a-rough-k-means-code-anyone-help-me-plzz-the-following-code-has-some-errors

评论： ali seraj 2019-11-8

function [rkm] = Km(I,K,wlower,threshold)
I    = imread('1.jpg'); 
I = im2double(I);                   % Load Image
F = reshape(I,size(I,1)*size(I,2),3);
wupper=1-wlower;
K    = 8;
F_ORG=F;
N=fix(size(F,1)/k);
for i=1:k
    lowap(:,:,i) = F( ceil(rand(N,1)*size(F,1)) ,:);
    F(ismember(F,lowap(:,:,i),'rows'),:)=[]
end
lowap(:,:,i)=(lowap(:,:,i);F);
    for i=1:k
    upap(:,:,i)=lowap(:,:,i);
    end
DAL   = zeros(size(F_ORG,1),K+2);                         % Distances and Labels
KMI   = 10;                                           % K-means Iteration
for n = 1:KMI
   for i = 1:size(F,1)   
upap(:,:,i)(ismember(upap(:,:,i),lowap(:,:,i),'rows'),:)=[];          
if (~isempty(lowap(:,:,i)))&& (isempty(upap(:,:,i))) 
    cents(k,size(F,1))=sum(lowapp(:,:,i),1)/size(lowap(:,:,i));
elseif (isempty(lowap(:,:,i))) && (~isempty(upap(:,:,i)))
    cents(k,size(F,1)=sum(upap(:,:,i)/size(upap(:,:,i);
else
    cents(k,size(F,1)=wlower*(sum(lowap(:,:,i))/size(lowap(:,:i)))+wupper*(sum(upap(:,:,i))/size(upap(:,:,i)));
         for j = 1:K  
            DAL(i,j) = norm(F(i,:) - cents(j,:));      
         end
[Distance CN] = min(DAL(i,1:K));                % 1:K are Distance from Cluster Centers 1:K 
DAL(i,K+1) = CN;                                % K+1 is Cluster Label
DAL(i,K+2) = Distance;                          % K+2 is Minimum Distance
   for a=1:k
       for b=1:k
       if DAL(:,a)/DAL(:,b)<=threshold
end
X = zeros(size(F_ORG));
for i = 1:K
idx = find(DAL(:,K+1) == i);
X(idx,:) = repmat(CENTS(i,:),size(idx,1),1); 
end
rkm = reshape(X,size(I,1),size(I,2),3);
end

4 个评论
显示 2更早的评论隐藏 2更早的评论

Jan 2019-2-20

Does this mean, that you have posted a completely flawed code? Why?

What's wrong with using one of the many published codes for a valid KNN-search? Take a look in the FileExchange.

ali seraj 2019-11-8

Hi

did you finally find Rough kmeans matlab code?

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

John D'Errico 2019-2-20

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/445804-i-need-a-rough-k-means-code-anyone-help-me-plzz-the-following-code-has-some-errors#answer_361830

编辑：John D'Errico 2019-2-20

Easy, peasy. It is right here:

https://www.mathworks.com/products/statistics.html

Writing your own code to do what others have done far better than you know how to do is a bad idea. And it costs what? Compared to the amount of your own time you have already wasted should say something. Is your own time worth nothing at all?

Or do a search for kmeans on the file exchange. I'd bet you will find at least one code. You would be better off getting the stats TB.

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

saradindu rana 2019-2-20

i have already wasted time..but i cant camplete..k means code in matlab is availabe but i have not found any codes for rough k means.

请先登录，再进行评论。

Answer 2

Image Analyst 2019-2-20

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/445804-i-need-a-rough-k-means-code-anyone-help-me-plzz-the-following-code-has-some-errors#answer_361831

Not sure what this confusing, uncommented code does. But if you want an example of how kmeans can do a segmentation of intensity in images, see the attached.

That said, I think that in most cases, kmeans is a lousy way to segment images, as you'll probably find out, so don't come back and say "Hey kmeans did a lousy job segmenting my images."

By the way, no sense in passing in the badly-named I when I is just overwritten by calling imread() in the function.