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;
???