我写的 matlab 二次差值 图像 旋转

[x,cmap]=imread('C:/Documents and Settings/Administrator/桌面/matlab/xiaoxiang6.jpg');
% image(x)
[m,n,p]=size(x);
im1=zeros(m,n)
im2=zeros(m,n)
im3=zeros(m,n)
x=double(x)
for i=1:m
    for j=1:n
       a=x(i,j,1);
        im1(i,j)=a;
        b=x(i,j,2);
        im2(i,j)=b;
        c=x(i,j,3);
        im3(i,j)=c;
    end
end
 
col=500;  %行   屏幕的 大小
row=500;   %列
im=zeros(col,row,3);
for i=1:m                % 将 原图像在新的画布里移动一小段距离，便于后面的 旋转，以免出现画面不完整的情况
    for j=1:n
        a=im1(i,j);
        im(i +150,j,1)=a;
        b=im2(i,j);
        im(i+150 ,j,2)=b;
        c=im3(i,j);
        im(i+150 ,j,3)=c;
    end;
end;
 newim=zeros(col,row,3);
for   ii=1:500
    for jj=1:500
         i=(sqrt(3)*ii+jj)/2;                               %  反向用双线性插值法  求出新的坐标在原来坐标系中的位置，正数的话就直接赋值，
         j=(6*jj-2*sqrt(3)*ii)/(4*sqrt(3));                 %不是的话要进行插值运算求出此点的rgb各个矩阵的元素值然后把此值赋给新的坐标
         if i<=150+m && i>=150 && j<=n && j>=0
               if i==round(i) && j==round(j)
        if i==0                                             %避免下标出现0的情况
            i=1;
        end;
        if j ==0
            j =1;
        end;
           
        a=im(i,j,1);
        newim(ii,jj,1)=a;                                 % 正数就直接赋值
        b=im(i,j,2);
        newim(ii,jj,2)=b;
        c=im(i,j,3);
        newim(ii,jj,3)=c;
    else
        a=i-round(i-0.5);                                %不是正数 要进行插值运算。
        b=j-round(j-0.5);                               
        c=round(i-0.5);                                   %这里要进行计算与这个不是正数坐标点相邻的点的坐标  还可以吧  哈哈 自己发明的
        d=round(j-0.5);
        if c==0
            c=1
        end;
        if d==0
            d=1;
        end;                                            %底下是求新的坐标的rgb各个矩阵的元素值大小  然后反向赋值
       
        a1=round(( im(c+1,d,1)- im(c,d,1))*a+( im(c,d+1)- im(c,d,1))*b+( im(c+1,d+1,1)+ im(c,d,1)- im(c,d+1,1)- im(c+1,d,1))*a*b+ im(c,d,1));
        newim(ii,jj,1)=a1;
        b1=round(( im(c+1,d,2)- im(c,d,2))*a+( im(c,d+2,2)- im(c,d,2))*b+( im(c+1,d+1,2)+ im(c,d,2)- im(c,d+1,2)- im(c+1,d,2))*a*b+ im(c,d,2));
        newim(ii,jj,2)=b1;
        c1=round(( im(c+1,d,3)- im(c,d,3))*a+( im(c,d+1,3)- im(c,d,3))*b+( im(c+1,d+1,3)+ im(c,d,3)- im(c,d+1,3)- im(c+1,d,3))*a*b+ im(c,d,3));
        newim(ii,jj,3)=c1;
    end;
end;
       
   
       
    end;
end;

   newim=uint8(newim);
   im=uint8(im)
  newim                                                                     % 一定要强制转化为  uint8类型  ，才能 用image() 函数
   image(newim);