图像处理课程实践论文 题目： 图像增强功能组件的GUI设计
院（系） 信息科学与工程学院 专 业 通信工程 届 别 班 级 学 号 姓 名 指导老师 目录 一、空域增强……………………………………………………………1 1.1 直方图均衡化…………………………………………………1 1.2 平滑滤波………………………………………………………2 1.3 锐化滤波………………………………………………………4 1.4 自适应滤波……………………………………………………5 二、频域增强……………………………………………………………6 2.1 低通滤波………………………………………………………6 2.2 高通滤波………………………………………………………7 三、其他增强方法………………………………………………………9 3.1 灰度变换………………………………………………………9 3.2 椒盐噪声加入…………………………………………………9 3.3 图像镜像………………………………………………………10 3.4 图像求反………………………………………………………10 3.5 边缘检测………………………………………………………11 3.6 图像旋转………………………………………………………11 3.7 打开 保存 撤销 退出 初始化………………………………12 四、GUI界面的使用说明………………………………………………14 4.1 简介……………………………………………………………14 4.2具体操作………………………………………………………14 五、功能组件测试报告………………………………………………15 六、参考文献……………………………………………………………16 七、附录…………………………………………………………………17 1 一、空域增强
1.1 直方图均衡化 1.1.1 知识点 直方图均衡是指将一个已知灰度分布的图像经过一种变换，使之变成一幅具有均匀灰度分布的新图像。 把原始图的直方图变换为均匀分布的形式，这样就增加了象素灰度值的动态范围，从而达到增强图像整体对比度的效果 直方图均衡化的增强函数需要满足两个条件： 1) EH(s) 在 0 ≤s≤ L-1 范围内是一个单值单增函数； 2) 对应 0 ≤ s ≤L-1 有 0 ≤ EH(s) ≤L-1。 上面第一个条件保证原图各灰度级在变换后仍保持从黑到白 (或从白到黑) 的排列次序。第二个条件保证变换前后灰度值动态范围的一致性。
1.1.2 程序设计思路 通过histeq（X）函数实现直方图均衡。因为此函数只能对灰度图像进行直方图均衡。故应先将彩图转为灰度图像。 直方图均衡实现程序段如下： global T axes(handles.axes2); T=getimage; k=histeq(handles.k); imshow(k); title('经过直方图均衡后的图像'); handles.k=k; guidata(hObject,handles);
1.1.3 结果与分析说明 1、结果 ①将彩图转为灰度图像：
图1 2
②在上一步的基础上对第二幅图进行直方图均衡： 图2 2、分析和说明 图1和图2的右侧图分别是左侧图的灰度变换图和经过直方图均衡化后的图像。由于直方图均衡化增加了图像灰度动态范围，所以也增加了图像的对比度，反映在图像上就是图像有较大的反差，许多细节可以看到比较清晰了。
1.2 平滑滤波 1.2.1 知识点 最简单的平滑滤波是用一个像素领域的平均值作为滤波结果，此时滤波模板的所有系数都取为1。为了保证输出图仍在原来的灰度值范围，在算得卷积值R后要将其除以系数总个数在进行赋值。例如对3X3的模板来说，在算得R后要将其除以系数9。领域平均的一般表达式为
),(),(2),(1),(yxNtstsfnyxg
其中，N(x,y)对应f(x,y)中(x,y)的n×n领域，与模板W所覆盖的范围对应。 1.2.2 程序设计思路 conv函数可计算起点为k=0的两个序列的卷积。利用此函数和卷积的性质，编写平滑滤波器的程序对图像进行领域平均平滑滤波。 平滑滤波实现的程序段如下： global T axes(handles.axes2); T=getimage; h=[1 1 1;1 1 1;1 1 1]; H=h/9; i=double(handles.k); k=convn(i,h); imshow(k,[]); 3
title('经过线性平滑后的图像'); handles.k=k; guidata(hObject,handles);
1.2.3 结果与分析说明 1、结果 ①将彩图转为灰度图像：
图3 ②在上一步的基础上对第二幅图叠加参数为0.05的椒盐噪声：
图4 ③对图4中的灰度图像进行多次的平滑滤波
图5 2、分析说明 图5依次为用3x3,5×5,7×7,9×9平滑模板对图4中的灰度图像进行平滑滤波的结果。有图可见，当所用平滑模板尺寸增大时，对噪声的消除效果有所增强。不过同时所得到的图像变得更为模糊，可视的细节逐渐减小。 4
1.3 锐化滤波 1.3.1 知识点 图像锐化的目的是加强图像中景物的边缘和轮廓。因为边缘和轮廓都位于灰度突变的地方，所以锐化算法的实现基于微分运算。 拉普拉斯算子是各向同性的阶微分算子。数字图像f(x,y)的二阶偏导数可表示为：
22222
yfxff
)1,()1,(),(2),1(),1(),(22222yxfyxfyxfy
f
yxfyxfyxfx
f
拉普拉斯算子可增强图像中的灰度不连续边缘，减弱灰度值缓慢变化区域的对比度，将结果叠加到原始图像上，就可以得到边缘增强后的图像。
1.3.2 程序设计思路 fspecial函数用于建立预定义的滤波算子，k = fspecial(type，para)其中type指定算子的类型，para指定相应的参数。故利用该函数进行锐化滤波的程序编写。 锐化滤波实现的程序段如下： global T axes(handles.axes2); T=getimage; K=fspecial('laplacian',0.7); k=filter2(K,T)/100; imshow(k); title('经过线性锐化后的图像'); handles.k=k; guidata(hObject,handles);
1.3.3 结果与分析说明 1、结果 对灰度图像经过锐化滤波后：
2、分析说明 通过对左右图的比较，可以看出经过锐化之后的图像，其中景物的边缘和轮廓更加的清晰。故锐化滤波器可以用于景物边缘的检测和提取。 5
1.4 自适应滤波 1.4.1 知识点 自适应滤波法是能够在一定程度上根据环境的改变，使用自适应算法来改变滤波器的参数和结构的滤波法。
1.4.2 程序设计思路 wiener函数用于建立预定义的滤波算子，k = wiener2(I,[m n],noise)其中I是处理的数据，用[m,n]指定相应的参数。故利用该函数进行自适应滤波的程序编写。 自适应滤波实现的程序段如下： global T axes(handles.axes2); T=getimage; k=wiener2(handles.k,[5,5]); imshow(k); title('经过自适应滤波后的图像'); handles.k=k; guidata(hObject,handles);
1.3.3 结果
经过一次自适应滤波后 6
二、频域增强 2.1 低通滤波 2.2.1 知识点 物理上可以实现的一种低通滤波器是巴特沃斯低通滤波器。一个阶为
n，截断频率为0D的巴特沃斯高通滤波器的转移函数为
nvuDDvuH20)],(/[11),(
低通巴特沃斯滤波器在高低频率间的过渡比较光滑，所以用巴特沃斯低通滤波器得到的输出图其振铃现象不明显。具体来说，阶为1时没有振铃现象，而随着阶的增加振铃现象也增加。 一般情况下，常取使H最大值降到某个百分比的频率为截断频率。当
D(u，v）=0D时，H(u,v)=0.5。
2.1.2 程序设计思路 通过对离散傅立叶变换的学习，我们知道一个有限长序列可以用N个频域抽样值唯一的确定，即其中H(k)是h(n)的离散傅立叶变换，是H(z)在单位圆上的均匀抽样值，其列长N，即直接从频域出发，对理想频响进行抽样，然后，用离散傅立叶反变换求单位脉冲响应再用Z变换求系统函数。 低通滤波实现的程序段如下： axes(handles.axes2); k1=handles.k; f=double(k1);% 数据类型转换，matlab不支持图像的无符号整型的计算 g=fft2(f);% 傅里叶变换 g=fftshift(g);% 转换数据矩阵 [M,N]=size(g); nn=2;%二阶巴特沃斯低通滤波器 d0=50;%截止频率50 m=fix(M/2); n=fix(N/2); for i=1:M for j=1:N d=sqrt((i-m)^2+(j-n)^2); h=1/(1+0.414*(d/d0)^(2*nn));% 计算低通滤波器传递函数 result(i,j)=h*g(i,j); end end result=ifftshift(result);
