雙線性插值matlab實(shí)現(xiàn)

1、 張俊飛一、算法簡(jiǎn)介假設(shè)ca到ca+1之間是線性變化的，那么對(duì)于浮點(diǎn)數(shù)x( a <= x < a+1)c(x) = ca+1*( x - a) + ca*( 1 + a - x);c(x) = ca + (ca+1-ca)/(b-a)*( x - a);b = a+1;c(x) = ca + ca+1 - ca*( x - a);把這種插值方式擴(kuò)展到二維情況：對(duì)于一個(gè)二維數(shù)組c，我們假設(shè)對(duì)于任意一個(gè)浮點(diǎn)數(shù)i,c(a,i)到c(a+1,i)之間是線性變化的,c(i,b)到c(i,b+1)之間也是線性變化的(a,b都是整數(shù)) ，那么對(duì)于浮點(diǎn)數(shù)的坐標(biāo)(x,y)滿足(a <= x &

2、lt; a+1, b <= y < b+1)，可以先分別求出c(x,b)和c(x,b+1):c(x,b) = ca+1b*( x - a) + cab*( 1 + a - x);c(x,b+1) = ca+1b+1*( x - a) + cab+1*( 1 + a - x);現(xiàn)在已經(jīng)知道c(x,b)和c(x,b+1)了，而根據(jù)假設(shè)c(x,b)到c(x,b+1)也是線性變化的，所以:c(x,y) = c(x,b+1)*( y - b) + c(x,b)*( 1 + b - y)二、matlab實(shí)現(xiàn)在command窗口輸入chahzi('cameraman.tif',2

3、)，這里cameraman.tif為灰度圖片，得到結(jié)果如下：在command窗口輸入chahzi('a.jpg',)，這里a.jpg為彩色圖片，得到結(jié)果為：由以上結(jié)果可知，該算法的魯棒性較強(qiáng)，圖片伸縮插值保真度高。三、代碼function =chahzi(str,k)% str input image% k zoo in or out k times% clear all% clc% k=2;I=imread(str);% if size(I,3)=1% I=rgb2gray(I);% endfigure,imshow(I);title('origin image

4、9;); I=double(I); rows cols tongdao=size(I); x_new=rows*k; y_new=cols*k; %縮放至k倍。 I_new=zeros(x_new,y_new,tongdao);for rgb=1:tongdaofor i=1:x_newfor j=1:y_new x=i/k;a=floor(x); y=j/k;b=floor(y); % 雙線性插值算法if a>0&&b>0&&a<rows&&b<cols cxb=I(a+1,b,rgb)*(x-a)+I(a,b,rgb)*(1+a-x); cxb1=I(a+1,b+1,rgb)*(x-a)+I(a,b+1,rgb)*(1+a-x); I_new(i,j,rgb)=round(cxb1*(y-b)+cxb*(1+b-y);endendendend% iftongdao=1%如果輸入是灰度圖，則% I_new=I_n