用于图像像素改组的2D Baker图matlab代码

Question

我正在尝试实施离散化的面包师地图以进行像素的改组，我无法清楚地理解它的离散化版本。我也无法在matlab上正确实现它总是显示索引越界错误。 链接在这里 http://link.springer.com/chapter/10.1007/978-3-540-95972-4_16 请有人帮忙。我是初学者，这是我的努力，我知道这是可悲的但仍然。

clc; 
Image=imread('C:\Users\AAA\Documents\lena_gray.gif');
[y,ymap]=imread('C:\Users\AAA\Documents\lena_gray.gif');
im=ind2gray(y,ymap);
%im=imresize(im,0.25);
I=im;
%ite=input('enter the number of iterations >');
[row,col]=size(im);

par=2;
for i=row:-1:1
for j=1:col
     if((1<=j)&&(j<257))
         rnew=floor(0.5*(i-mod(i-1,2)))+1;
         c_new=2*(j-1)+mod(i-1,2)+1;
      else
          rnew=floor(0.5*(i-mod(i-1,2)))+257;
          c_new=2*(j-257)+mod(i-1,2)+1;
     end
     encim(i,j)=im(rnew,c_new);
end
end 
encim=uint8(encim)
figure,imshow(encim)
impixelinfo
figure,imhist(encim)
if(imhist(encim)==imhist(I));
g=1
else
g=0
end