数字图象处理实验报告--直方图应用..docVIP

齐鲁工业大学实验报告 课程名称 数字图像处理 指导教师 任向河 实验日期 2014.6.5 院（系） 理学院 专业班级 光信11-1 实验地点 三机房 学生姓名 董萌 学号 20111021028 同组人 实验项目名称 直方图应用 一、 实验目的和要求 编程绘制数字图像的直方图。 直方图均衡处理。 二、 实验原理 设计思想或者流程图。 灰度直方图是将数字图像的所有像素，按照灰度值的大小，统计其所出现的频度。通常，灰度直方图的横坐标表示灰度值，纵坐标为半个像素个数，也可以采用某一灰度值的像素数占全图像素数的百分比作为纵坐标。 直方图均衡方法的基本原理是：对在图像中像素个数多的灰度值（即对画面起主要作用的灰度值）进行展宽，而对像素个数少的灰度值（即对画面不起主要作用的灰度值）进行归并。从而达到清晰图像的目的。 编辑matalab程序如下 首先将要处理的图片拖到window command 窗口 imshow(PS) %显示出来 title(输入的彩色JPG图像) imwrite(rgb2gray(PS),PicSampleGray.bmp); %将彩色图片灰度化并保存 PS=rgb2gray(PS); %灰度化后的数据存入数组 %二，绘制直方图 [m,n]=size(PS); %测量图像尺寸参数 GP=zeros(1,256); %预创建存放灰度出现概率的向量 for k=0:255 GP(k+1)=length(find(PS==k))/(m*n); %计算每级灰度出现的概率，将其存入GP中相应位置 end figure,bar(0:255,GP,g) %绘制直方图 title(原图像直方图) xlabel(灰度值) ylabel(出现概率) %三，直方图均衡化 S1=zeros(1,256); for i=1:256 for j=1:i S1(i)=GP(j)+S1(i); %计算Sk end end S2=round((S1*256)+0.5); %将Sk归到相近级的灰度 for i=1:256 GPeq(i)=sum(GP(find(S2==i))); %计算现有每个灰度级出现的概率 end figure,bar(0:255,GPeq,b) %显示均衡化后的直方图 title(均衡化后的直方图) xlabel(灰度值) ylabel(出现概率) %四，图像均衡化 PA=PS; for i=0:255 PA(find(PS==i))=S2(i+1); %将各个像素归一化后的灰度值赋给这个像素 end figure,imshow(PA) %显示均衡化后的图像 title(均衡化后图像) imwrite(PA,PicEqual.bmp); 三、 主要仪器设备 装有matalab 的计算机1台 四、 实验方法与步骤 1.将所有编写好的程序放在复制下来 2.打开在window command 命令窗口 3.将图片拖动到该窗口下 该窗口显示 uiopen(C:\Users\TOSHIBA\Desktop\22.jpg,1) 4.将程序复制粘贴在该命令窗口下 clear all %一，图像的预处理，读入彩色图像将其灰度化 PS=imread(C:\Users\TOSHIBA\Desktop\22.jpg); %读入JPG彩色图像文件 imshow(PS) %显示出来 title(输入的彩色JPG图像) imwrite(rgb2gray(