[工学]数字图像处理.doc

（程序设计类课程） 实验报告 课程名称： 数字图像处理 姓 名： 系： 信息与机电 专 业： 电子信息工程 年 级： 学 号： 指导教师： 职 称： 2011年 月 日 实验项目列表 序号 实验项目名称 成绩 指导教师 1 数字图像的放大与缩小 2 直方图均衡化 3 空间域数字图像的平滑与锐化 4 二维离散傅立叶变换的计算与显示 5 彩色图像处理 6 7 8 9 10 11 12 实验报告 系： 专业： 年级： 姓名： 学号： 实验室号_ _ 计算机号 实验时间： 指导教师签字： 成绩： 报告退发 (订正 、 重做) 实验一 数字图像的放大与缩小 实验目的和要求 掌握数字图像放大缩小算法。 实验内容和原理 利用最近邻域内插算法和双线性内插算法实现数字图像的放大和缩小。 双线性内插计算公式为：f(x,y)=ax’+by’+cx’y’+d; a=f(i+1,j)-f(i,j) ;b= f(i,j+1)-f(i,j) ; d=f(i,j); c= f(i+1,j+1)+f(i,j)-f(i+1,j)-f(i,j+1) 其中(i+1,j+1)、(i,j)、(i+1,j)、(i,j+1)是(i+x’,j+y’)的4个最近邻像素；(i+x’,j+y’)为(x,y)变换之前的坐标A=imread(flower.bmp); n=6; [t1,t2]=size(A); imshow(A); x=t1*n; y=t2*n; for i=1:x for j=1:y x1=round(i/n); y1=round(j/n); if(x11) x1=1; end if(y11) y1=1; end B(i,j)=A(x1,y1); end end figure; imshow(B); 双线性内插： A=imread(flower.bmp); n=3; [t1,t2]=size(A); imshow(A); x=t1*n; y=t2*n; for i=1:x for j=1:y x1=floor(i/n); y1=floor(j/n); if(x11)x1=1; end if(y11)y1=1; end if(x1=t1)x1=t1-1; end if(y1=t2)y1=t2-1; end a=A(x1+1,y1)-A(x1,y1); b=A(x1,y1+1)-A(x1,y1); c=A(x1+1,y1+1)+A(x1,y1)-A(x1+1,y1)-A(x1,y1+1); d=A(x1,y1); xx=i/n-x1; yy=j/n-y1; B(i,j)=a*xx+b*yy+c*xx*yy+d; end end figure; imshow(B); 调试过程 最近邻域内插： 刚开始实验时，因为没注意到文件的路径。出现了错误： 于是更改current directory。改为matlab安装路径下的work工作文件夹下。 更改过后再次运行。发现还是一样的错误。 经过仔细研究程序，发现是程序中的错误。应把B(i,j)=A(x,y)改为B(i,j)=A(x1,y1) 结果发现两张图重叠了。原因是没有进行清屏。修改过后。运行正确。 双线性内插： 由于有了前一次最近邻与内插的经验总结。在双线性内插中就避免了同样的错误。 6.实验结果 最近邻域内插： 双线性内插： 6.总结 通过比较最近邻域内插和双线性内插两种方法，发现放大一定倍数后会有块效应。但是双线性内插没有最近邻域内插的效果那么明显。 使用matlab软件时要注意一定要将.m文件和图片放在matlab安装路径下的work文件夹中。在给文件命名时也要注意文件一定要以字母为开头，可以有数字和下划线。 对于使用两种方法达到图像的放大和缩小，首先要理解程序中