欢迎来到天天文库
浏览记录
ID:37183365
大小:28.01 KB
页数:3页
时间:2019-05-21
《《(最新)直方图均衡报告》》由会员上传分享,免费在线阅读,更多相关内容在行业资料-天天文库。
1、这几次上机编程,由于我基础太差,都没有做出自己的程序,因此我改变策略,先读懂老师给的程序,这也是算是一点收获吧。下面就老师的程序,我解读一下。%直方图均衡clcpicture=input('请输入图片文件名:','s');输入一个黑白图像文件fid=fopen(picture);打开该文件dat=fread(fid);将该文件中的像素点读出来生成一个单列矩阵。d=dat';转置成一个一行的行矩阵。h=length(d);newname=input('输入新的文件名:','s');为后面生成的新图像文件命名N=zero
2、s(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)=p(
3、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);将均衡后的像素存入M中fo
4、ri=1:hk=d(i);M(i)=r(k+1);这一句不理解,希望得到指导。endjht=fopen(newname,'wb');创建新图片文件g=fwrite(jht,M);fclose('all');我对直方图均衡的理解:在数学上通过对图像中各灰度级像素依次往后的累加概率进行量化,使得原先靠近的灰度级合并,从而使得合并后该灰度级像素的增多,该灰度级得到加强。实质是牺牲了灰度级,拉大了图像对比度。电信科0803班吕鹏辉0807080324我的邮箱P616942659@163.com
此文档下载收益归作者所有