Retinex matlab实现RGB图像MSRCR.doc

%clc;clear all;I=imread(dark.jpg);%Ir=I(:,:,1); %R,GBIg=I(:,:,2);Ib=I(:,:,3);Ir_double=double(Ir);%RGBdoubleIg_double=double(Ig);Ib_double=double(Ib);m,n=size(Ir_double);%,mn%F1F2F3%L1L2L3%Rr,Rg,RbRGB%u,s,mi,ma%low,highif mnh=m;elseh=n;endsigma1=5;%double(uint8(h*0.03);%3%sigma2=40;%double(uint8(h*0.13);%13%sigma3=240;%double(uint8(h*0.40);%40%aa=125;%F1=fspecial(gaussian,3*sigma1,3*sigma1,sigma1);%F2=fspecial(gaussian,3*sigma2,3*sigma2,sigma2);%F3=fspecial(gaussian,3*sigma3,3*sigma3,sigma3);%R%L1=imfilter(Ir_double,F1,replicate,conv);%RL2=imfilter(Ir_double,F2,replicate,conv);L3=imfilter(Ir_double,F3,replicate,conv);for i=1:mfor j=1:nC(i,j)=(Ir_double(i,j)/(Ir_double(i,j)+Ig_double(i,j)+Ib_double(i,j);%endendfor i=1:mfor j=1:nG(i,j)=1/3*(log(Ir_double(i,j)+1)-log(L1(i,j)+1);%G(i,j)=1/3*(log(Ir_double(i,j)+1)-log(L2(i,j)+1)+G(i,j);G(i,j)=(1/3*(log(Ir_double(i,j)+1)-log(L3(i,j)+1)+G(i,j);%G(i,j)=exp(G(i,j);%endendu=mean2(G);%,s=std2(G);mi=min(min(G);ma=max(max(G);Rr=(G-mi)*255/(ma-mi);%G%L1=imfilter(Ig_double,F1,replicate,conv);%GL2=imfilter(Ig_double,F2,replicate,conv);L3=imfilter(Ig_double,F3,replicate,conv);for i=1:mfor j=1:nC(i,j)=(Ig_double(i,j)/(Ir_double(i,j)+Ig_double(i,j)+Ib_double(i,j);%endendfor i=1:mfor j=1:nG(i,j)=1/3*(log(Ig_double(i,j)+1)-log(L1(i,j)+1);%G(i,j)=1/3*(log(Ig_double(i,j)+1)-log(L2(i,j)+1)+G(i,j);G(i,j)=(1/3*(log(Ig_double(i,j)+1)-log(L3(i,j)+1)+G(i,j);%G(i,j)=exp(G(i,j);%endendu=mean2(G);%,s=std2(G);mi=min(min(G);ma=max(max(G);Rg=(G-mi)*255/(ma-mi);%B%L1=imfilter(Ib_double,F1,replicate,conv);%BL2=imfilter(Ib_double,F2,replicate,conv);L3=imfilter(Ib_double,F3,replicate,conv);for i=1:mfor j=1:nC(i,j)=(Ib_double(i,j)/(Ir_double(i,j)+Ig_double(i,j)+Ib_double(i,j);%endendfor i=1:mfor j=1:nG(i,j)=1/3*(log(Ib_double(i,j)+1)-log(L1(i,j)+1);%G(i,j)=1/3*(log(Ib_double(i,j)+1)-log(L2(i,j)+1)+G(i,j);G(i,j)=(1/3*(log(Ib_double(i,j)+1)-log(L3(i,j)+1)+G(i,j);%G(i,j)=exp(G(i,j);%endendu=mean2(G);%,s=std2(G);mi=min(min(G);ma=max(max(G);Rb=(G-mi)*255/(ma-mi);%Rr_01=Rr/255;Rg_01=Rg/255;R