BBHE/DSIHE/MMBEBHE代码 matlab版

BBHE/DSIHE/MMBEBHE代码 matlab版

clc
clear all
image = imread（‘f:\27.bmp‘）;  
hsvImg = rgb2hsv（image）;  
gf=hsvImg（::3）;  
[heightwidth]=size（gf）;  
gf = uint8（gf*255）;

gfMin=min（gf（:））;
gfMax=max（gf（:））;

%测量图像尺寸参数
[st]=size（gf）;

%计算每个灰度级出现的次数
h1=imhist（gf）;

%计算图像的平均灰度gray_mean
gray_mean=round（sum（gf（:））/（s*t））

    
%计算nlnu
nl=sum（h1（gfMin+1:gray_mean+1））;
nu=sum（h1（gray_mean+2:gfMax+1））;

% 计算BBH直方图均衡化后的输出Y（ij）
for i=1:s
    for j=1:t
        x=gf（ij）;
        if x<=gray_mean
            cl=sum（h1（gfMin+1:x+1））/nl;
            y（ij）=gfMin+（gray_mean-gfMin）*cl;
        else
            cu=sum（h1（gray_mean+2:x+1））/nu;
            y（ij）=（gray_mean+1）+（gfMax-（gray_mean+1））*cu;
        end  
    end
end
y=mat2gray（y）;
%显示BBH均衡化后的图像及其直方图
figureimshow（y）
title（‘BBH均衡化后的图像‘）;

y = im2double（y）; 
hsvImg（::3） = y;  
outputImg = hsv2rgb（hsvImg）;
% subplot（121）imshow（image）title（‘原始图片‘）; 
% subplot（122）imshow（outputImg）title（‘HE-HSV‘）; 
% imwrite（outputImg‘f:\27xxbbhe.bmp‘）;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       1089  2017-08-10 21:34  bbhe.m

     文件        989  2017-08-11 17:01  DSIHE_1999.m

     文件       1307  2017-08-10 21:34  MMBEBHE.m

----------- ---------  ---------- -----  ----

                 3385                    3