图像增强与MATLAB实现.doc

学习 参考 数字图像处理数字图像处理 课程设计课程设计 课设题目课设题目 图像增强与 MATLAB 实现 学校学院学校学院 华东交通大学理学院 学生班级学生班级 13 级信息计算 2 班 学生姓名学生姓名 魏超 学生学号学生学号 20130810010216 指导老师 指导老师 范自柱 学习 参考 图像增强与 MATLAB 实现 摘摘 要要 数字图像处理是指将图像信号转换成数字格式并利用计算机对其进行处理 的过程 图像增强是数字图像处理的过程中经常采用的一种方法 它对提高图 像质量起着重要的作用 本文先对图像增强的原理进行概述 然后对图像增强 的方法分类并给出直方图增强 对比度增强 平滑和锐化等几种常用的增强方 法的理论基础 通过 Matlab 实验得出的实际处理效果来对比各种算法的优缺点 讨论不同的增强算法的技术要点 并对其图像增强方法进行性能评价 关键字 关键字 图像 图像增强 算法 学习 参考 目 录 一 一 MATLABMATLAB 的简介的简介 1 1 1 1 MATLAB 主要功能 1 二 二 MATLABMATLAB 的主要功能的主要功能 1 1 2 1 数字增强技术概述 1 2 2 数字图像的表示 1 三 直方图的均衡化三 直方图的均衡化 2 2 3 1 图像的灰度 2 3 2 灰度直方图 2 3 3 直方图均衡化 3 四 图像二值化四 图像二值化 5 5 4 1 图像二值化 5 五 对比度增强五 对比度增强 6 6 5 1 对比度增强 6 5 2 灰度调整 6 5 3 对数变换 7 六 滤波六 滤波 8 8 6 1 平滑滤波 8 6 2 线性平滑滤波程序 8 6 3 非线性滤波 9 七 锐化七 锐化 1313 八 参考文献八 参考文献 1414 九 自我评价九 自我评价 1414 学习 参考 一 一 MatlabMatlab 的简介的简介 1 1 MATLAB 主要功能 MATLAB 是建立在向量 数组和矩阵基础上的一种分析和仿真工具软件包 包含各种能够进行常规运算的 工具箱 如常用的矩阵代数运算 数组运算 方程求根 优化计算及函数求导积分符号运算等 同时还提供了编程计算的编 程特性 通过编程可以解决一些复杂的工程问题 也可绘制二维 三维图形 输出结果可视化 目前 已成为工程领域中较常用的软件工具包之一 二 二 MATLABMATLAB 的主要功能的主要功能 2 1 数字增强技术概述 图像增强是按特定的需要突出一幅图像中的某些信息 同时 消弱或去除 某些信息使得图像更加实用 图像增强技术主要包含直方图修改处理 图像平 滑处理 图像尖锐化处理等 图像增强技术主要包括 直方图修改处理 图像平滑处理 图像尖锐化处 理 彩色图像处理 从纯技术上讲主要有两类 频域处理法和空域处理法 频域处理法主要是卷积定理 采用修改图像傅立叶变换的方法实现对图像 的增强处理技术 空域处理法 是直接对图像中的像素进行处理 基本上是以 灰度映射变换为基础的 2 2 数字图像的表示 图像并不能直接用计算机来处理 处理前必须先转化成数字图像 由于从 外界得到的图像多是二维 2 D 的 一幅图像可以用一个 2 D 数组表示 yxf 这里 x 和 y 表示二维空间 XY 中一个坐标点的位置 而 f 则代表图像在点的 yx 某种性质数值 为了能够用计算机对图像进行处理 需要坐标空间和性质空间 都离散化 学习 参考 三 直方图的均衡化三 直方图的均衡化 3 1 图像的灰度 常用的图像一般是灰度图 这时 f 表示灰度值 反映了图像上对应点的亮 度 亮度是观察者对所看到的物体表面反射光强的量度 作为图像灰度的量度 函数应大于零 人们日常看到的图像一般是从目标上反射出来的光组成 yxf 的 所以可看成由两部分构成 入射到可见场景上光的量 场景中目标 yxf 对反射光反射的比率 确切地说它们分别称为照度成分和反射成分 yxi 与和都成正比 可表示成 yxr yxf yxi yxr yxf yxi yxr 3 2 灰度直方图 灰度变换是图像增强的一种重要手段 使图像对比度扩展 图像更加清晰 特征更加明显 灰度级的直方图给出了一幅图像概貌的描述 通过修改灰度直 方图来得到图像增强 程序如下 灰度直方图 I imread D image rice bmp 读取图像 subplot 2 1 1 Imshow I 显示图像 title a 原图 subplot 2 1 2 imhist I 绘制图像的灰度直方图 title b 原图的灰度直方图 学习 参考 图 4 1 3 3 直方图均衡化 直方图均衡化过程如下 1 计算原图像的灰度直方图 Kr rP 2 计算原图像的灰度累积分布函数 进一步求出灰度变换表 k s 3 根据灰度变换表 将原图像各灰度级映射为新的灰度级 灰度直方图均衡化程序 灰度直方图均衡化 I imread D image bubbles bmp 读取图像 subplot 2 2 1 Imshow I 显示图像 title 原图 subplot 2 2 2 imhist I 绘制图像的灰度直方图 title 原图的灰度直方图 subplot 2 2 3 J histeq I 64 对图像进行均衡化处理 返回有 64 级灰度的图像 J Imshow J 显示图像 title 原图直方图均衡化 subplot 2 2 4 imhist J 绘制图像的灰度直方图 title 均衡后的灰度直方图 学习 参考 以下展示了直方图均衡化的效果 图4 2 优势 能够使得处理后图像的概率密度函数近似服从均匀分布 其结果扩 张了像素值的动态范围 是一种常用的图像增强算法 不足 不能抑制噪声 四 图像二值化四 图像二值化 4 1 图像二值化 图像的二值化处理就是将图像上的像素点的灰度值设置为 0 或 255 也就 是讲整个图像呈现出明显的黑白效果 首先 图像的二值化有利于图像的进一 步处理 使图像变得简单 而且数据量减小 能凸显出感兴趣的目标的轮廓 其次 要进行二值图像的处理与分析 首先要把灰度图像二值化 得到二值化 图像 图像二值化程序 图像二值化 选取一个域值 5 将图像变为黑白图像 I imread beauty tif bw im2bw I 0 5 选取阈值为 0 5 subplot 1 3 1 学习 参考 imshow I title 原图 subplot 1 3 2 imshow bw title 显示二值图像 J find I 150 I J 255 subplot 1 3 3 imshow I title 图像二值化 域值为 150 下图为图像二值化的效果 图 4 3 五 对比度增强五 对比度增强 5 1 对比度增强 对比度增强是按一定的规则修改输入图像每一个像素的灰度 从而改变图 像灰度的动态范围 例如 观察图 4 4 可以发现 该图的对比度不高其灰度直 方图没有低于 35 或高于 210 的值 如果将图像数据映射到整个灰度范围内 则 图像的对比度将大大增大 学习 参考 图 4 4 5 2 灰度调整 灰度调整程序 imadjust 函数 I imread D image rice bmp 读取图像 subplot 2 2 1 Imshow I 显示图像 title 原图 subplot 2 2 2 imhist I 绘制图像的灰度直方图 title 原图的灰度直方图 subplot 2 2 3 J imadjust I 0 3 0 7 对图像进行灰度变换 Imshow J 显示图像 title 原图直方图均衡化 subplot 2 2 4 imhist J 绘制图像的灰度直方图 title 均衡后的灰度直方图 以下展示了常用对比度扩展法的结果 学习 参考 图 4 5 从图 4 5 原图 可以看出原始图像动态范围较小 整体较暗 反映在直 方图上像素主要集中在低灰度的一侧 经过对比度调整 图像变亮 5 3 对数变换 对数变化常用来扩展低值灰度 压缩高值灰度 这样可以使低值灰度的图 像细节更容易看清 对数变换程序 对数变换程序 对数变换 I imread beauty tif 读取图像 I mat2gray I 对数变换不支持 uint8 类型数据 将一个矩阵转化为灰度 图像的数据格式 double J log I 1 subplot 1 2 1 Imshow I 显示图像 title 原图 subplot 1 2 2 Imshow J title 对数变换后的图像 学习 参考 图 4 6 对数变换确实能够扩展低值灰度 而压缩高值灰度 使低值灰度的图像细 节更容易看清 六 滤波六 滤波 6 1 平滑滤波 平滑技术用于平滑图像中的噪声 平滑噪声可以在空间域中进行 基本 方法是求像素灰度的平均或中值 为了既平滑噪声又保护图像信号 6 2 线性滤波 输出图像的值等于输入图像滤波后值的局部平均 各个项具有相同的权 下面是平滑窗口分别为矩形和圆形的情况 11111 11111 11111 11111 11111 25 1 kjhrect 01110 11111 11111 11111 01110 21 1 kjhcirc 对一些图像进行线性滤波可以去除图像中某些类型的噪声 如采用邻域平 均法的均值滤波器就非常适用于去除通过扫描得到的图像中的颗粒噪声 6 2 线性平滑滤波程序 线性平滑滤波 I imread beauty tif J imnoise I salt subplot 221 imshow I title 原图像 学习 参考 subplot 222 imshow J title 添加椒盐噪声图像 K1 filter2 fspecial average 3 J 255 应用 3 3 邻域窗口法 subplot 223 imshow K1 title 3x3 窗的邻域平均滤波图像 K2 filter2 fspecial average 7 J 255 应用 7 7 邻域窗口法 subplot 224 imshow K2 title 7x7 窗的邻域平均滤波图像 如图 4 8 a 为线性平滑滤波的例子 6 3 非线性滤波 中值滤波是一种最常用的图像增强技术 是非线性滤波 对椒盐噪声有 很好的去噪效果 下图是加高斯噪声后 中值滤波和平均滤波的滤波效果 程序 中值滤波和平均滤波 I imread girl bmp I rgb2gray I J imnoise I gaussian 0 0 01 subplot 2 2 1 imshow I title 原图 subplot 2 2 2 imshow J title noise 学习 参考 K fspecial average 5 K1 filter2 K J 255 subplot 2 2 3 imshow K1 title 平均滤波 L medfilt2 J 3 5 subplot 2 2 4 imshow L title 中值滤波 对比 图 4 9 a 但对于椒盐噪声 中值滤波效果比均值滤波效果好 中值滤波程序 二维中值滤波 hood 3 P map imread leno bmp I rgb2gray P imshow I map noisy imnoise I salt 加入强度为 0 05 的椒盐噪声 学习 参考 subplot 221 imshow noisy map title 加入椒盐噪声 filtered1 medfilt2 noisy hood hood 3 3 窗口二维中值滤波 subplot 222 imshow filtered1 map title 3 3 窗口 hood 5 filtered2 medfilt2 noisy hood hood 5 5 窗口二维中值滤波 subplot 223 imshow filtered2 map title 5 5 窗口 hood 7 filtered3 medfilt2 noisy hood hood 7 7 窗口二维中值滤波 subplot 224 imshow filtered3 map title 7 7 窗口 I imread leno bmp figure imshow I title 原图 实现中值滤波的效果图如下 图 4 9 b 比较发现 图像的噪声点被去除 图像的边缘稍微的变得模糊 高通滤波边缘增强程序 高通滤波边缘增强 I imread girl bmp 学习 参考 I rgb2gray I subplot 2 2 1 imshow I title original pic J fspecial average 3 J1 conv2 I J 255 subplot 2 2 2 imshow J1 title 3 3lowpass K fspecial prewitt K1 filter2 K J1 5 subplot 2 2 3 imshow K1 title prewitt L fspecial sobel L1 filter2 L J1 5 subplot 2 2 4 imshow L1 title sibel 下图为高通滤波边缘增强的例子 图 4 9 c 优势 去噪效果明显 并且能够较好的保持图像边缘位置和细节 不足 非线性滤波算法的实现相对线性滤波比较困难 学习 参考 七 锐化七 锐化 图像锐化处理的作用是使灰度反差增强 从而使模糊图像变得更加清晰 图像模糊的实质就是图像受到平均运算或积分运算 因此可以对图像进行逆运 算 如微分运算以突出图像细节使图像变得更为清晰 锐化程序 锐化 a imread dowels tif subplot 131 imshow a title 原图 b double a 将图像矩阵转化为 double 类型 s size b c zeros s 1 1 s 1 2 for x 2 s 1 1 1 for y 2 s 1 2 1 c x y b x 1 y b x 1 y b x y 1 b x y 1 4 b x y end end 用拉氏算子对图像进行滤波 这个过程相当于运用了一个 3 3 的掩膜 0 1 0 1 4 1 0 1 0 subplot 132 imshow c title Laplace 锐化滤波图像 d b c 当拉普拉斯掩膜中心系数为正时 增强图像为原图像于拉氏算子滤 波图像之和 d uint8 d 将图像矩阵变回 uint8 格式 subplot 133 imshow d title Laplace 锐化滤波结果 图 4 10 展示了在 MATLAB 中进行图像锐化的处理结果 图4 10 比较发现 图像模糊的部分得到了锐化 特别是模糊的边缘部分得到了增 学习 参考 强 边界更加明显 但是 图像显示清楚的地方 经过滤波却发生了失真 八 参考文献八 参考文献 1 汪志云 黄梦为 基本直方图的图像增强及其 MATLAB 实现 J 计算机工程 与科学 2006 28 2 54 56 2 刘榴娣 刘明奇 党长民 实用数字图像处理 M 北京理工大学出版社 1998 3 冈萨雷斯 RC 著 李叔梁译 数字图像处理 北京 科学出版社 1983 4 李月景编著 图像识别技术及其应用 北京 机械工业出版社 1985 九 自我评价九 自我评价 这次课程设计不仅使我对数字图像处理的知识有了进一步的了