xianxincz.m

双线性插值算法的实现

function  lkj
X=imread('lena.bmp');


 I=X(1:256,1:256);
 [n,m]=size(I);
 figure(1);

 imshow(I);
k=2;
B=zeros(k*n,k*m);
for i=1:n-1
    for j=1:m-1
        for t1=(k*(i-1)+1):(k*i+1)
            for t2=(k*(j-1)+1):(k*j+1)
                I1=I(i,j+1)-I(i,j);
                I2=I(i+1,j+1)-I(i+1,j);
                p1=I1*(t2-k*(j-1)-1)/k+I(i,j);
                p2=I2*(t2-k*(j-1)-1)/k+I(i+1,j);
                p3=(p2-p1)*(t1-k*(i-1)-1)/k+p1;
                B(t1,t2)=p3;
            end
        end
    end
end
B=uint8(B);
figure(2);
imshow(B);