效率很高的retinex方法 用MATLAB代码实现

这是效率很高的retinex方法 用MATLAB代码实现 分享给大家

I = imread（‘111.bmp‘）;
subplot（121）;
imshow（I）;


I=double（I）;
f=I（::1）;
ff=I（::2）;
fff=I（::3）;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
k1=4;
k2=5;
r=161;
alf=1458;
nn=floor（（r+1）/2）;
for i=1:r
    for j=1:r
        b（ij） =exp（-（（i-nn）^2+（j-nn）^2）/（k1*alf））/（k2*pi*alf*10000）;        %高斯函数1
   end
end

k1=8;                                                                  %高斯函数2
k2=8;
r=161;
alf=1458;
nn=floor（（r+1）/2）;
for i=1:r
    for j=1:r
        bb（ij） =exp（-（（i-nn）^2+（j-nn）^2）/（k1*alf））/（k2*pi*alf*10000）;     
   end
end

k1=0.5;
k2=0.5;
r=161;
alf=1458;
nn=floor（（r+1）/2）;
for i=1:r
    for j=1:r
        bbb（ij） =exp（-（（i-nn）^2+（j-nn）^2）/（k1*alf））/（k2*pi*alf*10000）;        %高斯函数23
end
%%%%%%%%%%%对R分量的处理%%%%%%%%%%%%%
Img = double（f）;
[mn]=size（f）;

aa=125;
for i=1:m
    for j=1:n 
        C（ij）=log（1+aa*（Img（ij）/I（ij）））;
    end
end

K=imfilter（Imgb）;
KK=imfilter（Imgbb）;
KKK=imfilter（Imgbbb）;

for i

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

     文件       2775  2009-03-08 18:33  retinex_frankle_mccann.m

     文件       2612  2009-09-06 20:55  MSRCR.m

     文件       3183  2009-09-07 10:07  MSR_new.m

     文件       2693  2010-04-23 21:00  retinex.m

     文件       5392  2009-03-08 18:33  retinex_mccann99.m

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

                16655                    5