[I,map] = imread('bw_L.gif');
J = ind2gray(I,map);
imhist(J)
K = J>.45;
imshow(K)