双三次插值放大2倍

clear all
close all
clc

filename='1.bmp';

data1=imread(filename);
imshow(data1);

[m,n]=size(data1);
%将待插值图像矩阵前后各扩展两行两列,共扩展四行四列  
data2=[data1(1,:);data1(1,:);data1;data1(end,:);data1(end,:)];
data3=[data2(:,1),data2(:,1),data2,data2(:,end),data2(:,end)];
f1=double(data3);

g1 =zeros(2*m,2*n);  

for i=1:2*m                 %利用双三次插值公式对新图象所有像素赋值  
   i1=floor(i/2)+2;  

   if rem(i,2)==1
       A=[-0.125 0.625 0.625 -0.125];
   else
       A=[0 1 0 0];
   end

    for j=1:2*n  
        j1=floor(j/2)+2;  

        if rem(j,2)==1
            C=[-0.125;0.625;0.625;-0.125];
        else
            C=[0;1;0;0];
        end     

        B=f1(i1-1:i1+2,j1-1:j1+2);
        g1(i,j)=(A*B*C);
    end  
end  
g=uint8(g1); 

figure();
imshow(g);