理论知识参考: https://blog.youkuaiyun.com/huang1024rui/article/details/46545329
程序:这个程序于我参考的文章的程序差不多,不一样的地方在于边界处理,
上边链接里采用方法1,扩成256*256,我的是方法2,扩成257*257
%%%%%双线性插值
clc
clear
%function [image1] = Bilinear_(image,x,y)
I=imread('E:\科研专用\散斑图平移\Image00.bmp');
[IH,IW,ID]=size(I);
zmf=2;
ZIH=round(IH*zmf);
ZIW=round(IW*zmf);
I1=zeros(ZIH,ZIW,ID);
IT=zeros(IH+1,IW+1,ID);%为什么+2
IT(1:IH,1:IW,:)=I;
% IT(1,2:IW+1,:)=I(1,:,:);
% IT(IH+2,2:IW+1,:)=I(IH,:,:);
% IT(2:IH+1,1,:)=I(:,1,:);
% IT(2:IH+1,IW+2,:)=I(:,IW,:);
% IT(1,1,:) = I(1,1,:);
% IT(1,IW+2,:) = I(1,IW,:);
% IT(IH+2,1,:) = I(IH,1,:);
% IT(IH+2,IW+2,:) = I(IH,IW,:);
IT(IH+1,1:IW)=I(IW,:);
IT(1:IH,IW+1)=I(:,IH);
IT(IH+1,IW+1)=I(IH,IW);
for zj = 1:ZIW % 对图像进行按列逐元素扫描
for zi = 1:ZIH
ii = (zi+1)/zmf;
jj = (zj+1)/zmf;
i = floor(ii);
j = floor(jj); % 向下取整
u = ii - i; v = jj - j;
% i=i+1;
% j=j+1;
ZI(zi,zj,:) = (1-u)*(1-v)*IT(i,j,:) +(1-u)*v*IT(i,j+1,:)+ u*(1-v)*IT(i+1,j,:) +u*v*IT(i+1,j+1,:);
end
end
ZI = uint8(ZI);
figure
imshow(I);
axis on
title(['原图像(大小: ',num2str(IH),'*',num2str(IW),'*',num2str(ID),')']);
figure
imshow(ZI);
axis on
title(['缩放后的图像(大小: ',num2str(ZIH),'*',num2str(ZIW),'*',num2str(ID)',')']);