Bit Plane Slicing

Aim:- Bit Plane Slicing
MSc Computer Science Image Processing Practical No. 13
Index of all Practicals ~ Click Here

Code:-

% **** 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:-