k_gray_histeq.rar_K._histeq的 _matlab histeq

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %函数功能，对图像进行直方图均衡 %直接读图像Img10_1.bmp，读到tu中 %graydis是原始直方图各灰度级像素个数 %原始直方图graydispro，利用原始直方图计算原始累计直方图graydispro %t[]计算和原始灰度对应的新的灰度t[]，建立映射关系，t坐标代表原始的灰度，t[]代表对应原始坐标的新坐标 %new_graydis是统计新直方图各灰度级像素个数 %计算新的灰度直方图new_graydispro,利用新的直方图计算新的累计直方图new_graydispro %计算直方图均衡后的新图new_tu %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% clear all close all tu=imread('lena.bmp'); graydis=zeros(1,256); graydispro=zeros(1,256); new_graydis=zeros(1,256); new_graydispro=zeros(1,256); [h w]=size(tu); new_tu=zeros(h,w); %计算原始直方图各灰度级像素个数graydis for x=1:h for y=1:w graydis(1,tu(x,y))=graydis(1,tu(x,y))+1; end end %计算原始直方图graydispro graydispro=graydis./sum(graydis); %计算原始累计直方图 for i=2:256 graydispro(1,i)=graydispro(1,i)+graydispro(1,i-1); end %计算和原始灰度对应的新的灰度t[]，建立映射关系 for i=1:256 t(1,i)=floor(254*graydispro(1,i)+0.5); end %统计新直方图各灰度级像素个数new_graydis for i=1:256 new_graydis(1,t(1,i)+1)=new_graydis(1,t(1,i)+1)+graydis(1,i); end %计算新的灰度直方图new_graydispro new_graydispro=new_graydis./sum(new_graydis); %计算直方图均衡后的新图new_tu for x=1:h for y=1:w new_tu(x,y)=t(1,tu(x,y)); end end figure,imshow(tu,[]) figure,imshow(new_tu,[]);