btplane slicing of image

how to decompose a image into 16planes?? is it possible??

Answer 1

[Omega]

Hi Surya,

Bit plane slicing involves decomposition of an image into its constituent bit planes which represent the bits of each pixel value. Standard 8-bit images have 8 bit per pixels and this results in 8 bit planes.

However, it is possible to decompose an image into 16 planes if the image has a higher bit depth, such as 16 bits per pixel.

Here is a sample code for decomposing an image into 16 bit planes:

it = imread('sample_16bit_image.png');   % read the image

% Assuming the image is stored in a 16-bit format
itemp = it(:,:,1);
[r, c] = size(itemp);         % get the dimensions of image
s = zeros(r, c, 16);          % pre allocate a variable to store 16 bit planes of the image

for k = 1:16
    for i = 1:r
        for j = 1:c
            s(i, j, k) = bitget(itemp(i, j), k);    % get kth bit from each pixel 
        end
    end
end

figure, imshow(uint16(itemp)); title('Original Image');    % display original image

figure;                                                    % display all the 16 bit planes

% Adjust the subplot layout to accommodate 16 planes
for k = 1:16
    subplot(4, 4, k);
    imshow(s(:,:,17-k));  % Display from MSB to LSB
    title(sprintf('%d Plane', 17-k));
end

% Reconstruct the original image from generated bit planes
rec = zeros(r, c);
for k = 1:16
    rec = rec + s(:,:,k) * 2^(k-1);
end

figure, imshow(uint16(itemp)); title('Recovered Image'); % display the reconstructed image

The code is inspired from the "Bit-Plane Slicing" by Akshay Bhosale.

https://www.mathworks.com/matlabcentral/fileexchange/48844-bit-plane-slicing