直方图均衡报告

《直方图均衡报告》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、这几次上机编程，由于我基础太差，都没有做出自己的程序，因此我改变策略，先读懂老师给的程序，这也是算是一点收获吧。下面就老师的程序，我解读一下。%直方图均衡clcpicture=input('请输入图片文件名：','s');输入一个黑白图像文件fid=fopen(picture);打开该文件dat=fread(fid);将该文件中的像素点读出来生成一个单列矩阵。d=dat';转置成一个一行的行矩阵。h=length(d);newname=input('输入新的文件名：','s');为后面生成的新图像文件命名N=ze

2、ros(1,256);该循环将h这个行矩阵中各灰度级的像素点加权统计存入N这个零矩阵中。N矩阵存放的就是各灰度级的像素数。fori=0:255forj=1:hifi==d(j)N(i+1)=N(i+1)+1;endendendfigure(1)通过对N矩阵的循环访问，各灰度级像素数除以总像素数h，就得到各灰度级的概率。也就是原始直方图。i=0:255;p=N/h;stem(i,p,'.')title('原始直方图')m=zeros(1,256);fork=1:256该循环是求Sk的过程，很简单。ifk==1m(k

3、)=p(k);elsem(k)=m(k-1)+p(k);endendfork=1:256r=round(m*255);量化Sk得到Sk‘ends=zeros(1,256);fora=1:256forb=1:256ifr(a)==r(b)s(a)=N(a)+N(b);该循环统计量化后个灰度级像素数。endendendq=s/h;量化后各灰度级概率figure(2)forj=1:256k(j)=r(j);endstem(k,q,'.');均衡直方图title('均衡直方图');M=zeros(1,h);将均衡后的像素

4、存入M中fori=1:hk=d(i);M(i)=r(k+1);这一句不理解，希望得到指导。endjht=fopen(newname,'wb');创建新图片文件g=fwrite(jht,M);fclose('all');我对直方图均衡的理解：在数学上通过对图像中各灰度级像素依次往后的累加概率进行量化，使得原先靠近的灰度级合并，从而使得合并后该灰度级像素的增多，该灰度级得到加强。实质是牺牲了灰度级，拉大了图像对比度。电信科0803班吕鹏辉0807080324我的邮箱P616942659@163.com