用MATLABR2016b的bitget将灰度图像Lena.bmp图像 分解成8个位平面，分别去掉从1到7个位平面，与剩余的位平面图像比较的源代码

好的，下面是MATLAB R2016b的代码，可以实现将灰度图像Lena.bmp图像分解成8个位平面，并分别去掉从1到7个位平面，并与剩余的位平面图像进行比较：

% 读取Lena.bmp图像
lena = imread('Lena.bmp');

% 将图像转换为二进制图像
lena_bin = de2bi(lena, 8);

% 分解成8个位平面
bit_planes = cell(1, 8);
for i = 1:8
    bit_planes{i} = lena_bin(:, :, i);
end

% 分别去掉从1到7个位平面，并与剩余的位平面图像进行比较
for i = 1:7
    % 去掉第1到第i个位平面
    removed_planes = bit_planes;
    for j = 1:i
        removed_planes{j} = zeros(size(removed_planes{j}));
    end
    
    % 重构图像
    removed_lena = bi2de(cat(3, removed_planes{:}));
    
    % 与剩余的位平面图像进行比较
    remaining_planes = bit_planes;
    for j = i+1:8
        remaining_planes{j} = zeros(size(remaining_planes{j}));
    end
    remaining_lena = bi2de(cat(3, remaining_planes{:}));
    diff = abs(removed_lena - remaining_lena);
    
    % 显示去掉第1到第i个位平面后的图像和与剩余的位平面图像比较的结果
    figure;
    subplot(1, 3, 1);
    imshow(removed_lena);
    title(sprintf('Removed bit planes %d-%d', 1, i));
    subplot(1, 3, 2);
    imshow(remaining_lena);
    title(sprintf('Remaining bit planes %d-%d', i+1, 8));
    subplot(1, 3, 3);
    imshow(diff, []);
    title('Difference');
end

这段代码将灰度图像Lena.bmp图像分解成8个位平面，然后分别去掉从1到7个位平面，并与剩余的位平面图像进行比较。最后显示去掉第1到第i个位平面后的图像和与剩余的位平面图像比较的结果。