clc;
clear;
X=imread('1.jpg');

Ib=rgb2gray(X);

G = fspecial('gaussian',[11 11],3);

filt = imfilter(Ib,G,'symmetric');

F1=im2double(filt);
imshow(filt)

I=Ib-filt;
imshow(I);

k=0.55;
p=k*max(max(I));
[x y]=size(I);
num=numel(I);
for n=1:num
   if I(n)>p F(n)=255;
       else F(n)=0;
   end
end
F;
bw=reshape(F,[x,y]);
figure, imshow(bw);