Matlab: seuillage avec méthode de Sauvola

I=imread('68JP.jpg');
figure(1),imshow(I);
I1=imcrop; %sélectionner une portion de l'image
ss=size(I1);%calculer la taille de l'image
Im=I1;
mf=zeros(ss(1),ss(2));

%parcours de l'image à seuiller
for i=1:15:ss(1)
 for j=1:15:ss(2)
     
     h=i;hi=i+15;
     s=j;sj=j+15;
%si la taile de la dernière portion de l'imge est inf à 15
     if(hi>ss(1))
         break;
     end
     if(sj>ss(2))
         break;
     end
%création d'une matrice vide de taille qui vaut à 15*15
     MM=zeros(hi-h,sj-s);

    %copier la portion 15*15 dans une matrice vide
    o=1;
     p=1;
     for m=h:hi
         for n=s:sj
              
                   MM(o,p)=Im(m,n);
                   o=o+1;
                   p=p+1;
               
         end
     end
     
    
     moyen=mean(mean(MM));
     ecart=std(std(MM));
    
     MM=zeros(hi-h+1,sj-s+1);

%calcul du seuil pour la portion 15*15
     seuil(m,n) = moyen *(1 + (1/10)*( 1-( ecart / 256) ));

%appliquer le seuillage sur l'image
k=1;
l=1;
     for k=h:hi
         for l=s:sj
                      disp(k);
             if(Im(k,l)=seuil(m,n))  %si le pixel sup ou egal à seuil, il devient blanc
                 mf(k,l)=256;
             end
         end
     end
 end
end

   
figure(2),imshow(mf);

M=MM(:)
moyen = mean(M)
ecart=std(M);