Retinex--matlab实现RGB图像MSRCR

clc clear all I imread dark jpg Ir I 1 R G B Ig I 2 Ib I 3 Ir double double Ir RGB double Ig double double Ig Ib double double Ib m n size Ir double m n F1 F2 F3 L1 L2 L3 Rr Rg Rb RGB u s mi ma low high if m n h m else h n end sigma1 5 double uint8 h 0 03 3 sigma2 40 double uint8 h 0 13 13 sigma3 240 double uint8 h 0 40 40 aa 125 F1 fspecial gaussian 3 sigma1 3 sigma1 sigma1 F2 fspecial gaussian 3 sigma2 3 sigma2 sigma2 F3 fspecial gaussian 3 sigma3 3 sigma3 sigma3 R L1 imfilter Ir double F1 replicate conv R L2 imfilter Ir double F2 replicate conv L3 imfilter Ir double F3 replicate conv for i 1 m for j 1 n C i j Ir double i j Ir double i j Ig double i j Ib do uble i j end end for i 1 m for j 1 n G i j 1 3 log Ir double i j 1 log L1 i j 1 G i j 1 3 log Ir double i j 1 log L2 i j 1 G i j G i j 1 3 log Ir double i j 1 log L3 i j 1 G i j G i j exp G i j end end u mean2 G s std2 G mi min min G ma max max G Rr G mi 255 ma mi G L1 imfilter Ig double F1 replicate conv G L2 imfilter Ig double F2 replicate conv L3 imfilter Ig double F3 replicate conv for i 1 m for j 1 n C i j Ig double i j Ir double i j Ig double i j Ib do uble i j end end for i 1 m for j 1 n G i j 1 3 log Ig double i j 1 log L1 i j 1 G i j 1 3 log Ig double i j 1 log L2 i j 1 G i j G i j 1 3 log Ig double i j 1 log L3 i j 1 G i j G i j exp G i j end end u mean2 G s std2 G mi min min G ma max max G Rg G mi 255 ma mi B L1 imfilter Ib double F1 replicate conv B L2 imfilter Ib double F2 replicate conv L3 imfilter Ib double F3 replicate conv for i 1 m for j 1 n C i j Ib double i j Ir double i j Ig double i j Ib do uble i j end end for i 1 m for j 1 n G i j 1 3 log Ib double i j 1 log L1 i j 1 G i j 1 3 log Ib double i j 1 log L2 i j 1 G i j G i j 1 3 log Ib double i j 1 log L3 i j 1 G i j G i j exp G i j end end u mean2 G s std2 G mi min min G ma max max G Rb G mi 255 ma mi Rr 01