% Fig. 2.27
% Setting the least significant bit to zero, then scale
% the difference to [0,255] for clarity

I = imread('Fig0227(a)(washington_infrared).tif');
J = bitset(I, 1, 0);
Diff = I - J;
D = 255*(I - J);
imshow(D)