位面分层实验记录

x=double(imread('C:\Users\Administrator\Desktop\Fig0314(a)(100-dollars).tif'));

[r c]=size(x);

subplot(3,3,1);

imshow(x,[]);

title('原始图像');

for k=1:8

     for i=1:r

        for j=1:c

           a(i,j)=bitget(x(i,j),k);

         end

     end

subplot(3,3,k+1);

imshow(a,[]);

title(['第',num2str(k),'个位平面']);

end