retinex算法代码MSR

搜集到的几个基于retinex算法的代码，用的是matlab，测试通过

f=imread（‘01_test.tif‘）;
fr=f（::1）; fg=f（::2）; fb=f（::3）;
mr=im2double（fr）; mg=im2double（fg）; mb=im2double（fb）;
n=141;%定义模板大小。 kid=141；
n1=floor（（n+1）/2）;%确定中心
a1=60; %定义标准差（尺度） kid=60；
for i=1:n
for j=1:n
b（ij） =exp（-（（i-n1）^2+（j-n1）^2）/（a1*a1））/（pi*a1*a1）; %高斯函数。
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+eps）;tg1=log（mg+eps）;tb1=log（mb+eps）;
yr1=（tr1-ur1）/3;yg1=（tg1-ug1）/3;yb1=（tb1-ub1）/3;
a2=10; %定义标准差（尺度）
for i=1:n
33
for j=1:n
a（ij） =exp（-（（i-n1）^2+（j-n1）^2）/（a2*a2））/（pi*a2*a2）; %高斯函数。
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+eps）;tg2=log（mg+eps）;tb2=log（mb+eps）;
yr2=（tr2-ur2）/3;yg2=（tg2-ug2）/3;yb2=（tb2-ub2）/3;
a3=150; %定义标准差（尺度）kid=150；
for i=1:n
for j=1:n
e（ij） =exp（-（（i-n1）^2+（j-n1）^2）/（a3*a3））/（pi*a3*a3）; %高斯函数。
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+eps）;tg3=log（mg+eps）;tb3=log（mb+eps）;
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（3crcgcb）;
subplot（121）;imshow（f）;title（‘原图‘）;
subplot（122）;imshow（z）;title（‘MSR 处理后‘）;%显示处理后的图像

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件        1627  2019-04-14 23:03  MSR.m
     文件        2301  2019-04-13 23:07  MSRCR.m
     文件        3823  2019-04-13 02:55  SSR.m