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

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

[Jan]

What is your question? "has some errors" is not helpful.

isempty(lowap(:,:,i)) cannot be true, because isempty checks, if teh input has elements.

Some code is strange:

    for i=1:k
    upap(:,:,i)=lowap(:,:,i);
    end

What about:

    upap = lowap;

[saradindu rana]

actually the whole code is wrong..so i need the complete code for rough k means can you plzz provide.

[Jan]

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]

Hi

did you finally find Rough kmeans matlab code?

Answer 1

[John D'Errico]

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.

[saradindu rana]

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]

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.

[saradindu rana]

sir i need rough k means code for image segmentation..not k means ..i need rough k means