精选文档图像处理实验3

昆明理工大学信息工程与自动化学院学生实验报告（2016—2017学年第一学期）课程名称：数字图像基础开课实验室：2016年月日年级、专业、班学号姓名成绩实验项目名称图像增强（1）指导教师教师评语该同学是否了解实验原理： A.了解□ B.基本了解□ C.不了解□该同学的实验能力： A.强□ B.中等□ C.差□该同学的实验是否达到要求： A.达到□ B.基本达到□ C.未达到□实验报告是否规范： A.规范□ B.基本规范□ C.不规范□实验过程是否详细记录： A.详细□ B.一般□ C.没有□教师签名：2016年11一、实验目的及内容目的：掌握和熟悉Matlab编程环境及语言；掌握直方图统计的算法和用途。内容：调试教材P25页例2.1输出类似教材图2.3的结果。调试教材P33页例2.4，编写一个程序，分别使用imhist、bar、stem、plot四种方式显示一幅灰度图像的直方图调试教材P37页例2.5。直方图均衡化的公式如下所示：

根据上式及课堂所讲直方图均衡化原理及方法，自己写一个Matlab函数实现对灰度图像的直方图均衡化功能（类似于Matlab提供的histeq函数）。（提示：实现中使用Matlab函数cumsum（P38）可能会使程序简单些）。二、要求描述直方图的概念并解释直方图均衡化原理。程序结构清晰，运行结果正确。对于第1、2、3小题在实验报告中给出所调试的程序，及其运行结果，对第4小题描述程序的设计、实现和结果，并对结果进行分析。

一、描述直方图的概念并解释直方图均衡化原理。直方图的概念：直方图是用来量化曝光量，能够使我们真实、直观地看出照片曝光情况的一张二维坐标系，其横轴代表的是图像中的亮度，从左向右，从全黑逐渐过渡到全白；纵轴代表的则是图像中处于这个亮度范围的像素的相对数量。当直方图中的黑色色块偏向于左边时，说明这张照片的整体色调偏暗，也可以理解为照片欠曝。而当直方图中的黑色色块集中在右边时，，说明这张照片整体色调偏亮，除非特殊构图需要，，否则我们可以理解为照片过曝。直方图均衡化原理：直方图均衡化是图像处理领域中利用图像直方图对对比度进行调整的方法。这种方法通常用来增加许多图像的局部对比度，尤其是当图像的有用数据的对比度相当接近的时候。通过这种方法，亮度可以更好地在直方图上分布。这样就可以用于增强局部的对比度而不影响整体的对比度，直方图均衡化通过有效地扩展常用的亮度来实现这种功能。二、程序代码以及运行结果第一题：代码：f=imread('breast.tif');>>imshow(f)>>g1=imadjust(f,[01],[10]);>>imshow(g1)>>g2=imadjust(f,[0.50.75],[01]);>>imshow(g2)>>g3=imadjust(f,[],[],2);>>imshow(g3)>>g=imadjust(f,stretchlim(f),[]);>>imshow(g)>>g=imadjust(f,stretchlim(f),[10]);>>imshow(g)第二题：代码：f=imread('pollen.tif');>>imhist(f);>>h=imhist(f,25);>>horz=linspace(0,255,25);>>bar(horz,h)>>axis([0255060000])>>set(gca,'xtick',0:50:255)>>set(gca,'ytick',0:20000:60000)>>stem(horz,h,'fill')>>axis([0255060000])>>set(gca,'xtick',0:50:255)>>set(gca,'ytick',0:20000:60000)>>hc=imhist(f);>>plot(hc)%Usethedefaultvalues.>>axis([0255015000])>>set(gca,'xtick',0:50:255)>>set(gca,'ytick',0:2000:15000)截图:第三题：代码：f=imread('pollen.tif');>>imshow(f);>>figure,imhist(f)>>ylim('auto')>>g=histeq(f,256);>>figure,imshow(g);>>figure,imhist(g);>>ylim('auto')>>截图：第四题：代码：functionE=Example(I)%UNTITLED此处显示有关此函数的摘要%此处显示详细说明[m,n]=size(I);GP=zeros(1,256);fork=0:255GP(k+1)=length(find(I==k))/(m*n);endS1=zeros(1,256);fori=1:256forj=1:iS1(i)=GP(j)+S1(i);endendS2=round((S1*256)+0.5);fori=1:256GPeq(i)=sum(GP(find(S2==i)));endfigure(1);bar(0:255,GP,'b');figure(2);bar(0:255,GPeq,'b');X=I;fori=0:255X(find(I==i))=S2(i+1);endfigure(3);imshow(I);fig