图像灰度的双三次插值的MATLAB实现

相比C/C++实现，图像灰度的双三次插值的MATLAB实现要方便的多，下面是MATLAB语言实现

clc,clear;
ff=imread('C:\Program Files\MATLAB\R2013a\bin\work\lena.bmp');
[mm,nn]=size(ff);%将图像隔行隔列抽取元素，得到缩小的图像f
m=mm/2;n=nn/2;
f=zeros(m,n);
for i=1:m
    for j=1:n
        f(i,j)=ff(2*i,2*j);
    end
end
k=3;%设置放大倍数
a=f(1,:);%取f的第1行
c=f(m,:);%取f的第m行
%将待插值图像矩阵前后各扩展两行两列,共扩展四行四列到f1
b=[f(1,1),f(1,1),f(:,1)',f(m,1),f(m,1)];
d=[f(1,n),f(1,n),f(:,n)',f(m,n),f(m,n)];
a1=[a;a;f;c;c];
b1=[b;b;a1';d;d];
ffff=b1';f1=double(ffff);
g1=zeros(k*m,k*n);
for i=1:k*m%利用双三次插值公式对新图象所有像素赋值
    u=rem(i,k)/k;
    i1=floor(i/k)+2;%floor()向左取整，floor(1.3)=floor(1.7)=1
    A=[sw(1+u) sw(u) sw(1-u) sw(2-u)];
    for j=1:k*n
        v=rem(j,k)/k;j1=floor(j/k)+2;
        C=[sw(1+v);sw(v);sw(1-v);sw(2-v)];
        B=[f1(i1-1,j1-1) f1(i1-1,j1) f1(i1-1,j1+1) f1(i1-1,j1+2); 
           f1(i1,j1-1) f1(i1,j1) f1(i1,j1+1) f1(i1,j1+2);
           f1(i1+1,j1-1) f1(i1+1,j1) f1(i1+1,j1+1) f1(i1+1,j1+2);
           f1(i1+2,j1-1) f1(i1+2,j1) f1(i1+2,j1+1) f1(i1+2,j1+2)];
        g1(i,j)=(A*B*C);
    end
end
g=uint8(g1);

imshow(uint8(f));

title('缩小的图像');%显示缩小的图像
figure,imshow(ff);title('原图');%显示原图像
figure,imshow(g);title('双三次插值放大的图像');%显示插值后的图像

运行结果