自己写的Retinex算法

S=imread(strcat('F:\人脸库\Yale_5G\face1\t\1_s7.bmp')); [m,n]=size(S); S=double(S); size=m*n; mean=sum(sum(S),2)/size; R=zeros(m,n); R=repmat(mean,m,n); %做高斯函数 c=250; for i=1:m for j=1:n N(i,j)=exp(-(i^2+j^2)/c^2); end end K=1/sum(N(:)); N1=K*N; G=conv2(S,N1,'same'); %计算SSR SS=zeros(m,n); for i=1:m for j=1:n SS(i,j)=log(S(i,j)); GG(i,j)=log(G(i,j)); V(i,j)= SS(i,j)-GG(i,j); end end Min=min(min(V)); Max=max(max(V)); for i=1:m for j=1:n VV(i,j)=(V(i,j)-Min)*255/(Max-Min); end end VV=uint8(VV); S=uint8(S); figure,subplot(1,2,1),imshow(S); subplot(1,2,2),imshow(VV);