retinex matlab实现带测试图

retinex的matlab实现。 打开retinex.m，点击运行，即可看到原图和处理后的图像。

clc
clear
f=imread（‘2.jpg‘）;
fr=f（: : 1）; fg=f（: : 2）; fb=f（: : 3）;%RGB通道
mr=mat2gray（im2double（fr））; mg=mat2gray（im2double（fg））; mb=mat2gray（im2double（fb））;%数据类型归一化
alf1=1458; %定义标准差alf=a^2/2  a=54
n=161;%定义模板大小 
n1=floor（（n+1）/2）;%计算中心 
for i=1:n 
    for j=1:n 
      b（ij） =exp（-（（i-n1）^2+（j-n1）^2）/（4*alf1））/（pi*alf1）; %高斯函数
    end 
end 
nr1 = imfilter（mrb‘conv‘ ‘replicate‘）;ng1 = imfilter（mgb‘conv‘ ‘replicate‘）;nb1 = imfilter（mbb‘conv‘ ‘replicate‘）;%卷积滤波
ur1=log（nr1）; ug1=log（ng1）; ub1=log（nb1）;
tr1=log（mr）;tg1=log（mg）;tb1=log（mb）;
yr1=（tr1-ur1）/3;yg1=（tg1-ug1）/3;yb1=（tb1-ub1）/3;
alf2=53.38; %定义标准差alf=a^2/2    a=10.3325
x=31;%定义模板大小 
x1=floor（（n+1）/2）;%计算中心 
for i=1:n 
    for j=1:n 
      a（ij） =exp（-（（i-n1）^2+（j-n1）^2）/（4*alf2））/（6*pi*alf2）; %高斯函数
    end 
end 
nr2 = imfilter（mra‘conv‘ ‘replicate‘）;ng2 = imfilter（mga‘conv‘ ‘replicate‘）;nb2 = imfilter（mba‘conv‘ ‘replicate‘）;%卷积滤波
ur2=log（nr2）; ug2=log（ng2）; ub2=log（nb2）;
tr2=log（mr）;tg2=log（mg）;tb2=log（mb）;
yr2=（tr2-ur2）/3;yg2=（tg2-ug2）/3;yb2=（tb2-ub2）/3;
alf3=13944.5; %定义标准差alf=a^2/2  a=167
l=501;%定义模板大小 
l1=floor（（n+1）/2）;%计算中心 
for i=1:n 
    for j=1:n 
      e（ij） =exp（-（（i-n1）^2+（j-n1）^2）/（4*alf3））/（4*pi*alf3）; %高斯函数
    end 
end 
nr3 = imfilter（mre‘conv‘ ‘replicate‘）;ng3 = imfilter（mge‘conv‘ ‘replicate‘）;nb3 = imfilter（mbe‘conv‘ ‘replicate‘）;%卷积滤波
ur3=log（nr3）; ug3=log（ng3）; ub3=log（nb3）;
tr3=log（mr）;tg3=log（mg）;tb3=log（mb）;
yr3=（tr3-ur3）/3;yg3=（tg3-ug3）/3;yb3=（tb3-ub3）/3;
dr=yr1+yr2+yr3;dg=yg1+yg2+yg3;db=yb1+yb2+yb3;
    
cr=im2uint8（dr）; cg=im2uint8（dg）; cb=im2uint8（db）;
z=cat（3 cr cg cb）;          
figure imshow（z）
 

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2017-11-15 09:14  retinex matlab代码\
     文件      267239  2017-10-23 15:15  retinex matlab代码\1.jpg
     文件       12564  2017-10-23 15:15  retinex matlab代码\10.jpg
     文件      305425  2017-10-23 15:15  retinex matlab代码\11.jpg
     文件       25001  2017-10-23 15:15  retinex matlab代码\12.jpg
     文件       14086  2017-10-23 15:15  retinex matlab代码\13.jpg
     文件      178720  2017-10-23 15:15  retinex matlab代码\14.jpg
     文件      326323  2017-10-23 15:15  retinex matlab代码\2.jpg
     文件      325289  2017-10-23 15:15  retinex matlab代码\3.jpg
     文件      544839  2017-10-23 15:15  retinex matlab代码\4.JPG
     文件      202104  2017-10-23 15:15  retinex matlab代码\5.jpg
     文件      312216  2017-10-23 15:15  retinex matlab代码\6.jpg
     文件       98152  2017-10-23 15:15  retinex matlab代码\7.jpg
     文件       43970  2017-10-23 15:15  retinex matlab代码\8.jpg
     文件        6635  2017-10-23 15:15  retinex matlab代码\9.jpg
     文件         664  2017-10-23 15:15  retinex matlab代码\Untitled.m
     文件        1415  2017-10-23 15:15  retinex matlab代码\Untitled2.m
     文件        1774  2017-11-15 09:14  retinex matlab代码\retinex.m
     文件        2009  2017-10-23 15:15  retinex matlab代码\retinex_matlab.txt
     文件        2736  2017-10-23 15:15  retinex matlab代码\tssrcolor20130517.m