Retinex算法Matlab代码——单尺度SSR

一句话简单概括单尺度SSR的思想：从原始图像中去除照明的影响来还原图像的本质特征，照明由高斯滤波与原图卷积得到。

代码如下：

clear all;close all;clc;
I = imread('image23.jpg');%读图
 
R = I(:, :, 1);
[N1, M1] = size(R);
R0 = double(R);
Rlog = log(R0+1);
Rfft2 = fft2(R0);
 
sigma = 80;%高斯核参数
F = fspecial('gaussian', [N1,M1], sigma);
Efft = fft2(double(F));
 
DR0 = Rfft2.* Efft;%卷积
DR = ifft2(DR0);%反变换到空域
 
DRlog = log(DR +1);
Rr = Rlog - DRlog;

EXPRr = Rr;
MIN = min(min(EXPRr));
MAX = max(max(EXPRr));
EXPRr = 255*(EXPRr - MIN)/(MAX - MIN);%线性拉伸

 
G = I(:, :, 2);
 
G0 = double(G);
Glog = log(G0+1);
Gfft2 = fft2(G0);
 
DG0 = Gfft2.* Efft;
DG = ifft2(DG0);
 
DGlog = log(DG +1);
Gg = Glog - DGlog;

EXPGg = Gg;
MIN = min(min(EXPGg));
MAX = max(max(EXPGg));
EXPGg = 255*(EXPGg - MIN)/(MAX - MIN);

 
B = I(:, :, 3);
 
B0 = double(B);
Blog = log(B0+1);
Bfft2 = fft2(B0);
 
DB0 = Bfft2.* Efft;
DB = ifft2(DB0);
 
DBlog = log(DB+1);
Bb = Blog - DBlog;

EXPBb = Bb;
MIN = min(min(EXPBb));
MAX = max(max(EXPBb));
EXPBb = 255*(EXPBb - MIN)/(MAX - MIN);

 
result = cat(3, EXPRr, EXPGg, EXPBb);
subplot(121), imshow(I);
subplot(122), imshow(uint8(result));

效果展示如图1所示：

图1 不同高斯核参数下的效果