Image segmentation by k-means algorithm
4 次查看(过去 30 天)
显示 更早的评论
I have an code for k-means segmentation but I have some problems when I applying it on my multispectral satellite image I facing error in dimensions (Subscripted assignment dimension mismatch), how can I solve this problem?
% Grayscale Image Segmentation Using K-Means Algorithm Function Kmeans segmentation
clc
close all
clear all
[im,map]=imread('pp1.bmp');
im=im2double(im);
[row,col]=size(im);
% number of clusters
nc=4;
% initial random cluster centroids
cs=rand(nc,1);
pcs=cs;
% number of iteration
T=50;
t=0;
D=zeros(row,col,nc);
tsmld=[];
eps=1.e-5;
cmx=1;
while (t<T && cmx>eps)
%Distance between centroids and image's pixel
for c=1:nc
D(:,:,c)=(im-cs(c)).^2;
end
%assign members (image pixels)to minimum distance clusters
[mv,ML]=min(D,[],3);
%updat cluster centroid
for c=1:nc
I=(ML==c);
cs(c)=mean(mean(im(I)));
end
%find maximum absolute difference between crrent and previous iteration
%cluster centroids
cmx=max(abs(cs-pcs));
pcs=cs;
t=t+1;
%sum difference between centroid and their members and store it for
%plotting energy minimization functions
tsmld=[tsmld; sum(mv(:))];
end
% assign a colour to each cluster
colors=hvs(nc);
sim=colors(ML,:);
sim=reshape(sim,row,col,3);
figure,subplot(1,2,1),imshow(im,map);
title('Input Image: pp1');
subplot(1,2,2);imshow(sim,map);
title('segmented Image:pp1')
figure;plot(tsmld,'*-b')
xlabel('Iteration');ylabel('Energy');
title('K-means energy minimization-pp1');
![](https://www.mathworks.com/matlabcentral/answers/uploaded_files/170206/image.bmp)
0 个评论
采纳的回答
Image Analyst
2018-2-11
I don't even see a call to kmeans() in your code. Adapt my attached example to however many spectral bands you want to use (it should be obvious how to do it).
10 个评论
Federico Alvarez
2024-1-21
Unrelated / related - thank you all (specially Image Analyst). I fell into a similar issue, and was suspecting that the way I was arranging my data was the problem. Looking at your code samples helped me confirm the problem.
更多回答(1 个)
Jaswinder Singh
2019-4-16
can anybody provide me the code to find coefficient of a hysteresis loop whose coordinates are given to me. data is attached
0 个评论
另请参阅
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!