# Nonuniform Illumination Correction Algorithm for Underwater Images Using Maximum Likelihood Estimation Method

Authors: Sonali Sachin Sankpal; Shraddha Sunil Deshpande
Journal: Journal of Engineering
 (2016)

---

## Abstract

Scattering and absorption of light is main reason for limited visibility in water. The suspended particles and dissolved chemical compounds in water are also responsible for scattering and absorption of light in water. The limited visibility in water results in degradation of underwater images. The visibility can be increased by using artificial light source in underwater imaging system. But the artificial light illuminates the scene in a nonuniform fashion. It produces bright spot at the center with the dark region at surroundings. In some cases imaging system itself creates dark region in the image by producing shadow on the objects. The problem of nonuniform illumination is neglected by the researchers in most of the image enhancement techniques of underwater images. Also very few methods are discussed showing the results on color images. This paper suggests a method for nonuniform illumination correction for underwater images. The method assumes that natural underwater images are Rayleigh distributed. This paper used maximum likelihood estimation of scale parameter to map distribution of image to Rayleigh distribution. The method is compared with traditional methods for nonuniform illumination correction using no-reference image quality metrics like average luminance, average information entropy, normalized neighborhood function, average contrast, and comprehensive assessment function.


---

## Body

## 1. Introduction

When light travels through the water it deviates from its path which is called scattering of the light. The scattering is caused by water molecules themselves, dissolved chemical compounds in water, and suspended particles. This scattering decreases visibility in water. As you go deeper in the water the scattering of natural light increases which results in poor visibility in deep water. This degrades the performance of the underwater imaging system and it goes worse for deep water imaging.Propagation of light in water is discussed by Jaffe in [1]. His image formation model divides reflected light in three components given by the following:(1)ET=Ed+Ef+Eb,where ET is total irradiance, Ed is direct component, Ef is forward scattered component, and Eb is backscattered component.Forward and backward scattering of light degrade the quality of underwater images. The degradation includes limited contrast, addition of blur, and diminished colors.The image enhancement techniques are used to compensate for the degradation like contrast enhancement and color correction. Underwater visibility can be increased by using artificial light source. But artificial light adds nonuniform illumination to image. This problem is ignored by many researchers. Very few researchers proposed a correction method for nonuniform illumination in underwater images.The nonuniform illumination problem in underwater images is shown in Figure1. An artificial light source is used in underwater imaging system to increase visibility in water. The light rays from artificial light source are reflected in backward direction without reaching the objects in the scene. This is called backward scattering. This reflected light produces bright spot in the center and is surrounded by dark region [2]. This illuminates the scene in nonuniform fashion. It is assumed that pixel intensity values of underwater images are dominated by Rayleigh scattering [3]. The scattering of light is wavelength dependent so for correction of nonuniform illumination each color component (R,G,B) should be processed separately. While providing solution for nonuniform illumination, the method assumes that the underwater images are Rayleigh distributed [4–6]. Consider(2)
f
x
x
=
x
α
e
-
x
2
/
2
α
2
,
x
≥
0
,
0
,
o
t
h
e
r
w
i
s
e
.Figure 1
Nonuniform illumination created in underwater image.The Rayleigh distribution function is bell shaped with major amount of pixels being concentrated in middle intensity levels. The probability distribution function of Rayleigh distributed random variablex with parameter α is given by (2), where α is mean square value of x.The method maps image to Rayleigh distribution by estimation of parameters using maximum likelihood method. The principle of maximum likelihood method [7] says that given observations x1,x2,…,xn, fx(x1,x2,…,xn;θ) are function of θ alone, and the value of θ that maximizes the above probability density function is the most likely value for θ, and it is chosen as its maximum likelihood estimation θML.This paper proposed a method for nonuniform illumination correction for underwater images. The rest of the paper is organized as follows. Section2 discussed state of the art for the problem, in Section 3 a proposed method for nonuniform illumination correction is explained, the image quality metrics are discussed in Section 4, results are reported in Section 5, and conclusion is given in Section 6.


## 2. Literature Review

Arnold-Bos et al. [8] suggested that global histogram equalization is not suitable when illumination in image is unequal, and local methods are needed. They proposed a method of histogram clipping and then equalized contrast by division method.Bazeille et al. [9] suggested homomorphic filtering to correct nonuniform illumination. The maximum and minimum coefficient values selected by them are γH=2.5 and γL=0.5.Garcia et al. [10] compared four methods for nonuniform illumination correction. The first method uses illumination reflectance model and illumination correction is obtained as (3)r~x,y=fx,yfsx,yδ,where f(x,y) is original and fs(x,y) is smoothed version of image. δ is constant. Contrast Limited Adaptive Histogram Equalization (CLAHE) and homomorphic filter are second and third methods. In fourth method illumination field is subtracted from original image for nonuniform illumination correction.Borgetto et al. [11] also discussed two methods for illumination correction: one method used homomorphic filter and the other is based on CCD camera radiometric correction required for mosaicking.Prabhakar and Praveen Kumar [12] proposed algorithm for enhancement of underwater images. They used homomorphic filter to correct nonuniform illumination.


## 3. Proposed Method

The proposed method in this paper improves quality of underwater images by correcting nonuniform illumination in image. Like traditional methods proposed by Iqbal et al. in [13, 14], histogram of the image is modified for correction. The method in this paper applied histogram stretching. The detailed process is given in Figure 2.Figure 2
Process of proposed algorithm.Input image is first decomposed into three channels (red, green, and blue). Then histogram stretching is performed on individual channels. With the assumption that each of R, G, and B channel is Rayleigh distributed, histogram stretching is done with respect to Rayleigh distribution. The scale parameter is estimated from given image. So the histogram stretching is adaptive. The scale parameter(θ) is estimated using the maximum likelihood method given below. First find log likelihood function of image as given by(4)lθ=ln⁡∏i=1nfxi;θ.Heref(xi;θ) is a probability density function of xi with scale parameter θ. Consider(5)lθ=ln⁡∏i=1nxiθe-xi2/2θ,lθ=∑i=1nln⁡xiθe-xi2/2θ,lθ=∑i=1nln⁡xi-nln⁡θ-∑i=1nxi22θ.Then find the value ofθ where log likelihood function has maximum value. It can be determined by taking a derivative of log likelihood function with respect to θ and equate it to zero:(6)∂∂θlθ=0.Solve (6) for θ: (7)∂∂θ∑i=1nln⁡xi-nln⁡θ-∑i=1nxi22θ=0,(8)-nθ-∑i=1nxi22-1θ2=0,(9)∑i=1nxi22θ2=nθ,(10)θML=12n∑i=1nxi2.Obtainedθ in (10) is estimated maximum likelihood value θML. Thus maximum likelihood values are estimated for all three channels using same process. These estimated values are used in histogram stretching with respect to Rayleigh distribution of R, G, and B components. As image is nonuniformly illuminated there are dark and bright patches in the image, so stretching is performed locally on small patches of image, instead of global stretching. A limit of 1% is applied to all three color components at minimum and maximum values. The limits are necessary to avoid under- and oversaturation [15]. The stretching process is applied in the range of 1% to 99%. Histogram stretching performed on each color channel with respect to Rayleigh distribution is given by the following:(11)iout=imin+2∗α2∗ln⁡1pii1/2,where iout is pixel value in transformed image, imin is minimum pixel value in the transformed image, α is parameter value, and pi(i) is cumulative distribution function of pixel values of input image. The effect of histogram stretching is shown in Figure 3. Figure 3 gives histogram of red, green, and blue components after and before stretching along with images.Figure 3
Original and processed images with their histogram for three color components (Original Image Courtesy NIO, Goa, India).


## 4. Image Quality Metrics

The image quality metrics used in this paper are no-reference color image quality metrics. These metrics are based on primary and fundamental vision parameters perceived by human vision system, which are luminance, contrast, information content, and added noise. According to human vision perception a good quality image should have appropriate luminance, information content, and contrast. These four vision parameters can be measured using mathematical model given by Xie and Wang [16]. These are average contrast (AC), average information entropy (AIE), and average luminance (AL).Average contrast is computed as(12)AC=1M-1×N-1∑x=0M-2∑y=0N-2Cx,y,where C(x,y) is magnitude of mean of gradients of three color components R, G, and B given as (13)C=13GradR2+GradG2+GradB2.Information entropy for single color channel is computed as(14)IE=-∑i=0255pilog2⁡pi,where p(i) is probability of ith color level.Then using (14) IE is calculated for red, green, and blue channels as IER, IEG, and IEB, respectively. Total AIE is defined as(15)AIE=13IER2+IEG2+IEB2.The maximal value of AIE for color image is 8 bits.Average luminance is computed as follows:(16)AL=1M×N∑x=0M-1∑y=0N-1lx,y.M  and  N are rows and columns of the image and l(x,y) is luminance value at pixel (x,y). When histogram of luminance component (lx,y∈[0,255]) of an image is equalized then ideal value of AL should be 127.5.So 127.5 is considered as an optimum value of luminance (OL) [16].Comprehensive assessment function (CAF) is a general image quality assessment function given by Xie and Wang [16] computed as follows: (17)CAF=AIEα×ACβ×NNFγ,where α, β, and γ are weight parameters; through experiment and comparison with subjective assessment Xie and Wang decided its values as 1, 1/4, and 3, respectively [16]. In the above equation (17) AIE is average information entropy computed using (15), AC is average contrast computed using (12), and NNF is normalized neighborhood function defined as(18)NNF=OL-distAL,OLOL,where OL is optimum luminance value and here it is considered [16] 127.5 and dist is absolute value. When AL is equal to OL then dist(AL,OL) is zero and NNF has its optimum value of 1. CAF computed by (17) is a convex function [16] and if the CAF is larger then the image quality is better.


## 5. Results

The proposed method is applied to five nonuniform illuminated images. The results are compared with the results of traditional methods of nonuniform illumination correction, which are histogram equalization (HE), adaptive histogram equalization (AHE), and homomorphic filter.Comparison of results is performed using no-reference image quality metrics for color images given in previous section.The results are shown in Figure4. Table 1 shows quantitative comparison of image quality metrics of original and processed images. The superior results are indicated in bold fonts in the table.Table 1
Comparison of results.
Method
AC
AIE
AL
NNF
CAF
Image 1
Original
7.88
6.68
80.35
0.63
8.61
Proposedmethod
22.93
7.11
121.38
0.95
9.70
HE
16.13
5.32
127.78
0.99
8.32
AHE
20.19
6.91
103.55
0.81
9.57
Homomorphicfilter
5.69
5.49
47.39
0.37
7.09
Image 2
Original
7.3
6.37
199.46
0.44
7.08
Proposedmethod
9.97
7.53
148.61
0.83
9.85
HE
7.17
5.76
125.88
0.98
8.35
AHE
6.31
7.23
167.6
0.68
9.14
Homomorphicfilter
3.53
7.61
176.86
0.61
8.61
Image 3
Original
2.44
6.4
198
0.45
7.74
Proposedmethod
9.1
7.04
171.11
0.66
9.10
HE
9.92
5.25
130.69
0.97
7.96
AHE
6.18
6.81
180.71
0.58
8.59
Homomorphicfilter
3.22
6.44
180.52
0.58
7.98
Image 4
Original
5.03
7.75
104.38
0.81
9.8
Proposedmethod
15.71
7.79
131.34
0.97
10.69
HE
5.94
5.97
125.51
0.98
8.48
AHE
12.00
7.84
115.34
0.9
10.4
Homomorphicfilter
4.09
6.7
68.64
0.53
8.29
Image 5
Original
2.17
7.37
118.57
0.93
9.38
Proposedmethod
11.67
7.61
141.23
0.89
10.18
HE
6.85
7.45
124.49
0.97
9.99
AHE
4.85
6.69
116.5
0.91
9.2
Homomorphicfilter
1.84
6.94
63.59
0.49
8.23Figure 4
Results of all four methods.
(a)
Original image
(b)
Proposed method
(c)
Histogram equalization method
(d)
Adaptive histogram equalization method
(e)
Homomorphic filter methodThe parameters values for homomorphic filter are selected as given by Bazeille et al. [9].


## 6. Conclusion

The proposed method is used for nonuniform illumination correction of underwater images. Nonuniform illumination affects the overall contrast of the image. The proposed method shows improvement in all the quality metrics when compared with original image. Also there is improvement in quality parameters like average contrast, average information entropy, and comprehensive assessment function when the results are compared with results of traditional methods. But there is small degradation in average luminance and normalized neighborhood function compared to the results of traditional methods. However the degradation is by very small amount as compared to the improvements in other quality parameters. Also general image quality assessment function (CAF) which includes all other quality parameters (AC, AIE, and NNF) shows improvement for proposed method when compared with other methods. So the proposed method shows improvement in the results.



---

*Source: 5718297-2016-01-27.xml*