图像处理的几何变换

一.  设计一个程序，对一幅灰度图像的实现如下几何变换,1)放大为原图1.5倍;2)绕中心旋转30度(CW);3)采用偏移量插值实现一个透视变换.灰度差值用最近邻插值和双线性插值。

1.放大为原图的1.5倍：

 (1) 最近邻插值法：

原理：输出图像像素的灰度值等于离它所映射到的位置最近的输入像素的灰度值。

具体实现代码：

编写图像缩放的函数myresize.m :

function myresize( I , d )
[M,N]=size(I);
I1=uint8(zeros(floor(M*d),floor(N*d)));
for i=2:floor(M*d)
    for j=2:floor(N*d)
        x=floor(i/d);   y=floor(j/d);
        if((x<M)&&(y<N))
            I1(i,j)=I(x,y);
        end
    end
end
figure;imshow(I);
figure;imshow(I1);
end

Command窗口中输入：

I=imread('lena3.jpg')
myresize(I,1.5)

处理效果:

                         处理前256×256

                   处理后384×384

(2)双线性插值法:

  原理：对于一个目的像素，设置坐标通过反向变换得到的浮点坐标为 (i+u,j+v),其中i，j为非负整数，u，v为[0,1]区间的浮点数,则这个像素的值 f(i+u,j+v)可由原图像中的坐标为(i,j),(i+1,j),(i,j+1),(i+1,j+1)所对应的周围四个像素的值决定。

 具体实现：

clear,clc
I=imread('lena3.jpg');
[M,N]=size(I);
d=1.5;
I1=uint8(zeros(floor(M*d),floor(N*d)));
M1=floor(M*d);  N1=floor(N*d);
for i=2:M1-1
    for j=2:N1-1