function res=DKu(im1,im2)
%im1,im2 : 2 images à comparer
%nx,ny   : largeur et hauteur des images
%
% Utilitaires pour l'évaluation de la segmentation d'images
% Toolbox matlab (version 5.3)
%
% (c) Laboratoire de Vision et Robotique (UPRES EA 2078)
% ENSI de Bourges - Université d'Orléans
%
% Sébastien Chabrier : sebastien.chabrier@ensi-bourges.fr
%
%
% Utilitaires pour l'évaluation de la segmentation d'images
% Toolbox matlab (version 5.3)
%
% (c) Laboratoire de Vision et Robotique (UPRES EA 2078)
% ENSI de Bourges - Université d'Orléans
%
% Sébastien Chabrier : sebastien.chabrier@ensi-bourges.fr
%
% Si vous utilisez cette toolbox, veuillez citer ce papier svp.
%
%S. Chabrier, H. Laurent, C. Rosenberger, Y.J. Zhang, 
%"Supervised evaluation of synthetic and real contour segmentation results" , 
%European Signal Processing Conference (EUSIPCO) 2006.

contour=0;
fond=255;

[a,b]=find(im1==contour & im2==contour);
Nand=length(a);

[a,b]=find(im1==contour);
Ni=length(a);

[a,b]=find(im2==contour);
Na=length(a);

%texte=strcat('Nombre de pixels contour Image segmentée1: ',int2str(Ni));
%disp(texte);

%texte=strcat('Nombre de pixels contour Image segmentée2: ',int2str(Na));
%disp(texte);

%texte=strcat('Nombre de pixels communs: ',int2str(Nand));
%disp(texte);

%texte=strcat('Nombre de pixels sursegmentés: ',int2str(Na-Nand));
%disp(texte);

%texte=strcat('Nombre de pixels soussegmentés: ',int2str(Ni-Nand));
%disp(texte);

im1=double(im1)+1;
im2=double(im2)+1;
tmp=im1-im2;

tmp=tmp.*log10(im1./im2);
   
somme=sum(sum(tmp));
taille=size(im1);
res=(1/(taille(1)*taille(2)))*somme;
%disp(strcat('DKu=',num2str(res,10)));