test1.m

来自「calculatePXTheta---Calculate the probabi」· M 代码 · 共 22 行

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