MEBHE.rar_MEBHE_depend499

clc %读入图片并显示 I=imread('C:\Users\HJ\Desktop\1.bmp'); I=rgb2gray(I); %将rgb图像转化为灰度图像 h1=imhist(I); %返回每个灰度的像素数目和对应灰度 figure,imshow(I); %显示原图像 title('原始图像') figure,imhist(I); %显示原图像直方图 title('原始图像直方图'); %直方图均值归一化 [H,W]=size(I); %计算图像大小 h1=h1/H/W*256; %直方图均值归一化 figure,set(stem(h1),'marker','none'); title('均值归一化直方图'); %直方图预处理 h2=zeros(size(h1)); %产生一全0数组 gamma=0.125; %设置gamma校正值 T=4; %设置钳位值 for i=1:256 %进行预处理操作：小于1进行gamma校正 if(h1(i,1)<=1) h2(i,1)=h1(i,1).^gamma; else h2(i,1)=h1(i,1); end end for i=1:256 %进行预处理操作：大于T钳位 if(h2(i,1)<=T) h2(i,1)=h2(i,1); else h2(i,1)=T; end end figure,set(stem(h2),'marker','none'); title('预处理后直方图'); %遍历双直方图均衡寻找Grs GMin=min(I(:)); %获得最大最小灰度值 GMax=max(I(:)); Gds=graythresh(I)*255; %通过大津法Ostu获得Gds B=zeros(size(h1)); %产生全0数组，存放熵值 for k=55:200 %从最小灰度值遍历到最大灰度值 Grs=k; %双直方图均衡操作?????????????????简易步骤直接调用函数？ nl=sum(h2(GMin+1:Grs+1)); %Grs左边求直方图总和(加1是因为0-255对应1-256) nr=sum(h2(Grs+2:GMax+1)); %Grs右边求和 M=zeros(size(I)); %产生一图像大小数组存放均衡后图像 for i=1:H %均衡操作 for j=1:W x=I(i,j); if x<=Grs cl=sum(h2(1:x+1))/nl; M(i,j)=GMin+(Grs-GMin)*cl; else cr=sum(h2(Grs+2:x+1))/nr; M(i,j)=(Grs+1)+(GMax-(Grs+1))*cr; end end end %计算对应熵值？调用函数？ M3=mat2gray(M); C=imhist(M3); A=C/H/W; S=0; for i=1:256 if A(i,1)~=0; S = -A(i,1)*log2(A(i,1))+S; end end B(k+1,1)=S; end [m,index]=max(B); %求得最大熵值及所属灰度 Gbrs=index-1; fprintf('Gbrs=%f',Gbrs); %利用LUT对照表进行最终双直方图均衡 h3=zeros(size(h2)); %nz=sum(h2(1:256)); nl=sum(h2(GMin+1:Gds+1)); %Grs左边求直方图总和(加1是因为0-255对应1-256) nr=sum(h2(Gds+2:GMax+1)); %Grs右边求和 h3(1,1)=h2(1,1)/nl; for i=2:Gds+1 h3(i,1)=h3(i-1,1)+h2(i,1)/nl; end h3(Gds+2,1)=h2(Gds+2,1)/nr; for i=Gds+3:256 %?范围到gmax或者256 h3(i,1)=h3(i-1,1)+h2(i,1)/nr; end h4=zeros(size(h2)); for i=1:Gds+1 h4(i,1)=(Gbrs-GMin)*h3(i,1)+GMin; end for i=Gds+2:256 h4(i,1)=(GMax-Gbrs)*h3(i,1)+Gbrs; end h4=uint8(h4); M2=intlut(I,h4); %M2=mat2gray(M2); %显示最终均衡化后的图像及其直方图 figure,imshow(M2); title('均衡化后的图像'); h5=imhist(M2); figure,set(stem(h5),'marker','none'); title('均衡化后的直方图');