retinex图像增强算法

Halcom

retinex图像增强算法:MATLAB代码如下：
单retinex图像增强算法---多retinex图像增强算法

1. clc,clear,close all
   
2. % single scale retinex
   
3. % I = imread('C:\Users\ysw\Desktop\视频讲解\图库\Tulips.jpg');
   
4. I = imread('C:\Users\ysw\Desktop\视频讲解\图库\shot.jpg');
   
5. [m,n,dim] = size(I);
   
6. c = 5e1;  % 可变参
   
7. [Y,X]=meshgrid(1:n,1:m);
   
8. 
   
9. %Fnok = exp(-((X.^2)+(Y.^2))./(c.^2));
   
10. Fnok = exp(-(((X-m/2).^2)+((Y-n/2).^2))./(c.^2));
    
11. K = 1/(sum(sum(Fnok)));
    
12. F = K.*Fnok;
    
13. 
    
14. IR = double(I(:,:,1));
    
15. IG = double(I(:,:,2));
    
16. IB = double(I(:,:,3));
    
17. FF  = fftshift(fft2(F));
    
18. IFR = fftshift(fft2(IR)); IFR=FF.*IFR; IFR=real(ifft2(ifftshift(IFR)));
    
19. IFG = fftshift(fft2(IG)); IFG=FF.*IFG; IFG=real(ifft2(ifftshift(IFG)));
    
20. IFB = fftshift(fft2(IB)); IFB=FF.*IFB; IFB=real(ifft2(ifftshift(IFB)));
    
21. 
    
22. RR = log10(double(IR))-log10(IFR);
    
23. RG = log10(double(IG))-log10(IFG);
    
24. RB = log10(double(IB))-log10(IFB);
    
25. 
    
26. OUT(:,:,1)=uint8(255*RR/(max(max([RR RG RB]))));
    
27. OUT(:,:,2)=uint8(255*RG/(max(max([RR RG RB]))));
    
28. OUT(:,:,3)=uint8(255*RB/(max(max([RR RG RB]))));
    
29. figure,
    
30. subplot(121),imshow(I,[])
    
31. subplot(122),imshow(OUT,[])
    
32. 
    
33. % multi scale retinex
    
34. % c = 45;  c = 100;  c=200;
    
35. % 加权和

复制代码

参考链接：http://zoi.utia.cas.cz/research