太原理工大学-数字图像处理实验.docxVIP

50 50 本科实验报告 课程名称： 数字图像处理 实验项目： 实验一空域图像增强技术 实验地点： 起点机房 专业班级： 通信1002 学号： 学生： 指导教师： 2013年 6 月 1 日 实验一空域图像增强技术 一、 实验目的 1结合实例学习如何在视频显示程序中增加图像处理算法； 2理解和掌握图像的线性变换和直方图均衡化的原理和应用； 3了解平滑处理的算法和用途，学习使用均值滤波、中值滤波和拉普拉斯锐 化进行图像增强处理的程序设计方法； 4 了解噪声模型及对图像添加噪声的基本方法。 实验原理 实验原理 1灰度线性变换就是将图像中所有点的灰度按照线性灰度变换函数进行变 换。 g(x,y) T[f(x, y)] f (x, y) 0 f(x, y) a g(x,y) [f (x, y) a] ga a f(x, y) b [f (x, y) b] gb b f (x, y) 255 x 1,2, ,m, y 1,2, n 2直方图均衡化通过点运算将输入图像转换为在每一级上都有相等像素点 数的输出图像。按照图像概率密度函数 PDF的定义： Pr（rk）匹 k 0,1,2,…丄 1 n 通过转换公式获得： k k山 Sk T（rQ Pr（rJ - k 0,1,2,...,L 1 j o j 0 n 3均值（中值）滤波是指在图像上，对待处理的像素给定一个模板，该模板 包括了其周围的临近像素。将模板中的全体像素的均值（中值）来代替原来像素 值的方法。 4拉普拉斯算子如下： 1 1 1 8 1 1 1 1 拉普拉斯算子首先将自身与周围的8个像素相减，表示自身与周围像素的 差异，再将这个差异加上自身作为新像素的灰度。 三、实验步骤 1启动MATLAB程序，对图像文件分别进行灰度线性变换（参考教材 页，例 4.1)、直方图均衡化(参考教材 56 页，例 4.6)、均值滤波(参考教材 60 页，例 4.9)、中值滤波(参考教材 64 页，例 4.11)和梯度锐化操作(参考教材 66 页，例 4.12)。添加噪声，重复上述过程观察处理结果。 记录和整理实验报告 四、 实验仪器 1 计算机； 2 MATLAB 程序； 记录用的笔、纸。 五、 实验报告容 1、对图像文件进行灰度线性变换 程序： I=imread(RX01.jpg); subplot(1,2,1 )； imshow(I);title( ‘原图 ); I=double(I); [M,N]=size(I); for i=1:M for j=1:N if I(i,j)=30 I(i,j)=I(i,j); else if I(i,j)=150 I(i,j)=(200-30)/(150-30)*(I(i,j)-30)+30; else I(i,j)=(255-200)/(255-150)*(I(i,j)-150)+200; end end end end subplot(122);imshow(uint8(I)); 结果： 加椒盐噪声再处理的结果: 加噪声原图 2、直方图均衡化 I l=imread( RX02.jpg ); subplot(221);imshow(l); subplot(222);imhist(l) I1=histeq(l); subplot(223);imshow(l1); subplot(224);imhist(l1); 结果： 10000 0 100 2CE 加椒盐噪声再处理的结果: 0 100 200 3、均值滤波 l=imread( 45.gif ); [M,N]=size(l); x=1; II1=zeros(M,N); for i=1:16 ll(:,:,i)=im noise(l, gaussia n ,0,0.01); Il1=ll1+double(ll(:,:,i)); if or(or(i==1,i==4),or(i==8,i==16)); subplot(2,2,x);imshow(ui nt8(ll1/i));title( num2str(x)); x=x+1; end ; end 结果: 加椒盐噪声再处理的结果: 1 2 4、中值滤波 l=imread(34.gif); 原始图像);添加椒盐噪声图像) 原始图像); 添加椒盐噪声图像) k2=medfilt2(J,[5,5]); k3=medfilt2(J,[7,7]); k4=medfilt2(J,[9,9]); 模板中值滤波); 模板中值滤波 模板中值滤波); 模板中值滤波); 模板中值滤波); 模板中值滤波); .障始图像 添加槻盐噪声囹僦 氛模板中值滤波 対模極中值滤滅 7x7^?中值滤液 旳g模板中值滤液 5、梯度锐化操作 l=imread(20.gif); subplot(