双直方图均衡化

function [J]=BHE（I）
%输入原图
%双直方图算法
%2012.11.12
%叶爱萍

%%求的整幅图的平均像素值
N=256;%65536
%%求的整幅图的平均像素值
m=round（mean（mean（I）））;
nLow=0;%统计低灰度的个数
nHigh=0;%统计高灰度的个数

%%分别统计两个范围内的灰度直方图
I=double（I）;%转换成double型，matlab不允许对uint8型进行计算%%I是原始图像
[hw]=size（I）;%得到图像高和宽
nhist = zeros（1N）;
ncdf=zeros（1N）;
for y=1:h
    for x=1:w
        if I（yx）<=m
            nhist（I（yx）+1）=nhist（I（yx）+1）+1;%统计低灰度部分直方图
            nLow=nLow+1;
        else
            nhist（I（yx）+1）=nhist（I（yx）+1）+1;%统计高灰度部分直方图
            nHigh=nHigh+1;