matlab图像几何变换

1、-. z.学 号光电图像处理实验报告 实验一：图像几何变换作 者肖剑洪专 业电子科学与技术学 院物理与电子学院指导教师王晓明完成时间实验一 图像几何变换一、实验目的图像几何变换是图像处理与分析的根底容之一。使学生掌握图像平移、镜像、缩放和旋转等几何变换的方法以及程序实现。二、实验设备联想图像处理工作站三、实验容及要求图像几何变换不改变图像的像素值，而是改变像素所在的几何位置，具有独特的方法。本实验的容如下：1学会使用Matlab的图像处理工具箱(Image Processing Toolbo*)。使学生初步具备使用该软件处理图像信息能力，并可以利用该软件完本钱课程规定的其他实验和作业。2编程实

2、现图像平移，要求平移后的图像大小不变； 3编程实现图像的镜像； 4编程实现图像的比例缩放，要求分别用最近邻插值、双线性插值和双三次插值三种方法来实现，并比拟缩放效果；5编程实现以任意角度对图像进展旋转变换。四、实验原理1图像平移从上到下、从左到右依次计算新图像中的每个像素*, y在源图像中对应的坐标值*0, y0； 即： *0 = * *，y0 = y y其中，*、y分别是水平和垂直偏移量。2图像镜像 从上到下、从左到右依次计算新图像中的每个像素*, y在源图像中对应的坐标值*0, y0； 即： *0 = fWidth *，y0 = y水平镜像 *0 = *，y0 = fHeight y垂直镜

3、像其中，fWidth、fHeight分别是图像的宽度和高度。3图像比例缩放图像比例缩放是指将给定的图像在*轴方向按比例缩放倍，在y轴方向按比例缩放倍，从而获得一幅新的图像。从上到下、从左到右依次计算新图像中的每个像素*, y在源图像中对应的坐标值*0, y0；即：*0 = * / f* + 0.5，y0 = y / fy + 0.5其中，f*、fy分别是水平和垂直方向的缩放比率，中括号表示对括号表达式取整。图像的比例缩放所产生的图像中的像素可能在源图像中找不到相应的像素点，这时就需要进展插值处理，通常可用最近邻插值和双线性插值算法来实现。(1) 最近邻插值对通过反变换得到的浮点坐标进展简单的取

4、整操作，得到一个整数型坐标，这个整数型坐标对应的原图像的像素值就是新图像的像素值。2双线性插值对通过反变换得到的浮点坐标，求出此浮点坐标对应原图像四个邻域整数型坐标的灰度值，然后经过双线性插值得到新图像的像素值。双三次插值对通过反变换得到的浮点坐标，求出此浮点坐标对应原图像16 个邻域整数型坐标的灰度值，然后经过双三次插值得到新图像的像素值。4图像旋转 一般图像的旋转是以图像的中心为原点，将图像上的所有像素都旋转一个一样的角度。如下图：新图像坐标与原图像坐标对应关系：两个常量：图像旋转之后会出现许多空洞点，须对这些空洞点进展插值处理。可用近邻插值，均值插值算法。五、实验步骤1根据实验容的要求在

5、Matlab软件中编写相应程序；2调试运行程序，并记录结果；3完成实验报告。六、实验程序清单和实验结果及分析1图像平移I = imread(f:/1.jpg);se = translate(strel(1), 25 25);J = imdilate(I,se);imshow(I), title(Original)figure, imshow(J), title(Translated);图像镜像A=imread(f:/1.jpg);height,width,dim=size(A);tform=maketform(affine,-1 0 0;0 1 0;width 0 1);B=imtransfo

6、rm(A,tform,nearest);tform2=maketform(affine,1 0 0;0 -1 0;0 height 1);C=imtransform(A,tform2,nearest);subplot(1,3,1),imshow(A);title(原图像);subplot(1,3,2),imshow(B);title(水平镜像);subplot(1,3,3),imshow(C);title(竖直镜像);图像比例缩放T1 = maketform(affine,2 0 0; 0 1 0; 0 0 1);T2 = maketform(affine,1 0 0; 0 2 0; 0 0

7、1);T3 = maketform(affine,4/5 0 0; 0 3/5 0; 0 0 1); I = imread(f:/1.jpg);I1 = imtransform(I,T1);I2 = imtransform(I,T2);I3 = imtransform(I,T3);imshow(I), figure, imshow(I1), figure, imshow(I2), figure, imshow(I3);图像旋转 Image=imread(f:/1.jpg);angle=30;pai=3.14;im1=(Image);subplot(211);imshow(im1);title(

8、1);roteang=angle/180*pai;im1*,im1Y=size(im1);rotmatri*=cos(roteang),-sin(roteang);sin(roteang),cos(roteang);rrotmatri*=cos(-roteang),-sin(-roteang);sin(-roteang),cos(-roteang);im2*Y1=im1*,im1Y*rotmatri*;im2*Y2=im1*,1*rotmatri*;im2*Y3=1,1*rotmatri*;im2*Y4=1,im1Y*rotmatri*;im2p11=min(im2*Y1(1),im2*Y2(

9、1),im2*Y3(1),im2*Y4(1),min(im2*Y1(2),im2*Y2(2),im2*Y3(2),im2*Y4(2);im2*=round(ma*(im2*Y1(1),im2*Y2(1),im2*Y3(1),im2*Y4(1)-im2p11(1);im2Y=round(ma*(im2*Y1(2),im2*Y2(2),im2*Y3(2),im2*Y4(2)-im2p11(2);im2=uint8(zeros(im2*,im2Y)+255;for i=1:im2* for j=1:im2Y pij=round(i+im2p11(1)-1,j+im2p11(2)-1*rrotmatri*); if (pij(1)=1)&(pij(2)=1)&(pij(1)=im1*)&(pij(2)=im1Y) im2(i,j)=im1(pij(1),pij(2); end endend;subp