苏州大学数字图像Matlab上机实验.docVIP

《数字图像处理与分析》 实验报告 专业 计算机科学与技术 年级 10级 姓名 学号 任课老师 龚声蓉 指导老师 蒋德茂 实验室 理工楼243机房 学期 2013（春） 苏州大学计算机科学与技术学院统一印制 二零一三年六月 实验一：空域及频域增强 实验目的 了解图像变换的意义和手段； ?熟悉傅里叶变换的基本性质； 通过本实验掌握利用MATLAB编程实现数字图像的傅立叶变换及滤波锐化和复原处理； 了解平滑处理的算法和用途，学习使用均值滤波、中值滤波和拉普拉斯锐化进行图像增强处理的程序设计方法。 实验内容 Matlab编程：实现Butterworth低通滤波和Butterworth高通滤波。 提示：预设Butterworth的阶等于2，d0=10。使用cameraman.tif作为原图像。 要求：使用subplot函数按3行2列分别显示（1）原始图（2）傅里叶频谱（3）Butterworth低通滤波（4）低通滤波图（5）Butterworth高通滤波（6）高通滤波图 实验分析 Roberts算子用来实现消除图像模糊地增强的方法,即“锐化”。此处理加强图像的边界和细节。Roberts算子提出的是在2*2的邻域上计算对角导数产生的偏导数在图像处理中心点上用均值或是绝对值求和的形式结合起来。但是此算法增强噪声的缺陷也在图像中体现了。 巴特沃思滤波器去的频率和通过的频率之间没有明显的不连续性，图像的模糊程度降低，而且也没有振铃效应，这是由于在低频和高频之间，滤波器平滑过渡的缘故。 程序及实验结果 程序代码： clc,clear all; I=imread(cameraman.tif); I=im2double(I); subplot(3,2,1);imshow(I,[]);%在一张图里依次显示,三行两列 title(原始图); F=fft2(I); G=fftshift(F); subplot(3,2,2);imshow(log(1+abs(G)),[]); title(傅里叶频谱); [m,n]=size(G); c1=floor(m/2); %中心行坐标 c2=floor(n/2); %中心列坐标 d0=10; for i=1:m for j=1:n d=sqrt((i-c1)^2+(j-c2)^2); h=1/(1+(d/d0)^4); A(i,j)=h*G(i,j); end end subplot(3,2,3);imshow(log(1+abs(A)),[]); title(Butterworth低通滤波); d0=10; for i=1:m for j=1:n d=sqrt((i-c1)^2+(j-c2)^2); h=1/(1+(d/d0)^4); A(i,j)=h*G(i,j); end end f=ifftshift(A); i=real(ifft2(f)); subplot(3,2,4);imshow(i,[]); title(低通滤波图); d0=10; for i=1:m for j=1:n d=sqrt((i-c1)^2+(j-c2)^2); h=1/(1+(d0/d)^4); A(i,j)=h*G(i,j); end end subplot(3,2,5);imshow(log(1+abs(A)),[]); title(Butterworth高通滤波); d0=40; for i=1:m for j=1:n d=sqrt((i-c1)^2+(j-c2)^2); h=1/(1+(d0/d)^4); B(i,j)=h*G(i,j); end end f=ifftshift(B); i=real(ifft2(f)); subplot(3,2,6);imshow(i,[]); title(高通滤波图); 实验结果： 实验二：图像编码与压缩 一．实验目的 了解图像编码与压缩过程； 实践算数编码及解码的Matlab编程实现； 二．实验内容 Matlab编程：算术编码及解码。 提示：已知字符事件出现的概率如下表 表：数据出现的概率和取