Research Article  Open Access
Haihua Wang, Xinxin Zhang, Shuli Mei, "ShannonCosine Wavelet Precise Integration Method for Locust Slice Image Mixed Denoising", Mathematical Problems in Engineering, vol. 2020, Article ID 4989735, 17 pages, 2020. https://doi.org/10.1155/2020/4989735
ShannonCosine Wavelet Precise Integration Method for Locust Slice Image Mixed Denoising
Abstract
A novel denoising method for removing mixed noise from locust slice images is proposed by means of Shannoncosine wavelet and the nonlinear variational model for the image processing. This method includes two parts that are the sparse representation of the slice images and the novel numerical algorithm for solving the variation model on image denoising based on the sparse representation. In the first part, a parametric Shannoncosine wavelet function is introduced to construct the multiscale wavelet transform matrix, which is applied to represent the slice images sparsely by adjusting the parameters adaptively based on the texture of the locust slice images. By multiplying the matrix with the signal, the multiscale wavelet transform coefficients of the signal can be obtained at one time, which can be used to identify the saltandpepper noises in the slice images. This ensures that the saltandpepper noise points are kept away from the sparse representation of the slice images. In the second part, a semianalytical method on solving the system of the nonlinear differential equations is constructed based on the sparse representation of the slice images, which is named the sparse wavelet precise integration method (SWPIM). Substituting the sparse representation of the slice images into the Perona–Malik model which is the famous edgepreserving smoothing model for removing the Gaussian noises of the biomedical images, a system of nonlinear differential equations is obtained, whose scale is far smaller than the one obtained by the difference method. The numerical experiments show that both the values of SSIM and PSNR of the denoised locust slice images are better than the classical methods besides the algorithm efficiency.
1. Introduction
The main purpose of locust microscopy experiments is to reveal the working mechanism of locust reproduction and pesticide extermination by studying the structure of locust cells and their neural tissues [1]. In addition, it is helpful for the research and development of biological pesticides and reduces the use of chemical pesticides [2]. Locust slice images are quite complex such as textures, edges, and structures, and slice tissue structures cannot be acquired repeatedly [3]. The locust slice images have the characteristics of biomedical image segmentation smoothing [4]. The saltandpepper noise and Gaussian noise can not only significantly deteriorate the quality of locust slice images but also bring some difficulty to the subsequent image analysis such as image segmentation, edge detection, and classification. How to effectively remove saltandpepper noise and Gaussian noise in locust slice images has been an important research topic in the field of image processing.
During the past several decades, many techniques have been developed to restore the corrupted image, but research on locust slice images is lacking. Most traditional filtering methods are designed for a certain type of noise, and the effect of removing saltandpepper and Gaussian mixed noises is not ideal. Among the traditional denoising methods, the standard median filter [5] is one of the most popular nonlinear filters for the removal of saltandpepper noise in terms of its good denoising capability and computational efficiency, but high noise density will lead to unsatisfactory denoising results with serious loss of image details. Wiener filtering [6] is a popular type of Gaussian noise removing technique with good denoising performance. However, it cannot remove the saltandpepper noise well due to a larger mean square error with the original image. Mean filtering [7] is not very effective in removing saltandpepper noise, and its application range is limited. Gaussian filtering [8] will damage the edge information of the image when denoising, and the denoising performance is poor. During the denoising process, the wavelet transform [9] is prone to oscillation, and the edge and texture details cannot be well preserved. Although the above algorithm can achieve good results when removing the special noise, it is not suitable for mixed noise to meet some situations. The detection and nonlocal mean denoising method based on robust singularity rate proposed in [10] had a large time complexity in the block search stage. In [11], combined adaptive median filtering (AMF) with improved sparse representation to remove mixed noise is presented. It uses the AMF to remove saltandpepper noise and then removes Gaussian noise by means of improved Ksingular value decomposition (KSVD) dictionary learning method and backtrackingbased adaptive orthogonal matching pursuit (BAOMP) sparse coding method. In [12], the new mixed noise removal scheme is based on a weighted averaging of pixels contained in a filtering block. The main advantage of the proposed scheme is that the new similarity measure is not influenced by the outliers injected into the image by the impulsive noise, and the averaging process ensures the effectiveness of the new filter in the reduction of Gaussian noise. The anisotropic dynamic diffusion model proposed in [13] could not only effectively remove noise but also retain more texture details. Generally speaking, the PM model is a nonlinear PDE that can identify textures and edge regions. After the denoising process, the texture and edges can be effectively retained to achieve the effect of edge preservation.
Image denoising is an extensively studied problem in the image processing community and continues to attract researchers who aim to perform better restoration in the presence of noise. Both of the Gaussian and saltandpepper noises are the common types of noises. Locust slice images are easy to be corrupted by saltandpepper noise due to bit errors in the process of image acquisition and transmission. The locust slice images may also be corrupted by Gaussian noise because the transmittance of different particles to light is Gaussian distribution. In this paper, the problem of image denoising with mixed noise is considered, in which the image is contaminated by both Gaussian noise and saltandpepper noise. However, the variation method cannot effectively remove saltandpepper noise, and all other existed filtering algorithms are all most with the limitations. So, a novel denoising method for removing mixed noise from locust slice images is proposed by means of Shannoncosine wavelet and the nonlinear variational model for the image processing. This method includes two parts that are the sparse representation of the slice images and the novel numerical algorithm for solving the variation model on image denoising based on the sparse representation. In order to accurately recognize Gaussian noise and image contour of the image, a method is designed to calculate the adaptive selection parameter using image contrast. With this new image removing method, the mixed noise of locust slice images can be processed precisely and efficiency.
2. Sparse Representation of the Slice Images
In the recent decades, many wavelets which have compact support, smoothness, and other properties have been constructed. Cattani studied the properties of the Shannon wavelet function, which possesses many advantages such as orthogonality, continuity, and differentiability [14, 15]. A perceived disadvantage of the Shannon scaling function is that it tends to zero quite slowly as [16]. Shannon–Gabor wavelet constructed by Aldroubi et al. [17, 18] is a Shannon wavelet with Gauss window, which is being called quasiwavelet with limited use. Recently, the Shannoncosine wavelet [19] was proposed with the characteristic of compact support, which inherits all the excellent characteristics of the Shannon wavelet and the normalized characteristics of the Shannon function. Furthermore, Shannoncosine wavelet has been used to solve the fractional calculus problems in the recent years [20]. Therefore, Shannoncosine wavelet is employed to construct the multiscale interpolator operator in our research, which can be expressed as follows:wherewhere is the matrix pulse function on and is a constant associated with the support interval. denotes the Heaviside function. The pulse function not only gives the Shannoncosine scale function tight support but also has no longer the continuity of at the endpoints. At the same time, the normalized property at the endpoint can be satisfied by selecting the parameters and [21]. So, the Shannoncosine wavelet is more suitable to locust slice images with excellent properties such as small support interval, compact support, and normalization.
2.1. Construction Schemes of Multiscale Interpolation Operators
Let be the scale basis function of the interpolation function, and the interpolation operator is defined as follows:
For function , large enough is found to make to approximate . Let the coefficients of wavelet function and scale function be and , respectively, thenwhere and is the characteristic point of the wavelet on the layer. And the wavelet coefficient is defined as follows:where and represents the wavelet interpolation operator of the layer. The calculation formula of the interpolated wavelet coefficients can be obtained by combining the above two formulas:
The wavelet transform coefficient has an intuitive geometric meaning, which is the error between the interpolation result and the exact result at the same collocation point. It reflects the local regularity of the approximated function. By multiplying the matrix with the signal, the multiscale wavelet transform coefficients of the signal can be obtained at one time, which can be used to identify the saltandpepper noises in the slice images. Then, after multiscale interpolation wavelet reconstruction, saltandpepper noise is effectively removed. According to the interpolation wavelet transform theory [22], the wavelet coefficients can be written aswhere denotes the multiscale wavelet function and the corresponding and derivatives with respect to and , respectively. is the position of the pixel after the discrete, and denotes the multilevel interpolation operator, where and are constants, which denote the wavelet scale number and the maximum of the scale number, respectively. The similar method can be used to obtain wavelet transform coefficients and at the points .
After wavelet transform, the useful signals in the image are mainly retained in the small wavelet coefficients. However, the signals energy of saltandpepper noise mainly exists in the large wavelet coefficients. So, signals between the different size thresholds are preserved by choosing the appropriate threshold. Then, the saltandpepper noise can be effectively removed. The description of image detail feature pixels at different scales can be achieved by setting wavelet coefficient thresholds at different scales, which iswhere is the wavelet threshold at and is the scaling parameter, which is usually taken as . During the iterative process in solving the nonlinear PDEs by means of the wavelet numerical method, most of the discrete points in which the wavelet transform coefficients are smaller than the threshold epsilon can be neglected, and this is helpful for improving the efficiency of the algorithm.
In order to obtain the multilevel interpolation operator, it is necessary to express the wavelet coefficients as a weighted sum of in all of the coefficient points in the level. So, the constraint operator is introduced. It can be seen from the meaning represented by interpolating wavelet operator. All the wavelet characteristic points on the layer can get the approximation function on the layer by interpolating the wavelet operator. The expression is as follows:
The interpolation wavelet operator is as follows:where .
2.2. Texture Features and Collocation Points
The locust microslice image is the one that possesses some properties such as texture, smooth area, and boundary, and its structure is relatively complicated, as shown in Figure 1.
The block inside the number 1 rectangular box represents the image texture. The number 2 is a smooth area of the locust image, and the number 3 is the image edge. Noise is mostly distributed in highfrequency information. But edges and texture details are also highfrequency information. So, the edge structure and texture information will be removed at the same time of denoising. The collocation points are image characteristic description operator based on multiscale space. Shannoncosine wavelet is applied to represent the slice images sparsely by adjusting the parameters adaptively based on the texture of the locust slice images. For contours with large gradient changes, more characteristic points are retained. However, in the detail texture of the image with small gradient change, the retained characteristic points are decreased; thereby, the denoising accuracy of the locust slice images is improved. Multiscale refinement can highlight the features of detail information in the image [23]. The collocation points for the locust slice image is obtained by the sparse interpolation operator, as shown in Figure 2.
Multiscale interpolation wavelet can be widely used in the field of image processing. The purpose of constructing the multilevel wavelet collocation method is to decrease the amount of the collocation points and then improve the efficiency of the algorithm. It can be seen from Figure 2 that the collocation points are adaptively collocationed by sparse interpolation wavelet, sparsely extracting points in smooth areas and densely extracting points in detailed texture areas. According to statistics, the total number of pixels in locust slice is 90000, and the number of adaptive collocation points is 12319. Obviously, the computation complexity is decreased greatly compared with processing all pixels. With the image represented sparsely, the collocation points are far less during image processing to save space and improve efficiency.
2.3. Sparse Representation and Noise Restoration
Multiscale interpolation operators are used to process the internal structure of the images. As shown in Figure 3(b), saltandpepper noises with a density of 0.03 dB are added. The denoising results of the sparse interpolation are shown in Figure 3(c). It can be seen from Figure 3 that the sparse representation wavelet can identify and remove saltandpepper noise through multiscale wavelet transform coefficients. The collocation points are adaptively captured by sparse interpolation wavelets, as shown in Figure 3(d). The total number of pixels in this image is 90000, and the number of adaptive collocation points is 13764. The efficiency is about 85% higher than the original efficiency. So, the multiscale interpolation wavelet can not only sparsely represent the image by the texture of the image but also identify and remove the saltandpepper noise of highfrequency information.
(a)
(b)
(c)
(d)
The noise points usually distribute in highfrequency information, while effective information distributes in lowfrequency information. Therefore, the noise point can be separated from the information point by setting threshold. The denoising effect of the image is affected by setting the size and number of thresholds. The wavelet transform coefficients are smaller than the threshold epsilon that can be neglected. In this paper, two thresholds are set to help improve the efficiency of the algorithm. With the increasing number of iterations, part of the texture at the edge of the images are treated as noise points, which decreases denoising performance. The effect of the proposed algorithm to remove saltandpepper noise is shown in Figure 4. The result shows that the number of iterations should not be too large when removing saltandpepper noise.
(a)
(b)
(c)
(d)
(e)
(f)
3. Sparse Wavelet Precise Integration Method (SWPIM)
3.1. Denoising Model and Problem Analysis of Nonlinear PDE
Perona–Malik adopted different diffusion speeds in the interior and edge of the feature area to achieve the effect of protecting the image boundary [24]. Perona–Malik edgepreserving smoothing model is defined as follows:where denotes the pixel point position and denotes the image to process. is the time parameter. denotes the image after processing, and denotes the initial value. is the diffusion operator. denotes the gradient of the image, and represents the diffusion coefficient. By substituting equations (10) and (11) into equation (12), the nonlinear PDE can be discretized into nonlinear ordinary differential equations. The ordinary differential equations are solved by precise integration. The discrete form of PM denoising equation is defined as follows:
Let the right side of formula (13) be noted as . Then, formula (13) can be written as follows:
Substituting the sparse representation of the slice images into the Perona–Malik model which is the famous edgepreserving smoothing model for removing the Gaussian noises of the biomedical images, a system of nonlinear differential equations is obtained. A semianalytical method on solving the system of the nonlinear differential equations is constructed based on the sparse representation of the slice images, which is named SWPIM. With the application of PDE in image processing, the wavelet precise integration method has also been widely used in image processing [25–27]. The essence of precise integration method is the solution of matrix, which is simple and has high calculation accuracy. As the calculation process is almost unlimited by the time step, the calculation efficiency is greatly improved. Perona–Malik introduced the following two diffusion coefficients within the model studied [28, 29]:where is used to control the degree of diffusion. With the increasing number of iterations, the gradient value of the image is smaller and smaller relative to . In order to ensure the relative proportion between and to prevent the edges of the image from being smoothed, it is unreasonable for k to be set as a constant. In the PM model, the image contour and noise are identified by the parameter . In order to accurately recognize Gaussian noise and image contour of the image, a method of adaptively selecting the parameter is proposed in this paper. Linear relationship between image contrast and value is expressed with the following formula:where denotes the contrast of the image and the and are constants.
The position of the PM model with a larger image gradient value is regarded as the image contour. When the diffusion coefficient decreases, the smoothing speed of the diffusion equation at the image contour will become slower. Since the difference value at the noise is smaller than the difference value at the contour, the attenuation speed will increase. That is, for the same number of iterations, the edge will be protected by PDE, and the saltandpepper noise will not be filtered out. There are many ways to solve partial differential equations, and the most typical method is the difference method [30]. This method uses the flat function to describe image, approximately the surface function. But, it is easy to cause artifacts phenomenon, affecting the accuracy of image denoising. Wavelet function has both smooth and compactly supported characteristics. Besides, performance of multiscale analysis can be used to construct the multiscale adaptive interpolation operator for solving nonlinear partial differential equations [31].
3.2. Determination of Parameter
As we all know, human vision is very sensitive to the contour and texture of the object. In the process of denoising, the edge texture of the images is usually treated as noise points, which makes the images blurred. So, it is necessary to effectively identify the boundary information by adjusting so that the boundary parts are not treated as noise points. The determination of value is closely related to the image contrast. The enhancement of image contrast will increase the gradient value of the contour pixels of the object in the image. for identifying contours and noise thresholds should also be adjusted accordingly. Based on this, a method for adaptively calculating the parameter is designed in this paper. The contrast of the locust slice images is calculated to determine the boundary texture of the image. At the same time, the Gaussian noise can be removed well by adaptively adjusting the value. For the locust slice images, the relationship formula between value and image contrast constructed in this paper is defined as follows:where denotes image contrast. Image contrast means difference of brightness in the image or the magnitude of gray contrast of the image [32]. The larger the range of difference, the greater the contrast. On the contrary, the smaller the difference range, the smaller the contrast. Early image contrast algorithms only considered image enhancement, not the measurement of image contrast. In this paper, this factor can be taken into account by using a contour detection operator in a window of appropriate size. Specifically, with the current pixel as the center, the number of pixels in the window should be odd rows and odd columns. The window should be large enough to ensure that noise is filtered out. At the same time, the window should be small enough to ensure that the visible details in the image can be preserved. Then, the boundary value of each pixel is calculated in each window . The boundary value is used to determine whether it is a weight of the boundary value. If the boundary value is 0, it cannot be determined as a boundary. If the boundary value is equal to 1 (after normalization), the pixel must be a point on the boundary. Laplacian operator and Sobel operator can be used to calculate the boundary value . These operators are not sensitive to noise compared with other operators. In this paper, the Laplacian operator is used to calculate the boundary weight :where denotes the average gray value of eight neighboring pixels of pixel . And is set to the gray value at pixel . In each window , the average boundary gray value (that is, the average gray value of the point at the target boundary) is calculated. The gray value of the target boundary can be obtained by calculating the weighted average of the gray value of the pixel. The definition of the average gray value of the boundary is defined as follows:
In this paper, a large window is used to ensure that these operators are not sensitive to noise. The image contrast at the center pixel of the window is
Figure 5 shows locust slice images with different contrasts. Figure 6 shows that Gaussian noise with a standard deviation of 0.002 is added to the locust slice image with different contrasts. Figure 7 shows that the locust slice images with Gaussian noise are denoised by the proposed algorithm. Table 1 shows the , PSNR, and SSIM values corresponding to the optimal denoising effect of locust slice images with different contrasts. From Figures 5–7 and Table 1, it can be seen that the corresponding values in the denoising process are different for locust slice images with different contrasts. And the effect of using the proposed algorithm for processing is also different. The image contrast is linearly correlated with the adaptive value within a certain range. The experimental results show that the boundary information can be effectively identified by the adaptive value so that the boundary part is not treated as a noise point. It provides convenience for subsequent denoising work and makes the image denoising effect more significant.
(a)
(b)
(c)
(d)
(e)
(f)
(g)
(h)
(i)
(j)
(k)
(l)
(m)
(n)
(o)
(p)
(a)
(b)
(c)
(d)
(e)
(f)
(g)
(h)
(i)
(j)
(k)
(l)
(m)
(n)
(o)
(p)
(a)
(b)
(c)
(d)
(e)
(f)
(g)
(h)
(i)
(j)
(k)
(l)
(m)
(n)
(o)
(p)

3.3. Sparse Interpolation Wavelet Precise Integration Method
A novel denoising method for removing mixed noise from locust slice images is proposed by means of Shannoncosine wavelet and the nonlinear partial differential equation variational model for the image processing. This method includes two parts that are the sparse representation of the slice images and the novel numerical algorithm for solving the variation model on image denoising based on the sparse representation. In the first part, a parametric Shannoncosine wavelet function is introduced to construct the multiscale wavelet transform matrix, which is applied to represent the slice images sparsely by adjusting the parameters adaptively based on the texture of the locust slice images. By multiplying the matrix with the signal, the multiscale wavelet transform coefficients of the signal can be obtained at one time, which can be used to identify the saltandpepper noises in the slice images. This ensures that the saltandpepper noises points are kept away from the sparse representation of the slice images. In the second part, a semianalytical method on solving the system of the nonlinear differential equations is constructed based on the sparse representation of the slice images, which is named SWPIM. Substituting the sparse representation of the slice images into the Perona–Malik model which is the famous edgepreserving smoothing model for removing the Gaussian noises of the biomedical images, a system of nonlinear differential equations is obtained, whose scale is far smaller than one obtained by the difference method. Figure 8 is the flow chart of the proposed algorithm.
4. Application of Coupling Algorithm in Biological Slice Images
According to the effect of noise on the image, the noise is divided into two categories: additive noise and multiplicative noise. This paper focuses on additive noise, and the mixed noise is considered to denoise. The mixed noise includes both Gaussian noise and saltandpepper noise. In order to verify the denoising performance of the proposed algorithm on locust slice images, the problem of image denoising with single noise (Gaussian noise or saltandpepper noise) is also considered. Then, these noise test images are restored by the proposed algorithm. The mathematical model of the noise image is as follows:where is the original image and is the image with noise . is the mixed noise. is the set of pixels. The density of saltandpepper noise pixels in the image is denoted by , while the standard deviation of the Gaussian noise is denoted by .
This article uses Matlab for simulation experiments. In order to demonstrate the effectiveness of the proposed algorithm, it is compared with the representative denoising methods. Locust slice images with different image characteristics are used as test images. In order to evaluate the performance, the denoising algorithm of this paper and other denoising algorithms are compared in terms of objective indicators and visual effects.
4.1. Evaluation and Analysis of Objective Indicators
4.1.1. Image Quality Evaluation Standard
The performance of image restoration results is quantitatively evaluated by two measures, PSNR (peak signaltonoise ratio) and SSIM (structural similarity index). The higher the value of PSNR is, the better the quality of the restored image is. PSNR can be expressed aswhere and denote the original pixel and the restored pixel, respectively. and are the sizes of an image for the width and the height, respectively. The denotes the mean square error between the original noisefree image and the restored image. The SSIM denotes the structural similarity between the original noisefree image and the restored image which is given bywhere and are the means of the input image and the restored image, respectively. is the covariance of the input image and the restored image. and are the variances of the input image and the restored image, respectively. and are the constants, which can be calculated according to the paper. The higher value of SSIM means that the structures in the restored image can be preserved more completely.
4.1.2. Objective Evaluation of Image Denoising Performance Contaminated by SaltandPepper Noise
In the first set of experiments, the problem of image denoising is considered with saltandpepper noise. The image in Figure 1 is used as the test image. It is added to the saltandpepper noise with different degrees: the density of saltandpepper noise varies with 0.01, 0.02, 0.03, 0.04, and 0.05. The proposed algorithm is compared with representative denoising methods including the improved median filtering (IMED), mean filtering (MEA), Wiener filtering (WI), wavelet transform (WA), improved nonlocal mean filtering (INLM), bilateral filtering (BI), median and Wiener combination (MWI), median and wavelet combination (MWA), and median and Gauss combination (MG). The quantitative indicators PSNR and SSIM data are shown in Tables 2 and 3, respectively. Table 2 presents the performance comparisons for the different denoising algorithms in terms of PSNR value. Table 3 presents the performance comparisons for the different denoising algorithms in terms of SSIM value.


As shown in Figure 9, the changes in the PSNR and SSIM evaluation indicators in Tables 2 and 3 are vividly and intuitively reflected through line charts. Figure 9 shows how PSNR and SSIM varies for the denoised images as the saltandpepper noise varies from 0 to 0.05. It can be seen that, in the case of low noise corruption ( less than 0.03), all the denoising algorithms can perform well.
(a)
(b)
It can be seen from Tables 2 and 3 that the proposed algorithm shows a good ability to deal with saltandpepper noise. The PSNR of the image restoration results do not decrease significantly. With the increase in saltandpepper noise from 0.01 dB to 0.03 dB, the PSNR value of the proposed denoising algorithm has always been the largest, and the maximum value is 29.1610 dB. The proposed algorithm reveals the superior denoising performance to the median filtering. The SSIM value of the proposed algorithm is also the largest. The similarity of the image processed using the proposed algorithm to the original image reaches 0.9612. It means that the proposed algorithm has a stronger capability than other compared denoising algorithm for preserving image structures. As to the cases of heavy noise corruptions (noise density higher than 0.03), although the proposed algorithm’s recovery results are not as good as the median filtering, the PSNR values are larger than other algorithms. Hence, Tables 2 and 3 indicate that the proposed denoising algorithm is able to more effectively remove saltandpepper noise than other compared denoising algorithm.
4.1.3. Objective Evaluation of Image Denoising Performance Contaminated by Gaussian Noise
In the second set of experiments, the problem of image denoising with Gaussian noise is considered. The image in Figure 10(a) is used as the test image. It is added to the Gaussian noise with different degrees: the standard deviation of Gaussian noise varies with 0.001, 0.002, 0.003, 0.004, and 0.005. The proposed algorithm is compared with representative denoising methods including the improved median filtering (IMED), mean filtering (MEA), Wiener filtering (WI), wavelet transform (WA), bilateral filtering (BI), improved nonlocal mean filtering (INLM), and median and Wiener combination (MWI). The quantitative indicators PSNR and SSIM data are shown in Tables 4 and 5, respectively. Table 4 presents the performance comparisons for the different denoising algorithms in terms of PSNR value. Table 5 presents the performance comparisons for the different denoising algorithms in terms of SSIM value.
(a)
(b)
(c)
(d)
(e)
(f)
(g)
(h)
(i)
(j)


In order to intuitively reflect the changes in the PSNR and SSIM evaluation indicators in Tables 4 and 5, data comparison line charts are given as shown in Figure 11. It shows how PSNR and SSIM varies for the denoised images as the Gaussian noise varies from 0.001 to 0.005.
(a)
(b)
The higher the value of PSNR is, the better the quality of the restored image is. It can be seen from Table 4 and Figure 11(a) that the proposed algorithm usually obtains the largest PSNR value, and the maximum value is 31.9292 dB. In the same noise corruption, the proposed algorithm can perform better than other algorithms. In the case of low noise corruption ( less than 0.005), all of the difference values between the proposed algorithm and the WI are greater than 1 dB and those between the proposed algorithm and the IMED are greater than 2 dB. At the same time, the PSNR value of INLM restoration image is about 3 dB higher than that of the proposed algorithm. From the quantitative comparison as shown in Table 5 and Figure 11(b), we can see that the SSIM of our image recovery results is always the highest in Gaussian lownoise environment, and the maximum value is 0.8966. As to the cases of heavy noise corruptions, the SSIM value of INLM restored image is about 0.2 higher than that of the proposed algorithm. So, in Gaussian noise environments, the proposed algorithm reveals the superior denoising performance than the other algorithms. In other words, the structures in the restoration results of the proposed algorithm can be preserved more completely than other algorithms. The proposed algorithm can achieve the larger values of PSNR and SSIM than other algorithms, which can indicate that the proposed algorithm performs better than the other algorithms in image denoising and detail preservation.
4.1.4. Objective Evaluation of Image Denoising Performance Contaminated by Mixed Noise
In the third set of experiments, image denoising with mixed noise is considered, which includes Gaussian noise and saltandpepper noise. Further quantitative results for the all algorithms are obtained using input images generated by adding saltandpepper noise with different levels ( = 0.01, 0.02, 0.03) mixed with Gaussian noise with different levels ( = 0.001, 0.002, 0.003, 0.004, 0.005). The locust slice image in Figure 12(a) is used as the test image. The proposed algorithm is compared with representative denoising methods including the improved median filtering (IMED), mean filtering (MEA), Wiener filtering (WI), wavelet transform (WA), bilateral filtering (BI), improved nonlocal mean filtering (INLM), median and Wiener combination (MWI), and median and wavelet combination (MWA). Tables 6 and 7 present the quantitative comparisons for the various algorithms in terms of the PSNR value and the SSIM value, respectively.
(a)
(b)
(c)
(d)
(e)
(f)
(g)
(h)
(i)
(j)
(k)


In order to visually illustrate performance, Figure 13 shows the trends of PSNR evaluation indicators for different mixed noise intensities. Figure 13(a) shows how PSNR of denoised images varies as the standard deviation of the Gaussian noise varies from 0.001 to 0.005 with a fixed saltandpepper noise density = 0.01. Figure 13(b) shows how PSNR varies for the denoised images as the varies from 0 to 0.03 with fixed Gaussian noise σ = 0.001.
(a)
(b)
From Figure 13(b) and Table 7, we can see that the ability of the proposed algorithm is better to process saltandpepper noise. In low saltandpepper noise environments, as the standard deviation of Gaussian noise increases from 0.001 to 0.005, the proposed algorithm is superior to other algorithms. And compared to other algorithms, the proposed algorithm can obtain the higher value of PSNR in the cases of various noise densities. The difference between PSNR value of MWI restored image and SSIM value of MEA is about 1 dB. The IMED and proposed algorithm denoising performance are almost similar. It can be seen from Table 7 that the SSIM values of the restoration results of the MEA, WI, and BI are significantly less than the MWI, IMED, INLM, and the proposed algorithm under the case of various noise densities. Compared with other algorithms, this proposed algorithm can achieve larger PSNR and SSIM values in a low noise densities environment. The proposed algorithm reveals the superior denoising performance than the IMED, MEA, INLM, and MWI. It means that the proposed algorithm has a stronger capability than other compared denoising algorithms for preserving image structures. Hence, Tables 6 and 7 and Figure 13 indicate that the proposed denoising algorithm is able to more effectively remove mixed noise than other compared denoising algorithms. The ability to remove the mixed noise of the locust slice image from the objective effect is very good. And the structures in the restoration results of the proposed algorithm can be preserved more completely than other algorithms. The proposed algorithm provides important features for identifying targets in locust bioslice images.
4.2. Visual Evaluation of Denoising Performance
Due to the large number of test images, partial noise recovery results are used to evaluate visual effects in each set of experiments. Figure 14 shows the recovery results of various denoising algorithms in the first set of experiments. And the test images are corrupted by saltandpepper noise with 0.01 noise density. The original images in these experiments are shown in Figure 1. It can be found that the WI, BI, and INLM fail to restore the noise corrupted image. MEA, WA, and MG suffer from a great quantity of blurred effects for edges and details. MWI suffers from an excessive smoothing effect for restored image. IMED and proposed algorithm can efficiently remove the noise and be free from the blurred effect in the restored images. But the proposed algorithm has no blurred effect in the restored image. So, the proposed algorithm achieves the better visual effect than other denoising algorithms. The proposed algorithm can retain more texture details than the other compared algorithms, especially for the marked area.
(a)
(b)
(c)
(d)
(e)
(f)
(g)
(h)
(i)
(j)
(k)
(l)
Figure 10 shows the recovery results of various denoising algorithms in the second set of experiments. And the test image is corrupted by Gaussian noise with 0.001 noise density. The original image in this experiment is shown in Figure 10(a). Although other algorithms can remove the Gaussian noise from the corrupted images, they still suffer from the blurred effect for the edges and details in the restored images. MEA has a certain effect of removing Gaussian noise, but it makes the image blurred. The INLM and BI can not only filter the Gaussian noise well but also protect the edge information of the image to the greatest extent and achieve better filtering effect. However, the proposed algorithm can retain more texture details than the INLM and BI, especially for the marked area. And, the proposed algorithm has a higher image clarity after restoration.
Figure 12 shows the recovery results of various denoising algorithms in the third set of experiments. And the test image is corrupted by saltandpepper noise with a noise intensity of 0.01 and Gaussian noise with standard deviation of 0.002. Original image in this experiment is shown in Figure 12(a). It can be seen from the figure that MEA, WI, WA, and BI are well restored by single noise corruption. These algorithms suffer from blurring effect in the restored image. And there are still some noise points after filtering. The BI can restore image edges and details from the noise corrupted image, and plenty of residual noise still exists in the restored images. The INLM has a certain effect of removing noise and can retain the edge details of the image to the greatest extent. However, the filtered image appears blurred and has some noise points. Although MWI can remove the mixed noise from the corrupted images, they still suffer from the blurred effect for the edges and details in the restored images. The proposed algorithm can efficiently remove the noise and be free from the blurred effect in the restored images. Besides, the proposed algorithm achieves the better visual effect than other algorithms because they can preserve more edges and contain less noise in the restoration images. Therefore, the above restoration results clearly indicate that the proposed algorithm obtains the best visual effect in terms of noise suppression and detail preservation, especially for the marked area.
5. Conclusions
Shannoncosine wavelet precise integration method is a new image denoising method based on Shannoncosine wavelet and the nonlinear PDE variational model, which was used to construct adaptive wavelet interpolation operator due to multiscale characteristics of wavelet transform. The adaptive interpolation operator in the Shannoncosine wavelet precision integration method can reduce the amount of the collocation points and improve the calculation efficiency. The efficiency is about 86% higher than the original efficiency. Multiscale wavelet transform coefficients are used to identify and remove saltandpepper noise in slice images. The proposed algorithm shows superior denoising performance to representative denoising algorithms.
The method makes full use of the multiscale characteristics and the high precise performance of precise integration method. Based on the edge retention smoothing model, SWPIM is used to remove Gaussian noise in biomedical images. The precise integration method is used to calculate the nonlinear ordinary differential equations, which improves the accuracy of the calculation. Compared to the other methods of image denoising, object boundary of locust slice images obtained by SWPIM denoising method is clear. With the increase in noise, the proposed algorithm can achieve the larger values of PSNR and SSIM than other algorithms. Experimental results show that the proposed calculation method of adaptive selection parameter k can accurately identify Gaussian noise and image contours in the image. SWPIM method is very suitable for measurement of image denoising, such as locust slice images.
Data Availability
The locust slice data used to support the findings of this study have been deposited in the Locustsslicedata repository (https://github.com/feifei123feifei/Locustsslicedata.git). All data generated or analyzed during this study are included in these files.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
The authors would like to thank their colleagues for all their help and support during this work. This work was supported by the Beijing Natural Science Foundation (Grant no. 4172034) and National Natural Science Foundation of China (Grant no. 61871380).
References
 L. Li, N. Zhang, S. Mei, and X. Li, “Image denoising of locust sections based on adaptive wavelet and partial differential equation method,” Transactions of the Chinese Society of Agricultural Engineering, vol. 31, no. 20, pp. 172–177, 2015. View at: Google Scholar
 L. Li, S. Guo, S. Mei, and N. Zhang, “Image restoration of locust slices based on nearest unit matching,” Transactions of the Chinese Society of Agricultural Machinery, vol. 46, no. 8, pp. 15–19, 2015. View at: Google Scholar
 S. Mei, “Denoising for locust slice image with texture preserving based on coupling technology of variational method and shearlet transform,” Transactions of the Chinese Society of Agricultural Engineering, vol. 32, no. 17, pp. 152–159, 2016. View at: Google Scholar
 W. Tian, T. Ma, Y. Zheng et al., “Weighted curvaturepreserving PDE image filtering method,” Computers & Mathematics with Applications, vol. 70, no. 6, pp. 1336–1344, 2015. View at: Publisher Site  Google Scholar
 G. Liu and W. Guo, “Application of improved arithmetic of median filtering denoising,” Computer Engineering and Applications, vol. 46, no. 10, pp. 187–189, 2010. View at: Google Scholar
 S. Yi and J. He, “Image denoising method based on bemd and adaptive wiener filter,” Computer Engineering and Applications, vol. 49, no. 10, pp. 156–158, 2013. View at: Google Scholar
 X. Liu, Y. Tian, H. Hui, and Y. Zhong, “Improved nonlocal means algorithm for image denoising,” Computer Engineering, vol. 38, no. 4, pp. 199–201, 2012. View at: Google Scholar
 X. Qian and L. Guo, “Adaptive Gaussian filter based on object scale,” Computer Engineering and Applications, vol. 46, no. 12, pp. 14–20, 2010. View at: Google Scholar
 R. Xing and S. Mei, “Engineering drawings denoising based on orthogonal wavelet packets transormation,” Journal of China Agricultural University, vol. 8, no. 6, pp. 62–65, 2003. View at: Google Scholar
 B. Xiong and Z. Yin, “A universal denoising framework with a new impulse detector and nonlocal means,” IEEE Transactions on Image Processing, vol. 21, no. 4, pp. 1663–1675, 2011. View at: Google Scholar
 J. Zhang, L. Pei, X. Li, and H. Wu, “Hybrid denoising algorithm combining median filtering with sparse representation,” Computer Engineering, vol. 42, no. 9, pp. 240–245, 2016. View at: Google Scholar
 B. Smolka and D. Kusnik, “Robust local similarity filter for the reduction of mixed Gaussian and impulsive noise in color digital images,” Signal, Image and Video Processing, vol. 9, no. 1, pp. 49–56, 2015. View at: Publisher Site  Google Scholar
 H. Chen, Y. Zhang, M. K. Kalra et al., “Lowdose CT with a residual encoderdecoder convolutional neural network (REDCNN),” IEEE Transactions on Medical Imaging, vol. 36, 2017. View at: Google Scholar
 L. Li, N. Zhang, and S. Mei, “Anisotropic dynamic diffusion model for texture preserving denoising of tomato images,” Transactions of the Chinese Society for Agricultural Machinery, vol. 47, no. 11, pp. 18–24, 2016. View at: Google Scholar
 C. Cattani, “Shannon wavelets theory,” Mathematical Problems in Engineering, vol. 2008, Article ID 16480, 24 pages, 2008. View at: Publisher Site  Google Scholar
 C. Cattani, A. Ciancio, and B. Lods, “On a mathematical model of immune competition,” Applied Mathematics Letters, vol. 19, no. 7, pp. 678–683, 2006. View at: Publisher Site  Google Scholar
 A. Aldroubi, M. Unser, and A. Aldroubi, “Sampling procedures in function spaces and asymptotic equivalence with shannon’s sampling theory,” Numerical Functional Analysis and Optimization, vol. 15, no. 1–2, pp. 1–21, 1994. View at: Publisher Site  Google Scholar
 D. K. Hoffman, G. W. Wei, D. S. Zhang, and D. J. Kouri, “Shannon–gabor wavelet distributed approximating functional,” Chemical Physics Letters, vol. 287, no. 1–2, pp. 119–124, 1998. View at: Publisher Site  Google Scholar
 G. W. Wei, “Quasi wavelets and quasi interpolating wavelets,” Chemical Physics Letters, vol. 296, no. 3–4, pp. 215–222, 1998. View at: Publisher Site  Google Scholar
 S.L. Mei, H.L. Lv, and Q. Ma, “Construction of interval wavelet based on restricted variational principle and its application for solving differential equations,” Mathematical Problems in Engineering, vol. 2008, Article ID 629253, 14 pages, 2008. View at: Publisher Site  Google Scholar
 H. Wang and S.L. Mei, “Shannon wavelet precision integration method for pathologic onion image segmentation based on homotopy perturbation technology,” Mathematical Problems in Engineering, vol. 2014, no. 3, Article ID 601841, 10 pages, 2014. View at: Publisher Site  Google Scholar
 S. Mei and W. Gao, “Shannoncosine wavelet spectral method for solving fractional FokkerPlanck equations,” International Journal of Wavelets, Multiresolution and Information Processing, vol. 16, no. 3, Article ID 1850021, 2018. View at: Publisher Site  Google Scholar
 Y. Wei, Y. Zhang, S. Mei, and S. Wei, “Image dehazing method based on dark channel prior and interval interpolation wavelet transform,” Transactions of the Chinese Society of Agricultural Engineering, vol. 33, no. 1, pp. 281–287, 2017. View at: Google Scholar
 S. Guo, S. Mei, and L. Li, “Texture extension method for farmland remote sensing image based on shannoncosine wavelet,” Transactions of the Chinese Society for Agricultural Machinery, vol. 48, no. S1, pp. 142–146, 2017. View at: Google Scholar
 Y. Zhang, Y. Wei, S. Mei, and M. Zhu, “Application of multiscale interval interpolation wavelet in beef image of marbling segmentation,” Transactions of the Chinese Society of Agricultural Engineering, vol. 32, no. 21, pp. 296–304, 2016. View at: Google Scholar
 J. L. Wu, “A wavelet operational method for solving fractional partial differential equations numerically,” Applied Mathematics and Computation, vol. 214, no. 1, pp. 31–40, 2009. View at: Publisher Site  Google Scholar
 S. Mei, “Adaptive wavelet precise integration method on nonlinear BlackScholes equations,” Journal of Quantitative Economics, vol. 29, no. 04, pp. 8–14, 2012. View at: Google Scholar
 S. Mei, X. Li, H. Zhao, Li Li, and S. Guo, “Method of denoising and removing artifacts for farm remote sensing image based on shearlet and total variation,” Transactions of the Chinese Society of Agricultural Engineering, vol. 33, no. 1, pp. 274–280, 2017. View at: Google Scholar
 H. Wang, M. Zhu, Li Li, L. Wang, H. Zhao, and S. Mei, “Regional weed identification method from wheat field based on unmanned aerial vehicle image and Shearlets,” Transactions of the Chinese Society of Agricultural Engineering, vol. 33, no. 1, pp. 99–106, 2017. View at: Google Scholar
 S.L. Mei and D. Zhu, “HPMbased dynamic sparse grid approach for peronamalik equation,” The Scientific World Journal, vol. 2014, Article ID 417486, 15 pages, 2014. View at: Publisher Site  Google Scholar
 S. Mei, “HPMbased dynamic wavelet transform and its application in image denoising,” Mathematical Problems in Engineering, vol. 2013, Article ID 309418, 10 pages, 2013. View at: Publisher Site  Google Scholar
 X. Li, S. Mei, and L. Li, “Image menoising of locust slices based on Meyer window function shearlet transform,” Transactions of the Chinese Society of Agricultural Engineering, vol. 47, no. z1, pp. 449–456, 2016. View at: Google Scholar
Copyright
Copyright © 2020 Haihua Wang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.