实验二数字图像处理孙博文.doc

昆明理工大学信息工程与自动化学院学生实验报告 （ 2016 — 2017 学年 第 一 学期 ） 年级、专业、班 2014级计科143班 学号 201410405325 姓名 孙博文 成绩 实验项目名称 图像灰度变换 指导教师 张 云 教 师 评 语 该同学是否了解实验原理： A.了解□　　B.基本了解□　C.不了解□ 该同学的实验能力： A.强　□ B.中等 □ C.差 □ 该同学的实验是否达到要求： A.达到□　 B.基本达到□　C.未达到□ 实验报告是否规范： A.规范□　 B.基本规范□　C.不规范□ 实验过程是否详细记录： A.详细□　 B.一般 □　C.没有 □ 教师签名： 年 月 日 实验二、图像灰度变换 实验目的 理解图像增强的概念； 掌握灰度线性变换和非线性变换的图像增强方法及应用； 掌握灰度直方图均衡化的图像增强方法及应用。 实验原理 术语“空域”指的是图像平面本身。在空域上，图像处理方法是通过直接对图像像素的处理来实施的。空域处理方法可分为两种：灰度级变换与空域滤波。空域技术所使用的对像素的直接操作可用下式表示： ， 其中f(x,y)为输入图像，g(x,y)为输出图像，T是对图像f进行处理的操作符，它定义在像素点(x,y)所指定的空间邻近像区（简称为邻域）内。? 定义像素点(x,y)的邻域的主要方法是，使用以(x,y)为中心的方形或矩形像区。当这样的邻域的中心从左上角原点的开始遍历像素点时，将覆盖图像中不同的像区。当T应用于每个被遍历的像素点f(x,y)时，便得到在该像素点的输出图像g(x,y)。在计算g(x,y)时，只使用在(x,y)邻域中的像素。? 在处理灰度图像并且当邻域大小为1×1时，T的形式最简单，它成为一个灰度（或亮度）级的变化函数（即，灰度变换）。此时,g(x,y)的值仅由在(x,y)处的亮度f(x,y)来决定。由于灰度变换T仅取决于灰度的值，与(x,y)无关，所以T可写成如下的简单形式： ， 其中，r表示图像f中相应点(x,y)的灰度，s表示图像g中相应点(x,y)的灰度。 灰度变换T通常包括灰度变换增强与直方图增强。灰度变换增强可以通过改变图象的灰度范围及分布来实现。在实际应用中,为了突出图像中感兴趣的研究对象，常常要求局部扩展拉伸某一范围的灰度值，或对不同范围的灰度值进行不同的拉伸处理，即分段拉伸。另一方面，直方图增强是以累积分布函数变换法为基础的直方图修正法。 所用仪器、材料 微型计算机一台、Matlab软件。 实验过程、内容 实验内容 输入图像读入。 灰度变换：请自己编程和调用Matlab中的图像调整函数（imadjust）对输入图像进行分段线性变换和非线性变换，以达到理想的图像视觉效果，同屏显示处理前后的图像及其灰度直方图，比较处理效果的异同，并解释图像灰度调整后其直方图的分布情况。 直方图均衡化：请自己编程和调用Matlab中函数来统计原图像的灰度直方图，并利用直方图均衡化处理进行图像增强，同屏显示处理前后的图像及其灰度直方图，比较与原图的异同，并解释图像均衡化后其直方图分布情况。 附：灰度变换的基本代码 （A）灰度线性变换（函数imadjust(f,[low_in; high_in],[low_out; high_out],gamma)） clc,clear; %清空命令窗口，清空工作空间的变量 f=imread(blob.tif); subplot(2,2,1) imshow(f),title(灰度图); subplot(2,2,2) imhist(f,256); %显示图像直方图（256个分组） title(灰度图的直方图); g1=imadjust(f,[0 1],[1 0]); %灰度转换，实现明暗转换(负片图像) subplot(2,2,3) imshow(g1),title(灰度图的负片图像); g2=imadjust(f,[0.5;0.75],[0;1],1); %将f的[0.5,0.75]线性映射到g2的[0,1] subplot(2,2,4) imshow(g2),title(灰度级的线性变换); （B）灰度非线性变换（对数映射） clc,clear; %清空命令窗口，清空工作空间的变量 f=imread(point.tif); subplot(2,2,1) imshow(f),title(灰度图