retinex相关代码汇总

混合方法 SSR.m

matlab代码，本来是RGB，改成了处理灰度图像的。

%%%%%%%%%%%%%%%RGB normalisation%%%%%%%%%%%%%%%%%%%%%%
%its cascaded implementain of  section of paper "A FAST SKIN REGION DETECTOR" by
%Phil Chen, Dr.Christos Greecos
%and
%section 2.1 of paper "Simple and accurate face detection in color images" by
%YUI TING PAI et al
% Coding by Madhava.S.Bhat, Dept of Electronics an communication
%Dr.Ambedkar Institute of Technology, Bangalore
%madhava.s@dr-ait.org
function[]= SSR(IMG)
IMG=rgb2gray(IMG);
figure,imshow(IMG)
title('original')
R=double(IMG);
[H,W]=size(R);
% minR=;
% minG=;
% minB=;
% [srow,scol]=find(R== & G== & B==);
% if(isempty(srow) && isempty(scol))
    minR=min(min(R));
% end
R=R-minR;  

S=zeros(H,W);
[srow,scol]=find(R==);
[sm,sn]=size(srow);  

for i=:sm
     S(srow(i),scol(i))=;
end
mstd=sum(sum(S))
Nstd=(H*W)-mstd;  

Cst=;
Cst=double(Cst);
for i=:H
    for j=:W
         a=R(i,j);
         b=R(i,j);
         Cst=a+b+Cst;
    end
end
%%%%sum of black pixels%%%%%%%%%%%
blacksumR=;
for i=:sm
    blacksumR=blacksumR+R(srow(i),scol(i));
end
Cstd = Cst/(*Nstd)
CavgR=sum(sum(R))./(H*W)
Rsc=Cstd./CavgR
R=R.*Rsc;
C=R;
C=C/;
Y=C;
figure,imshow(C)
title('aft 1st stage of compensation')
%normalize Y
minY=min(min(Y));
maxY=max(max(Y));
Y=255.0*(Y-minY)./(maxY-minY);
YEye=Y;
Yavg=sum(sum(Y))/(W*H)  

T=;
if (Yavg<)
    T=1.4
elseif (Yavg>)
    T=0.6
end
T
if (T~=)
    RI=R.^T;
else
    RI=R;
end  

Cfinal=uint8(RI);
figure,imshow(Cfinal)
title('Light intensity compensated')
  

clc;
clear all;
%%%%%%sourceImage%%%%%%%%%
sourceImage = imread('F:\Users\*****\Desktop\**图像\test.bmp');
%imshow(sourceImage); 

SSR(sourceImage);

效果不对，不知道哪有问题：

看灰度图的处理效果还可以。

【转载自】

图像增强算法四种，图示与源码，包括retinex（ssr、msr、msrcr）和一种混合算法 - onezeros的专栏 - CSDN博客 https://blog.csdn.net/onezeros/article/details/6342661

【其他】

SSR matlab实现单尺度retinex算法，实现图像增强。 联合开发网 - pudn.com http://www.pudn.com/Download/item/id/2776026.html

%retinex%
clc;
clear all;
%%%%%%sourceImage%%%%%%%%%
sourceImage = imread('C:\Users\huang\Desktop\����1.bmp');
imshow(sourceImage); 

%%%R,G,B%%%%
image_r=sourceImage(:,:,);
image_g=sourceImage(:,:,);
image_b=sourceImage(:,:,); 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

r=;
cParam=;
block=*r+;
sum=0.0;
for i=:block
    for j=:block
        gauss(i,j)=exp((-((i-r-)^+(j-r-)^))/cParam^);
        sum=sum+gauss(i,j);
    end
end
for i=:block
    for j=:block
        gauss(i,j)=gauss(i,j)/sum;
    end
end 

%%%%%%%%%%%R%%%%%%%%%%%%%
Img = double(image_r);
K=imfilter(Img,gauss,'replicate','conv');
G=log(Img+)-log(K+);  

mi=min(min(G));
ma=max(max(G));
destImage_r=(G-mi)*/(ma-mi);  

%%%%%%%%%%%%%%G%%%%%%%%%%%%%
Img = double(image_g);
K=imfilter(Img,gauss,'replicate','conv');
G=log(Img+)-log(K+); 

mi=min(min(G));
ma=max(max(G));
destImage_g=(G-mi)*/(ma-mi);
%%%%%%%%%%%%%B%%%%%%%%%%%%%%%%%%
Img = double(image_b);
K=imfilter(Img,gauss,'replicate','conv');
G=log(Img+)-log(K+); 

mi=min(min(G));
ma=max(max(G));
destImage_b=(G-mi)*/(ma-mi);
%%%%%%%%%%%%%%%%%%
destImageSSR=cat(,destImage_r,destImage_g,destImage_b);
figure;
imshow(uint8(destImageSSR)); 

%etinex
clc
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
r=;
cParam1=;
cParam2=;
cParam3=;
block=*r+;
sum1=0.0;
sum2=0.0;
sum3=0.0;
for i=:block
    for j=:block
        gauss1(i,j)=exp((-((i-r-)^+(j-r-)^))/cParam1^);
        sum1=sum1+gauss1(i,j);
        gauss2(i,j)=exp((-((i-r-)^+(j-r-)^))/cParam2^);
        sum2=sum2+gauss2(i,j);
        gauss3(i,j)=exp((-((i-r-)^+(j-r-)^))/cParam3^);
        sum3=sum3+gauss3(i,j);
    end
end
for i=:block
    for j=: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+)-(/)*(log(K1+)+log(K2+)+log(K3+));  

mi=min(min(G));
ma=max(max(G));
destImage_r=(G-mi)*/(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+)-(/)*(log(K1+)+log(K2+)+log(K3+));  

mi=min(min(G));
ma=max(max(G));
destImage_g=(G-mi)*/(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+)-(/)*(log(K1+)+log(K2+)+log(K3+));  

mi=min(min(G));
ma=max(max(G));
destImage_b=(G-mi)*/(ma-mi); 

%%%%%%%%%%%%%%%%%%
destImageMSR=cat(,destImage_r,destImage_g,destImage_b);
figure;
imshow(uint8(destImageMSR));
 

【别人博客里的一些效果图】

retinex图像增强算法的研究 - zmshy2128 - 博客园 https://www.cnblogs.com/zmshy2128/p/6126987.html

关于Retinex图像增强算法的一些新学习。 - Imageshop - 博客园 https://www.cnblogs.com/Imageshop/p/3810402.html

Retinex图像增强算法代码 - liyingjiang22的博客 - CSDN博客 https://blog.csdn.net/liyingjiang22/article/details/52903127

图像增强处理之：同态滤波与Retinex算法(二)McCann Retinex和McCann99 Retinex迭代算法 - 满城风絮 - CSDN博客 https://blog.csdn.net/piaoxuezhong/article/details/78248219