matlab練習程序（圖像放大/縮小，雙線性插值）

　　半年前寫過matlab最鄰近插值的圖像縮放，沒怎么考慮邊界問題。更早之前用Opencv寫過雙線性插值圖像放大，不過寫的比較混亂。所以這里用matlab重新再清楚的寫一遍。

close all;
clear all;
clc;

m=1.8;              %放大或縮小的高度
n=2.3;              %放大或縮小的寬度
img=imread('lena.jpg');
imshow(img);
[h w]=size(img);
imgn=zeros(h*m,w*n);
rot=[m 0 0;0 n 0;0 0 1];                                   %變換矩陣

for i=1:h*m
    for j=1:w*n
        pix=[i j 1]/rot;   
        
        float_Y=pix(1)-floor(pix(1)); 
        float_X=pix(2)-floor(pix(2));
       
        if pix(1) < 1　　　　　　　　%邊界處理
            pix(1) = 1;
        end
        
        if pix(1) > h
            pix(1) = h;
        end
        
        if pix(2) < 1
            pix(2) =1;
        end
        
        if pix(2) > w
            pix(2) =w;
        end
        
        pix_up_left=[floor(pix(1)) floor(pix(2))];　　　　%四個相鄰的點
        pix_up_right=[floor(pix(1)) ceil(pix(2))];
        pix_down_left=[ceil(pix(1)) floor(pix(2))];
        pix_down_right=[ceil(pix(1)) ceil(pix(2))];     
    
        value_up_left=(1-float_X)*(1-float_Y);　　　　　　%計算臨近四個點的權重
        value_up_right=float_X*(1-float_Y);
        value_down_left=(1-float_X)*float_Y;
        value_down_right=float_X*float_Y;
　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　%按權重進行雙線性插值
        imgn(i,j)=value_up_left*img(pix_up_left(1),pix_up_left(2))+ ...
                  value_up_right*img(pix_up_right(1),pix_up_right(2))+ ...
                  value_down_left*img(pix_down_left(1),pix_down_left(2))+ ...
                  value_down_right*img(pix_down_right(1),pix_down_right(2));        
    end
end

figure,imshow(uint8(imgn))

原圖

 放大后的。