% cells.m generates Poisson patches of different colors

m = 100;                                      % avg number of points
s = 0;
N = -1;
while (s<m)
  s = s - log(rand);
  N = N + 1;                                  % N ~ Poisson(m)
end

a = rand(N,1);
b = rand(N,1);
c = rand(N,1);                                % grayscale levels of points

clf
plot(a,b,'*')
axis([0 1 0 1]);
hold on
drawnow
pause

voronoi(a,b)
axis([0 1 0 1]);
drawnow

gridsize = [100 100];                           % number of grid points in x, y
gridrange= [0 1 0 1];                         % default [xmin xmax ymin ymax]
[x,y]    = easygrid(gridrange, gridsize);     % matrices of x and y values
Z        = zeros(gridsize);

tri = delaunay(a,b);
for i=1:gridsize(1),
  for j=1:gridsize(2),
    k = dsearch(a,b,tri,x(i,j),y(i,j));
    Z(i,j) = c(k);
  end
end
 
clf
colormap(gray);
axis([0 1 0 1]);
hold on
pcolor(x,y,Z)
shading flat