%双三次插值具体实现
clc,clear;
fff=imread('E:\Documents\BUPT\DIP\图片\lena.bmp');
ff =
rgb2gray(fff);%转化为灰度图像
[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=5; %设置放大倍数
bijiao1 =
imresize(f,k,'bilinear');%双线性插值结果比较
bijiao =
uint8(bijiao1);
a=f(1,:);c=f(m,:); %将待插值图像矩阵前后各扩展两行两列,共扩展四行四列,为什么要扩展
%出四行四列,是因为,后面插值的时候为了保证插值的数保证在i1的后面,同时符合插值算法而设计的。
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);