光线补偿的matlab源码

function im=LightCompensate（I）

[m0n0l]=size（I）;
%figure（1）imshow（I）
thresholdco=0.05;   %比例系数
thresholdnum=100;   %像素个数的临界常数
histogram=zeros（1256）;  %灰度级 数组
if m0*n0*thresholdco    disp（‘输入图像太小，请换一张！‘）;
    return
end
gray=0;
index0=0;
for i=1:m0        %scan image
     for j=1:n0
          gray=round（I（ij1）*.299+I（ij2）*.587+I（ij3）*.114）;
          index0=gray+1;
          histogram（1index0）= histogram（1index0）+1;
     end
 end
 calnum=0;
 total=m0*n0;
 num=0;
 %next获得满足系数thresholdco的临界灰度级
 index1=0;
 for i=1:256
     if calnum/total         index1=256-i+1;
         calnum=calnum+histogram（1index1）;
         num=i;
     e