Next Article in Journal
Low-Complexity QRD-M with Path Eliminations in MIMO-OFDM Systems
Next Article in Special Issue
A Thermographic Measurement Approach to Assess Supercapacitor Electrical Performances
Previous Article in Journal
A Comparison between Horizontal and Vertical Interchannel Decorrelation
Previous Article in Special Issue
Pulsed Thermography Applied to the Study of Cultural Heritage
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Non-Reference Image Denoising Method for Infrared Thermal Image Based on Enhanced Dual-Tree Complex Wavelet Optimized by Fruit Fly Algorithm and Bilateral Filter

1
School of Mechatronic Engineering, China University of Mining & Technology, No.1 Daxue Road, Xuzhou 221116, China
2
Jiangsu Key Laboratory of Mine Mechanical and Electrical Equipment, China University of Mining & Technology, No. 1 Daxue Road, Xuzhou 221116, China
*
Authors to whom correspondence should be addressed.
Submission received: 19 October 2017 / Revised: 13 November 2017 / Accepted: 13 November 2017 / Published: 22 November 2017

Abstract

:
To eliminate the noise of infrared thermal image without reference and noise model, an improved dual-tree complex wavelet transform (DTCWT), optimized by an improved fruit-fly optimization algorithm (IFOA) and bilateral filter (BF), is proposed in this paper. Firstly, the noisy image is transformed by DTCWT, and the noise variance threshold is optimized by the IFOA, which is enhanced through a fly step range with inertia weight. Then, the denoised image will be re-processed using bilateral filter to improve the denoising performance and enhance the edge information. In the experiment, the proposed method is applied to eliminate both addictive noise and multiplicative noise, and the denoising results are compared with other representative methods, such as DTCWT, block-matching and 3D filtering (BM3D), median filter, wiener filter, wavelet decomposition filter (WDF) and bilateral filter. Moreover, the proposed method is applied as pre-processing utilization for infrared thermal images in a coal mining working face.

1. Introduction

Image denoising is a very important part of image processing, and it is the basis of all subsequent processing, such as image recognition [1], object detection and tracking [2]. The denoising methods are mainly divided into spatial domain and transform domain filters. Spatial domain filters have high processing speed since they operate the pixels on the original images, such as median filter [3], wiener filter [4] and bilateral filter (BF) [5]. These algorithms can remove noise quickly, but the denoising abilities are mediocre. Transform domain filters are more robust since the images are transformed into another domain. Two major transform domain filters are Fourier transform (FT) [6] and wavelet transform (WT) [7]. WT overcomes the problem of FT that the window size is fixed, and does not change with frequency. Thus, WT gets better effect than FT. Discrete wavelet transform (DWT), proposed by Shensa [8], is widely applied in the image processing field. Dual-tree complex wavelet transform (DTCWT) proposed by N. Kingsbury in 1998 [9] can improve the DWT in several aspects. It owns three additional properties: approximate shift invariance, good directional selectivity, and limited data redundancy [10]. Recently, DTCWT has been used to eliminate the noise of images, and the method has been improved by many researchers. Shahdoosti et al. [11] proposed a denoising method using the dual contourlet transform (DCT), which combined the directional filter bank (DFB) and DTCWT. The bivariate shrinkage function (BSF) with local variance estimation was adopted in [12] to process wavelet coefficients of DTCWT and the denoised images were reconstructed by the shrunk wavelet coefficients. Through changing the window size and shape of DTCWT, Liu et al. [13] acquired good quality of denoised images. However, these methods cannot reach the best denoising performance since the parameters are fixed. Therefore, intelligent methods that search the best parameters through the optimization algorithm have been proposed by some researchers.
In recent years, optimization problems are more and more extensive and many meta-heuristic algorithms has been proposed, such as genetic algorithm [14], ant colony algorithm [15], fish swarm algorithm [16], particle swarm algorithm [17], Virus colony search (VCS) algorithm [18], fire-fly algorithm [19], and so on. Among these algorithms, fruit-fly optimization algorithm (FOA), proposed by Professor Pan Wenchao in 2011 [20], has the advantages of simple structure, easy parameter setting, and fast convergence rate. Thus, FOA is widely used in economic management [21], forecasting [22], planning [23], etc. Nevertheless, similar to other optimal evolutionary algorithms, the original FOA may sink into local optima due to fixed fly step size. Meanwhile, the fitness function is a key problem that affects the denoising performance. Some researchers have put forward some examples. However, Most of them are full-reference. In practice, reference images cannot be acquired and these methods are not capable for practical usage.
Moreover, the original thermal images cannot be used for analyzing because of the low sharpness, noise problem and inhomogeneous. Therefore, it is a significant task to create a method for thermal images that is capable of acquiring best image denoised performance without reference images. In this work, a non-reference image denoising method based on enhanced dual-tree complex wavelet optimized by fruit fly algorithm and bilateral filter (IFOA-DTCWT-BF) is proposed and the method is compared with DTCWT and other denoising algorithms. Some common objective indexes are used to analyze the performance of these algorithms.
The rest of this paper is organized as follows. In Section 2, some related works are outlined based on the literature. In Section 3, the basic theories of DTCWT, FOA and bilateral filter are presented. In Section 4, the image denoising method based on IFOA-DTCWT-BF is elaborated, and a fly step range improved by inertia weight is performed. In Section 5, experiments are conducted and the proposed denoising method is tested, and an application is performed. Some conclusions and outlooks are summarized in Section 6.

2. Literature Review

This paper is mainly concerned with three research streams: image denoising method, dual-tree wavelet transform, and fruit fly optimization algorithm. In this section, we will summarize the relevant literature.

2.1. Image Denoising Method

In recent years, image denoising becomes an urgent need in image processing field. Many creative methods are put forward to wipe out the image noise. HR. Shahdoosti and SM. Hazavei [11] applied the hidden Markov tree (HMT) with mixtures of one-sided exponential densities to denoise images. To excavate the spatial information in hyperspectral image, a novel denoising method which integrates superpixel segmentation (SS) into low-rank representation (LRR), is presented in [24]. Because of the advantage in preserving texture in image, total variation is widely applied in the problems of image denoising [25]. Jia Zhi-Gang et al. [26] proposed a new TV-Stokes model for image deblurring with a good geometry explanation. In reality, the noise models of noisy images are unknown beforehand. Zhu et al. [27] consider the problem and proposes a novel blind image denoising algorithm to recover the clean image from noisy one with the unknown noise model. Kernel greedy algorithm for robust denoising (KGARD), put forward in [28], employs sparse modeling arguments to explicitly model and estimate the outliers, adopting a greedy approach. Recently, intelligent algorithms are applied in the image denoising field. Zhang Lin et al. [29] proposed a non-model dual-tree wavelet thresholding for image denoising that uses the improved chaotic drosophila optimization algorithm (CDOA) to estimate the noise variance.

2.2. Dual-Tree Complex Wavelet Transform

DTCWT is suitable for image processing applications [9]. In [30], DTCWT is used to create an effective contrast enhancement method and achieves a good performance in low illumination image enhancement. A new method for the interpolation of a full high-definition (HD) image which applied DTCWT to decompose the low-resolution image into different sub-bands was presented in [31]. To acquire good performance, the Wiener filter is used to remove the noise component of the dual-tree complex wavelet transform (DTCWT) coefficients [32]. In [33], the DTCWT coefficients are preprocessed by median filtering. Zhao et al. [34] presented a novel biometric image encryption algorithm based on DTCWT and compressed sensing. To produce better fusion effects, Q-shift DTCWT was proposed to simplify the construction of filters in DTCWT [35].

2.3. Fruit Fly Optimization Algorithm

FOA is easy to understand and has lower computational complexity. That is why FOA can be successfully applied in many areas such as the financial distress model solving [20], annual power load forecasting [36], the multidimensional knapsack problem solving [22], tuning of PID controller [37], neural network training [38] and so on. In [39], W Lijuan et al. proposed a novel approach, namely SFOASVR, which hybridizes SVR model with FOA and the seasonal index adjustment to forecast monthly tourist flow. Niu jinwei et al. [40] applied FOA to optimize the operation of the Texaco gasification process. However, the original FOA also has the possibility of falling into local extreme due to the fixed fly range and limited search space. To cover these problems, many improved methods have been proposed. Niu et al. [41] divided fruit group of FOA into two parts and used Cauchy mutation process to make fruit fly individuals variant. Xu et al. [42] improved the fly distance range (FR) of the FOA by using normal distribution F R ~ N ( 0 , L 2 ) . A novel multi scale cooperative mutation FOA was proposed by Zhang et al. [43] to tackle the limitation of local optimum. In [22], the parallel search is employed to balance exploitation and exploration of FOA. Wang Lin et al. [44] applied a level probability policy and a new mutation parameter to balance the population diversity and stability of FOA.

2.4. Discussion

Although the above research works are valuable, there are still some drawbacks remained. First, the denoising methods were improved by optimization algorithms and acquired better performances in recent decades, at the same time they were dragged down by the great amount of calculation and slow convergence rate of the optimization algorithms. Second, improved denoising algorithms always need reference images to get good performances, while, in practice, reference images do not exist. Third, despite the fine iteration rate and easy encoding, the original FOA is still likely to fall into local extreme due to the settled parameters. Moreover, the images processed by DTCWT always have small speckles which influence the quality of denoised images.
Therefore, an adaptive dual-tree complex wavelet threshold denoising method for non-reference image denoising based on improved fruit fly algorithm and bilateral filter is proposed in this paper, and a series of simulations and applications prove the effectiveness of the method.

3. Basic Theory

3.1. Dual-Tree Complex Wavelet Transform Based on Bivariate Shrinkage Function

When processing image with DTCWT, the image is resolved into a desired level by two separable 2D DWT branches, branch a (real part) and branch b (imagery part), whose filters are explicitly designed to meet the Hilbert pair requirement.
As shown in Figure 1, the h 0 ( n ) ,   h 1 ( n ) are the low pass filter and high pass filter of the filter bank above which produce the scale coefficient and wavelet coefficient of the real part, while the g 0 ( n ) ,   g 1 ( n ) are the low pass filter and high pass filter of the filter bank below that create the scale coefficient and wavelet coefficient of the imaginary part. In the first stage, the image has a line transformation through filter { g 0 ( n ) , g 1 ( n ) } , then it has a rank transformation through filter { h 0 ( n ) ,   h 1 ( n ) } . In the second stage, the image has a line transformation through filter { h 0 ( n ) ,   h 1 ( n ) } , then has a rank transformation through filter { g 0 ( n ) ,   g 1 ( n ) } . The decomposition at each stage produces one low-pass sub band LL and three high-pass sub bands HL, LH, and HH. The sum or difference of each pair of sub bands forms low-frequency wavelet coefficients and six wavelet transform coefficients in six directions. Then, bivariate shrinkage function is applied to process the wavelet coefficients.
Bivariate shrinkage function was proposed by Sendur L and Selesnick I W [45] in 2002. The DTCWT with bivariate shrinkage function can be described as follows.
Assuming that the image is contaminated by noise, the noised image can be denoted as:
f = x n ,
where the f is a noisy image; is the operator represents corresponding to different noise models, e.g., Additive noise is +, and multiplicative noise is × ; x is the original image; and n is the noise. It is transformed into dual-tree complex wavelet domain:
y = ω + n ,
where y stands for the wavelet coefficient of the observed image, ω represents the wavelet coefficient of the original image (to be estimated), and n is the noise wavelet coefficient which is independent and identically distributed and obeys the N(0, σ 2 ).
Using the maximum a posteriori probability estimate (MAP), the ω is estimated from the noise factor y:
ω ^ ( y ) = a r g   max ω   p ω | y ( ω | y ) ,
where p ω | y ( ω | y ) is the posterior probability function.
According to Bayesian estimation theory, Formula (3) can be written:
ω ^ ( y ) = a r g   max ω [ p n ( n ) p ω ( ω ) ] ,
p n ( n ) = 1 2 π σ 2 e x p ( n 1 2 + n 2 2 2 σ 2 ) ,
p ω ( ω ) is the bivariate probability distribution function that proposed in [45]:
p ω ( ω ) = 3 2 π σ ω 2 · e x p ( 3 σ ω 1 2 + ω 2 2 ) ,
ω 1 is the wavelet coefficient, ω 2 is the coefficient of the parent node at the next level of ω 1 , and σ ω 2 is the variance of the signal and is related to the location of the wavelet coefficient. The maximum posterior probability estimate of ω 1 can be obtained:
ω 1 ^ = ( y 1 2 + y 2 2 3 σ 2 σ ω ) + y 1 2 + y 2 2 y 1 ,
( g ) + = { 0 ,   g < 0 g ,   o t h e r w i s e ,
This is the bivariate shrinkage function.
The rectangular window region N(k) centered on the first k wavelet coefficients is defined. The signal variance σ ω 2 of this region is estimated by the wavelet coefficients in the N(k) region, and the noise variance σ 2 is estimated by the median filtering method:
σ 2 ^ = M e d i a n ( | y i | ) 0.6745 ,   y i s u b b a n d   H H ,
where subband HH is an aggregate of wavelet coefficients in the N(k) region.
From the observation model, we can obtain: σ y 2 = σ ω 2 + σ 2 , the empirical estimation of σ y 2 is:
σ y 2 ^ = 1 M y i N ( k ) y i 2 ,
M is the size of the rectangle region N(k), and the σ ω is estimated to be:
σ ω ^ = ( σ y 2 ^ σ 2 ^ ) + ,
After the noisy image processed by DTCWT, the wavelet coefficients of DTCWT is handled with the local adaptive bivariate shrinkage model which has local variance estimation, and then the denoised image can be obtained via inverse DTCWT of new wavelet coefficients.

3.2. Fruit Fly Optimization Algorithm

Fruit fly Optimization Algorithm (FOA), similar to other biomimetic interactive evolutionary computation algorithm, simulates the foraging process of fruit fly. The procedure of basic FOA can be summarized as follows and the schematic diagram is showed in Figure 2.
Step 1.1: Initialize the parameters: population size PS, maximum iteration number I N m a x , location range LR and the fly distance range FR.
The initial position of the fruit fly group is represented as follows:
{ X _ a x i s = r a n d ( L R ) Y _ a x i s = r a n d ( L R ) ,
Step 1.2: An individual fruit fly uses its olfactory sensation to find food, the random searching direction and distance is given as follows:
{ X i = X _ a x i s + r a n d ( F R ) Y i = Y _ a x i s + r a n d ( F R )
Step 1.3: Because of the unknown location of the food, the distance to the origin of coordinates ( D i s t i ) and the judgment value of smell concentration ( S i ) are calculated as follows:
D i s t i = X i 2 + Y i 2 ,   S i = 1 / D i s t i
Step 1.4: The fruit fly with maximum smell concentration among the fruit fly group searches food according to the smell concentration judgment function (or called fitness function), which represents as follows:
S m e l l ( i ) = f u n c t i o n ( S i ) ,   [ b e s t s m e l l ,   b e s t i n d e x ] = m a x ( S m e l l ( i ) ) ,
where bestsmell donates the maximal smell concentration, bestindex is the corresponding fruit fly number and Smell is the smell concentration of the group.
Step 1.5: The new smell concentration is compared with the previous best smell concentration; if it is better than the last generation, the current best location and smell concentration can be presented as follows:
s m e l l b e s t = b e s t s m e l l ,   { X _ a x i s = X ( b e s t i n d e x ) Y _ a x i s = Y ( b e s t i n d e x ) ,
Otherwise, repeated Step 1.2 to Step 1.4 to find other better smell concentrations.
Step 1.6: When the smell concentration reaches the preset precision value or the iteration number reaches the maximal IN, the circulation stops. Otherwise, repeat the process from Step 1.2 to Step 1.5.

3.3. Bilateral Filter Algorithm

As a kind of nonlinear spatial filter, BF is able to achieve uniform spatial filtering and reserve the edge information of an image. BF considers both the grey value of the point in the neighborhood scope, and the distance from center of geometric. Thus, the formula of grey value after filtering is presented as follows:
h x ( x ) = k 1 f ( ξ ) c ( ξ x ) s ( f ( ξ ) f ( x ) ) d ξ ,
where k is the normalization coefficient, the expression is as follows:
k ( x ) = c ( ξ x ) s ( f ( ξ ) f ( x ) ) d ξ ,
in which h x ( · ) is the corresponding point of the grey value before filtering, and x is the value after filtering. c ( · ) is the space similarity value of the center point with its neighborhood points, s ( · ) is the gray similarity value of the center point with its neighborhood points.
The functions c ( · ) and s ( · ) can be implemented using Gaussian function, which can be defined as follows:
c ( ξ x ) = e 1 2 ( ξ x σ ) 2 ,
s ( f ( ξ ) f ( x ) ) = e 1 2 ( f ( ξ ) f ( x ) σ ) 2 ,

4. The Proposed Method

In this section, the architecture of proposed method is presented and the flow of the proposed method is described as well. To enhance the threshold optimization ability of proposed method, an improved FOA is proposed.

4.1. Architecture of the Proposed Method

This algorithm is a new method that deals with noisy images without reference images. The architecture of algorithm can be divided into two parts. In the first part, DTCWT enhanced by IFOA (IFOA-DTCWT) is applied to remove noise. Since there are no reference images, the denoised images can only compare with the noisy images to estimate the denoising performance. The IFOA searches the noise variance randomly, and a new wavelet coefficient is gained through the bivariate shrinkage function with the noise variance. Then, the denoised image is obtained through the inverse DTCWT of the new wavelet coefficient. The process continues until the optimized denoised image is acquired. After being denoised by IFOA-DTCWT, many processing areas of the square shape are remained in the denoised image. Then in the second part, bilateral filter is applied to smooth the area and enhance the edge information of the image at the same time. The combination of the two methods acquired good performance. The architecture is shown in Figure 3.

4.2. Improved FOA

The fly step range affects convergence speed and search ability of fruit fly optimization algorithm. The global search ability of the algorithm will be obvious improved when fly step range is large enough. The problem is that its speed will slow down. In the condition of small fly step range, the algorithm convergences quickly but it will be easy to fall into local optimum.
In this article, an inertia weight is applied for modification of fly step range to balance the global search ability and search speed of FOA. While in initial stage of the iteration, large inertia weight is applied to make FOA have large step range. In the late stage of the iteration, the inertia weight diminishes and makes the step range small. In this way, FOA is intended to have better global search ability and be more precise in local development. Some inertia weight functions are presented as follows:
ω ( g ) = 1 ( 1 × ( g 1 ) ) / g max ,
ω ( g ) = 1 ( ( g 1 ) / g max ) 2 ,
ω ( g ) = cos ( π 2 ( g 1 ) / g max ) ,
where g is the current number of iteration, and g m a x is the maximum iteration number. The three inertia weights above are applied to improve FOA and formed IFOA1, IFOA2 and IFOA3. The flowchart of improved FOA is presented in Figure 4.

4.3. Flow of the Proposed Method

The proposed IFOA-DTCWT-BF can be summarized as follows:
Step 2.1: The noisy image is first transformed into a gray image. The parameters of IFOA are initialized, like population size PS, maximum iteration number I N m a x , location range LR and the fly step range FR.
Step 2.2: The location of each fly individual is gained through the fly group and FR. The distance and smell concentration of each fly individual are calculated according to Equation (14). Each smell concentration is considered as a potential threshold. A new threshold is just a new noise variance. Then, the noisy image is denoised through DTCWT and bivariate shrinkage function of the new variance. To evaluate the performance, the fitness function f is calculated as follows:
f = 1 / 10 l o g 255 2 M S E d B ,
where MSE is the mean square error between the noisy image and denoised image with size m × n , and it can be calculated as follows:
M S E = 1 m × n x = 1 m y = 1 n [ f ( x , y ) f ^ ( x , y ) ] ,
where f ( x , y ) is the noisy image before processing, f ^ ( x , y ) is the denoised image using proposed IFOA-DTCWT. The fitness function f reaches the maximal when MSE is the maximum.
Step 2.3: Fruit fly with maximal fitness is selected as bestsmell. If the present bestsmell is bigger than the former value, smellbest, the corresponding coordinates are updated. Otherwise, smellbest, X _ a x i s and Y _ a x i s are reserved.
Step 2.4: If the ending conditions are researched, smellbest, X _ a x i s and Y _ a x i s are treated as the optimum. Otherwise, Steps 2.2 and 2.3 are repeated.
Step 2.5: The optimal threshold is regarded as the noise variance of noisy image and the noisy image is denoised by the DTCWT based on bivariate shrinkage function.
Step 2.6: After being denoised by the DTCWT, the image will be processed with bilateral filter. Then the final denoised image is obtained.
The Pseudocode code of the proposed method is demonstrated as Algorithm 1:
Algorithm 1: IFOA-DTCWT-BF
  Inputs: Noisy image y
  Outputs: Denoised image y*
  % Initialization
  Set the parameters of the IFOA: population size PS, maximum iteration number INmax, location range LR, the fly step range FR and inertia weight ω(g).
  % Randomly search the noise variance σ
  For i = 1, 2, …, PS
   { X i = X _ a x i s + r a n d ( F R ) Y i = Y _ a x i s + r a n d ( F R )
   σ ( i ) = ( 1 / X ( i ) + Y ( i ) )
  % obtain denoised image through DTCWT
  K=DTCWT(y, σ(i))
  % obtain smell concentration
  Smell(i) = Objective function(y, K) through Equations (24) and (25)
  End
  % Set individual with best smell concentration
  Sbest max i = 1 , 2 P S ( S m e l l ( i ) )
  For g = 1: INmax
  For i = 1, 2, … PS
   { X i = X _ a x i s + r a n d ( ω ( g ) F R ) Y i = Y _ a x i s + r a n d ( ω ( g ) F R )
   σ ( i ) = ( 1 / X ( i ) + Y ( i ) )
  % obtain denoised image through DTCWT
  K=DTCWT(y,σ(i))
  % obtain smell concentration
  Smell(i) = Objective function(y, K)through Equations (24) and (25)
  End
  bestS max i = 1 , 2 P S ( S m e l l ( i ) )
  If bestS >Sbest than Sbest = bestS; X_axis = X(bestindex);Y_axis = Y(bestindex);
  Xbest(g) = X_axis;
  Ybest(g) = Y_axis;
   σ b e s t = ( 1 / X b e s t ( g ) + Y b e s t ( g ) )
  End
  y = DTCWT(y,σbest)
  % obtain final denoised image through bilateral filter
  y* = bilateral filter(y)

5. Experiment and Application

In this part, the improved FOA was tested and the proposed method IFOA-DTCWT-BF was used to remove noise of the images which were polluted by Gaussian noise and Speckle noise. Moreover, the denoising method was applied to denoise the real images from the surface experiment. To test the performance of improved fruit fly algorithm, a Schaffer function that has only one minimum point min f(0, 0) = 1 was used to test the search ability of algorithms. The Schaffer function is always used to test the algorithms’ performance, since there are many local extremum points around the minimum point, forming two loop ditches. The optimization algorithms are easy to fall into the local minimum point. In addition, the typical image denoising methods, e.g., DTCWT, block-matching and 3D filtering (BM3D), median filtering, wiener filtering, wavelet decomposition filtering (DWF) and bilateral filtering (BF), were applied to compare with the proposed method. Experiments and algorithm implementation were performed using Matlab, the experiment configuration was showed in Table 1, and the testing gray images were showed in Figure 5.
Some key parameters of IFOA are initialized: maximum number of generation is 30, population size is 20, location range is [−10, 10], etc.

5.1. Algorithm Performance Testing

To compare the improved FOA with previous FOA and test the performances of three different inertia weights, the four algorithms were tested by searching the minimum points of Schaffer function to compare their performances. The Schaffer function has only one minimum point min (0, 0) = −1, and many local extremum points around the minimum point form two loop ditches. The function values are 0.990283 and 0.962776 respectively, so it is easy to fall into these tiny points.
Some parameters of FOA: maximum number of generation was 90, population size was 20, the initial location range was [−1.0, 1.0], fly step range was [−2.0, 2.0]. The three IFOAs changed the original fly step range with different inertia weights. Because of the inertia weight, the fly step range of improved FOA could be set two times larger than normal step or more. The improved fly range can be expressed as ω ( g ) [−5.0, 5.0]. In Figure 6, IFOA had faster speed in the early stage of the global search. Then, in the later local optimization, IFOA had higher precision of search capabilities and found the extremum more accurately.
Figure 6 demonstrates that improved FOAs were faster than original FOA in search for extremum at the beginning. When in the later of the process, the local search became more detailed to find accurate extremum. Moreover, different inertia weights had different performances. The IFOA1 and IFOA2 searched quickly, and the IFOA3 was getting faster and faster and surpassed the two algorithms at the fifth generation. In Figure 7, Figure 8, Figure 9 and Figure 10, the search route of IFOA was much wider than original FOA. The IFOA2 and IFOA3 had almost same search route. Overall, the fly step improved by inertia weights had good performance.
Then, the FOA and IFOAs were applied in denoising algorithm. Firstly adding Gaussian noise on the original image formed the noisy image. It could be described as g(x, y) = f (x, y) + v (x, y). To simulate the Gaussian noise on the original image, the test image Lena.bmp (size: 512 × 512) was polluted with zero-mean Gaussian noise with the noise variance σ = 0.1. To compare the optimization ability of the FOA and IFOAs, four algorithms were applied to optimize the threshold of DTCWT to remove Gaussian noise on Lena image. The window size of the DTCWT was 5 × 5 which was convenient for bilateral filter to remove the small noise in the small areas. Some parameters of FOA include: maximum number of generation was 30, population size was 20, the flies’ initial location range was [−10.0, 10.0], and fly step range of [−1.0, 1.0]. IFOA parameters include: the fly step length of ω ( g ) [−2.0, 2.0], and other parameters were same as FOA.
In Figure 11, the IFOAs convergences faster than FOA, and searched more meticulous than FOA. It showed that IFOA possessed the ability of fast global search and precise local search. The IFOA2 with the inertia weight of Formula (22) had the best performance in searching the optimal threshold of the denoising method. As a result, we used IFOA2 to improve the proposed method.

5.2. Experiment with Gaussian Noise

To compare denoising performance, Gaussian noise with different noise variance was added to different original images, and other denoising methods such as DTCWT, median filter, wiener filter, bilateral filter, and wavelet decomposition filter were used to remove the noise of the image. The image quality evaluation strategies were peak signal to noise ratio (PSNR) and two-dimensional correlation coefficient of image matrix r. The PSNR and r here were calculated from denoised images and original images. Denoising results are listed in Table 2, Table 3 and Table 4.
The image denoising results showed that IFOA-DTCWT-BF had powerful performance dealing with Gaussian noise. In Table 3, we could see that although DTCWT acquired good data of PSNR, the r was not always better than the proposed method. Moreover, the proposed method got clearer edge information. Median filter, Wiener filter and WDF had general performance. BM3D was incapable of dealing with heavy noise. Bilateral filter was not good at denoising Gaussian noise and acquired worst results, while, in the proposed method, it played an important role, as shown in Figure 12. The edge information of denoised Lena images processed by DTCWT and the proposed method is shown in Figure 13. The denoised images are shown in Figure 14, Figure 15 and Figure 16.

5.3. Experiment with Speckle Noise

Speckle noise, which is known as multiplicative noise, added to original images forms noisy images. It can be described as: g(x, y) = f(x, yn(x, y), in which n(x,y) is the multiplicative noise. The original images were polluted by Speckle noise using Matlab function IMNOISE in the experiment. Different noise variances were set to pollute the original image, and different methods were utilized to compare the denoising performances. Denoising results are listed in Table 5, Table 6 and Table 7.
As shown in Table 5, Table 6 and Table 7, when the noise variance of Speckle noise was 0.01, the DTCWT, BM3D, Wiener and bilateral filter did good jobs. However, the proposed method performed much better than other methods when removing the Speckle noise of big noise variances. The results indicated that the proposed method was insensitive to the Speckle noise. In other words, it can be adopted in any condition. All denoised images are shown in Figure 17, Figure 18 and Figure 19.

5.4. Application

At present, the video monitoring of coal winning machine under the mine is getting more and more attention. However, using the camera for video monitoring is challenged by low illumination and coal dust environment in the mining face. When cutting the coal wall, the cutting unit of the mining machine produces great heat and radiates infrared ray. Thus, in 2016, the infrared thermal imager was applied to monitor the coal winning machine in Zhangjiakou Coal Mining Machinery Co., Ltd., Zhangjiakou, China. The infrared thermal imager taking into use in mining face could acquire clearer images than camera, as shown in Figure 20C, and it was more likely to realize the video monitoring and tracking of coal mining machine.
However, the infrared thermal images have same unavoidable weaknesses. The low sharpness, noise problem and inhomogeneous are three major problems of the thermal images. Because of the low detectivity and spatial resolution of the thermal imaging system, the thermal images owns blurred visual effect. The random disturbances of the external environment like lamp light, dust, etc. cause the noise problem and the inconsistent of the response characteristics of each pixel in infrared detector results in the inhomogeneous of the thermal images. Before the infrared thermal image can be used for tracking coal mining machine, it needs to be denoised. In this work, the original infrared thermal images were denoised by the proposed method and other methods, and the denoised images are shown in Figure 21.
As the images were denoised, the processed images would have less information than before. Therefore, information entropy was used to represent the amount of information contained in the image: the greater the value of information entropy, the more information the image contains. As the noise in the images was wiped by denoising methods, the denoised images were smoother than before. The information entropy (Ie) and smoothness (Smo) were introduced to evaluate the images in Figure 21, which could be calculated as
I e = k = 1 K H k log H k .
S m o = 1 M × N i = 1 m j = 1 n S b e s t ( i , j )
where K was the sum of gray scales in an image, Hk is the proportion of the pixels of gray value k in an image, and S b e s t ( i , j ) was the average gray value of the pixels in the best active window.
The information entropy and smoothness of original images and processed images are plotted in Figure 22. As shown in Figure 22, we could find that the information entropies of images processed by different methods are almost the same. It showed the denoised images also had good quality and the major information of images was retained. In Figure 21 and Figure 23, the images processed by proposed method had best smoothness. The smoothness of the images processed by IFOA-DTCWT-BF was 24.07, 24.28 and 23.41, which could not be acquired by other methods. It was beneficial for the following work of cutting unit tracking, as there were less unnecessary textures. The analysis results indicated that the images denoised by IFOA-DTCW-BF were better than other methods, which had better image textures and quality.

6. Conclusions and Future Work

To eliminate noise from a noisy image without reference image and obtain good performance, this paper proposed a non-reference image denoising method based on enhanced dual-tree complex wavelet optimized by fruit fly algorithm and bilateral filter. In the absence of reference image, a hypothesis was proposed that the MSE of the denoising image and noisy image was supposed to increase as the noise is denoised. The IFOA searched the optimal threshold to gain best denoising performance. The bilateral filter was applied to process images in the late work, which can significantly improve the performance of images and enhance the edge information. To verify the performance of the proposed IFOA, it was compared with original FOA. Moreover, different inertia weights were tested to get the best inertia weight. As comparisons, six representative methods, DTCWT, BM3D, Median filter, Wiener filter, DWT and BF, were implemented. The results showed that the proposed denoising method outperformed other methods. Finally, the proposed IFOA-DTDWT-BF was applied to infrared thermal image of roller in running to test the actual effect.
This work is valuable for image processing, which provides a non-reference denoising method. However, there are also some shortcomings. On the one hand, the calculation duration is still a problem, although it is faster than other adaptive methods. On the other hand, the proposed algorithm adopts the MSE as the denoising evaluation index, which may have imperfections. To surmount these shortcomings, the authors intend to improve the method with more effective intelligence algorithms and form an evaluation function with several evaluation indexes to assess the quality of the denoised image. Besides, the authors plan to combine other efficient denoising methods with the proposed methods to create a better method and apply the method in the aspect of industrial vision in the future.

Acknowledgments

The authors would like to thank all of the reviewers for their constructive comments. The supports of National Natural Science Foundation of China (No. U1510117 and 51605477), National Key Basic Research Program of China (No. 2014CB046301), and the Priority Academic Program Development (PAPD) of Jiangsu Higher Education Institutions in carrying out this research are gratefully acknowledged.

Author Contributions

Zhongbin Wang and Yiwen Liu conceived and designed the experiments; Lin Zhang and Lei Si performed the experiments; Jing Xu and Chao Tan analyzed the data; and Yiwen Liu wrote the paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Fan, J.; Yi, H.; Xu, L.; Zhao, T. A Histogram-Based Denoising Algorithm in a Joint-Fourier Transform Correlator for Image Recognition. Photonic Optoelectron. 2012, 1–3. [Google Scholar] [CrossRef]
  2. Walia, G.S.; Raza, S.; Gupta, A.; Asthana, R.; Singh, K. A novel approach of multi-stage tracking for precise localization of target in video sequences. Expert Syst. Appl. Int. J. 2017, 78, 208–224. [Google Scholar] [CrossRef]
  3. Gan, S.; Wang, S.; Chen, Y.; Chen, X.; Xiang, K. Separation of simultaneous sources using a structural-oriented median filter in the flattened dimension. Comput. Geosci. 2016, 86, 46–54. [Google Scholar] [CrossRef]
  4. Bhamre, T.; Zhang, T.; Singer, A. Denoising and covariance estimation of single particle cryo-em images. J. Struct. Biol. 2016, 195, 72–81. [Google Scholar] [CrossRef] [PubMed]
  5. Dai, L.; Yuan, M.; Zhang, X. Speeding up the bilateral filter: A joint acceleration way. IEEE Trans. Image Process. 2016, 25, 2657–2672. [Google Scholar] [CrossRef] [PubMed]
  6. Ran, Q.; Yuan, L.; Zhao, T. Image encryption based on nonseparable fractional fourier transform and chaotic map. Opt. Commun. 2015, 348, 43–49. [Google Scholar] [CrossRef]
  7. Cui, S.S.; Chen, G. Terahertz digital holography image denoising using stationary wavelet transform. In Proceedings of the SPIE Conferences of the Photoelectronic Technology Committee of the Chinese Society of Astronautics, Changchun, China, 13–15 August 2014; Volume 9522. [Google Scholar]
  8. Shensa, M.J. The discrete wavelet transform: Wedding the a trous and mallat algorithms. IEEE Trans. Signal Process. 1992, 40, 2464–2482. [Google Scholar] [CrossRef]
  9. Kingsbury, N. The dual-tree complex wavelet transform: A new technique for shift invariance and directional filters. In Proceedings of the IEEE Digital Signal Processing Workshop, Bryce Canyon, UT, USA, 9–12 August 1998; pp. 2543–2560. [Google Scholar]
  10. Kingsbury, N. Complex wavelets for shift invariant analysis and filtering of signals. Appl. Comput. Harmonic Anal. 2001, 10, 234–253. [Google Scholar] [CrossRef]
  11. Shahdoosti, H.R.; Hazavei, S.M. Image denoising in dual contourlet domain using hidden markov tree models. Dig. Signal Process. 2017, 67, 17–29. [Google Scholar] [CrossRef]
  12. Gong, W.G.; Liu, X.Y.; Li, W.H.; Li, J.F. Local adaptive image denoising based on double-density dual-tree complex wavelet transform. Opt. Precis. Eng. 2009, 17, 1171–1180. [Google Scholar] [CrossRef]
  13. Liu, Y.; Jiang, M. Image denoising algorithm based on dtcwt and adaptive windows. In Proceedings of the 3rd IEEE International Conference on Computer Science and Information Technology, Chengdu, China, 9–11 July 2010; pp. 207–210. [Google Scholar]
  14. Lee, S.H.; Lee, J. Optimization of three-dimensional wings in ground effect using multiobjective genetic algorithm. J. Aircr. 2011, 48, 1633–1645. [Google Scholar] [CrossRef]
  15. Bououden, S.; Chadli, M.; Karimi, H.R. An ant colony optimization-based fuzzy predictive control approach for nonlinear processes. Inf. Sci. 2015, 299, 143–158. [Google Scholar] [CrossRef]
  16. Wang, H.B.; Fan, C.C.; Tu, X.Y. Afsaocp: A novel artificial fish swarm optimization algorithm aided by ocean current power. Appl. Intell. 2016, 45, 1–16. [Google Scholar] [CrossRef]
  17. Cheng, R.; Jin, Y. A social learning particle swarm optimization algorithm for scalable optimization. Inf. Sci. 2015, 291, 43–60. [Google Scholar] [CrossRef]
  18. Li, M.D.; Zhao, H.; Weng, X.W.; Han, T. A novel nature-inspired algorithm for optimization. Adv. Eng. Softw. 2016, 92, 65–88. [Google Scholar] [CrossRef]
  19. Yang, X.S. Nature-Inspired Metaheuristic Algorithms, 2nd ed.; Luniver Press: Beckington, UK, 2010; pp. 92–116. ISBN 9781905986286. [Google Scholar]
  20. Pan, W.T. A new fruit fly optimization algorithm: Taking the financial distress model as an example. Knowl. Based Syst. 2012, 26, 69–74. [Google Scholar] [CrossRef]
  21. Lin, S.M. Analysis of service satisfaction in web auction logistics service using a combination of fruit fly optimization algorithm and general regression neural network. Neural Comput. Appl. 2013, 22, 783–791. [Google Scholar] [CrossRef]
  22. Li, H.Z.; Guo, S.; Li, C.J.; Sun, J.Q. A hybrid annual power load forecasting model based on generalized regression neural network with fruit fly optimization algorithm. Knowl. Based Syst. 2013, 37, 378–387. [Google Scholar] [CrossRef]
  23. Wang, L.; Shi, Y.; Liu, S. An improved fruit fly optimization algorithm and its application to joint replenishment problems. Expert Syst. Appl. 2015, 42, 4310–4323. [Google Scholar] [CrossRef]
  24. Fan, F.; Ma, Y.; Li, C.; Mei, X.; Huang, J.; Ma, J. Hyperspectral image denoising with superpixel segmentation and low-rank representation. Inf. Sci. Int. J. 2017, 397–398, 48–68. [Google Scholar] [CrossRef]
  25. Wu, Q.; Li, Y.; Lin, Y. The application of nonlocal total variation in image denoising for mobile transmission. Multimedia Tools Appl. 2017, 76, 1–13. [Google Scholar] [CrossRef]
  26. Jia, Z.G.; Wei, M. A new TV-Stokes model for image deblurring and denoising with fast algorithms. J. Sci. Comput. 2017, 72, 1–20. [Google Scholar] [CrossRef]
  27. Zhu, F.; Chen, G.; Hao, J.; Heng, P.A. Blind image denoising via dependent dirichlet process tree. IEEE Trans. Pattern Anal. Mach. Intell. 2017, 39, 1518–1531. [Google Scholar] [CrossRef] [PubMed]
  28. Papageorgiou, G.; Bouboulis, P.; Theodoridis, S. Robust non-linear regression: A greedy approach employing kernels with application to image denoising. IEEE Trans. Signal Process. 2016, 99. [Google Scholar] [CrossRef]
  29. Zhang, L.; Zhou, X.; Wang, Z.; Tan, C.; Liu, X. A nonmodel dual-tree wavelet thresholding for image denoising through noise variance optimization based on improved chaotic drosophila algorithm. Int. J. Pattern Recognit. Artif. Intell. 2017, 31. [Google Scholar] [CrossRef]
  30. Jung, C.; Yang, Q.; Sun, T.; Fu, Q.; Song, H. Low light image enhancement with dual-tree complex wavelet transform. J. Vis. Commun. Image Represent. 2016, 42. [Google Scholar] [CrossRef]
  31. Lama, R.K.; Choi, M.R.; Kwon, G.R. Image interpolation for high-resolution display based on the complex dual-tree wavelet transform and hidden markov model. Multimedia Tools Appl. 2016, 75, 1–12. [Google Scholar] [CrossRef]
  32. Zhang, X. Image denoising using dual-tree complex wavelet transform and wiener filter with modified thresholding. J. Sci. Ind. Res. 2016, 75, 687–690. [Google Scholar]
  33. Wei, T.; Gao, Q.; Ma, N.; Li, N.; Wang, J.; Lei, P.; Ji, X. Feature-level image fusion through consistent region segmentation and dual-tree complex wavelet transform. J. Imaging Sci. Technol. 2016. [Google Scholar] [CrossRef]
  34. Zhao, Z.; Dong, J.; Li, H. A novel biometric image encryption algorithm based on compressed sensing and dual-tree complex wavelet transform. In Proceedings of the Eighth International Conference on Digital Image Processing, Chengdu, China, 20–22 May 2016. [Google Scholar]
  35. Tian, Y.; Luo, J.; Zhang, W.; Jia, T.; Wang, A.; Li, L. Multifocus image fusion in q-shift dtcwt domain using various fusion rules. Math. Probl. Eng. 2016, 2016, 1–12. [Google Scholar] [CrossRef]
  36. Meng, T.; Pan, Q.K. An improved fruit fly optimization algorithm for solving the multidimensional knapsack problem. Appl. Soft Comput. 2017, 50, 79–93. [Google Scholar] [CrossRef]
  37. Xu, C. An improved fruit fly optimization algorithm and its application of pid parameters tuning. J. Inf. Comput. Sci. 2015, 12, 3647–3654. [Google Scholar] [CrossRef]
  38. Yang, Y.; Luo, Z.; Xie, C. Study and application on fruit fly optimization algorithm optimized general regression neural network in mined-out area stability analysis. World Sci. Tech. R D 2015, 3, 230–234. [Google Scholar] [CrossRef]
  39. Wu, L.; Cao, G. Seasonal svr with foa algorithm for single-step and multi-step ahead forecasting in monthly inbound tourist flow. Knowl. Based Syst. 2016, 110, 157–166. [Google Scholar] [CrossRef]
  40. Niu, J.; Zhong, W.; Liang, Y.; Luo, N.; Qian, F. Fruit fly optimization algorithm based on differential evolution and its application on gasification process operation optimization. Knowl. Based Syst. 2015, 88, 253–263. [Google Scholar] [CrossRef]
  41. Niu, D.; Ma, T.; Liu, B. Power load forecasting by wavelet least squares support vector machine with improved fruit fly optimization algorithm. J. Comb. Optim. 2016, 33, 1–22. [Google Scholar] [CrossRef]
  42. Xu, J.; Wang, Z.; Tan, C.; Si, L.; Liu, X. A novel denoising method for an acoustic-based system through empirical mode decomposition and an improved fruit fly optimization algorithm. Appl. Sci. 2017, 7, 215. [Google Scholar] [CrossRef]
  43. Zhang, Y.; Cui, G.; Wu, J.; Pan, W.T.; He, Q. A novel multi-scale cooperative mutation fruit fly optimization algorithm. Knowl. Based Syst. 2016, 114, 24–35. [Google Scholar] [CrossRef]
  44. Wang, L.; Liu, R.; Liu, S. An effective and efficient fruit fly optimization algorithm with level probability policy and its applications. Knowl. Based Syst. 2016, 97, 158–174. [Google Scholar] [CrossRef]
  45. Şendur, L.; Selesnick, I.W. Bivariate shrinkage functions for wavelet-based denoising exploiting interscale dependency. IEEE Trans. Signal Process. 2002, 50, 2744–2756. [Google Scholar] [CrossRef]
Figure 1. Dual-tree complex wavelet transform: (A) analysis filter bank; and (B) synthesis filter bank.
Figure 1. Dual-tree complex wavelet transform: (A) analysis filter bank; and (B) synthesis filter bank.
Applsci 07 01190 g001
Figure 2. Fruit fly optimization algorithm.
Figure 2. Fruit fly optimization algorithm.
Applsci 07 01190 g002
Figure 3. Architecture of enhanced dual-tree complex wavelet optimized by fruit fly algorithm and bilateral filter (IFOA-DTCWT-BF).
Figure 3. Architecture of enhanced dual-tree complex wavelet optimized by fruit fly algorithm and bilateral filter (IFOA-DTCWT-BF).
Applsci 07 01190 g003
Figure 4. Flowchart of the improved fruit fly optimization algorithm.
Figure 4. Flowchart of the improved fruit fly optimization algorithm.
Applsci 07 01190 g004
Figure 5. Original testing gray images (size: 512 × 512). (A) Lena; (B) Columbia; (C) Peppers.
Figure 5. Original testing gray images (size: 512 × 512). (A) Lena; (B) Columbia; (C) Peppers.
Applsci 07 01190 g005
Figure 6. Optimization result of FOA and IFOAs with different inertia weights.
Figure 6. Optimization result of FOA and IFOAs with different inertia weights.
Applsci 07 01190 g006
Figure 7. (A) Optimization process of FOA; and (B) fruit fly flying route of FOA.
Figure 7. (A) Optimization process of FOA; and (B) fruit fly flying route of FOA.
Applsci 07 01190 g007
Figure 8. (A) Optimization process of FOA1; and (B) fruit fly flying route of FOA1.
Figure 8. (A) Optimization process of FOA1; and (B) fruit fly flying route of FOA1.
Applsci 07 01190 g008
Figure 9. (A) Optimization process of FOA2; and (B) fruit fly flying route of FOA2.
Figure 9. (A) Optimization process of FOA2; and (B) fruit fly flying route of FOA2.
Applsci 07 01190 g009
Figure 10. (A) Optimization process of FOA3; and (B) fruit fly flying route of FOA3.
Figure 10. (A) Optimization process of FOA3; and (B) fruit fly flying route of FOA3.
Applsci 07 01190 g010
Figure 11. Optimization result using FOA and IFOAs to remove Gaussian noise from noisy image (Lena).
Figure 11. Optimization result using FOA and IFOAs to remove Gaussian noise from noisy image (Lena).
Applsci 07 01190 g011
Figure 12. The effect of the bilateral filter: the left is no BF method and right is proposed method. (A) Lena of noise σ = 0.01; (B) Lena of noise σ = 0.05; (C) Lena of noise σ = 0.1; (D) Lena of noise σ = 0.2.
Figure 12. The effect of the bilateral filter: the left is no BF method and right is proposed method. (A) Lena of noise σ = 0.01; (B) Lena of noise σ = 0.05; (C) Lena of noise σ = 0.1; (D) Lena of noise σ = 0.2.
Applsci 07 01190 g012
Figure 13. The edge information comparison of DTCWT and proposed method: left image is denoised by DTCWT, and right is proposed method. (A) Part of lena of noise σ = 0.01; (B) Part of lena of noise σ = 0.05; (C) Part of lena of noise σ = 0.1; (D) Part of lena of noise σ = 0.2.
Figure 13. The edge information comparison of DTCWT and proposed method: left image is denoised by DTCWT, and right is proposed method. (A) Part of lena of noise σ = 0.01; (B) Part of lena of noise σ = 0.05; (C) Part of lena of noise σ = 0.1; (D) Part of lena of noise σ = 0.2.
Applsci 07 01190 g013
Figure 14. Denosing results of Lena with Gaussian noise using IFOA-DTCWT-BF. (A) Lena of noise σ   = 0.01; (B) Lena of noise σ = 0.05; (C) Lena of noise σ = 0.1; (D) Lena of noise σ = 0.2; (E) Denoised image of σ = 0.01; (F) Denoised image of σ = 0.05; (G) Denoised image of σ = 0.1; (H) Denoised image of σ = 0.2.
Figure 14. Denosing results of Lena with Gaussian noise using IFOA-DTCWT-BF. (A) Lena of noise σ   = 0.01; (B) Lena of noise σ = 0.05; (C) Lena of noise σ = 0.1; (D) Lena of noise σ = 0.2; (E) Denoised image of σ = 0.01; (F) Denoised image of σ = 0.05; (G) Denoised image of σ = 0.1; (H) Denoised image of σ = 0.2.
Applsci 07 01190 g014
Figure 15. Denosing results of Columbia with Gaussian noise using IFOA-DTCWT-BF. (A) Columbia of noise σ = 0.01; (B) Columbia of noise σ = 0.05; (C) Columbia of noise σ = 0.1; (D) Columbia of noise σ = 0.2; (E) Denoised image of σ = 0.01; (F) Denoised image of σ = 0.05; (G) Denoised image of σ = 0.1; (H) Denoised image of σ = 0.2.
Figure 15. Denosing results of Columbia with Gaussian noise using IFOA-DTCWT-BF. (A) Columbia of noise σ = 0.01; (B) Columbia of noise σ = 0.05; (C) Columbia of noise σ = 0.1; (D) Columbia of noise σ = 0.2; (E) Denoised image of σ = 0.01; (F) Denoised image of σ = 0.05; (G) Denoised image of σ = 0.1; (H) Denoised image of σ = 0.2.
Applsci 07 01190 g015
Figure 16. Denoising results of Peppers with Gaussian noise using IFOA-DTCWT-BF. (A) Peppers of noise σ = 0.01; (B) Peppers of noise σ = 0.05; (C) Peppers of noise σ = 0.1; (D) Peppers of noise σ = 0.2; (E) Denoised image of σ = 0.01; (F) Denoised image of σ = 0.05; (G) Denoised image of σ = 0.1; (H) Denoised image of σ = 0.2.
Figure 16. Denoising results of Peppers with Gaussian noise using IFOA-DTCWT-BF. (A) Peppers of noise σ = 0.01; (B) Peppers of noise σ = 0.05; (C) Peppers of noise σ = 0.1; (D) Peppers of noise σ = 0.2; (E) Denoised image of σ = 0.01; (F) Denoised image of σ = 0.05; (G) Denoised image of σ = 0.1; (H) Denoised image of σ = 0.2.
Applsci 07 01190 g016
Figure 17. Denoising results of Lena with Speckle noise using IFOA-DTCWT-BF. (A) Lena of noise σ = 0.01; (B) Lena of noise σ = 0.05; (C) Lena of noise σ = 0.1; (D) Lena of noise σ = 0.2; (E) Denoised image of σ = 0.01; (F) Denoised image of σ = 0.05; (G) Denoised image of σ = 0.1; (H) Denoised image of σ = 0.2.
Figure 17. Denoising results of Lena with Speckle noise using IFOA-DTCWT-BF. (A) Lena of noise σ = 0.01; (B) Lena of noise σ = 0.05; (C) Lena of noise σ = 0.1; (D) Lena of noise σ = 0.2; (E) Denoised image of σ = 0.01; (F) Denoised image of σ = 0.05; (G) Denoised image of σ = 0.1; (H) Denoised image of σ = 0.2.
Applsci 07 01190 g017
Figure 18. Denoising results of Columbia with Speckle noise using IFOA-DTCWT-BF. (A) Columbia of noise σ = 0.01; (B) Columbia of noise σ = 0.05; (C) Columbia of noise σ = 0.1; (D) Columbia of noise σ = 0.2; (E) Denoised image of σ = 0.01; (F) Denoised image of σ = 0.05; (G) Denoised image of σ = 0.1; (H) Denoised image of σ = 0.2.
Figure 18. Denoising results of Columbia with Speckle noise using IFOA-DTCWT-BF. (A) Columbia of noise σ = 0.01; (B) Columbia of noise σ = 0.05; (C) Columbia of noise σ = 0.1; (D) Columbia of noise σ = 0.2; (E) Denoised image of σ = 0.01; (F) Denoised image of σ = 0.05; (G) Denoised image of σ = 0.1; (H) Denoised image of σ = 0.2.
Applsci 07 01190 g018
Figure 19. Denoising results of Peppers with Speckle noise using IFOA-DTCWT-BF. (A) Peppers of noise σ = 0.01; (B) Peppers of noise σ = 0.05; (C) Peppers of noise σ = 0.1; (D) Peppers of noise σ = 0.2; (E) Denoised image of σ = 0.01; (F) Denoised image of σ = 0.05; (G) Denoised image of σ = 0.1; (H) Denoised image of σ = 0.2.
Figure 19. Denoising results of Peppers with Speckle noise using IFOA-DTCWT-BF. (A) Peppers of noise σ = 0.01; (B) Peppers of noise σ = 0.05; (C) Peppers of noise σ = 0.1; (D) Peppers of noise σ = 0.2; (E) Denoised image of σ = 0.01; (F) Denoised image of σ = 0.05; (G) Denoised image of σ = 0.1; (H) Denoised image of σ = 0.2.
Applsci 07 01190 g019
Figure 20. Original image form the surface experiment. (A) Image form surface mining engineering; (B) Cutting unit of mining machine; (C) Infrared thermal image of cutting unit.
Figure 20. Original image form the surface experiment. (A) Image form surface mining engineering; (B) Cutting unit of mining machine; (C) Infrared thermal image of cutting unit.
Applsci 07 01190 g020aApplsci 07 01190 g020b
Figure 21. Original images and denoised images using different methods. (A) Infrared thermal image of cutting unit from the side; (B) Infrared thermal image of cutting unit from close distance; (C) Infrared thermal image of cutting unit from the front.
Figure 21. Original images and denoised images using different methods. (A) Infrared thermal image of cutting unit from the side; (B) Infrared thermal image of cutting unit from close distance; (C) Infrared thermal image of cutting unit from the front.
Applsci 07 01190 g021
Figure 22. Information entropy of images using different methods.
Figure 22. Information entropy of images using different methods.
Applsci 07 01190 g022
Figure 23. Smoothness of images using different methods.
Figure 23. Smoothness of images using different methods.
Applsci 07 01190 g023
Table 1. Experiment configurations. RAM: random-access memory.
Table 1. Experiment configurations. RAM: random-access memory.
IndexConfiguration
ProcessorIntel Xeon E5506, 2.13 GHz
RAM8 G
SystemWin7, ×64
SoftwareMatlab
Table 2. Comparison of the five methods (Lena with Gaussian noise). WDF: wavelet decomposition filter; BM3D: block-matching and 3D filtering; DTCWT: dual-tree complex wavelet transform; IFOA-DTCWT-BF: enhanced dual-tree complex wavelet optimized by fruit fly algorithm and bilateral filter.
Table 2. Comparison of the five methods (Lena with Gaussian noise). WDF: wavelet decomposition filter; BM3D: block-matching and 3D filtering; DTCWT: dual-tree complex wavelet transform; IFOA-DTCWT-BF: enhanced dual-tree complex wavelet optimized by fruit fly algorithm and bilateral filter.
Parametersσ = 0.01σ = 0.05σ = 0.1σ = 0.2
Methods PSNRrPSNRrPSNRrPSNRr
Noised20.070.883713.640.644111.320.50539.5180.3731
Median27.860.976623.410.936220.930.894518.400.8281
Wiener27.700.972620.900.888118.660.810516.930.7062
WDF27.480.974523.940.941121.770.900819.570.8289
Bilateral25.440.960014.930.697111.980.53459.8630.3877
BM3D31.960.990916.020.741912.520.554310.260.3975
DTCWT31.010.988726.510.968724.300.952521.700.9229
IFOA-DTCWT-BF28.660.980726.340.968524.150.951221.570.9155
Table 3. Comparison of the five methods (Columbia with Gaussian noise).
Table 3. Comparison of the five methods (Columbia with Gaussian noise).
Parametersσ = 0.01σ = 0.05σ = 0.1σ = 0.2
Methods PSNRrPSNRrPSNRrPSNRr
Noised20.180.871413.940.626611.640.49159.6820.3583
Median27.220.968923.280.924221.090.881418.710.8127
Wiener26.950.967520.950.875118.580.792716.570.6828
WDF26.810.965723.570.929621.300.887418.820.8125
Bilateral25.120.950615.240.679612.310.520010.020.3719
BM3D31.400.988216.560.734912.910.541410.430.3824
DTCWT30.120.984125.510.957622.890.934520.130.8964
IFOA-DTCWT-BF27.700.972925.260.957422.870.937220.090.8973
Table 4. Comparison of the five methods (Peppers with Gaussian noise).
Table 4. Comparison of the five methods (Peppers with Gaussian noise).
Parametersσ = 0.01σ = 0.05σ = 0.1σ = 0.2
Methods PSNRrPSNRrPSNRrPSNRr
Noised20.130.903913.750.686911.430.54899.5760.4135
Median28.290.983223.670.952221.180.918118.510.8607
Wiener27.100.978020.850.907218.600.840616.710.7449
WDF26.930.977023.570.949821.460.918219.120.8554
Bilateral25.510.968615.050.738712.090.57929.9210.4294
BM3D31.670.992416.190.781212.610.598210.340.4420
DTCWT30.510.990025.940.972923.550.958421.000.9369
IFOA-DTCWT-BF28.120.983025.610.972223.490.959120.740.9279
Table 5. Comparison of the six methods (Lena with Speckle noise).
Table 5. Comparison of the six methods (Lena with Speckle noise).
Parametersσ = 0.01σ = 0.05σ = 0.1σ = 0.2
Methods PSNRrPSNRrPSNRrPSNRr
Noised25.680.963718.830.852015.990.752613.250.6226
Median29.440.983925.820.963423.660.941021.270.9030
Wiener31.490.989924.920.955322.220.918019.600.8547
WDF28.760.981026.980.971325.400.958623.290.9328
Bilateral31.510.990022.050.918017.750.810014.190.6610
BM3D33.400.993526.790.970419.790.868315.080.6978
DTCWT31.880.990825.830.963323.650.939221.280.8954
IFOA-DTCWT-BF29.140.982928.410.979627.390.974925.550.9633
Table 6. Comparison of the six methods (Columbia with Speckle noise).
Table 6. Comparison of the six methods (Columbia with Speckle noise).
Parametersσ = 0.01σ = 0.05σ = 0.1σ = 0.2
Methods PSNRrPSNRrPSNRrPSNRr
Noised28.570.978221.680.904618.750.831315.930.7250
Median29.270.981226.950.967825.350.953123.370.9260
Wiener32.090.989926.160.961523.340.928320.720.8742
WDF28.050.974327.220.968826.260.961024.930.9468
Bilateral31.200.987925.010.949521.090.886117.470.7777
BM3D33.080.992029.700.982524.130.939718.990.8292
DTCWT31.550.988825.100.951722.400.914219.950.8554
IFOA-DTCWT-BF28.210.975927.810.973527.290.970126.350.9630
Table 7. Comparison of the six methods (Peppers with Speckle noise).
Table 7. Comparison of the six methods (Peppers with Speckle noise).
Parametersσ = 0.01σ = 0.05σ = 0.1σ = 0.2
Methods PSNRrPSNRrPSNRrPSNRr
Noised25.760.971518.890.879516.090.796513.420.6788
Median30.140.989226.190.967823.810.954821.240.9212
Wiener31.510.992124.740.963222.030.932519.430.8792
WDF28.040.982226.520.974825.220.965923.180.9459
Bilateral31.950.992822.010.933817.810.846414.340.7147
BM3D33.120.994526.200.973419.680.891815.280.7518
DTCWT30.730.990524.400.960422.230.935120.130.8948
IFOA-DTCWT-BF28.660.984828.010.982327.140.978625.450.9703

Share and Cite

MDPI and ACS Style

Liu, Y.; Wang, Z.; Si, L.; Zhang, L.; Tan, C.; Xu, J. A Non-Reference Image Denoising Method for Infrared Thermal Image Based on Enhanced Dual-Tree Complex Wavelet Optimized by Fruit Fly Algorithm and Bilateral Filter. Appl. Sci. 2017, 7, 1190. https://0-doi-org.brum.beds.ac.uk/10.3390/app7111190

AMA Style

Liu Y, Wang Z, Si L, Zhang L, Tan C, Xu J. A Non-Reference Image Denoising Method for Infrared Thermal Image Based on Enhanced Dual-Tree Complex Wavelet Optimized by Fruit Fly Algorithm and Bilateral Filter. Applied Sciences. 2017; 7(11):1190. https://0-doi-org.brum.beds.ac.uk/10.3390/app7111190

Chicago/Turabian Style

Liu, Yiwen, Zhongbin Wang, Lei Si, Lin Zhang, Chao Tan, and Jing Xu. 2017. "A Non-Reference Image Denoising Method for Infrared Thermal Image Based on Enhanced Dual-Tree Complex Wavelet Optimized by Fruit Fly Algorithm and Bilateral Filter" Applied Sciences 7, no. 11: 1190. https://0-doi-org.brum.beds.ac.uk/10.3390/app7111190

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