%retinex%
clc;
clear all;
%%%%%%sourceImage%%%%%%%%%
sourceImage = imread('C:\Users\huang\Desktop\����1.bmp');
imshow(sourceImage); %%%R,G,B%%%%
image_r=sourceImage(:,:,1);
image_g=sourceImage(:,:,2);
image_b=sourceImage(:,:,3); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% r=100;
cParam=200;
block=2*r+1;
sum=0.0;
for i=1:block for j=1:block gauss(i,j)=exp((-((i-r-1)^2+(j-r-1)^2))/cParam^2); sum=sum+gauss(i,j); end
end
for i=1:block for j=1:block gauss(i,j)=gauss(i,j)/sum; end
end %%%%%%%%%%%R%%%%%%%%%%%%%
Img = double(image_r);
K=imfilter(Img,gauss,'replicate','conv');
G=log(Img+1)-log(K+1);  mi=min(min(G));
ma=max(max(G));
destImage_r=(G-mi)*255/(ma-mi);  %%%%%%%%%%%%%%G%%%%%%%%%%%%%
Img = double(image_g);
K=imfilter(Img,gauss,'replicate','conv');
G=log(Img+1)-log(K+1); mi=min(min(G));
ma=max(max(G));
destImage_g=(G-mi)*255/(ma-mi);
%%%%%%%%%%%%%B%%%%%%%%%%%%%%%%%%
Img = double(image_b);
K=imfilter(Img,gauss,'replicate','conv');
G=log(Img+1)-log(K+1); mi=min(min(G));
ma=max(max(G));
destImage_b=(G-mi)*255/(ma-mi);
%%%%%%%%%%%%%%%%%%
destImageSSR=cat(3,destImage_r,destImage_g,destImage_b);
figure;
imshow(uint8(destImageSSR)); %etinex
clc
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
r=100;
cParam1=30;
cParam2=100;
cParam3=300;
block=2*r+1;
sum1=0.0;
sum2=0.0;
sum3=0.0;
for i=1:block for j=1:block gauss1(i,j)=exp((-((i-r-1)^2+(j-r-1)^2))/cParam1^2); sum1=sum1+gauss1(i,j); gauss2(i,j)=exp((-((i-r-1)^2+(j-r-1)^2))/cParam2^2); sum2=sum2+gauss2(i,j); gauss3(i,j)=exp((-((i-r-1)^2+(j-r-1)^2))/cParam3^2); sum3=sum3+gauss3(i,j); end
end
for i=1:block for j=1:block gauss1(i,j)=gauss1(i,j)/sum1; gauss2(i,j)=gauss2(i,j)/sum2; gauss3(i,j)=gauss3(i,j)/sum3; end
end %%%%%%%%%%%R%%%%%%%%%%%%%
Img = double(image_r);
K1=imfilter(Img,gauss1,'replicate','conv');
K2=imfilter(Img,gauss2,'replicate','conv');
K3=imfilter(Img,gauss3,'replicate','conv'); G=log(Img+1)-(1/3)*(log(K1+1)+log(K2+1)+log(K3+1));  mi=min(min(G));
ma=max(max(G));
destImage_r=(G-mi)*255/(ma-mi);  %%%%%%%%%%%%%%G%%%%%%%%%%%%%
Img = double(image_g);
K1=imfilter(Img,gauss1,'replicate','conv');
K2=imfilter(Img,gauss2,'replicate','conv');
K3=imfilter(Img,gauss3,'replicate','conv'); G=log(Img+1)-(1/3)*(log(K1+1)+log(K2+1)+log(K3+1));  mi=min(min(G));
ma=max(max(G));
destImage_g=(G-mi)*255/(ma-mi);
%%%%%%%%%%%%%B%%%%%%%%%%%%%%%%%%
Img = double(image_b);
K1=imfilter(Img,gauss1,'replicate','conv');
K2=imfilter(Img,gauss2,'replicate','conv');
K3=imfilter(Img,gauss3,'replicate','conv'); G=log(Img+1)-(1/3)*(log(K1+1)+log(K2+1)+log(K3+1));  mi=min(min(G));
ma=max(max(G));
destImage_b=(G-mi)*255/(ma-mi); %%%%%%%%%%%%%%%%%%
destImageMSR=cat(3,destImage_r,destImage_g,destImage_b);
figure;
imshow(uint8(destImageMSR));