function G = granulometrie (I,p,alpha,K,KK); 

% function G = granulometrie(I,p,alpha,K,KK); 
%
% Calcule (et affiche) la granulometrie d'une image de taille (N1,N2) (pairs) 
% pour les boules de bi-rayons (N1/alpha^(K+1-k),N2/alpha^(K+1-k), k=1,2,...,KK. 
% (KK inferieur ou egal a K) lorsque 1<alpha <=2. 

[N1,N2] = size(I); 
G = zeros (N1,N2,KK); 
for k=1:KK; 
    ak = N1/alpha^(K+1-k); 
    bk = N2/alpha^(K+1-k); 
    G(:,:,k) = ouverture (I,ak,bk,p); 
    figure
    imshow(G(:,:,k));
end 
end

