hw3main.m

这是一个用Kmean处理图像的matlab文件

clear all;
close all;
img=imread('hw3p3_im.jpg');
i=rgb2gray(img);
i=double(i);

c=kmeansimage(i,1);
c=sort(c);
newimg1=i;
newimg1(find(i<=c(1)))=0;
newimg1(find(i>c(1)))=255;
imwrite(newimg1,'c1.jpg');

newimg=i;
c=kmeansimage(i,2);
c=sort(c);
newimg(find(i<=c(1)))=0;
for t=1:2
    newimg(find(i>c(t)))=ceil(255/2)*t;
end
imwrite(newimg,'c2.jpg');

newimg=i;
c=kmeansimage(i,3);
c=sort(c);
newimg(find(i<=c(1)))=0;
for t=1:3
    newimg(find(i>c(t)))=ceil(255/3)*t;
end
imwrite(newimg,'c3.jpg');

newimg=i;
c=kmeansimage(i,4);
c=sort(c);
newimg(find(i<=c(1)))=0;
for t=1:4
    newimg(find(i>c(t)))=ceil(255/4)*t;
end
imwrite(newimg,'c4.jpg');

newimg=i;
c=kmeansimage(i,5);
c=sort(c);
newimg(find(i<=c(1)))=0;
for t=1:5
    newimg(find(i>c(t)))=ceil(255/5)*t;
end
imwrite(newimg,'c5.jpg');

newimg=i;
c=kmeansimage(i,6);
c=sort(c);
newimg(find(i<=c(1)))=0;
for t=1:6
    newimg(find(i>c(t)))=ceil(255/6)*t;
end
imwrite(newimg,'c6.jpg');

newimg=i;
c=kmeansimage(i,7);
c=sort(c);
newimg(find(i<=c(1)))=0;
for t=1:7
    newimg(find(i>c(t)))=ceil(255/7)*t;
end
imwrite(newimg,'c7.jpg');

newimg=i;
c=kmeansimage(i,8);
c=sort(c);
newimg(find(i<=c(1)))=0;
for t=1:8
    newimg(find(i>c(t)))=ceil(255/8)*t;
end
imwrite(newimg,'c8.jpg');

newimg=i;
c=kmeansimage(i,9);
c=sort(c);
newimg(find(i<=c(1)))=0;
for t=1:9
    newimg(find(i>c(t)))=ceil(255/9)*t;
end
imwrite(newimg,'c9.jpg');

newimg=i;
c=kmeansimage(i,10);
c=sort(c);
newimg(find(i<=c(1)))=0;
for t=1:10
    newimg(find(i>c(t)))=ceil(255/10)*t;
end
imwrite(newimg,'c10.jpg');