图像缩放程序

一个是双线性插值的,SX,SY是缩放倍数:

function transImage=imscale(orgImage,sx,sy)
[orgH,orgW,layer]=size(orgImage);
transH=round(orgH*sy);
transW=round(orgW*sx);
transImage=zeros(transH,transW);
for m=1:transH
    for n=1:transW
        x=(n-1)/sx+1;   y=(m-1)/sy+1;
        xi=fix(x); yi=fix(y);        xm=x-xi;   ym=y-yi;
        
        if x<orgW & y<orgH
            tmp=orgImage(yi,xi)*(1-xm)*(1-ym)+orgImage(yi+1,xi)*(1-xm)*ym+orgImage(yi,xi+1)*xm*(1-ym)+orgImage(yi+1,xi+1)*xm*ym;
        elseif x==orgW & y<orgH
            tmp=orgImage(yi,xi)*(1-xm)*(1-ym)+orgImage(yi+1,xi)*(1-xm)*ym;
        elseif x<orgW & y==orgH
            tmp=orgImage(yi,xi)*(1-xm)*(1-ym)+orgImage(yi,xi+1)*xm*(1-ym);
        elseif x==orgW & y==orgH
            tmp=orgImage(yi,xi)*(1-xm)*(1-ym);
        end;            
        transImage(m,n)=round(tmp);
    end;
end;

另外一个:

function I=ImageZoom(I1,x,y)
%功能：实现图像的任意倍数缩放
%x---水平缩放系数
%y---竖直缩放系数

if length(size(I1))>2
I1=rgb2gray(I1);
end

figure,imshow(I1);
[m,n]=size(I1);

newWidth=round(x*m);
newHeight=round(x*n);

T=[x 0 0;0 y 0;0 0 1];
tform=maketform('affine',T);

tx=zeros(newWidth,newHeight);
ty=zeros(newWidth,newHeight);
for i=1:newWidth
for j=1:newHeight
tx(i,j)=i;
ty(i,j)=j;
end
end

[w z]=tforminv(tform,tx,ty); %反向坐标值

I=uint8(zeros(newWidth,newHeight));

%给新图像各像素点赋值
for i=1:newWidth
for j=1:newHeight
S_x=w(i,j);
S_y=z(i,j);
if(S_x>=m-1||S_y>=n-1||double(uint16(S_x))<=0||double(uint16(S_y))<=0) 
I(i,j)=0; %不在原图像上
else
if (S_x/double(uint16(S_x))==1.0&S_y/double(uint16(S_y))==1.0)
I(i,j)=I1(uint16(S_x),uint16(S_y));%整数点
else
%不是整数点
a=double(uint16(S_x));
b=double(uint16(S_y));
u=S_x-a;
v=S_y-b;
x11=double(I1(a,b));
x12=double(I1(a,b+1));
x21=double(I1(a+1,b));
x22=double(I1(a+1,b+1));
I(i,j)=uint8((1-u)*(1-v)*x11+(1-u)*v*x12+u*(1-v)*x21+u*v*x22);
end
end
end
end

figure,imshow(I)