Adaptive Multi-Scale Entropy Fusion De-Hazing Based on Fractional Order

This paper describes a proposed fractional filter-based multi-scale underwater and hazy image enhancement algorithm. The proposed system combines a modified global contrast operator with fractional order-based multi-scale filters used to generate several images, which are fused based on entropy and standard deviation. The multi-scale-global enhancement technique enables fully adaptive and controlled color correction and contrast enhancement without over exposure of highlights when processing hazy and underwater images. This in addition to the illumination/reflectance estimation coupled with global and local contrast enhancement. The proposed algorithm is also compared with the most recent available state-of-the-art multi-scale fusion de-hazing algorithm. Experimental comparisons indicate that the proposed approach yields a better edge and contrast enhancement results without a halo effect, without color degradation, and is faster and more adaptive than all other algorithms from the literature.


Introduction
Hazy and underwater images share similar characteristics in terms of reduced visibility and low contrast due to the nature of image formation [1,2].Several single image-based enhancement and restoration models and algorithms have been proposed to solve this problem [1,2].However, they work with varying degrees of success at the cost of increased structural and computational complexity.Furthermore, color correction combined with these highly complex de-hazing algorithms have been used to restore underwater images.However, there are relatively few digital hardware realizations and reduced real-time prospects for such schemes due to a high computational cost.
In this work, we propose a fractional order-based algorithm for the enhancement of hazy and underwater images.The algorithm performs color correction and multiscale spatial filter-based localized enhancement.We compare results with other algorithms from the literature and show that the proposed system is effective with the fastest execution time.
The paper is outlined as follows.The Section 2 provides the background, motivation, and key contributions of the proposed system.Section 3 presents the proposed algorithms for both underwater and hazy image enhancement in addition to solutions to problems and modifications.Section 4 presents and compares the results (obtained using the proposed system) to other algorithms from the literature.The Section 5 explicitly compares the proposed approach against a recent algorithm from the literature, which further strengthens the justification of the proposed scheme.The Section 6 presents the conclusion.

Hazy Image Processing Algorithms
Hazy image processing also deals with visibility restoration of image scenes degraded by weather conditions and can be multi-image or single-image based solutions [31].Furthermore, hazy image processing algorithms can also be classified as either restoration or enhancement-based schemes.In the restoration-based hazy image processing, the de-hazing process is based on the hazy image formation model [31].The objective is, therefore, meant to obtain the de-hazed image from the input hazy image.The algorithms in this class include the popular DCP method by He et al. [32], which has been adopted and modified in various forms and a review of several DCP-based methods can be found in Reference [1].
Other schemes include works based on segmentation [33][34][35], fusion [36,37], geometry [38], Weighted Least Squares [39], variational [37,[40][41][42], and regularization approaches [34] using sparse priors [43], other boundary constraints [44], biological retina-based model [45], and multi-scale convolutional neural networks [46].The enhancement-based hazy image processing method is based on directly obtaining the by-product of radiance scene recovery through visibility restoration by contrast enhancement/maximization.The algorithms in this category utilize contrast limited adaptive histogram equalization (CLAHE), histogram specification (HS) [47], and Retinex [48][49][50].Additionally, some of these algorithms combine dark channel priors and transmission map extraction with contrast enhancement for refinement.However, consistently good results are not guaranteed since some images will depict color fading/distortion and darkening of regions in addition to over-enhancement of sky/homogeneous regions.Thus, some threshold and segmentation-based algorithms [33][34][35]51] have been developed to solve the peculiar problems of these algorithms.Furthermore, recently developed algorithms using partial differential equations (PDEs) and gradient metric-based optimization were developed [52,53] to avoid the usage of DCP-based stages and multiple (and manual adjustment of) parameters.Recently, an Artificial Multiple-Exposure Image Fusion (AMEF) de-hazing algorithm was proposed by Galdran [54], which represents the current state-of-the-art technology.
Physical methods depend on prior image information obtained by capturing the image scenes at different times under varying conditions using physical hardware/optical equipment such as cameras and lighting rigs [2].They may also incorporate multi-image processing schemes for either hazy or underwater images.However, consistently good results are not assured due to the unpredictable nature of weather and aquatic medium conditions.In addition, the cost of such hardware imaging systems is prohibitive and are usually not universally applicable.Such schemes are fully listed and described in work by Li et al. [5].Single-image-based software implementations offer the best outcome when factors such as cost, time, replicability, and convenience are considered since they do not necessarily require prior knowledge of the environment or an image acquisition process for operation [2,5].Thus, the scope of this work is limited to single-image-based enhancement of both hazy and underwater images.
The primary motivation for this work is to develop fast, practical, and effective algorithms for underwater and hazy image enhancement that are amenable to hardware implementation for real-time operation.

Key Contributions and Features of a Proposed Scheme
The key contributions and features of this work include:

•
A modified global contrast enhancement and a multi-scale illumination/reflectance model-based algorithm using fractional order calculus-based kernels.

•
Relatively low-complexity underwater image enhancement algorithm utilizing color correction and contrast operators.

•
Frequency-based approach to image de-hazing and underwater image enhancement using successive, simultaneous high frequency component augmentation and low frequency component reduction.

•
Feasible hazy and underwater image enhancement algorithm for relatively easier hardware architecture implementation.

•
Avoidance of the dark channel prior based stages and iterative schemes by utilizing combined multi-level convolution using fractional derivatives.

Proposed Algorithms
Underwater image enhancement usually involves some color correction/white balancing in addition to contrast enhancement process, which is usually a local/global operation.The first step for reducing the need for such involved local processing was to avoid the over-exposure of bright regions while enhancing the dark regions.Initial logarithmic solutions were ineffective and flattened the images in addition to fading colors.Thus, a new formulation for the global contrast operator had to be devised to achieve this objective.We present the modification and realization of the improved global contrast operator and spatial filter based system for processing underwater and hazy images.Furthermore, the simplified scheme using fractional calculus is presented in the form of spatial masks based on the Grunwald-Letnikov definition [74].

Selection and Modification of Global Contrast Operator
Previously, extensive experiments were conducted to determine effectiveness on several contrast stretching algorithms [26].Due to the adjustable nature of the high and low values by modifying the percentiles, the contrast stretching (CS) algorithm appeared to be much more versatile than the other algorithms.However, it works best for faded low-contrast images but not so well for underwater images since it does not perform adequate color correction unless applied iteratively.Conversely, some of the other algorithms were too harsh and had no effect or minimal impact on most underwater images while others resulted in color bleeding.The selected algorithms such as the piecewise linear transform (PWL) [75] and the gain offset correction (GOC2) [76] were selected for incorporation into effective PDE-based formulations [26,27].This was because some underwater images responded better to GOC2 (due to its mainly color correction ability) than to PWL (due to its generality) and vice versa.Thus, there is the need to develop a global contrast operator that would merge the advantages of both GOC2 and PWL while mitigating their weaknesses.
Since the linear contrast stretch (similar to the PWL and GOC2) does not utilize any edge enhancement features or region-based methods, it does not enhance noisy edge artefacts.However, several of these contrast stretching algorithms lead to oversaturation of already bright regions of the image (whitening out/over-exposure).This is in addition to a threshold effect when applied to images with bimodal histogram.The linear contrast stretching can be applied to both greyscale and color images with excellent results similar to the PWL.However, the PWL method also suffers from the threshold of images when there are distinct regions of dark and light intensity, which leads to whiting out of bright areas.This is because it truncates values at upper and lower limits to maximum and minimum possible pixel values in the image without taking into account pixels in those regions.The linear contrast stretch seeks to expand the range based on the surrounding pixels in the distribution.
Underwater image enhancement usually involves some color correction/white balancing in addition to contrast enhancement process, which is usually a local/global operation.The GOC2 algorithm adequately processed underwater images, which required mild color correction and contrast enhancement and, thus, avoided overexposure of highlights unlike most other tested contrast enhancement algorithms [26].This necessitated the incorporation of a local contrast operator such as the CLAHE, which, even though effective, further added to the computational complexity of the algorithms and introduced additional parameters.The first step to reducing the need for such involved local processing was to avoid the over-exposure of bright regions while enhancing the dark regions.Initial logarithmic solutions were not effective and flattened the images in addition to fading colors.Thus, a new formulation for the global contrast operator had to be devised to achieve this objective.

Gain Offset Correction-Based Stretching (GOCS)
The expression for the GOC2 mapping function [76] is given as shown in Equation (1).
The contrast stretching function is given below.
In Equations (1) and ( 2), I o GOC2 and I o CS are the enhanced images using GOC2 and CS respectively, I max , I min are maximum and minimum pixel intensities in the input image, I i , L is the number of grey intensity levels (L = 256 for unsigned integer, eight-bit-per-pixel (uint8, 8 bpp) image format) while I low and I high are the lower and upper percentiles of the image pixel intensity distribution normally set at 5% and 95%, respectively.
The faults of the GOC2 lie in the statistics such as maximum and minimum pixel intensity values utilized in its computation.Since an image, which is already utilizing its full dynamic range, will not be affected by such statistics, we needed to realize a more influential statistic.The contrast stretching operator utilizes lower and upper percentiles of the image intensity distribution for its computation.Consequently, the contrast stretching operator does not suffer over-exposure effects and performs adequate contrast enhancement.Conversely, the GOC2 performs sufficient color correction but minimal contrast enhancement.Thus, by replacing the maximum and minimum pixel intensity values with the upper and lower percentiles in the formulation, we can realize a new formula for the global contrast operation as seen in the equation below.
Initial experiments using the 5th and 95th percentiles led to some pixels being over-exposed and, as we widened the range between the percentiles, the results improved and, in some cases, settled on the 1st and 99th percentiles for the best results.Increasing the range to its maximum yields a result similar to GOC2, which is expected since the high and low percentiles now become the maximum and minimum pixel intensity values.The GOCS is related to the CS in the following form.

Proposed Multi-Scale Local Contrast Operator
We present the development of the multi-scale algorithm for local contrast enhancement, which replaces the CLAHE used in previous work by drastically reducing complexity and run-time.The initial derivations of the illumination-reflectance estimation part of the algorithm using integer-based low-pass and high-pass filters can be found in previous work [2].We focus on new aspects such as multi-scale, fractional order-based entropy and standard deviation guided enhancement.

Multi-Scale Fractional Order-Based Illumination/Reflectance Contrast Enhancement (Multi-fractional-IRCES)
The fractional derivative-based re-definitions for high-pass and low-pass filtering of arbitrary order α are obtained using the equations below. and Leading to the expression using fractional order as seen in Equation ( 7) below.
We further extend the application to hazy image enhancement as seen in the equations below.
In the expressions of Equations ( 5)- (10), k is the power factor, U(x, y) is the inverted image, I max is the maximum pixel intensity of the input image, I(x, y), ∇ α U(x, y) is the fractional derivative of the inverted image, and Ω ∇ α U(x, y)dΩ denotes the fractional order integral.Additionally, U α e (x, y) is the enhanced inverted image using fractional order-based operations and U α e max is the maximum pixel intensity of U α e (x, y) while I α e (x, y) is the de-hazed image using fractional order-based operations.Additionally, we wish to reduce the computational load of computing both the derivative and the integral especially in the fractional order-based version.Thus, we simple obtain the fractional integral of the input image and subtract it from the original image and multiply by the appropriate factor to obtain the fractional order derivative.This saves resources especially on digital hardware implementations since only one operator is utilized and re-used.This is easily expressed as the equations below.
I HPF (x, y) = I(x, y) − I LPF (x, y) (11) which gives the expressions using fractional order calculus as seen below.
The scheme for hazy image enhancement can also be updated accordingly without much effort.The central idea is that, by further decomposing a low-pass filtered image and enhancing the details at each level and recombining the results, we would obtain much finer local enhancement.Additionally, using the fractional order reduces or minimizes the issue of noise enhancement as high frequency components are amplified at each stage, which further reduces or minimizes the low frequency components at each stage.Since the haze is a low frequency phenomenon, we expect that such effects would be greatly reduced after processing without enhancing noise.The entropy and standard deviation measures are utilized to select the best outcome for the processed image in terms of the value of the exponent, k.The mathematical expressions for the algorithm are shown in Equations ( 14)- (20).
In Equation ( 14), I i (x, y) is the enhanced image at level i and N is the number of decomposition levels while I HPFi (x, y) and I LPFi (x, y) are high-pass and low-pass filtered images obtained at level i.Based on experiments, we set N = 5.The obtained level images are then aggregated to obtain the final images I A k or I B k for the different values of the power factor and k in Equations ( 15) and ( 16).The values for the power factor are chosen to be multiples of two (2) due to hardware design considerations to enable fast computation by bit shifting.This is explained by the fact that the value 0.5 is a multiple of 2 since 2 n for n = −1; 2 −1 = 1 2 = 0.5 or 000.1.This is a right bit shift by one-bit position in binary logic for the binary representation of one (1) as 001.0.Conversely 2 n for n = 1; 2 1 = 2, or 010.0.This makes the operation synthesizable for digital hardware realization since exponential operations become bit shifting operations, which are much faster.
The respective entropies (e A k , e B k ) and standard deviations (σ A k , σ B k ) of the aggregated images are computed in Equations ( 17) and ( 18) and used to decide the best image outcome, f (x, y) in Equation (19), which is then processed with the modified global contrast enhancement algorithm to obtain the final output image, f o (x, y), in Equation (20).This is based on the simultaneous multi-level high frequency component (edges and details) enhancement and multi-level low frequency component (haze) attenuation.
The diagram of the proposed algorithm for enhancement of both hazy and underwater images is shown in Figure 1.All processing operations are achieved with spatial filter kernels using fractional order-based calculus, which yields better results than the integer order in terms of balanced edge enhancement.

Preliminary Results
After testing several images, it was discovered that some images were better enhanced when using the 5th and 95th percentiles rather than the 1st and 99th percentiles.The representative images of these two groups include those unaffected by wide ranges while the other exhibits over-exposure for narrow ranges.This was partly the reason that the PWL approach was utilized in previous work [27].Thus, one approach would be to devise a means of selecting the appropriate percentiles for these two groups of images.A simple compromise was to set the range between the 2nd and 98th percentiles.However, we would still be faced with the issue of outlier images, which resist color correction attempts.Thus, the need for the localized operator to aid in the detail recovery in the otherwise over-exposed regions when global contrast operations are performed.

Problems and Solutions
The initial developed scheme worked extremely well for underwater images and several hazy images.However, problems were observed in other hazy images.These issues included color fading, distortion, discoloration, image darkening, inadequate haze removal, and over-enhanced edges.Thus, we devised solutions to some of these problems.The color correction routine was omitted and the output, (, ), was reformulated as Equation ( 21) below.
This improved results and resolved color distortion in the affected hazy images.However, there was some color fading in RGB and HSI/HSV versions.Thus, we utilized the red-green-blue-intensity/value (RGB-IV) formulation [77] to improve the color rendition, which resulted in color enhancement but with dark images.We also investigated the use of CLAHE to improve local contrast, which results in drastic improvements.However, enhanced images also exhibited halo effects and color distortion, which persisted despite a combination with the multi-scale IRCES algorithm.Furthermore, there was drastic color loss/fading using CLAHE in addition to increased computational complexity, which defeats the initial objective of the proposed

Preliminary Results
After testing several images, it was discovered that some images were better enhanced when using the 5th and 95th percentiles rather than the 1st and 99th percentiles.The representative images of these two groups include those unaffected by wide ranges while the other exhibits over-exposure for narrow ranges.This was partly the reason that the PWL approach was utilized in previous work [27].Thus, one approach would be to devise a means of selecting the appropriate percentiles for these two groups of images.A simple compromise was to set the range between the 2nd and 98th percentiles.However, we would still be faced with the issue of outlier images, which resist color correction attempts.Thus, the need for the localized operator to aid in the detail recovery in the otherwise over-exposed regions when global contrast operations are performed.

Problems and Solutions
The initial developed scheme worked extremely well for underwater images and several hazy images.However, problems were observed in other hazy images.These issues included color fading, distortion, discoloration, image darkening, inadequate haze removal, and over-enhanced edges.Thus, we devised solutions to some of these problems.The color correction routine was omitted and the output, f (x, y), was reformulated as Equation ( 21) below.
This improved results and resolved color distortion in the affected hazy images.However, there was some color fading in RGB and HSI/HSV versions.Thus, we utilized the red-green-blueintensity/value (RGB-IV) formulation [77] to improve the color rendition, which resulted in color enhancement but with dark images.We also investigated the use of CLAHE to improve local contrast, which results in drastic improvements.However, enhanced images also exhibited halo effects and color distortion, which persisted despite a combination with the multi-scale IRCES algorithm.Furthermore, there was drastic color loss/fading using CLAHE in addition to increased computational complexity, which defeats the initial objective of the proposed approach.Thus, alternatives were considered to resolve these issues.
Wavelet-based fusion of I A k (x, y) and I B k (x, y) using mean, minimum, or maximum configurations was implemented.Good results were observed in images with mostly a uniform haze.Conversely, sky regions were degraded in hazy images with uneven haze or considerable sky regions.Furthermore, dark bands and outlines were observed around edges in some processed images.Overall, image results were inconsistent using this scheme.Thus, we reformulated the multi-scale algorithm after extensive analysis.
Redundant frequencies, which were unnecessary in hazy image enhancement results, were observed.This was due to the nature of the generation of the two combined images.I A k (x, y) and I B k (x, y) lead to unbalanced contributions of frequency components.Constant variance of weights for both images and corresponding results led to inconsistent results.Thus, a more formalized, systematic approach was required.Based on the analysis of the Fourier Transform of the images, we require subtle enhancement of the high frequency components and a drastic reduction of the contributions of the low frequency components.This informed the reformulation of the multi-scale algorithm for hazy images as seen in Equation ( 22) below.
In Equations ( 22)-( 27), U(x, y) and I(x, y) are the original and reversed hazy image, respectively, while U max is the maximum pixel intensity value of the image.I LPF i (x, y), I HPF i (x, y), S LPF i and S HPF i are the low-pass and high-pass filtered images of level (or scale) including i and their respective summations.The terms S total , p LPF i , p HPF i are the total sum and the percentage of low and high frequency components, respectively.In order to balance the high and low frequency components, we create new constants, c 1 and c 2 , to be dependent on each other using the percentages.
After evaluation of the two constants, we use the expression to obtain the enhanced level image as seen in Equation ( 29) below.
J. Imaging 2018, 4, 108 The level images are subsequently added to obtain the enhanced image, which is shown in Equation (30).
The de-hazed image, U (x, y) is obtained by inverting the image as shown Equation (31).
Based on experiments, we set c 1 and c 2 as 1.21 and 0.8264, respectively, since they are always constant.These are the default values for balanced enhancement of high and low frequency components to avoid visual artefacts.However, the values may be increased or decreased gradually for maximum visual effect in certain images.This new formulation solves the edge over-enhancement, color distortion, and halo effect problem.The results are shown in Figure 2 for processed images using previous and improved configurations of PA (test images used were obtained from the University of Texas at Austin Laboratory for Image & Video Engineering (LIVE) dataset, which can be downloaded from the website: http://live.ece.utexas.edu/research/fog/fade_defade.html)[78].Note the elimination of the color distortion and reduced degree of noise enhancement for images in Figure 2b when compared to Figure 2a.The estimated computational complexity of the proposed approach is given as ( ) for D levels using the spatial window size and w of the fractional order-based filter for an image with N rows and M columns.Additionally, the algorithm can be speeded up by exploiting symmetric convolutional structures to reduce the number of multiplications and additions.

Results
We present the result comparisons of the proposed approach (PA) with other algorithms from the literature.We utilize metrics such as entropy (E), (relative) average gradient (RAG) [79], global contrast factor (GCF) [80], and colorfulness or color enhancement factor (CEF) [81] for underwater images.For hazy images, we utilize the RAG, ratio of visible edges, Qe [1], and saturation parameter/percentage of black or white pixels, and σ [1] to evaluate results.Higher values indicate better results for all the used metrics except for the saturation parameter metric where lower values imply improvement.
The hardware specifications of the computing platform are PC with Intel ® Core i7-6500U x64-based processor (Intel, Santa Clara, CA, USA) at 2.5 GHz/2.59GHz, 12 GB RAM running 64-bit OS (Microsoft ® Windows™ 10 Home, Redmond, Washington, WA, USA) and NVIDIA ® GeForce™ 940M GPU with a compute capability of 5.0 (NVIDIA, Santa Clara, CA, USA).The estimated computational complexity of the proposed approach is given as O N Mw 2 D for D levels using the spatial window size and w of the fractional order-based filter for an image with N rows and M columns.Additionally, the algorithm can be speeded up by exploiting symmetric convolutional structures to reduce the number of multiplications and additions.

Results
We present the result comparisons of the proposed approach (PA) with other algorithms from the literature.We utilize metrics such as entropy (E), (relative) average gradient (RAG) [79], global contrast factor (GCF) [80], and colorfulness or color enhancement factor (CEF) [81] for underwater images.For hazy images, we utilize the RAG, ratio of visible edges, Qe [1], and saturation parameter/percentage of black or white pixels, and σ [1] to evaluate results.Higher values indicate better results for all the used metrics except for the saturation parameter metric where lower values imply improvement.

Subjective Evaluation
Even though objective metrics have been used for underwater image enhancement and restoration evaluation, they were initially developed for images acquired on land.Thus, some measures have been developed for subjective evaluation of underwater images.These include the Underwater Color Image Quality Evaluation (UCIQE) metric proposed by Yang and Sowmya [82] and the underwater image quality metric (UIQM) [83].However, the colorfulness parameter, which is utilized in this scenario, is incorporated into the UCIQE metric.The colorfulness and the color enhancement factor have high correlations with the subjective, human-based mean opinion score (MOS) [81].Additionally, according to Reference [82], "CIELAB color space-based subjective evaluation indicates that sharpness and colorful factors are highly correlated with the subjective image quality perception."Additionally, they note that "subjective quality metrics give the most reliable outcomes but are expensive, time-consuming, and impractical for real-time implementation and system integration" [82].These real-time and integration concerns are the very objectives of the realization of the proposed algorithm.
As noted by El Khoury et al. [84], no single algorithm can yield the best or worst performance for all images.In addition, they state that color and sharpness are always utilized in de-hazing evaluation.Thus, the use of colorfulness, gradient, and sharpness measures are justified.More than 50 underwater images from the literature and the Internet were used in the experiment.Some duplications were permitted due to the varying resolutions and formats of some images since this is another challenge of underwater image evaluation.The images used are shown in Figure 3a while the enhanced images are shown in Figure 3b.The original images were obtained from the website: https://github.com/IsaacChanghau/ImageEnhanceViaFusion [85].Even though objective metrics have been used for underwater image enhancement and restoration evaluation, they were initially developed for images acquired on land.Thus, some measures have been developed for subjective evaluation of underwater images.These include the Underwater Color Image Quality Evaluation (UCIQE) metric proposed by Yang and Sowmya [82] and the underwater image quality metric (UIQM) [83].However, the colorfulness parameter, which is utilized in this scenario, is incorporated into the UCIQE metric.The colorfulness and the color enhancement factor have high correlations with the subjective, human-based mean opinion score (MOS) [81].Additionally, according to Reference [82], "CIELAB color space-based subjective evaluation indicates that sharpness and colorful factors are highly correlated with the subjective image quality perception."Additionally, they note that "subjective quality metrics give the most reliable outcomes but are expensive, time-consuming, and impractical for real-time implementation and system integration" [82].These real-time and integration concerns are the very objectives of the realization of the proposed algorithm.
As noted by El Khoury et al. [84], no single algorithm can yield the best or worst performance for all images.In addition, they state that color and sharpness are always utilized in de-hazing evaluation.Thus, the use of colorfulness, gradient, and sharpness measures are justified.More than 50 underwater images from the literature and the Internet were used in the experiment.Some duplications were permitted due to the varying resolutions and formats of some images since this is another challenge of underwater image evaluation.The images used are shown in Figure 3a while the enhanced images are shown in Figure 3b.The original images were obtained from the website: https://github.com/IsaacChanghau/ImageEnhanceViaFusion[85].The high-boost filter configuration results in brighter images as opposed to using the high-pass filter configuration (HBFC).However, the contrast and edge enhancement using the high-boost filter version (HPFC) is slightly less when compared to the high-pass filter version.The global contrast operator was unmodified for all experiments.The PA performs well for a majority of the images except for images with non-overlapping color channel histograms.

Comparison with Other Algorithms
In Figure 4, we use visual results from Reference [86], which is appended with the visual results The high-boost filter configuration results in brighter images as opposed to using the high-pass filter configuration (HBFC).However, the contrast and edge enhancement using the high-boost filter version (HPFC) is slightly less when compared to the high-pass filter version.The global contrast operator was unmodified for all experiments.The PA performs well for a majority of the images except for images with non-overlapping color channel histograms.

Comparison with Other Algorithms
In Figure 4, we use visual results from Reference [86], which is appended with the visual results of PA using both HPFC and HBFC versions.For the video frames in Figure 4(1) from Reference [86], PA yields the sharpest detail, contrast enhancement, and color correction with the least color distortion when compared to the other methods apart from the method by Ancuti et al. [87] unlike Li et al. [88].In Figure 4(2), PA yields the sharpest and best color corrected image when compared to results obtained by Emberton et al. [86,89] and Drews-Jr et al. [90].This is verified in terms of the colorfulness metric used to compute the CEF, which closely correlates with the mean opinion score (MOS) based on the human subjective evaluation [81].PA yields the sharpest detail, contrast enhancement, and color correction with the least color distortion when compared to the other methods apart from the method by Ancuti et al. [87] unlike Li et al. [88].In Figure 4(2), PA yields the sharpest and best color corrected image when compared to results obtained by Emberton et al. [86,89] and Drews-Jr et al. [90].This is verified in terms of the colorfulness metric used to compute the CEF, which closely correlates with the mean opinion score (MOS) based on the human subjective evaluation [81].
(1) In Table 1, we verify the quality of the results of PA compared with the other algorithms using CEF, UCIQE, and UIQM.The bolded values indicate best results in all tables.The table shows that PA has the most consistent and highest values for the previously mentioned metrics.In Emberton et al. [86], it was claimed that the standard deviation of chroma, σc, correlates well with subjective evaluation.Thus, we replicate the table from Reference [86] and only select the σc and UCIQE values, which are also computed for the results obtained with PA.The results are shown in Table 1a.The In Table 1, we verify the quality of the results of PA compared with the other algorithms using CEF, UCIQE, and UIQM.The bolded values indicate best results in all tables.The table shows that PA has the most consistent and highest values for the previously mentioned metrics.In Emberton et al. [86], it was claimed that the standard deviation of chroma, σ c , correlates well with subjective evaluation.Thus, we replicate the table from Reference [86] and only select the σ c and UCIQE values, which are also computed for the results obtained with PA.The results are shown in Table 1a.The contrast of the luminance (not shown) is consistently highest for PA even though it is not correlated with subjective evaluation and was omitted.Once more, PA yields the most consistent high scores for both metrics.In Table 1b, PA yields the most colorful and enhanced image, which is shown by the CEF, UIQM, and UCIQE values.Though PA is not optimized as previous or contemporary methods, it is easily the fastest algorithm among the available implementations with acceptable visual outcomes.Furthermore, its multiple applications (alluded to in previous work) and real-time amenability compensates for its minimal impact on the underwater images with non-overlapping histograms.Furthermore, additional experiments showed that the HBFC and HPFC configurations work differently for various images.The percentiles of the global contrast enhancement function are fixed for all experiments to maintain consistency.Only the fractional order and gain are used to switch from HBFC to HPFC.However, the fractional order is kept constant to also maintain consistency.Thus, only the filter gain is the tuning parameter to switch between the two filter configurations.

Hazy Image Enhancement Results
We compare PA against available de-hazing implementations with respective parameters, which is stated by the authors in their papers [53].

•
Artificial Multi-Exposure-based Fusion (AMEF) by Galdran [54]: clip limit value, c as specified.The original method by He et al. took about a minimum of 30 s on the current platform for small images.Thus, we focused on the fast version using guided filters.We also present results and comparisons for hazy image contrast enhancement with algorithms from the literature using 53 real benchmark images employed in de-hazing experiments.In addition, the FRIDA3 dataset [92,93] consisting of left and right views of 66 synthetic images was also tested.The algorithms include Tarel and Hautiere [94], Dai et al. [95], Nishino et al. [96], He et al. [32], Galdran et al. [41], Wang and He [97], Zhu et al. [91], Ren et al. [46], Ju et al. [98], partial differential equation-based single scale Retinex GOC-CLAHE (PDE-GOC-SSR-CLAHE or PDE-Retinex) [52], PDE-IRCES [53], and PA.All test images (apart from FRIDA3) used were obtained from the University of Texas at Austin Laboratory for Image & Video Engineering (LIVE) dataset, which can be downloaded from the website: http: //live.ece.utexas.edu/research/fog/fade_defade.html)[78].More detailed and extensive results can be found in additional material [99,100].
Additionally, we present the numerical results for the available algorithm implementations compared with PA in Table 2.The bolded values indicate the best results.

Visual Comparison of AMEF and PA and Discussion
The key components of the enhancement capability of AMEF are the CLAHE and Gamma Correction (GC) algorithm.Unlike the PDE-GOC-CLAHE, which included the CLAHE and minimized its negative effects [52], the AMEF does not possess such features.We directly compare

Visual Comparison of AMEF and PA and Discussion
The key components of the enhancement capability of AMEF are the CLAHE and Gamma Correction (GC) algorithm.Unlike the PDE-GOC-CLAHE, which included the CLAHE and minimized its negative effects [52], the AMEF does not possess such features.We directly compare and present a sample of visual results of the state-of-the-art AMEF with PA in Figures 6-9.Based on visual observation, AMEF generally yields poor results without constant tuning of the clip limit, Figure 6c shows that PA can replicate the results of the AMEF by utilizing a high-boost filter with slightly better contrast than AMEF without CLAHE.Adding CLAHE to PA yields better results than AMEF with CLAHE.The AMEF de-hazing algorithm yields images with halos and color distortion similar to or worse than the CLAHE-based or Retinex-based de-hazing algorithms, which is seen in the Brick house image in Figure 7b.The AMEF is mainly suited to images with thick haze as seen in the Train image in Figure 7b even though there is color fading.The Horses image was processed using c = 0.03 for AMEF and PA was processed using both filter settings, which is shown in Figure 8.This is one of the images where AMEF performs adequately even though any slight increase in c leads to heavy color distortion.Increasing the clip limit of the CLAHE in the AMEF leads to increased color distortion.
Additionally, the AMEF algorithm is neither optimized nor adaptive and requires constant tuning of this clip limit parameter to obtain the best results for each hazy image.This makes the AMEF algorithm impractical for effective batch (or real-time image de-hazing) processing since The AMEF de-hazing algorithm yields images with halos and color distortion similar to or worse than the CLAHE-based or Retinex-based de-hazing algorithms, which is seen in the Brick house image in Figure 7b.The AMEF is mainly suited to images with thick haze as seen in the Train image in Figure 7b even though there is color fading.The Horses image was processed using c = 0.03 for AMEF and PA was processed using both filter settings, which is shown in Figure 8.This is one of the images where AMEF performs adequately even though any slight increase in c leads to heavy color distortion.Increasing the clip limit of the CLAHE in the AMEF leads to increased color distortion.
Additionally, the AMEF algorithm is neither optimized nor adaptive and requires constant tuning of this clip limit parameter to obtain the best results for each hazy image.This makes the AMEF algorithm impractical for effective batch (or real-time image de-hazing) processing since these issues were consistently observed using several benchmark hazy images.Ultimately, PA is much faster than the AMEF algorithm while yielding good enhancement results without halos, color degradation, or the need to constantly adjust parameters.Additionally, the AMEF is unable to enhance underwater images while PA effortlessly performs this operation, which is seen in Figure 9.

Conclusions
A fast, adaptive, and versatile multi-scale, fractional order-based hazy and underwater image enhancement algorithm with a relatively simplified structure suitable for hardware implementation has been proposed and developed.The earlier problems of the algorithm were addressed by automated balanced weighting of the filtered images used in the fusion process.The adherence to image entropy and standard deviation features coupled with global and local contrast enhancement ensures that visibility is greatly improved in the final result.Furthermore, comparisons with a recent state-of-the-art, multi-scale algorithm shows that the proposed approach is unmatched in several aspects such as speed, consistency, versatility, adaptability, and flexibility.Results show that the proposed scheme achieves the stated objectives and can be easily realized in hardware systems for fast image processing in challenging imaging environments.
Funding: This research received no external funding.

Figure 1 .
Figure 1.Proposed algorithm (PA) for enhancing hazy and underwater images.

Figure 1 .
Figure 1.Proposed algorithm (PA) for enhancing hazy and underwater images.
processed images using previous and improved configurations of PA (test images used were obtained from the University of Texas at Austin Laboratory for Image & Video Engineering (LIVE) dataset, which can be downloaded from the website: http://live.ece.utexas.edu/research/fog/fade_defade.html)[78].Note the elimination of the color distortion and reduced degree of noise enhancement for images in Figure2bwhen compared to Figure2a.

Figure 2 .
Figure 2. Processed images using (a) previous configuration (b) improved configuration of PA.

Figure 2 .
Figure 2. Processed images using (a) previous configuration (b) improved configuration of PA.

Figure 5 .
Figure 5. Runtime comparison of various algorithms using (a) 53 real and 66 synthetic (b) left and (c) right view hazy images.

Table 2 .
[53] ratio of visible edges, and saturation parameter values for images processed with He et al.[32], Zhu et al.[90], Ren et al., PDE-GOC-SSR-CLAHE[52], PDE-IRCES[53], and PA.indicate the best results.Only the method by Ren et al. is fully optimized for GPU computation with PA and other algorithms using parallel computation where possible.Results indicate that PA is the fastest algorithm of all the compared ones.Furthermore, the revised formulation combined with the RGB-IV does not increase run-time considerably except for images with very large dimensions.Nevertheless, the run-time is still much less than the algorithms by He et al., Zhu et al., and Ren et al.The revised scheme is also much easier to implement in digital hardware than the earlier version due to spatial filter structures. values