4 Pages
4 Pages
Usually there are differences between a recorded color image and the direct
observation of the same scene. One of the big differences is the fact that the human
visual system is able to distinguish details and vivid colors in shadows and in scenes
that contain illuminant shifts.
In this project we examine the performance of an algorithm called A Multiscale
Retinex with Color Restoration (MSRCR), presented in [1]. This algorithm tries to
imitate human visual “computation” while observing scenes that contains lighting
variations. MSRCR is an extension of a former algorithm called Single Scale
center/surround Retinex (SSR) [2]. MSRCR achieves simultaneous dynamic range
compression, color consistency and lightness rendition. In this project we look for the
optimal value of each parameter of the SSR and the MSRCR. The MSRCR algorithm
was compared to other two, well known, image enhancement methods – gamma
correction and histogram equalization. Results indicated that the performance of the
Multiscale Retinex is consistently good, while the performance of the others is quite
variable.
Basic Approach
The general mathematical formulation of the center/surround Retinex is
R̂ i ( x , y) = log I i ( x , y) − log[ F( x , y) ∗ I i ( x , y)]
where R̂ i (x, y) denotes the Retinex output, Ii(x,y) the image distribution in the i'th
color spectral band, “*” the convolution operation, and F(x,y) the surround function
−r 2 / c 2
F( x , y) = Ke
where c is the Gaussian surround space constant and K is selected such that
∫∫F( x, y)dxdy =1
The MSR output is then simply a weighted sum of the outputs of several different
SSR output.
N
R̂ MSR i = ∑ w n R̂ n ,i
n =1
where N is the number of scales, R̂ n ,i the i'th component of the n’th scale, R̂ MSR the i
i'th spectral component of the MSR output, and wn the weight associated with the
n’th scale. The only difference between R̂(x, y) and R̂ n (x, y) is that the surround
function is now given by
2
/ c 2n
Fn ( x , y) = K n e − r
The color restoration method for the MSR is given by
S
C i ( x , y) = log[ αI i ( x , y)] − log ∑ I i ( x , y)
i =1
where α is a constant parameter of the color restoration function. The MSRCR is
given by
R̂ MSRCRi ( x , y) = C i ( x , y) R̂ MSRi ( x , y)
2. Daniel J. Jobson, Zia-ur Rahman, Glenn A. Woodell; “Properties and Performance of a Center/ Surround Retinex”, IEEE
Transaction on Image Processing, vol. 6, no. 3, March 1997.
Project Goal
In this project we examine several aspects concerning the SSR:
• the form of the surround function
• the placement of the log function
• the scale constant for the surround function
• the treatment of the Retinex result before display
Tools
The algorithm was realized as a MATLAB program. In order to ease the research
and to make it more convenient to get the visual results, we developed a graphic user
interface (GUI), see Figure 1. The GUI allows the user to change the MSRCR
constants before running the algorithm and to compare the resulting performance to
the performance of a standard brightness correction (gamma correction).
Results
Acknowledgments
We would like to thank our supervisor Dr. Doron Shaked for his support and
guidance throughout this project. Also we would like to thank Johanan Erez for his
attention and help during the whole project. This project was supported by the
Ollendorff Research Center Fund.