Next Article in Journal
A Power-Efficient Clustering Protocol for Coal Mine Face Monitoring with Wireless Sensor Networks Under Channel Fading Conditions
Previous Article in Journal
Estimating Forest Aboveground Biomass by Combining Optical and SAR Data: A Case Study in Genhe, Inner Mongolia, China
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Realistic Image Rendition Using a Variable Exponent Functional Model for Retinex

1
Key Laboratory of Photoelectronic Imaging Technology and System, Ministry of Education of China, Beijing Institute of Technology, Beijing 100081, China
2
School of Science, Communication University of China, Beijing 100024, China
*
Author to whom correspondence should be addressed.
Sensors 2016, 16(6), 832; https://doi.org/10.3390/s16060832
Submission received: 17 February 2016 / Revised: 11 May 2016 / Accepted: 16 May 2016 / Published: 7 June 2016
(This article belongs to the Section Physical Sensors)

Abstract

:
The goal of realistic image rendition is to recover the acquired image under imperfect illuminant conditions, where non–uniform illumination may degrade image quality with high contrast and low SNR. In this paper, the assumption regarding illumination is modified and a variable exponent functional model for Retinex is proposed to remove non–uniform illumination and reduce halo artifacts. The theoretical derivation is provided and experimental results are presented to illustrate the effectiveness of the proposed model.

1. Introduction

Realistic image rendition aims to represent human perception of natural scenes. The meaning of “realistic” is to provide machine vision with ideal images according to the human visual system. A complete visual pathway includes the optic nerve, retina, optic tract, optic chiasm, superior colliculus, lateral geniculate nucleus, optic radiation, and cortex, as shown diagrammatically in Figure 1 [1]. The main features of realistic image rendition include color constancy, image enhancement, high dynamic range compression, etc. The physiological basis for color constancy involves specialized neurons in the primary visual cortex that compute local ratios of cone activity [2], which is the same calculation as Land’s Retinex algorithm [3,4] used to achieve color constancy. The existence of these specialized cells, double–opponent cells, has been proven using receptive field mapping. Receptive field [5,6] is the basic unit of visual information processing, and can be separated into two types: On–Center and Off–Center ganglion cells. Figure 2 shows the receptive field in the retina.
Algorithms of realistic image rendition based on visual characteristics generally include Retinex algorithms for color constancy. The word “Retinex” is a combination of “retina” and “cortex”. The aim of Retinex theory is to tell whether human eyes can determine reflectance when both the illumination and reflectance are unknown. Land and McCann [6] first proposed the Retinex theory, a path-based algorithm, as a model of color perception of the human visual system (HSV). Many algorithms [7,8,9] are based on this approach, which differ in how the path is selected. However, these methods have high computation complexity and require numerous parameters. McCann [10,11,12] replaced the path calculation by a recursive matrix computation which greatly improved computational efficiency. However, the terminal criterion is not clear and can strongly influence the result. In PDE based models [13], the Retinex principles are often translated into a physical form. These algorithms are developed based on solving a Poisson equation which can yield fast and exact implementation using only two fast Fourier transforms. The main assumption in this algorithm type is that the reflectance performs as the sharp details in the image, while illumination varies smoothly. Based on the assumptions used in PDE formulations, Kimmel et al. [14] proposed a general variational model for the Retinex problem that unified previous methods. Ma and Osher [15,16] proposed a total variation and nonlocal total variation(TV) regularized model using the same assumptions. Ng el at. [17] investigated the TV model with more constraints. Recently, Liang and Zhang [18] established a new higher order total variation L1 decomposition model (HoTVL1) which can correct the piecewise linear shadows. Zosso [19,20] proposed a unifying Retinex model based on non-local differential operators.
To the best of our knowledge, almost all of the important assumptions about illumination in existing Retinex models require spatial smoothness. However, many images with non-uniform illuminations have non smooth illumination, actually. In this paper, we assume:
(a)
The reflecting object a Lambertian reflector and reflectance corresponds to sharp details in the image;
(b)
Illumination is smooth in most regions, but may contain non-smooth part(s).
Based on these assumptions, we propose a new Retinex model using a variable exponent functional. We assume that the illumination function belongs to some Sobolev space with variable exponents. The proposed model solution existence is proved here. Although the proposed model is developed for specific cases, it can also be applied to general degraded images and significantly reduces the halo artifact.
In Section 2, we argue the reasonability of the assumption and present the proposed model. We also present a proof of solution existence for the proposed model, and introduce an efficient iterative solution method. In Section 3, we present several numerical examples to demonstrate the effectiveness of the proposed model. Concluding remarks are presented in Section 4.

2. New Assumption and Proposed Model

2.1. New Assumption

To illustrate the proposed assumption, let us consider the images with different illumination conditions and their corresponding surfaces in Figure 3. The corresponding surfaces illustrate the shadow shapes. The illumination of the text image varies smoothly, whereas that of the book image has an apparent non-smooth component. Figure 4 shows a single row from the two images, illustrating the text image curve changes relatively smooth in the shadow area, while that of the book image changes dramatically at the edge of the shadow and relatively smooth in the interior of the shadow.
The above examples support the proposed assumption. Indeed, every severe non-uniform illumination case is likely to have a non-smooth part. Our aim is to extract illumination images and recover realistic images.

2.2. Proposed Model

First, we introduce the variable exponent functional and some related models. Blomgren et al. [21] proposed the variable exponent functional for image denoising problems. They tried to minimize:
E ( u ) = Ω | u | p ( | u | ) d x
where u is the image function and p is a monotonically decreasing function with lim s 0 p ( s ) = 2 , lims→∞p(s) = 1. Choosing p = 1 produces the widely used Rudin-Osher-Fatemi (ROF) model [22] which preserves edge sharpness, but often causes the “staircase” effect; p = 2 produces isotropic diffusion, which avoids the “staircase” effect but smears edges. Thus, it is natural to combine their benefits with a variable exponent. However, because p relies on u , it is difficult to establish the lower semi continuity property of the functional. Chen et al. [23] proposed a variable exponent linear growth functional model for image denoising, enhancement and restoration, which is extended by Li et al. [24], using variable exponent functionals in image restoration.
For simplicity, we formulate and discuss our model based on grayscale images. For color images, we simply map the color into HSV(hue, saturation, value) color space, process only the V channel, then transform it back to the RGB domain. This method is called HSV Retinex [14,17].
Let I be an image defined in image domain Ω . The primary goal of Retinex theory is to decompose I into the reflectance image, R, and the illumination image, L, as shown in Figure 5, such that, at each point in the image domain [25]:
I = R L
and following [14,17], we may further assume that:
L I > 0
We first convert Equation (2) into the logarithmic domain:
i = log ( I ) ,   l = log ( L ) ,   r = log ( R )
so that:
i = l + r
Based on our new assumption, the illumination image may contain non-smooth parts. Weuse a total variation like regularizer near non-smooth parts and a Tikhonov like regularizer for smooth parts. We minimize the objective function as follows:
E ( l ) = Ω | l i | 2 d x + λ Ω | l | p ( x ) d x
where λ is a positive number, and p ( x ) = 1 + 1 1 + w | d | 2 , where d is the ideal illumination image, discussed in Section 2.4.2. The first fidelity term on the right side of model (3) measures the similarity of the gradient between the illumination and the original image, and the second is the regularization term. Clearly, p 1 near the edges of d where the gradient is large, and so the regularizer is similar to a TV regularizer which can preserve edges; p 2 in the homogeneous regions where the gradient is small, and here the regularizer is similar to a Tikhonov like regularizer, which is superior to total variation. In other regions, the penalty is adjusted by p(x).
The classical Retinex algorithm uses a Gaussian filter, equivalent to a Tikhonov regularizer, to obtain the illumination image. However, a Gaussian filter smears edges, which is the main cause of halo artifacts [26]. Using the adaptive TV like regularizer for the high contrast edges in the image, our model not only prevents halo artifacts but also extracts the edges of non-uniform illumination from the image.

2.3. Solution Existence

Let us recall some definitions and basic properties of variable exponent Lebesgue and Sobolev spaces, following [24,27].
Definition 
(variable exponent spaces): Let Ω be a bounded open set with Lipchitz boundary and p ( x ) : Ω [ 1 , + ) be a measurable function, with the family of all measurable functions on Ω being P ( Ω ) . We define a functional, which is also called modular:
Q p ( x ) ( u ) = Ω | u | p ( x ) d x
and a norm:
| | u | | p ( x ) = inf { λ > 0 : Q p ( x ) ( u / λ ) 1 }
Then the variable exponent Lebesgue and Sobolev spaces are, respectively:
L p ( x ) ( Ω ) = { u : Ω R |   | | u | | p ( x ) < }
and:
W 1 , p ( x ) ( Ω ) = { u : Ω R |   u L p ( x ) ( Ω ) , u L p ( x ) ( Ω ) }
With the norm | | u | | 1 , p ( x ) = | | u | | p ( x ) + | | u | | p ( x ) , W 1 , p ( x ) ( Ω ) becomes a Banach space.
Lemma 1. 
(relationship between modular and norm [27]): Let Q p ( x ) be a modular on X and u X , then | | u | | p ( x ) Q p ( x ) ( u ) +1.
Lemma 2. 
(embedding theorem [24]): Let p ( x ) , q ( x ) P ( Ω ) , and p ( x ) q ( x ) for a.e. x Ω . Then L q ( x ) ( Ω ) is continuously embedded in L p ( x ) ( Ω ) .
Lemma 3. 
(convexity [27]): Let F ( l , x ) = | l | p ( x ) , with p ( x ) = 1 + 1 1 + w | d | 2 as in model (3). Then for each x , F ( ξ , x ) is convex in ξ .
Lemma 4. 
(weak lower semi continuity [27]) Let F ( ξ , x ) be bounded from below, and the map ξ F ( ξ , x ) is convex in each x Ω . Then the energy functional, I = Ω F ( l , x ) d x , is weak lower semi-continuous in W 1 , p ( x ) .
Theorem 1. 
Let Ω R 2 be a bounded open set with Lipchitz boundary, i W 1 , p ( x ) ( Ω ) L 2 ( Ω ) , then the minimization problem:
min l W 1 , p ( x ) ( Ω ) L 2 ( Ω ) { E ( l ) = Ω | l i | 2 d x + λ Ω | l | p ( x ) d x }
has a minimizer l W 1 , p ( x ) ( Ω ) L 2 ( Ω ) .
Proof of Theorem 1. 
Let { l k } k = 1 , l k W 1 , p ( x ) ( Ω ) L 2 ( Ω ) be the minimizing sequence for E ( l ) . Then:
Ω | l k i | 2 d x < M  and  Ω | l k | p ( x ) d x < M
where M denotes a universal positive constant that may differ from line to line. Hence Ω | l k | 2 d x < M . Thanks to Poincare inequality, we have Ω l k 2 d x < M , and from Lemma 2, L 2 ( Ω ) L p ( x ) ( Ω ) . Therefore, Ω | l k | p ( x ) d x < M , and together with the inequality Ω | l k | p ( x ) d x < M , we obtain Q p ( x ) ( l k ) + Q p ( x ) ( l k ) < M . This implies that { l k } k = 1 is a uniformly bounded sequence in W 1 , p ( x ) ( Ω ) due to Lemma 1, and { l k } k = 1 is also uniformly bounded in L 2 ( Ω ) . Since W 1 , p ( x ) ( Ω ) L 2 ( Ω ) is a reflexive Banach space, up to a subsequence, there exists l * W 1 , p ( x ) ( Ω ) L 2 ( Ω ) such that { l k j } converges weakly to l * in W 1 , p ( x ) ( Ω ) L 2 ( Ω ) . From Lemma 4, E ( l ) is lower semi continuous in W 1 , p ( x ) ( Ω ) L 2 ( Ω ) . Thus:
E ( l * ) lim k E ( l k ) = inf l W 1 , p ( x ) ( Ω ) L 2 ( Ω ) E ( l )
Therefore, l * is the minimum point of E ( l ) . ☐

2.4. Implementation

We formulate the basic procedure for solving problem Equation (3) following the split Bregman [28,29,30,31,32,33] technique. We solve the minimization by introducing an auxiliary variable b:
min { Ω | l i | 2 d x + λ Ω | b | p ( x ) d x }  subject to  b = l
By adding one quadratic penalty function term, we convert Equation (4) to an unconstrained splitting formulation:
min { Ω | l i | 2 d x + λ Ω | b | p ( x ) d x + γ Ω | b l | 2 d x }
where γ is a positive parameter which controls the weight of the penalty term. Similar to the split Bregman iteration, we propose the scheme:
{ ( l k + 1 , b k + 1 ) = arg min l , b Ω | l i | 2 d x + λ Ω | b | p ( x ) d x + γ Ω | b l - t k | 2 d x t k + 1 = t k + l k + 1 b k + 1
Alternatively, this joint minimization problem can be solved by decomposing into several subproblems.

2.4.1. Subproblem l with Fixed b and t

Given the fixed variable bk and tk , our aim is to find the solution of the problem:
l k + 1 = arg min l Ω | l i | 2 d x + γ Ω | b k l - t k | 2 d x
which has the optimality condition:
( γ + 1 ) Δ l = γ ( b k t k ) + Δ i
where b = ( b x , b y ) and t = ( t x , t y ) . Since the discrete system is strictly diagonally dominant with Neumann boundary condition, the most natural choice is the Gauss-Seidel method. The Gauss-Seidel solution to this subproblem can be written componentwise as:
l i , j k + 1 = γ 4 ( γ + 1 ) ( b x , i 1 , j k + b y , i , j 1 k b x , i , j k b y , i , j k + t x , i 1 , j k + t y , i , j 1 k t x , i , j k t y , i , j k ) + 1 4 ( γ + 1 ) ( 4 i i , j k i i + 1 , j k i i 1 , j k i i , j 1 k i i , j + 1 k ) + γ + 1 4 ( γ + 1 ) ( l i + 1 , j k + l i 1 , j k + l i , j + 1 k + l i , j 1 k )
Note that this subproblem can also be solved by FFT with periodic boundary condition.

2.4.2. Subproblem b with Fixed l and t

Similarly, we solve:
b k + 1 = arg min b λ Ω | b | p ( x ) d x + γ Ω | b l k + 1 t k | 2 d x
which has the optimality condition:
{ λ p ( x ) | b | p ( x ) 2 b x + 2 γ ( b x x l k + 1 t x k ) = 0 λ p ( x ) | b | p ( x ) 2 b y + 2 γ ( b y y l k + 1 t y k ) = 0
where l = ( x l , y l ) . If b x and b y are not zero, then:
b x = x l k + 1 + t x k y l k + 1 + t y k b y
Substituting Equation (11) into Equation (10):
sgn ( b y ) T b y p ( x ) 1 + 2 γ ( b y y l k + 1 t y k ) = 0
where T = λ p ( x ) ( ( x l k + 1 + t x k y l k + 1 + t y k ) 2 + 1 ) p ( x ) 2 2 . Note that:
sgn ( b x ) = sgn ( x l k + 1 + t x k )
sgn ( b y ) = sgn ( y l k + 1 + t y k )
So Equation (12) can be expressed as:
sgn ( y l k + 1 + t y k ) T b y p ( x ) 1 + 2 γ ( b y y l k + 1 t y k ) = 0
Unfortunately, we cannot obtain the explicit solution of the Equation (15). We can use the Newton method to get an approximate solution. If b y is solved, b x can be easily determined using Equations (11) and (13). The process is shown as Algorithm 1.
Algorithm 1: (Newton’s method)
Input: b y 0 = y l k + 1
 If y l k + 1 + t y k = 0
  Output b y k + 1 = 0
 Else
  while not converged
b y j + 1 = sign ( y l k + 1 + t y k ) max { b y j T ( b x j ) p ( x ) 1 + 2 γ ( b y j | y l k + 1 + t y k | ) ( p ( x ) 1 ) T ( b x j ) p ( x ) 2 + 2 γ b y j , 0 }

  End
 Output b y k + 1
 If x l k + 1 + t x k = 0
  Output b x k + 1 = 0
 Else
b x k + 1 = x l k + 1 + t x k y l k + 1 + t y k b y k + 1
 End
End
Another problem is that in practice we don't know d in p(x). We have tested two ways to approximate d. One way is to use edge preserving filter (e.g., bilateral filter) to give an approximation of d and keep the exponent fixed during the iteration; Another way is to replace d with G ( l k + 1 ) during the iteration and G ( ) represents the Gauss convolution operator, In most cases, both methods can generate similar prominent results. However, in some cases, dynamic approximation would give better results than fixed approximation because dynamic approximation can give a more accurate approximation of d along with the iteration. To illustrate this, consider the associated heat flow to problem Equation (3):
l t = 2 ( Δ l + Δ i ) + λ p ( x ) ( | l | p ( x ) 2 l ) = ( 2 + λ p ( x ) | l | p ( x ) 2 ) l T T + ( 2 + λ p ( x ) ( p ( x ) 1 ) | l | p ( x ) 2 ) l N N
where l T T and l N N are the second derivatives of l in the tangent and normal direction to the isophote lines respectively. From Equation (16), we have two critical conclusions:
  • The illumination image, l k + 1 , becomes increasingly smooth over time.
  • Diffusion speed in the tangent direction is always faster than that in the normal direction.
The first conclusion conforms to the smooth assumption of the illumination image. If the illumination image has non smooth parts, then the second conclusion guarantees that the solution can preserve these parts. Thus, l k + 1 continuously gets closer to d with calculation. However, the convergence proof of the algorithm is difficult since the exponent is changing during the iterations. If the exponent is fixed, the convergence proof can be directly obtained because the objective function is fixed and the iteration of split Bregman is monotone decreasing in the function values. The strict proof can be found in references [34,35]. If the exponent changes during the iterations, then the objective function changes as well. The convergence proof in references [34,35] cannot be applied here. However, we have tested numerous experiments and our algorithm did converges in all the tests. We leave it for further work.

2.4.3. Update:

t k + 1 = t k + l k + 1 b k + 1

2.4.4. Update l k + 1 :

l k + 1 = max { l k + 1 , i }
which corresponds to the constraint L I > 0 .
The process is shown as Algorithm 2.
Algorithm 2: (Variable Exponent Functional Retinex)
Input: Image I
Transform into log domain i = log ( I + 1 ) ;
Initialization: l 0 = i , b 0 = 0 , t 0 = 0 and k = 0
While | | l k l k - 1 | | | | l k | | ε
  (1) Given b k and t k , update l k + 1 by solving Equation (8).
  (2)Given l k + 1 and t k , update b k + 1 by using algorithm 1.
  (3)Update t k + 1 = t k + l k + 1 b k + 1 .
  (4) l k + 1 = max { l k + 1 , i } ;
  (5) k = k + 1 ;
End
Output: Image R = exp ( i ) exp ( l )

2.5. Relation to Previous Methods

Let us revisit the model in Section 2.2. If we set p ( x ) = 2 and remove the constraint li, our model is equivalent to homomorphic filtering [36]. Retaining the constraint li and fixing p(x) = 2, it is similar to a random walk, Ng’s model [17] and McCann algorithm [12]. Thus, our proposed model generalizes previous models.

3. Numerical Results

We present numerical results to demonstrate the efficiency of the proposed model and algorithm. For color images examples, we use HSV Retinex. We compare our proposed model with three state of the art methods, HoTVL1 [18], Ng’s method [17] and multiscale Retinex [37].
For all the tests, the recovered reflectance of our model is:
R = I L
where L = exp(l) is the illumination function obtained from Section 2.4, and I = exp(i) is the original image. Note that the reflectance image is sometimes over enhanced, and we add the Gamma correction illumination to the reflectance image after decomposition. The Gamma correction of L with parameter s is:
L = W ( L W ) 1 s
where W is the value of the white pixel. Parameter s was set to Section 2.2 in the tests. Thus:
I = L R
and the global framework of our proposed method is illustrated in Figure 6.

3.1. Synthetic Images

In this subsection, we set λ = 80, γ = 103 and w = 109. Simulated illumination is added to the original texture images, as shown in Figure 7, with the numerical results shown in Figure 8. The recovered image following our proposed method is visually superior. We use signal to noise ratio (SNR) to measure the similarities between the original and recovered images, as shown in Figure 9. SNR from our proposed method is significantly superior to the other methods. We further use structural similarity index (SSIM) and CIEDE2000 color difference to measure the texture similarities and perceptual difference between the original and recovered images respectively, as shown in Table 1 and Table 2. We can see from tables that our proposed method is superior to the other methods. We note that HoTVL1 failed in these tests. The main reason is that the assumption in HoTVL1 is piecewise constant and piecewise linear, which means that the shadow should be piecewise linear. However, this is not the case of these tests. The shadow part is almost a constant and also has sharp edges. Hence the result of HoTVL1 is not satisfactory.

3.2. Natural Images

For all tests, we set λ = 80, γ = 103 and w = 103. We begin with Andelson’s checkerboard shadow image, as shown in Figure 10a. Region A looks darker than region B, although they have the same values. Figure 11 shows the reconstructed illumination and reflection images using Ng’s, HoTVL1, multiscale Retinex and our proposed model. HoTVL1 and the proposed method produce superior results to Ng’s method and multiscale Retinex. The recovered illumination using our proposed method contains less reflectance information than HoTVL1, e.g., the outline of the cylinder. Our proposed method also contains less shadow information in the reflectance image than other methods. Table 3 compares the recovered intensity values of the two regions for the four methods. The contrast of the marked areas using our proposed method is superior to the other three methods.
Consider the degraded image shown in Figure 10b. Figure 12 shows the reconstruction for the four methods. Note that in this example, we adopted a Gamma correction step, as discussed above. Our proposed method has superior visual outcome. Ng’s method suffers halo artifacts, e.g., near the edges of the tower and the roof of the building, which rarely appear in HoTVL1, multiscale Retinex and our proposed method. However, many fine structures lost in the HoTVL1 and multiscale Retinex reproduced image, which are retained in our proposed method.
The next illustrative example is recovery of non-uniform degraded images. The two images in Figure 10c and Figure 3b suffer from the strong shadow areas. Figure 13 shows the comparison for the considered methods. The shadow is almost entirely removed by our proposed method, whereas the other methods retain partly shadowed regions.
In the end, we test the effect of different approximations of d. We use bilateral filter to approximate d and keep the exponent fixed during the iterations. Figure 14 shows the numerical results. We see in Figure 14 that the illumination and the reflectance of the results are not as good as those in Figure 11. This experiment supports our discussion in Section 2.4.2.

4. Conclusions

We proposed a variable exponent functional model for Retinex, proved the existence of the solution for the model and provided the theoretical derivation. The proposed method can be applied to general degraded cases as well. Experimental results validatethat our proposed method can remove non-uniform illumination and significantly reduce halo artifacts.

Acknowledgments

The authors would like to thank the reviewers of this manuscript for their helpful comments and suggestions, Stanley Osher for his introduction of Retinex during his visit to Harbin Institute of Technology in 2015, and Jingwei Liang for sharing his code and helpful discussions. This work was supported by the National High Technology Research and Development Program of China (863 Program, grant No.2014AA7026082), the Natural Science Foundation of Beijing, China (grant No. 4152045), the National Natural Science Foundation of China (grant No. 61071148) and the Science Research Project of CUC (grant No. 3132014XNG1417, 3132016XNG1612).

Author Contributions

Zeyang Dou and Kun Gao proposed model, Bin Zhang derived the optimization scheme of the model, Xinyan Yu proved the solution existence of the model, Lu Han and Zhenyu Zhu made the numerical comparisions.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Visual System. Available online: https://en.wikipedia.org/wiki/Visual_system (accessed on 15 February 2016).
  2. Daw, N.; Goldfish, W. Retina: Organization for simultaneous color contrast. Science 1967, 158, 942–944. [Google Scholar] [CrossRef] [PubMed]
  3. Land, E.H. An alternative technique for the computation of the designator in the retinex theory of color vision. Proc. Natl. Acad. Sci. USA 1986, 83, 3078–3080. [Google Scholar] [CrossRef] [PubMed]
  4. Land, E.H.; McCann, J. Lightness and retinex theory. J. Opt. Soc. Am. 1971, 61, 1–11. [Google Scholar] [CrossRef] [PubMed]
  5. Alonso, J.M.; Chen, Y. Receptive field. Scholarpedia 2009, 4, 5393. [Google Scholar] [CrossRef]
  6. Lindeberg, T. A computational theory of visual receptive fields. Biol. Cybern. 2013, 107, 589–635. [Google Scholar] [CrossRef] [PubMed]
  7. Marini, D.; Rizzi, A. A computational approach to color adaptation effects. Image Vis. Comput. 2000, 18, 1005–1014. [Google Scholar] [CrossRef]
  8. Provenzi, E.; Marini, D.; de Carli, L.; Rizzi, A. Mathematical definition and analysis of the Retinex algorithm. JOSA A 2005, 22, 2613–2621. [Google Scholar] [CrossRef] [PubMed]
  9. Cooper, T.J.; Baqai, F.A. Analysis and extensions of the Frankle-McCann Retinex algorithm. J. Electron. Imaging 2004, 13, 85–92. [Google Scholar] [CrossRef]
  10. Funt, B.; Ciurea, F.; McCann, J. Retinex in MATLAB. J. Electron. Imaging 2004, 13, 48–57. [Google Scholar] [CrossRef]
  11. McCann, J. Lessons learned from mondrians applied to real images and color gamuts. In Proceedings of the Color and Imaging Conference, Scottsdale, AZ, USA, 16–19 November 1999; pp. 1–8.
  12. McCann, J.; Sobel, I. Experiments with Retinex; Technical Report for HPL Color Summit: Hewlett Packard Laboratories, Israel, 1998. [Google Scholar]
  13. Morel, J.M.; Petro, A.B.; Sbert, C. A PDE formalization of retinex theory. IEEE Trans. Image Process. 2010, 19, 2825–2837. [Google Scholar] [CrossRef] [PubMed]
  14. Kimmel, R.; Elad, M.; Shaked, D. A variational framework for retinex. Int. J. Comput. Vis. 2003, 52, 7–23. [Google Scholar] [CrossRef]
  15. Ma, W.; Morel, J.M.; Osher, S. An L 1-based variational model for Retinex theory and its application to medical images. In Proceedings of the IEEE Conference on Computer Vision and Pattern Recognition, Colorado Springs, CO, USA, 21–23 June 2011; pp. 153–160.
  16. Ma, W.; Osher, S. A TV Bregman iterative model of Retinex theory. UCLA CAM Rep. 2010, 10–13. [Google Scholar] [CrossRef]
  17. Ng, M.K.; Wang, W. A total variation model for Retinex. SIAM J. Imaging Sci. 2011, 4, 345–365. [Google Scholar] [CrossRef]
  18. Liang, J.; Zhang, X. Retinex by Higher Order Total Variation L1Decomposition. J. Math. Imaging Vis. 2015, 52, 345–355. [Google Scholar] [CrossRef]
  19. Zosso, D.; Tran, G.; Osher, S. Non-Local Retinex-A Unifying Framework and Beyond. J. SIAM Imaging Sci. 2015, 8, 787–826. [Google Scholar] [CrossRef]
  20. Dominique, Z.; Tran, G.; Osher, S. A unifying retinex model based on non-local differential operators. In Proceedings of the SPIE—The International Society for Optical Engineering, San Diego, CA, USA, 25–29 August 2013; pp. 1–12.
  21. Blomgren, P.; Chan, T.F.; Mulet, P. Total variation image restoration: numerical methods and extensions. In Proceedings of the IEEE Conference on Image Processing, Santa Barbara, CA, USA, 26–29 October 1997; pp. 384–389.
  22. Rudin, L.I.; Osher, S.; Fatemi, E. Nonlinear total variation based noise removal algorithms. Phys. D 1992, 60, 259–268. [Google Scholar] [CrossRef]
  23. Chen, Y.; Levine, S.; Rao, M. Variable exponent, linear growth functionals in image restoration. SIAM J. Appl. Math. 2006, 66, 1383–1406. [Google Scholar] [CrossRef]
  24. Li, F.; Li, Z.; Pi, L. Variable exponent functionals in image restoration. Appl. Math. Comput. 2010, 216, 870–882. [Google Scholar] [CrossRef]
  25. Xie, S.J.; Lu, Y.; Yoon, S. Intensity Variation Normalization for Finger Vein Recognition Using Guided Filter Based Singe Scale Retinex. Sensors 2015, 15, 17089–17105. [Google Scholar] [CrossRef] [PubMed]
  26. Meylan, L.; Süsstrunk, S. High dynamic range image rendering with a retinex-based adaptive filter. IEEE Trans. Image Process. 2006, 15, 2820–2830. [Google Scholar] [CrossRef] [PubMed]
  27. Diening, L.; Harjulehto, P.; Hästö, P. Lebesgue and Sobolev Spaces with Variable Exponents; Springer: Berlin, Germany, 2011. [Google Scholar]
  28. Goldstein, T.; Osher, S. The split Bregman method for L1-regularized problems. SIAM J. Imaging Sci. 2009, 2, 323–343. [Google Scholar] [CrossRef]
  29. Cai, J.F.; Osher, S.; Shen, Z. Split Bregman methods and frame based image restoration. Multiscale Model. Simul. 2009, 8, 337–369. [Google Scholar] [CrossRef]
  30. Nien, H.; Fessler, J.A. A convergence proof of the split Bregman method for regularized least-squares problems. SIAM J. Imaging Sci. 2014, 4371, 1402. [Google Scholar]
  31. Elvetun, O.L.; Nielsen, B.F. The split Bregman algorithm applied to PDE-constrained optimization problems with total variation regularization. Comput. Optim. Appl. 2016, 1–26. [Google Scholar] [CrossRef]
  32. Zou, J.; Fu, Y. Split Bregman algorithms for sparse group Lasso with application to MRI reconstruction. Multidimens. Syst. Signal. Process. 2015, 26, 787–802. [Google Scholar] [CrossRef]
  33. Chen, S.; Liu, H.; Hu, Z. Simultaneous Reconstruction and Segmentation of Dynamic PET via Low-Rank and Sparse Matrix Decomposition. IEEE Trans. Biomed. Eng. 2015, 62, 1784–1795. [Google Scholar] [CrossRef] [PubMed]
  34. Esser, E. Applications of Lagrangian-based alternating direction methods and connections to split Bregman. CAM Rep. 2009, 9, 31. [Google Scholar]
  35. Stephen, B.; Neal, P.; Eric, C. Distributed optimization and statistical learning via the alternating direction method of multipliers. Found. Trends Mach. Learn. 2011, 3, 1–122. [Google Scholar]
  36. Petro, A.B.; Sbert, C.; Morel, J.M. Multiscale retinex. IPO 2014, 71–88. [Google Scholar] [CrossRef]
  37. Buttkus, B. Homomorphic filtering—Theory and practice. Geophys. Prospect. 1975, 23, 712–748. [Google Scholar] [CrossRef]
Figure 1. Elementary structure of the retina, lateral geniculate nucleus, and cortex.
Figure 1. Elementary structure of the retina, lateral geniculate nucleus, and cortex.
Sensors 16 00832 g001
Figure 2. Receptive field in the retina. (a) on-center ganglion cell; (b) off-center ganglion cell.
Figure 2. Receptive field in the retina. (a) on-center ganglion cell; (b) off-center ganglion cell.
Sensors 16 00832 g002
Figure 3. Images with different illumination and the corresponding surfaces: (a) text image; (b) book image; (c) surface of (a); (d) surface of (b).
Figure 3. Images with different illumination and the corresponding surfaces: (a) text image; (b) book image; (c) surface of (a); (d) surface of (b).
Sensors 16 00832 g003
Figure 4. A single row extracted from: (a) text image; (b) book image.
Figure 4. A single row extracted from: (a) text image; (b) book image.
Sensors 16 00832 g004
Figure 5. Schematicdiagram for Retinex.
Figure 5. Schematicdiagram for Retinex.
Sensors 16 00832 g005
Figure 6. Global framework of the proposed method.
Figure 6. Global framework of the proposed method.
Sensors 16 00832 g006
Figure 7. Synthetic example. (a,d) original image; (b,e) simulated illumination; (c,f) synthetic image.
Figure 7. Synthetic example. (a,d) original image; (b,e) simulated illumination; (c,f) synthetic image.
Sensors 16 00832 g007
Figure 8. Image reconstruction by different methods. (a,e) Ng; (b,f) HoTVL1; (c,g) proposed; (d,h) mutiscale Retinex.
Figure 8. Image reconstruction by different methods. (a,e) Ng; (b,f) HoTVL1; (c,g) proposed; (d,h) mutiscale Retinex.
Sensors 16 00832 g008
Figure 9. SNR curves. (a) For Figure 7a; (b) For Figure 7d.
Figure 9. SNR curves. (a) For Figure 7a; (b) For Figure 7d.
Sensors 16 00832 g009
Figure 10. Test images. (a) Checkerboard image; (b) Tower image; (c) Text image.
Figure 10. Test images. (a) Checkerboard image; (b) Tower image; (c) Text image.
Sensors 16 00832 g010
Figure 11. Illumination (ad) and reflection (eh) images. (a,e) Ng’s method; (b,f) HoTVL1; (c,g) proposed method; (d,h) multiscale Retinex.
Figure 11. Illumination (ad) and reflection (eh) images. (a,e) Ng’s method; (b,f) HoTVL1; (c,g) proposed method; (d,h) multiscale Retinex.
Sensors 16 00832 g011
Figure 12. Reproduced degraded image (Figure 10b). (a) Ng’s model; (b) HoTVL1; (c) proposed model; (d) multiscale Retinex.
Figure 12. Reproduced degraded image (Figure 10b). (a) Ng’s model; (b) HoTVL1; (c) proposed model; (d) multiscale Retinex.
Sensors 16 00832 g012
Figure 13. Reproduced shadowed images (Figure 10c and Figure 3b): (a,e) Ng’s model; (b,f) HoTVL1; (c,g) proposed model; (d,h) multiscale Retinex.
Figure 13. Reproduced shadowed images (Figure 10c and Figure 3b): (a,e) Ng’s model; (b,f) HoTVL1; (c,g) proposed model; (d,h) multiscale Retinex.
Sensors 16 00832 g013
Figure 14. Reproduced illumination and reflectance (Figure 10a) (a) Illumination with fixed approximation; (b) Reflectance with fixed approximation. Note the circled part by the red square.
Figure 14. Reproduced illumination and reflectance (Figure 10a) (a) Illumination with fixed approximation; (b) Reflectance with fixed approximation. Note the circled part by the red square.
Sensors 16 00832 g014
Table 1. SSIM of the four methods.
Table 1. SSIM of the four methods.
ImageNgHoTVL1Proposed MethodMultiscale Retinex
Figure 7a0.90760.63420.92890.8764
Figure 7d0.76030.66330.82910.7238
Table 2. CIEDE2000colordifference of the four methods.
Table 2. CIEDE2000colordifference of the four methods.
ImageNgHoTVL1Proposed MethodMultiscale Retinex
Figure 7a26.892729.672926.618227.9038
Figure 7d24.547426.011521.997225.4528
Table 3. Recovered intensity values for regions A and B of Figure 10a.
Table 3. Recovered intensity values for regions A and B of Figure 10a.
Image OriginalNgHoTVL1Proposed MethodMultiscale Retinex
CheckerboardA12014085135109
B120180174230149

Share and Cite

MDPI and ACS Style

Dou, Z.; Gao, K.; Zhang, B.; Yu, X.; Han, L.; Zhu, Z. Realistic Image Rendition Using a Variable Exponent Functional Model for Retinex. Sensors 2016, 16, 832. https://doi.org/10.3390/s16060832

AMA Style

Dou Z, Gao K, Zhang B, Yu X, Han L, Zhu Z. Realistic Image Rendition Using a Variable Exponent Functional Model for Retinex. Sensors. 2016; 16(6):832. https://doi.org/10.3390/s16060832

Chicago/Turabian Style

Dou, Zeyang, Kun Gao, Bin Zhang, Xinyan Yu, Lu Han, and Zhenyu Zhu. 2016. "Realistic Image Rendition Using a Variable Exponent Functional Model for Retinex" Sensors 16, no. 6: 832. https://doi.org/10.3390/s16060832

APA Style

Dou, Z., Gao, K., Zhang, B., Yu, X., Han, L., & Zhu, Z. (2016). Realistic Image Rendition Using a Variable Exponent Functional Model for Retinex. Sensors, 16(6), 832. https://doi.org/10.3390/s16060832

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop