BIT PLANE SLICING

% **** BIT PLANE SLICING ****

a=imread('Moon.tif');

b=double(a);

[row col]=size(b);

%r=input('enter bitval :  ');

for x=1:row

    for y=1:col

        c=dec2bin(b(x,y),8);

        d=c(1);

        w1(x,y)=double(d);

        if w1(x,y)==49

            w1(x,y)=255;

        else w1(x,y)=0;

        end

    end

end

for x=1:row

    for y=1:col

        c=dec2bin(b(x,y),8);

        d=c(2);

        w2(x,y)=double(d);

        if w2(x,y)==49

            w2(x,y)=255;

        else w2(x,y)=0;

        end

    end

end

for x=1:row

    for y=1:col

        c=dec2bin(b(x,y),8);

        d=c(3);

        w3(x,y)=double(d);

        if w3(x,y)==49

            w3(x,y)=255;

        else w3(x,y)=0;

        end

    end

end

for x=1:row

    for y=1:col

        c=dec2bin(b(x,y),8);

        d=c(4);

        w4(x,y)=double(d);

        if w4(x,y)==49

            w4(x,y)=255;

        else w4(x,y)=0;

        end

    end

end

for x=1:row

    for y=1:col

        c=dec2bin(b(x,y),8);

        d=c(5);

        w5(x,y)=double(d);

        if w5(x,y)==49

            w5(x,y)=255;

        else w5(x,y)=0;

        end

    end

end

for x=1:row

    for y=1:col

        c=dec2bin(b(x,y),8);

        d=c(6);

        w6(x,y)=double(d);

        if w6(x,y)==49

            w6(x,y)=255;

        else w6(x,y)=0;

        end

    end

end

for x=1:row

    for y=1:col

        c=dec2bin(b(x,y),8);

        d=c(7);

        w7(x,y)=double(d);

        if w7(x,y)==49

            w7(x,y)=255;

        else w7(x,y)=0;

        end

    end

end

for x=1:row

    for y=1:col

        c=dec2bin(b(x,y),8);

        d=c(8);

        w8(x,y)=double(d);

        if w8(x,y)==49

            w8(x,y)=255;

        else w8(x,y)=0;

        end

    end

end

subplot(3,3,1);

imshow(uint8(w1));

subplot(3,3,2);

imshow(uint8(w2));

subplot(3,3,3);

imshow(uint8(w3));

subplot(3,3,4);

imshow(uint8(w4));

subplot(3,3,5);

imshow(uint8(w5));

subplot(3,3,6);

imshow(uint8(w6));

subplot(3,3,7);

imshow(uint8(w7));

subplot(3,3,8);

imshow(uint8(w8));

subplot(3,3,9);

imshow(a);

Output

Bit Plane Slicing