msr

Matlab用Retinex的MSR实现图像增强的代码，可运行，f=imread('c.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/2a=54n=161;%定义模板大小n1=floor((n+1)/2);%计算中心fori=1:nforj=1:nb(i,j)=exp(-((i-n1)^2+(j-n1)^2)/(4*alf1))/(pi*alf1);%高斯函数endendnr1=imfilter(mr,b,'conv','replicate');ng1=imfilter(mg,b,'conv','replicate');nb1=imfilter(mb,b,'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/2a=10.3325x=31;%定义模板大小x1=floor((n+1)/2);%计算中心fori=1:nforj=1:na(i,j)=exp(-((i-n1)^2+(j-n1)^2)/(4*alf2))/(6*pi*alf2);%高斯函数endendnr2=imfilter(mr,a,'conv','replicate');ng2=imfilter(mg,a,'conv','replicate');nb2=imfilter(mb,a,'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/2a=167l=501;%定义模板大小l1=floor((n+1)/2);%计算中心fori=1:nforj=1:ne(i,j)=exp(-((i-n1)^2+(j-n1)^2)/(4*alf3))/(4*pi*alf3);%高斯函数endendnr3=imfilter(mr,e,'conv','replicate');ng3=imfilter(mg,e,'conv','replicate');nb3=imfilter(mb,e,'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)