Bit plane decomposition to bit plane decoder

I am having code of bit plane decomposition.Please give me some ideas to generate codes for the decoder .

     if true
        % X=double(imread('lenna.gif'));
    coef=X;
    [nRow, nColumn] = size(X)
    M=size(coef,1);
    WW=zeros(M,M,8);
    for l=1:8 %λƽÃæµÄ²ãºÅ
      for m=1:M  %MΪϵÊý¾ØÕóµÄ³ß´ç
        for n=1:M
            if    mod(fix(coef(m,n)/2^(l-1)),2)==1
                  WW(m,n,l)=1;
            else  WW(m,n,l)=0;
            end

        end
      end
    end
    subplot(3,3,1);imshow(coef,[]);
    Bit1=WW(:,:,1);
    subplot(3,3,2);imshow(Bit1,[]);

    Bit2=WW(:,:,2);
    subplot(3,3,3);imshow(Bit2,[]);

    Bit3=WW(:,:,3);
    subplot(3,3,4);imshow(Bit3,[]);
    Bit4=WW(:,:,4);
    subplot(3,3,5);imshow(Bit4,[]);
    Bit5=WW(:,:,5);
    subplot(3,3,6);imshow(Bit5,[]);
    Bit6=WW(:,:,6);
    subplot(3,3,7);imshow(Bit6,[]);
    Bit7=WW(:,:,7);
    subplot(3,3,8);imshow(Bit7,[]);

    Bit8=WW(:,:,8);
    subplot(3,3,9);imshow(Bit8,[]);
    MSE = sum(sum((X).^2))/nRow / nColumn;
    Q = 255;
    fprintf('The psnr performance is %.2f dB\n', 10*log10(Q*Q/MSE));
      end

Will the reverse process helps in bit plane decoder?

[Walter Roberson]

Sorry, your code has only two comments, both of which appear to be in the 1512 character set. This is not sufficient for us to determine what variety of bitplane decomposition you are trying to do or what algorithm you think you are following.