clear all;
clc;
 
InPic=imread('C:\test.jpg');
[MM,NN,ZZ] = size(InPic);
if(ZZ>1)
  pic=uint8(rgb2gray(InPic));
else
  pic=uint8(InPic);
end;
%========================================================
% Preprocessing: Binarized input image
%========================================================
for i=1:MM-40
for j=1:NN-40
if(pic(i+20,j+20)<128)
          imgBinary(i,j)=0;
end;
if(pic(i+20,j+20)>=128)
          imgBinary(i,j)=255;
end;
end;
end;
%========================================================
% Computing fractal dimension with 5 scales
%========================================================
imgBinary=mat2gray(imgBinary);
imshow(imgBinary); 
[M,N]=size(imgBinary);
for k=1:5   % scale
  Sz=2*k;
  sum=0;
for i=1:Sz:M-Sz
for j=1:Sz:N-Sz
          flag=0;
          mask((1:Sz),(1:Sz))=imgBinary((1:Sz)+i-1,(1:Sz)+j-1);
for p=1:Sz
for q=1:Sz
if(mask(p,q)==0)                    
                      flag=flag+1;
end;
end;
end;
if(flag>=Sz)
              sum=sum+1;
end;
end;
end;
  w(k,1)=Sz;
  w(k,2)=sum;
  x(k,1)=log10(1/Sz);
  y(k,1)=log10(sum);
end;
 
m = robustfit(x,y) ;
%figure,plot(x,y,'o');
%figure,plot(x,y,'-');
result=m(2,1);
disp(result);