绘制数字图像灰度直方图实验报告MATLAB实现.pdfVIP

数字图像处理 实验报告 实验一 绘制直方图 学号 姓名 日期 实验一 绘制直方图 一、实验内容 1、编程绘制数字图像的直方图。 2、直方图均衡处理。 二、实验步骤 1、设计思想或者流程图。 灰度直方图是将数字图像的所有像素，按照灰度值的大小，统计其所出 现的频度。通常，灰度直方图的横坐标表示灰度值，纵坐标为半个像素个数， 也可以采用某一灰度值的像素数占全图像素数的百分比作为纵坐标。 直方图均衡方法的基本原理是：对在图像中像素个数多的灰度值 （即对 画面起主要作用的灰度值）进行展宽，而对像素个数少的灰度值 （即对画面 不起主要作用的灰度值）进行归并。从而达到清晰图像的目的。 2、源程序并附上注释。 clear all %一，图像的预处理，读入彩色图像将其灰度化 PS imread(1.jpg); %读入JPG彩色图像文件 imshow(PS) %显示出来 title(输入的彩色JPG 图像) imwrite(rgb2gray(PS),PicSampleGray.bmp);%将彩色图片灰度化并保存 PS rgb2gray(PS); %灰度化后的数据存入数组 %二，绘制直方图 [m,n] size(PS); %测量图像尺寸参数 GP zeros(1,256); %预创建存放灰度出现概率的向量 fork 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 forj 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); 3、程序运行结果比较。（包括结论） 原始图像 均衡化后图象 原图像直方图 均衡化后直方图 通过原始图像和均衡化后的图像对比，均衡化后图像由彩色图像变为了 灰度图像。原始图像直方图与均衡化后直方图对比，均衡化后的直方图灰度