RGB=imread('miasto_NY.jpg') ;
I=rgb2gray(RGB);
I=double(I);
step=32;
partition = [step:step:256];
r = [0:step:256];
[m,n] = size(I);
I = reshape(I,m*n,1);
min=10000;
res=1:9;
for i=1:11
    [partition, codebook] = lloyds (I,r);
    [index,quant,distor] = quantiz(I,partition, codebook);
    res(i) = distor;
    if distor<min
        min=distor;
        bestP=partition;
        bestC=codebook;
        bestR=r;
    end
    r =[0:step:256];
    r=round(sort(randn(1,9)*step+r));
end
[index,quant,distor]=quantiz(I,bestP,bestC);
quant= uint8(reshape(quant,m,n));
imshow(quant);