Retinex图像增强方法

Halcom

Retinex图像增强方法： 百度网盘链接：http://pan.baidu.com/s/1pL9H0QR 具体链接在halcom.cn论坛，联系人QQ：3283892722 该论坛是一个学习交流平台，我会逐一的和大家分享学习。 欢迎大家录制视频，并提交给我，我来设置视频，你可在论坛进行打赏分享。 视频专用播放器：http://halcom.cn/forum.php?mod=viewthread&tid=258&extra=page%3D1

可参考链接：http://blog.csdn.net/carson2005/article/details/9502053
http://www.cnblogs.com/Imageshop/p/3810402.html
MATLAB retinex：http://halcom.cn/forum.php?mod=v ... ge=1&extra=#pid1127

1. read_image (Image, 'C:/Users/ysw/Desktop/视频讲解/图库/Tulips.jpg')
   
2. count_channels(Image, Channels)
   
3. if(Channels==1)
   
4.     return()
   
5. else
   
6.     decompose3(Image,R,G,B)
   
7.     c :=500
   
8.     get_image_size(R,Width, Height)
   
9.     tuple_gen_sequence(0, Width*Height-1, 1, Sequence)
   
10.     tuple_mod(Sequence,Width, CF)
    
11.     tuple_div(Sequence,Width,RF)
    
12.     Fnok := exp( - ( (RF-Height/2)*(RF-Height/2)+ (CF-Width/2)*(CF-Width/2) )*1.0/c/c )
    
13.     K:=1.0/sum(Fnok)
    
14.     F := K*Fnok
    
15.    
    
16.     convert_image_type(R,IR,'real')
    
17.     convert_image_type(G,IG,'real')   
    
18.     convert_image_type(B,IB,'real')
    
19.    
    
20.     gen_image_const(Filters, 'real', Width, Height)
    
21.     set_grayval(Filters, RF,CF,F)
    
22.     *fft_generic(Filters, FF,'to_freq', -1, 'sqrt', 'dc_center', 'complex')
    
23.     fft_image(Filters,FF)
    
24.    
    
25.     fft_generic(IR, IFR, 'to_freq', -1, 'sqrt', 'dc_center', 'complex')
    
26.     mult_image(IFR,FF,IFR,1.0,0)
    
27.     fft_generic(IFR,IFR,'from_freq',1,'sqrt','dc_center','real')
    
28.    
    
29.     fft_generic(IG, IFG, 'to_freq', -1, 'sqrt', 'dc_center', 'complex')
    
30.     mult_image(IFG,FF,IFG,1.0,0)
    
31.     fft_generic(IFG,IFG,'from_freq',1,'sqrt','dc_center','real')
    
32.    
    
33.     fft_generic(IB, IFB, 'to_freq', -1, 'sqrt', 'dc_center', 'complex')
    
34.     mult_image(IFB,FF,IFB,1.0,0)
    
35.     fft_generic(IFB,IFB,'from_freq',1,'sqrt','dc_center','real')
    
36.    
    
37.     log_image(IR, LIR, 10)
    
38.     log_image(IG, LIG, 10)
    
39.     log_image(IB, LIB, 10)
    
40.    
    
41.     log_image(IFR, LIFR, 10)
    
42.     log_image(IFG, LIFG, 10)
    
43.     log_image(IFB, LIFB, 10)
    
44.    
    
45.    sub_image(LIR, LIFR, RR,  1.0, 0)
    
46.    sub_image(LIG, LIFG, GG,  1.0, 0)
    
47.    sub_image(LIB, LIFB, BB,  1.0, 0)
    
48.    
    
49.    scale_image_max(RR,RR)
    
50.    scale_image_max(GG,GG)
    
51.    scale_image_max(BB,BB)
    
52.    
    
53.    compose3(RR,GG,BB,ImageRetinex)
    
54. endif
    
55. 
    
56. 
    

复制代码