图像的放大与缩小(MATLAB 代码)

版权声明:本文为博主原创文章,未经博主允许不得转载。 https://blog.csdn.net/u013146742/article/details/52997878

首先第一种:

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)

 

然后再输入I=imread('lena.jpg');

       B=ImageZoom(I,0.5,0.5);

结果显示:

第二种 双线性插值

close all;
clear all;
clc;

m=1;              %放大或缩小的高度
n=2;              %放大或缩小的宽度
img=imread('lena.jpg');
img=rgb2gray(img);
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(2),imshow(uint8(imgn)) ;

结果显示

注意 不过因为zeros 只能输入整数小数无法显示出来,有一定的缺陷

 

阅读更多
想对作者说点什么?

博主推荐

换一批

没有更多推荐了,返回首页