数字图像处理.doc

数 字 图 像 处 理 专业班级：测绘1103班 姓 名： 学 号： 指导老师： 梅 小 明 实验一、空域图像处理 实验目的 1熟悉CCS集成开发环境的操作和基本功能； 2熟悉MATLAB基本图像操作； 3结合实例学习如何在程序中增加图像处理算法； 4理解和掌握图像的线性变换和直方图均衡化的原理和应用； 5了解平滑处理的算法和用途，学习使用均值滤波、中值滤波和拉普拉斯锐化进行图像增强处理的程序设计方法； 6 了解噪声模型及对图像添加噪声的基本方法。 实验原理 1 灰度线性变换就是将图像中所有点的灰度按照线性灰度变换函数进行变换。 2 直方图均衡化通过点运算将输入图像转换为在每一级上都有相等像素点数的输出图像。按照图像概率密度函数PDF的定义： 通过转换公式获得： 3 均值（中值）滤波是指在图像上，对待处理的像素给定一个模板，该模板包括了其周围的临近像素。将模板中的全体像素的均值（中值）来代替原来像素值的方法。 4 拉普拉斯算子如下： 拉普拉斯算子首先将自身与周围的8个像素相减，表示自身与周围像素的差异，再将这个差异加上自身作为新像素的灰度。 实验步骤 1 启动MATLAB程序，对图像文件分别进行灰度线性变换、直方图均衡化、均值滤波、中值滤波和拉普拉斯锐化操作；添加噪声，重复上述过程观察处理结果。 2 记录和整理实验报告 实验仪器 1计算机； 2 MATLAB程序； 3记录用的笔、纸。 五 实验用到的代码及结果 I=imread(数字图像.bmp); B=imnoise(I,salt pepper); subplot(3,2,4);imshow(B); H=fspecial(sobel); C=imfilter(I,H); subplot(3,2,5);imshow(C); subplot(3,2,5);imshow(H)； D=medfilt2(B); subplot(3,2,6);imshow(D); %显示中值滤波效果 实验二、图像的傅立叶变换 实验目的 1了解图像变换的意义和手段； 2熟悉傅里叶变换的基本性质； 3熟练掌握FFT方法及应用； 4通过实验了解二维频谱的分布特点； 5通过本实验掌握利用MATLAB编程实现数字图像的傅立叶变换及滤波锐化和复原处理； 6 了解理想、巴特沃兹、高斯等不同滤波器的结构及滤波效果。 实验原理 1应用傅立叶变换进行图像处理 傅里叶变换是线性系统分析的一个有力工具，它能够定量地分析诸如数字化系统、采样点、电子放大器、卷积滤波器、噪音和显示点等的作用。通过实验培养这项技能，将有助于解决大多数图像处理问题。对任何想在工作中有效应用数字图像处理技术的人来说，把时间用在学习和掌握博里叶变换上是很有必要的。 2傅立叶（Fourier）变换的定义 对于二维信号，二维Fourier变换定义为：　　 　　二维离散傅立叶变换为： 图像的傅立叶变换与一维信号的傅立叶变换变换一样，有快速算法，具体参见参考书目，有关傅立叶变换的快速算法的程序不难找到。实际上，现在有实现傅立叶变换的芯片，可以实时实现傅立叶变换。 3利用MATLAB软件实现数字图像傅立叶变换的程序： I=imread(‘数字图像.gif’); %读入原图像文件 imshow(I); %显示原图像 fftI=fft2(I); %二维离散傅立叶变换 sfftI=fftshift(fftI); %直流分量移到频谱中心 RR=real(sfftI); %取傅立叶变换的实部 II=imag(sfftI); %取傅立叶变换的虚部 A=sqrt(RR.^2+II.^2);%计算频谱幅值 A=（A-min(min(A))）/(max(max(A))-min(min(A)))*225; %归一化 figure; %设定窗口 imshow(A); %显示原图像的频谱 实验步骤 1打开计算机，安装和启动MATLAB程序；程序组中“work”文件夹中应有待处理的图像文件； 2利用MatLab工具箱中的函数编制FFT频谱显示的函数; 3 a).调入、显示图像；图像存储格式应为“.gif”; b)对图像做FFT并利用自编的函数显示其频谱; c)讨论不同的图像内容与FFT频谱之间的对应关系。 4 对频谱分别进行高斯、巴特沃兹、高斯高通和低通滤波，比较其锐化和平滑效果 5记录和整理实验报告。 实验仪器 1计算机； 2 MATLAB程序； 3移动式存储器（软盘、U盘等