数字图像处理实验报告-直方图规定化0528

1、数字图像处理实验报告-直方 图规定化数字图像处理实验报告Sanz1直方图规定化实验名称: 姓名：班级学号专业电子信息工程（2+2）指导教师：陈华华实验日期：2012年5月24日直方图匹配(规定化)直方图均衡化能够自动增强整个图像的对比度，但它的具体增强效果 不容易控制，处理的结果总是得到全局均匀化的直方图。实际上有时需要 变换直方图，使之成为某个特定的形状，从而有选择地增强某个灰度值范 围内的对比度。这时可以采用比较灵活的直方图规定化。一般来说正确地选择规定化的函数可以获得比直方图均衡化更好的效果。所谓直方图规定化，就是通过一个灰度映像函数，将原灰度直方图改 造成所希望的直方图。所以，直方图修

2、正的关键就是灰度映像函数。直方图匹配方法主要有 3个步骤(这里设 M和N分别为原始图和规定 图中的灰度级数，且只考虑NKM的情况)：(1)如同均衡化方法中，对原始图的直方图进行灰度均衡化：(2)规定需要的直方图，并计算能使规定的直方图均衡化的变换：(3)将第1个步骤得到的变换反转过来，即将原始直方图对应映射到 规定的直方图，也就是将所有pf (fi )对应到pu( uj)去。A图直方图规定B图Fil Cd Vic Inj( Toe Desk, Wind fto *膈人的屈像尸FJ FigureF住髭立工壬R旧 Edit事处r-ttll Tcok4 W ild口璃中口闻昌ll瞰国力 及 口国产:

3、-旧消Matlab 程序:%直方图规定化clear allA=imread(C:UsershpDesktopA.tif);%读入A图像imshow(A)%显示出来title(输入的A图像)%绘制直方图m,n=size(A);%测量图像尺寸B=zeros(1,256);%预创建存放灰度出现概率的向量for k=0:255B(k+1)=length(find(A=k)/(m*n);% 计算每级灰度出现的概率，将其存入B中相应位置 endfigure,bar(0:255,B,g);%绘制直方图title(A图像直方图)xlabel(灰度值)ylabel(出现概率)axis(0,260,0,0.015

4、)C=imread(C:UsershpDesktopB.tif);% 读入 B图像imshow(C)%显示出来title(输入的B图像)%绘制直方图m,n=size(C);%测量图像尺寸D=zeros(1,256);%预创建存放灰度出现概率的向量for k=0:255D(k+1)=length(find(C=k)/(m*n);% 计算每级灰度出现的概率，将其存入 D中相应位 置 endfigure,bar(0:255,D,g);%绘制直方图 title(B图像直方图) xlabel(灰度值) ylabel(出现概率) axis(0,260,0,0.015)S1=zeros(1,256);fo门

5、=1:256for j=1:iS1(i)=D(j)+S1(i);% 计算 B 灰度图累计直方图 end endcounts=Bfigure,bar(0:255,counts,r)title(A图像直方图) S2=zeros(1,256);fo门=1:256for j=1:iS2(i)=counts(j)+S2(i); endend;%累计”规定化直方图%对比直方图，找到相差最小的灰度级 fo门=1:256for j=1:255if S1(j)=S2(i) ifabs(S1(j)-S2(i)0for k=(T(i-1)+1):T(i)H(i)=H(i)+D(k);endelse H(i)=0;e

6、ndendfigure,bar(0:255,H,g)%显示规定化后的直方图title(A规定B后的直方图)xlabel(灰度值)ylabel(出现概率)axis(0,260,0,0.6)%显示规定图PA=C;%将各个像素归一化后的灰度值赋给这个像素for i=1:mfor j=1:nfor k=2:255if T(k-1)=C(i,j)&C(i,j)=T(k)PA(i,j尸k； break;endendendendfigure,imshow(PA)%显示规定化后的图像title(A规定B后图像)imwrite(PA,guidinghua.bmp);二、用已知直方图规定A图规定灰度为,4zero

7、s(1,49),03zeros(1,49),0.1,zeros(1,49),0.2,zeros(1,49),0.1,zeros(19),0.2,zeros(1,49),0.1. Figure- 2ini. Edrf ，；f ：i%e 。& *4ta& ehDIFie Ed I 5m wen TcqI& DeUitoc Window. Hlp，Ur. 盛*星羽E击 yw biCiri rode DhllepHlp d 口白k做戏心g 口目口嶷范化良方爸 。等Me idit 中电 kxrf Tgk DmL r&* 出履H臣Ie段反9启K二I我定七后井通方三D1Q5Matlab 程序：clear a

8、llA=imread(C:UsershpDesktopB.tif);%读入A图像imshow(A)%显示出来title(输入的A图像)%绘制直方图m,n=size(A);%测量图像尺寸B=zeros(1,256);%预创建存放灰度出现概率的向量%计B中相应位置for k=0:255B(k+1)=length(find(A=k)/(m*n);算每级灰度出现的概率，将其存入endfigure,bar(0:255,B,g);绘制直方图title(A图像直方图,)xlabel(灰度值)ylabel(出现概率)axis(0,260,0,0.015)S1=zeros(1,256);fo门=1:256for

9、 j=1:iS1(i尸B(j)+S1(i);%计算原灰度图累计直方图endendcounts=zeros(1,49),0.3,zeros(1,49),0.1,zeros(1,49 ),02zeros(1,49),0.1,zeros(1,49),02zeros(1,49),0. 1;%规定化直方图figure,bar(1:300,counts,r)title(规定化直方图)S2=zeros(1,256);for i=1:256for j=1:iS2(i)=counts(j)+S2(i);endend;%累计规定化直方图%对比直方图，找到相差最小的灰度级fo门=1:256for j=1:256if S1(j)=S2ifabs(S1(j)-S2(i)0for k=(T(i-1)+1):T(i)H(i尸H(i)+B(k);endelse H(i)=0;endendfigure,bar(0:255,H,g)%显示规定化后的直方图title(规定化后的直方图)xlabel(灰度值)ylabel(出现概率)axis(0,260,0,0.6)%显示规定图PA=A；%将各个像素归一化后的灰度值赋给这个像素for