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

巴特西

retinex相关代码汇总

最新文章

热门文章