www.pudn.com > colorseg.zip > test1.m, change:2003-03-03,size:430b

im = imread('images/berries.jpg','jpg');
im = double(im) ./ 256;

load clusters % clusters was created using the previous command

M = calculatePXTheta(im,clusters);

figure(1)
colormap gray;

for i=1:2:10
  subplot(5,2,i)
  imagesc(clusters(:,:,(i+1)/2));
  title(sprintf('Cluster %d',(i+1)/2));
end

for i=2:2:10
  subplot(5,2,i)
  imagesc(M(:,:,i/2)); colorbar
  title(sprintf('Probability of being in cluster %d',i/2));
end