MATLAB实现双线性内插算法

用MATLAB简易的实现双线性内插算法算法，实现图像的缩放。具体的细节，程序注释里面里面会有。

本文这里主要是实现双线性内插，如果需要仔细深究一下原理，请参考：

讲的还是清晰易懂 

本实现代码，也是参考上面这篇博文的讲解实现的

% 双线性插值实现
% f(x,y)=w1*p1+w2*p2+w3*p3+w4*p4 p1,p2,p3,p4是与(x,y)最近的四个像素点
% w1,w2,w3,w4分别是对应的权重
% src是原图像（M*N）,dst是目标图像（P*Q）
% 缩放因子是M/P,N/Q
%%双线性内插本质上就是从x方向和y方向分别作了一次线性插值
src = imread('E:\pythonfile_withpycharm\cluster.jpg');
% 得到原图像的大小
[M,N,Z] = size(src);
% 设置目标图像的大小是300*300*3
dst = zeros(300,300,3);
[P,Q,H] = size(dst);
scaler_x = M/P;
scaler_y = N/Q;
% 插值实现
for i=1:P
    for j=1:Q
        % 后向映射得到在原图上的大概坐标
        % 注意一般是浮点类型
        x = i*scaler_x;
        y = j*scaler_y;
        % 取得与(x,y)距离最近的4个像素点
        X1 = floor(x);
        X2 = ceil(x);
        Y1 = floor(y);
        Y2 = ceil(y);
        P1_1 = src(X1,Y1,:);
        P1_2 = src(X1,Y2,:);
        P2_1 = src(X2,Y1,:);
        P2_2 = src(X2,Y2,:);
        % 先在x轴方向进行二次线性内插
        fx1 = (X2-x)/(X2-X1).*P1_1+(x-X1)/(X2-X1).*P2_1;
        fx2 = (X2-x)/(X2-X1).*P1_2+(x-X1)/(X2-X1).*P2_2;
        % 在y方向进行一次线性内插
        f = (Y2-y)/(Y2-Y1).*fx1+(y-Y1)/(Y2-Y1).*fx2;
        dst(i,j,:) = f;
    end
end
dst = im2uint8(mat2gray(dst));
figure();
imshow(src)
title('原图像512*512*3')
figure();
imshow(dst)
title('目标图像300*300*3') 

结果如下图所示:

可以看到缩放后的图像出现了黑色的网格线，可能是由于坐标系的设置问题。