Research Article  Open Access
Bo Chen, Jinbin Zou, Weiqiang Zhang, "Speckle Noise Removal by Energy Models with New Regularization Setting", Journal of Function Spaces, vol. 2020, Article ID 3936975, 18 pages, 2020. https://doi.org/10.1155/2020/3936975
Speckle Noise Removal by Energy Models with New Regularization Setting
Abstract
In this paper, we introduce two novel total variation models to deal with speckle noise in ultrasound image in order to retain the fine details more effectively and to improve the speed of energy diffusion during the process. Firstly, two new convex functions are introduced as regularization term in the adaptive total variation model, and then, the diffusion performances of Hypersurface Total Variation (HYPTV) model and Logarithmic Total Variation (LOGTV) model are analyzed mathematically through the physical characteristics of local coordinates. We have shown that the larger positive parameter in the model is set, the greater energy diffusion speed appears to be, but it will cause the image to be too smooth that required adequate attention. Numerical experimental results show that our proposed LOGTV model for speckle noise removal is superior to traditional models, not only in visual effect but also in quantitative measures.
1. Introduction
With the development of digital image technology, a large number of digital images are transmitted and compressed through various channels. However, image corruption usually is unavoidable during transmission and storage, and the resultant noise quite often seriously affects the visual effect of the image. Clearly, highquality images are desirable in many areas, such as in medical imaging and pattern recognition. Thus, image denoising is a critical step in image processing and computer version, which plays an important role in various applied areas, especially in medical imaging, video processing, and remote sensing. On the other hand, it is also an important preprocessing process for other image processing that relies on subsequent processing.
Today, image denoising becomes research of focus, and many image denoising methods have been proposed such as Lee filter [1], Kuan filter [2], locally adaptive statistic filters [3–5], PDEbased and curvaturebased methods [6, 7], wavelet transform based thresholding methods [8], and total variational [9–11]. In addition, the method based on machine learning [12–15] has received wide attention in recent years, such as deep learning [12, 13], linear regression [14], and Bayesian learning [13, 15].
In 1992, Rudin et al. [10] proposed a denoising model based on the total variation: where represents the TV regularization term, is an original image that is without noise, is noisy image, and represents the Gaussian random noise with mean zero and standard deviation . represents the regularization parameter which can be used to balance the fidelity terms and regularized terms in TV model, and represents the norm of the image gradient.
It is well known that medical ultrasonic images may have many noises of speckle, which will bring a significant problem in terms of the quality of ultrasonic images and cover up the lesions of certain important tissues. Further, it brings great difficulties to the actuate diagnosis and identification of certain specific diseases and can create the potential risk of missed diagnosis and misdiagnosis. Thus, it is very desirable to eliminate the speckle noise in ultrasound image and simultaneously retain the important features in practices. As mentioned in article [3], the speckle noise in medical ultrasonic images can be formulated in the following form: where is a noisy image, and represents the Gaussian random noise with zeros mean and standard deviation .
In this paper, we focus on the image denoising form by using variation method, where the noisy image is ultrasound speckle noise. Based on the model (2) and the characteristics of the Gaussian distribution, Krissian et al. in the article [16] derived a date fidelity term:
Within the variational framework, the data fidelity function is derived from the degradation model (2). One of the technical approaches to solve the variational model is the regularization technique, which minimizes the cost function to obtain stable and accurate solutions. In general, the image denoising variation method is to consider: where is a regularization term that represents a prior information about the object to be restored, and is a fidelity term to ensure that the restoration is not far from the original observation . represents the regularization parameter which can balance the fidelity term and the regularization term.
In [11], motivated by the classical ROF model [10], the authors proposed a convex variational model for removing the speckle noise in ultrasound image. The convex variational model involves the TV regularization term and convex fidelity term (see Equation (5)): where represents the TV regularization term, and represents the directional gradient of . The existence and uniqueness of the solution of model (5) is proved in [11]. In this paper, we call the model in (5) the “JIN’s model.” In [17], the authors proposed a wellbalanced speckle noise reduction (WBSN) model that can detect edges.
Although TV regularization is effective for image denoising, it also leads to some staircase effects that is undesirable. In order to solve this problem, many methods based on improved TV regularization are proposed, such as highorder TV regularization [18–20], several hybrid TV regularization [21], the improved infimal convolution [22, 23], nonlocal TV model [24, 25], fractional order TV model [26], and anisotropic TV model [27, 28]. Fractional theory [29, 30], wavelet [31], and statistical information [32–34] are also employed to deal with intensity inhomogeneity or noise. Although these denoising methods can reduce the staircase effects in the restoration of additive noise images, however, there are more staircase effects that appeared in the restoration of speckle noise images. In this paper, we introduce HYPTV model and LOGTV model to reduce speckle noise and staircase effects in ultrasound images in an effective way.
In numerical algorithms, most of the energy function minimization problems can be transformed to an EulerLagrange equation and then be solved by using the finite difference method. However, the choice of adequate regularization terms is critical in terms of solution accuracy. Moreover, when solving the EulerLagrange equation, the energy diffusion form of the noise image in different regions is required to handle differently based on the physical characteristics of local coordinates, since the diffusion velocity of different parameters to the energy function is different in the two models. The new proposed HYPTV model and the LOGTV model, as shown in this paper, not only can preserve the edges of the restored images well when restoring the ultrasonic image with speckle noise but also better reduce the staircase effect generated during the recovery process.
The rest of this paper is as follows: in Section 2, we review some background knowledge. In Section 3, we propose two new models based on variation; meanwhile, we not only analyze diffusion performance of the proposed models but also give the corresponding numerical algorithms. Section 4 shows five different experiments and results. The paper ends with concluding remarks in Section 5.
2. Preliminary
2.1. Some Theoretical Background
The EulerLagrange equation was developed in the 1750s by Euler and Lagrange in connection with their studies of the tautochrone problem.
Theorem 1. A multidimensional generalization comes from considering a function on variables. If is open, bounded Lipschitz domain in , then is extremized only if satisfies the partial differential equation: For model (5), the corresponding EulerLagrange equation is as follow: where and , respectively, represent gradient operators and divergence operators. Using gradient descent method, we can get the model as follows: where is the unit out normal vector of .
Without losing generality, in the following, we consider the grayscale images as matrices.
Definition 2. Let , the gradient operators on the space and the divergence operators on the space are defined as: where are the firstorder forward and backward discrete derivation operators in the direction and direction, respectively, which are defined as: Applying the above gradient operators and divergence operators to model (9), we can obtain the equivalent minimization problem.
Definition 3. Let be a convex subset of , A function is called convex if
The following facts are easily checked:
Theorem 4. If functions and are convex and have the same domain definition, then is also convex.
Proof. Let functions and are convex. According to the Definition 3, for any and , we have Hence is convex.
Theorem 5. If a differentiable function satisfied then is convex.
2.2. The Condition of TV Regularization Term
Although TV regularization is very effective in image restoration, it usually generates some staircase effects. Thus, it is suggested in literature to use general variational methods, i.e., to consider: where represents a potential function, and the case leads to the total variation regularization term. In the literature [35], the author Costanzino chooses that leads to the wellknown harmonic model.
In practice, we prefer good smoothing in some domain where the intensity of variations is relatively weak. This can be achieved by requiring a function to satisfy the following conditions:
Near the edge of the image, the intensity of variations is strong. If we would like to preserve the edge, then the function should satisfy the following conditions:
With the conditions (16) and (17), the function is convex and nondecreasing function, such as:
There functions are convex and nondecreasing on , as shown in Figure 1.
In this paper, we will use two new functions and , which appear to be quite effective for image processing, in particular, for ultrasound image denoising. Obviously, the functions and are convex and nondecreasing.
3. The Proposed Restoration Model
In this section, we propose adaptive total variation model for image restoration. We use the finite difference method to solve the EulerLagrange equation directly, and then find the minimum value of the energy function.
3.1. The Adaptive Total Variation Model
Apply the selected function to model (15), we propose adaptive total variation model, where is a regularization term, and ; is a fidelity term; represents the regularization parameter which can balance fidelity term and regularization term.
Firstly, the energy functional is convex, which guarantees the existence of the minimal solution of the model (19).
Theorem 6. The energy functional is convex. That is to say, for any , and , we have: where .
Proof. Firstly, the function is convex; according to Definition 3, for any , and , we have
Meanwhile, we have (The proof step of the inequality is in the appendix.)
Therefore,
This proof is established.
Secondly, the uniqueness of the minimum solution of the model (19) can also be proved.
Theorem 7. If and are two minimize solutions of model (19), then we have .
Proof. According to Theorem 6, we have If , then the above assumption gives a contradiction that is not a minimize solution.
3.2. Diffusion Performance
In this subsection, we mainly analyze the diffusion performance and diffusion speed of the energy function of HYPTV model and LOGTV model.
3.2.1. Diffusion Performance of HYPTV Model
Firstly, we use the finite difference method to solve the HYPTV model, associated with the potential function .
From Definition 2, we can obtain the corresponding EulerLagrange equation HYPTV model that as follows:
Using gradient descent method, Equation (25) can be transformed to: where is the unit out normal vector of .
In order to analyze the diffusion performance, local image coordinate system is established. As shown in Figure 2, the axis represents the direction parallel to the image gradient at the pixel level, and the axis is the corresponding vertical direction.
According to Figure 2, we can know:
So, Equation (26) can be rewritten as: where
The and are control functions of the diffusion along the direction and direction, respectively. Now, we consider the diffusion of image restoration. Some test images are shown in Figure 3.
(a)
(b)
(c)
(d)
(e)
(f)
(g)
(1) Smooth area. When , . This shows that the diffusion form of the energy Equation (19) is isotropic. In other words, the energy diffusion rate along direction and direction is very close in the process of image restoration in the smooth region. And the rate of energy diffusion is obviously positively correlated with the parameter .
(2) Sharp area. When , we obtain . This shows that the diffusion form of the energy Equation (19) is anisotropic. In other words, the energy diffusion rate in direction in Equation (28) is much larger than that in the direction in the sharp region. But the gradient does not exceed 255, so . One can see that the larger the parameter is set, the smaller the limit becomes. And the rate of energy diffusion is obviously positively correlated with the parameter .
3.2.2. Diffusion Performance of LOGTV Model
Secondly, we use the finite difference method to solve the LOGTV model, which the potential function is .
From the Definition 2, we can obtain the corresponding EulerLagrange equation LOGTV model that as follows:
Using gradient descent method, Equation (30) can be transformed to: where is the unit out normal vector of .
Hence, Equation (31) can be rewritten as: where
and are control functions of the diffusion along the direction and direction, respectively. Now, we consider the diffusion of image restoration.
(1) Smooth area. When , and , we can obtain and . This shows that the diffusion form of the energy Equation (19) is isotropic. In other words, the energy diffusion rate along direction and direction is very close in the process of image restoration in the smooth region. When , and , we can obtain and . This shows that the diffusion form of the energy Equation (19) is anisotropic. However, the gradient of noise image is relatively large, so in the smooth region, whatever the value of , it has little effect on the model.
(2) Sharp area. When , we can obtain . This shows that the diffusion form of the energy Equation (19) is anisotropic. In other words, the energy diffusion rate in direction in Equation (28) is much larger than that in the direction in the sharp region. But the gradient does not exceed 255, so . One can see that the larger the parameter is set, the smaller the limit becomes. And the rate of energy diffusion is clearly positively correlated with the parameter .
3.3. Numerical Implementation
We will describe the corresponding numerical algorithm in this section. Firstly, the HYPTV model can be solved by discretization as follows: where , and represents time step. Furthermore, the iterative formula can approximate as: for ; , and represent the size of the image. Here: where , and is a positive parameter that is close to zero. With boundary conditions:
Now, note Equation (25), the two sides are multiplied by , and then, the integral on the domain can be obtained:
Because the Gaussian noise have mean 0 and variance , we can obtain:
In the process of iteration, we always use the previous solution to calculate the next solution. The optimization algorithm for HYPTV model is given in the following (Algorithm 1).
Secondly, the LOGTV model can be solve by discretization as follows: where . Furthermore, the iterative formula can approximate as:
Similar to the HYPTV model, the optimization algorithm for LOGTV model is given in the following (Algorithm 2).
4. Experimental Results
In this section, we present numerical results to demonstrate the effectiveness of the HYPTV and LOGTV model in image restoration. Firstly, to evaluate the quality of restored images, we use the peak signaltonoise ratio (PSNR) value and the structure similarity (SSIM) index, which are defined as follows: where is the clean image, and is the restored image. is the average of , is the standard deviation of , and and are some constants for stability. Secondly, we calculated the noise deviation reduction (NDR) at each iteration as a convergence condition; where represents the results of the iterations, respectively. Finally, the stopping condition (NDR) for the HYPTV and LOGTV models is as follows:
In the numerical experiment, we will use the noise image as the initial value, that is, . Moreover, the gray values of all original images are in range [0, 255].
4.1. Denoising Effect of Different Parameters of HYPTV Model and LOGTV Model
In this example, we use different parameter values in the algorithm to test the effect of HYPTV and LOGTV models. The test images is shown in Figure 3(a), and the noise levels are and . Figure 4 shows the different PSNR values and iteration numbers when different values are used in the HYPTV model algorithm. Figure 5 shows the different PSNR values and iteration numbers when different values are used in the LOGTV model algorithm. Firstly, we can also see that PSNR is the largest in . Secondly, with the increase of parameters, the speed of image restoration is faster. Based on the above analysis, when , the denoising performance of HPYTV and LOGTV models is close to the best. Therefore, in the following experiment, we choose in HPYTV and LOGTV models.
4.2. Denoising Effect of the HYPTV and LOGTV Model
In this subsection, we use the algorithm to test the effect of the HYPTV and LOGTV model on image denoising. In Figures 6 and 7, we show the original images, the noise images, and the restored images by HYPTV and LOGTV models. Figures 6(a) and 7(a) are original images, Figures 6(b) and 7(b) are noisy images, which the noise levels are and , respectively. Figures 6(c) and 7(c) are corresponding restored images by HYPTV model. Figures 6(d) and 7(d) are corresponding restored images by LOGTV model. Table 1 shows that the PSNR values, SSIM values, and the numbers of iteration for the different test images can be got by using the HYPTV model and the LOGTV model. “Noise image PSNR” is the peak signaltonoise ratio of noisy images and original images. “Denoising image PSNR” is the peak signaltonoise ratio of restored images and original images. “Iter” is the number of iterations of the algorithm. From the results, it is obvious that the HYPTV model and the LOGTV model are fairly effective in reducing the speckle noise and edgepreserving.
(a)
(b)
(c)
(d)
(a)
(b)
(c)
(d)

4.3. Reduction of Staircase
In this subsection, we test the reduction of the staircase effect in HYPTV model and LOGTV model by image in Figure 8. Figure 8(a) shows original images (“house”). Figure 8(b) shows detailed images of Figure 8(a), respectively. Figure 8(c) shows noise images which noise standard deviation . Figure 8(d) shows detailed images of Figure 8(c). The original images contain a lot of details information, such as textures edges and in homogeneous regions.
(a)
(b)
(c)
(d)
Figure 9 displays the restoration results of the noisy “house” image, respectively. Figures 9(a)–9(e) was restored ROF model [10], ATV model [36], JIN’s model [11], HYPTV model, and LOGTV model, respectively. Figures 9(f)–9(j) is corresponding details images.
(a)
(b)
(c)
(d)
(e)
(f)
(g)
(h)
(i)
(j)
According to the results, we can see that JIN’s model, HYPTV model, and LOGTV model had a good effect in removing image noise and preserving image edges. In addition, we can clearly see staircase effect in the detailed images obtained by ROF model, ATV model, and JIN’s model. At the same time, the staircase effect of HYPTV model and LOGTV model restored images is reduced. Although both the HYPTV model and the LOGTV model can reduce the ladder effect, the LOGTV model is better than the HYPTV model.
4.4. Comparison with ROF Model, ATV Model, and JIN’s Model
In this subsection, we compare the effect of the ROF model, ATV model, and JIN’s model with HYPTV and LOGTV models for some images. The test original images are shown in Figure 3(“house,” “peppers,” “boat,” “pirate,” and “bird”), with two sizes of and three sizes of . Figure 10 shows some noise images with different standard deviation.
(a)
(b)
(c)
(d)
(e)
Figures 11–15 display the restoration results for images (“peppers,” “boat,” “house,” “pirate,” and “bird”) through ROF model, ATV model, JIN’s model, HYPTV model, and LOGTV model. The noise versions of “peppers”, “boat” and “house”, and “pirate” and “bird” are obtained by model (2) with standard deviations 2, 3, and 4, respectively. In addition, the detailed images of the restored images are also displayed. Table 2 shows the PSNR and SSIM values for different test images by using the ROF model, ATV model, JIN’s model, HYPTV model, and LOGTV model. From the results of Figures 11–15 and Table 2, the noise standard deviation . We can observe that although the four models can effectively remove the noises while preserving the edges and details, the restored images by HYPTV model and LOGTV model have better visual effect with less staircase effects than by the ROF model, ATV model, and JIN’s model. The noise standard deviation . The visual effect of restored images by the ROF model and ATV model is particularly poor, but JIN’s model, HYPTV model, and LOGTV model can effectively remove the noises. Finally, Table 2 shows that LOGTV model has higher PSNR and SSIM values than other four models. This means that our proposed LOGTV model is available in reducing the speckle noise in some images.
(a)
(b)
(c)
(d)
(e)
(f)
(g)
(h)
(i)
(j)
(a)
(b)
(c)
(d)
(e)
(f)
(g)
(h)
(i)
(j)
(a)
(b)
(c)
(d)
(e)
(f)
(g)
(h)
(i)
(j)
(a)
(b)
(c)
(d)
(e)
(f)
(g)
(h)
(i)
(j)
(a)
(b)
(c)
(d)
(e)
(f)
(g)
(h)
(i)
(j)
 
Best denoising performance are given in bold. 
4.5. Denoising Results of Real Ultrasound Images
In this subsection, we test some real ultrasound images. Figure 16 shows the experimental results of real ultrasound images by applying JIN’s model, HYPTV model, and LOGTV model. Table 3 shows the different iteration for the different test images by using the JIN’s model, HYPTV model, and LOGTV model. We find that LOGTV model is much effective than JIN’s model and HYPTV model in obtaining the satisfactory restored images.
(a)
(b)
(c)
(d)
(e)
(f)
(g)
(h)
(i)
(j)
(k)
(l)
 
Best performance are given in bold. 
5. Concluding Remarks
In this paper, we propose a new speckle noise restoration model based on adaptive TV method. Two new convex functions are introduced as the TV regularization term. By analyzing the diffusion performance of the proposed two models, one can see that the LOGTV model has faster diffusion speed than the HYPTV model. Moreover, we introduced two iterative numerical algorithms to solve the proposed models. The experiment results show the affectivity of our proposed model and the similarity between JIN’s model and HYPTV model. In addition, we compared the effect of the ROF model, ATV model, and JIN’s model with the LOGTV model, and the experiment results show high efficiency of LOGTV model in image restoration.
Appendix
The Proof of the Inequality of Theorem 6
Proof. For any , and , .
Here, we know ,
so, we have . Therefore, the proof holds to be true.
Data Availability
The experimental data are obtained by MATLAB R2017a, 2.93 GHz cup, 4 G RAM, and Windows 7.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Authors’ Contributions
All authors typed, read, and approved the final manuscript.
Acknowledgments
This paper is partially supported by the Natural Science Foundation of Guangdong Province (2018A030313364), the Special Innovation Projects of Universities in Guangdong Province (2018KTSCX197), the Science and Technology Planning Project of Shenzhen City (JCYJ20180305125609379), the Natural Science Foundation of Shenzhen (JCYJ20170818091621856), and the China Scholarship Council Project (201508440370).
References
 H. H. Arsenault, “Speckle suppression and analysis for synthetic aperture radar images,” Optical Engineering, vol. 25, no. 5, pp. 636–643, 1986. View at: Google Scholar
 D. Kuan, A. Sawchuk, T. Strand, and P. Chavel, “Adaptive restoration of images with speckle,” IEEE Transactions on Acoustics, Speech, and Signal Processing, vol. 35, no. 3, pp. 373–383, 1987. View at: Publisher Site  Google Scholar
 T. Loupas, W. N. McDicken, and P. L. Allan, “An adaptive weighted median filter for speckle suppression in medical ultrasonic images,” IEEE Transactions on Circuits and Systems, vol. 36, no. 1, pp. 129–135, 1989. View at: Publisher Site  Google Scholar
 Y. J. Yu and S. T. Acton, “Speckle reducing anisotropic diffusion,” IEEE Transactions on Image Processing, vol. 11, no. 11, pp. 1260–1270, 2002. View at: Publisher Site  Google Scholar
 K. Krissian, C. F. Westin, R. Kikinis, and K. G. Vosburgh, “Oriented speckle reducing anisotropic diffusion,” IEEE Transactions on Image Processing, vol. 16, no. 5, pp. 1412–1424, 2007. View at: Publisher Site  Google Scholar
 J. Shen and T. F. Chan, “Mathematical models for local nontexture inpaintings,” SIAM Journal on Applied Mathematics, vol. 62, no. 3, pp. 1019–1043, 2002. View at: Publisher Site  Google Scholar
 J. Shen, S. H. Kang, and T. F. Chan, “Euler's elastica and curvaturebased inpainting,” SIAM Journal on Applied Mathematics, vol. 63, no. 2, pp. 564–592, 2003. View at: Publisher Site  Google Scholar
 J. Jin, Y. Liu, Q. Wang, and S. Yi, “Ultrasonic speckle reduction based on soft thresholding in quaternion wavelet domain,” in 2012 IEEE International Instrumentation and Measurement Technology Conference Proceedings, pp. 255–262, Graz, Austria, May 2012. View at: Publisher Site  Google Scholar
 M. Kang, M. Kang, and M. Jung, “Total generalized variation based denoising models for ultrasound images,” Journal of Scientific Computing, vol. 72, no. 1, pp. 172–197, 2017. View at: Publisher Site  Google Scholar
 L. I. Rudin, S. Osher, and E. Fatemi, “Nonlinear total variation based noise removal algorithms,” Physica D: Nonlinear Phenomena, vol. 60, no. 14, pp. 259–268, 1992. View at: Publisher Site  Google Scholar
 Z. Jin and X. Yang, “A variational model to remove the multiplicative noise in ultrasound images,” Journal of Mathematical Imaging and Vision, vol. 39, no. 1, pp. 62–74, 2011. View at: Publisher Site  Google Scholar
 K. G. Lore, A. Akintayo, and S. Sarkar, “LLNet: a deep autoencoder approach to natural lowlight image enhancement,” Pattern Recognition, vol. 61, pp. 650–662, 2017. View at: Publisher Site  Google Scholar
 N. Wang, D. Tao, X. Gao, X. Li, and J. Li, “A comprehensive survey to face hallucination,” International Journal of Computer Vision, vol. 106, no. 1, pp. 9–30, 2014. View at: Publisher Site  Google Scholar
 Y. Hu, N. Wang, D. Tao, X. Gao, and X. Li, “SERF: a simple, effective, robust, and fast image superresolver from cascaded linear regression,” IEEE Transactions on Image Processing, vol. 25, no. 9, pp. 4091–4102, 2016. View at: Publisher Site  Google Scholar
 N. Wang, X. Gao, L. Sun, and J. Li, “Bayesian face sketch synthesis,” IEEE Transactions on Image Processing, vol. 26, no. 3, pp. 1264–1274, 2017. View at: Publisher Site  Google Scholar
 K. Krissian, R. Kikinis, C. F. Westin, and K. Vosburgh, “Speckleconstrained filtering of ultrasound images,” in 2005 IEEE Computer Society Conference on Computer Vision and Pattern Recognition (CVPR'05), pp. 547–552, San Diego, CA, USA, June 2005. View at: Publisher Site  Google Scholar
 C. A. Z. Barcelos and L. E. S. R. Vieira, “Ultrasound speckle noise reduction via an adaptive edgecontrolled variational method,” in 2014 IEEE International Conference on Systems, Man, and Cybernetics (SMC), pp. 145–151, San Diego, CA, USA, October 2014. View at: Publisher Site  Google Scholar
 F. Knoll, K. Bredies, T. Pock, and R. Stollberger, “Second order total generalized variation (TGV) for MRI,” Magnetic Resonance in Medicine, vol. 65, no. 2, pp. 480–491, 2011. View at: Publisher Site  Google Scholar
 M. Lysaker, A. Lundervold, and XueCheng Tai, “Noise removal using fourthorder partial differential equation with applications to medical magnetic resonance images in space and time,” IEEE Transactions on Image Processing, vol. 12, no. 12, pp. 1579–1590, 2003. View at: Publisher Site  Google Scholar
 K. Bredies, K. Kunisch, and T. Pock, “Total generalized variation,” SIAM Journal on Imaging Sciences, vol. 3, no. 3, pp. 492–526, 2010. View at: Publisher Site  Google Scholar
 F. Li, C. Shen, J. Fan, and C. Shen, “Image restoration combining a total variational filter and a fourthorder filter,” Journal of Visual Communication and Image Representation, vol. 18, no. 4, pp. 322–330, 2007. View at: Publisher Site  Google Scholar
 A. Chambolle and P. L. Lions, “Image recovery via total variation minimization and related problems,” Numerische Mathematik, vol. 76, no. 2, pp. 167–188, 1997. View at: Publisher Site  Google Scholar
 S. Setzer, G. Steidl, and T. Teuber, “Deblurring Poissonian images by split Bregman techniques,” Journal of Visual Communication and Image Representation, vol. 21, no. 3, pp. 193–199, 2010. View at: Publisher Site  Google Scholar
 X. Liu and L. Huang, “A new nonlocal total variation regularization algorithm for image denoising,” Mathematics and Computers in Simulation, vol. 97, no. 1, pp. 224–233, 2014. View at: Publisher Site  Google Scholar
 S. Osher, “Nonlocal operators with applications in imaging,” SIAM Journal on Multiscale Modeling and Simulation, vol. 7, no. 3, pp. 1005–1028, 2008. View at: Google Scholar
 Z. Ren, C. He, and Q. Zhang, “Fractional order total variation regularization for image superresolution,” Signal Processing, vol. 93, no. 9, pp. 2408–2421, 2013. View at: Publisher Site  Google Scholar
 J. S. Moll, “The anisotropic total variation flow,” Mathematische Annalen, vol. 332, no. 1, pp. 177–218, 2005. View at: Publisher Site  Google Scholar
 H. Chen, C. Wang, Y. Song, and Z. Li, “Split Bregmanized anisotropic total variation model for image deblurring,” Journal of Visual Communication and Image Representation, vol. 31, pp. 282–293, 2015. View at: Publisher Site  Google Scholar
 B. Chen, S. Huang, Z. Liang, W. Chen, and B. Pan, “A fractional order derivative based active contour model for inhomogeneous image segmentation,” Applied Mathematical Modelling, vol. 65, no. 1, pp. 120–136, 2019. View at: Publisher Site  Google Scholar
 B. Chen, Q. H. Zou, W. S. Chen, and Y. Li, “A fast regionbased segmentation model with Gaussian kernel of fractional order,” Adv. Math. Phys., vol. 2013, article 501628, 7 pages, 2013. View at: Publisher Site  Google Scholar
 B. Chen and W. S. Chen, “Noisy image segmentation based on wavelet transform and active contour model,” Applicable Analysis, vol. 90, no. 8, pp. 1243–1255, 2010. View at: Publisher Site  Google Scholar
 B. Chen, Q. Zou, and Y. Li, “A new image segmentation model with local statistical characters based on variance minimization,” Applied Mathematical Modelling, vol. 39, no. 12, pp. 3227–3235, 2015. View at: Publisher Site  Google Scholar
 B. Chen, Y. Li, and J. Cai, “Noisy image segmentation based on nonlinear diffusion equation model,” Applied Mathematical Modelling, vol. 36, no. 3, pp. 1197–1208, 2012. View at: Publisher Site  Google Scholar
 B. Chen, J. L. Cai, W. S. Chen, and Y. Li, “A multiplicative noise removal approach based on partial differential equation model,” Mathematical Problems in Engineering, vol. 2012, Article ID 242043, 14 pages, 2012. View at: Publisher Site  Google Scholar
 N. Costanzino, “Structure inpainting via variational methods,” 2002, https://www.lems.brown.edu/nc. View at: Google Scholar
 J. Fehrenbach and J. M. Mirebeau, “Sparse nonnegative stencils for anisotropic diffusion,” Journal of Mathematical Imaging and Vision, vol. 49, no. 1, pp. 123–147, 2014. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2020 Bo Chen 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.