比特平面

originalImg = imread('Fig0314(a)(100-dollars).tif');
tempImg = originalImg;
figure;
subplot(3,3,1);
imshow(originalImg);
title('原始图像');
height = size(originalImg,1);
width = size(originalImg,2);
for n = 1:8
    for i=1:height
        for j=1:width
            gray =bitand( originalImg(i,j), 2^(n-1) );%位与操作判断
           if(gray==2^(n-1))
                tempImg(i,j) = 255;%二值化 突出比特平面            
            else
                tempImg(i,j) = 0;
           end                          
        end
    end
    subplot(3,3,n+1);
    imshow(tempImg);
    title(['第',num2str(n),'比特图像 ']);
end