box_cuont_dif.m

Code matlab pour l extraction de la signature de l iris.Nous avons appliquer la méthode de box cuont

%function BCD=Box_cuont_dif(L)
clear all
clc
L=2;
%-----------Image au niveaux de gris----------%
x=imread('Resultats4/001_1_2.bmp');


x(129:520,512:520)=0;

figure(5),imshow(x)
[lin col]=size(x);

% 
% h(1:256)=0;
%     for i=1:lin
%         for j=1:col
%             h(x(i,j)+1)=h(x(i,j)+1)+1;        
%         end
%     end
%    % h(1)=0;
% figure,plot(h)
% k=length(h); 

%_____________________Box_cuonting differentiel____________________%
%__________________La Taille de la boite est " L "_________________%

c=1;
res=1;
moy_T=0;



for f=1:33

Boite=0;

limite_L=round((lin/L)-0.50);
limite_C=round((col/L)-0.50);

c1=1; 
moy_T=[];

for i=1:limite_L
c2=1;

    for j=1:limite_C
      

comp=1;
moy_Box=0;
mi=0;ma=0;
 %__________________________________Parcourir la boite______________________________%
        for k=c1:c1+L-1
            for m=c2:c2+L-1
                              
                boxe(comp)=x(k,m);
                comp=comp+1;
                
                
            end

        end 
%______________________________________________________________________________________%
       
for n=1:size(boxe,2)
    if(boxe(n)~=0)
        Boite=Boite+1;
        break;
    end
end

boxe;
    mi=min(boxe)/2;
    ma=max(boxe)/2;
    moy_Box=round((ma+mi));
        
        
        moy_T(res)=moy_Box;
        
        res=res+1;
        
        c2=c2+L;
    end
   c1=c1+L;
end

Total(f)=sum(moy_T)/Boite;

Total

%_______Evolution de la taille de la boite________%
L=L+1;
%_________________________________________________%

%______________________   Matrice_resultat pour l'apparition  _____________________%
X(f)=log(Boite);
Y(f)=(log(1/L));

%_________La dimention fractal________%
DimF(f)=X(f)/Y(f);

%_____________________________________%
end

NG=sum(Total)/size(Total,2);
NG
% ________________________________________%

figure(5),
plot3(-Y,DimF,NG);
grid on
P=s(Y,DimF,NG,2);
P
%__________________________________________%



% 
% 
% L=L+1;
% 
% % figure(3);
% % plot(Y,X);
% 


%----------------------------------------------------%