matlab图像几何变换

1、最新资料推 姑斤希 段工学卩免学 光电图像处理实验报告 实验一：图像几何变换 作 者肖剑洪 专 业电弄I学皱术 学 院物理与电子学院指导老师 王晓明 完成时间 2013.12.2 实验一图像几何变换 一、实验目的 图像儿何变换是图像处理与分析的基础内容之一。使学生掌握图像平 移、 镜像、缩放和旋转等儿何变换的方法以及程序实现。 二、实验设备 联想图像处理工作站 三、实验内容及要求 图像儿何变换不改变图像的像素值，而是改变像素所在的几何位置，具有独 特的方法。本实验的内容如下： 1. 学会使用Mat lab的图像处理工具箱仃mage Processing Toolbox

2、) o使学 生初步具备使用该软件处理图像信息能力，并可以利用该软件完成本课程规定 的其他实验和作业。 式取整。 图像的比例缩放所产生的图像中的像素可能在源图像中找不到相应的像素 点，这时就需要进行插值处理，通常可用最近邻插值和双线性插值算法来实现。 (1) 最近邻插值 对通过反变换得到的浮点坐标进行简单的取整操作，得到一个整数型坐标， 这个整数型坐标对应的原图像的像素值就是新图像的像素值。 (2) 双线性插值 对通过反变换得到的浮点坐标，求出此浮点坐标对应原图像四个邻域整数型 坐标的灰度值，然后经过双线性插值得到新图像的像素值。 (3) 双三次插值 对通过反变换得到的浮点坐标，求出此浮点坐标对

3、应原图像16个邻域整数 型坐标的灰度值，然后经过双三次插值得到新图像的像素值。 4.图像旋转 一般图像的旋转是以图像的中心为原点，将图像上的所有像素都旋转 个相同的角度。如图所示： 新图像坐标与原图像坐标对应关系: 两个常量: 图像旋转之后会出现许多空洞点，须对这些空洞点进行插值处理。可用近邻插 值，均值插值算法。 五、实验步骤 1. 根据实验内容的要求在Matlab软件中编写相应程序； 2. 调试运行程序，并记录结果： 3. 完成实验报告。 六、实验程序清单和实验结果及分析 1.图像平移 I 二 imread(f:/l .jpg)； se = translate(strel(l)z 25 2

4、5); J = imdila te(l,se); imshow(l), title(Originar) figure, imshow(J)z title(Translated); UriqinalTranslated 2.图像镜像 A=imread(f:/1 .jpg)； height,width,dim二 size( A); tformmaketformfaffineJ-l 0 0;0 1 0;width 0 1); B=imtransform (Aztform/neoest); tform2=maketform(affine,l 0 0;0 -1 0;0 heigh11); C=imtr

5、on sform (Aztform2/neorest); subplotfl ,3,1 )Jmshow(A); title(原图像)； subplot(l ,3,2),imshow(B); title fTk 平镜像)； subplot(l ,3,3),imshow(C); title (*!g 直镜像)； 原團像 3.图像比例缩放 水平镜像 竖直镜像 T1 二 maketform(affine,2 0 0; 0 1 0; 0 0 1); T2 二 maketformfaffineJl 0 0; 0 2 0; 0 0 1); T3 = maketform(affine,4/5 0 0; 0 3

6、/5 0; 0 0 1); I 二 imread(f:/l .jpg); II = imtr on sform(IJl); 12 = imtr on sform(IJ2); 13 = imtransform(IJ3); imshowfl), figure, imshow(ll)z figure, imshow(l2), figure, imshow(l3); 4. 图像旋转 lmage=imread(,f:/1 .jpg1)； angle=30; pai=3.14; iml=(lmage);subplot(211); imshow(iml); titlefr); 最新资料推 示心3_K3_ r

7、oteang=angle/180*pai; imlXziml Y=size(iml); rotmatrix=cos(rote on g)z-sin(roteang);si n(oteong),cos(oteo ng)； rrotmatrix= cos(-roteang)z-sin(-roteang):sin (-rotea ng) ,cos (-ot eang); im2XY 1 = im 1 Xzim 1Y *rotmatrix; im2XY2= imlXzl*rotmatrix; im2XY3=l z 1 *rotmatrix; im2XY4= Liml Y*rotmatrix; im2p

8、l l=min(im2XYl (1 )zim2XY2(l )Jm2XY3(l)jm2XY4(l )zmin (im2XYl(2)zim2XY2(2)zim2XY3(2)zim2XY4(2)J; im2X=round(max(im2XYl (1 ),im2XY2(l) Jm2XY3(l) Jm2XY4(l) -im2pll(l); im2Y=round(max(im2XYl(2)jm2XY2(2)/im2XY3(2)zim2XY4(2) -im2pll(2); im2=uint8(zeros(im2Xjm2Y)+255; for i=l:im2X forj=l:im2 Y pij=round(i+im2pl 1 (l)-l/j+im2pl 1 (2)-l*rrotmatrix); if (pij(1)二 1 ) end end end; subplot(212); imshow(im2);