Next Article in Journal
Ensemble Learning for Multi-Label Classification with Unbalanced Classes: A Case Study of a Curing Oven in Glass Wool Production
Next Article in Special Issue
Modeling the Solution of the Pursuit–Evasion Problem Based on the Intelligent–Geometric Control Theory
Previous Article in Journal
Lipschitz Transformations and Maurey-Type Non-Homogeneous Integral Inequalities for Operators on Banach Function Spaces
Previous Article in Special Issue
Stochastic Time Complexity Surfaces of Computing Node
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Hybrid Lattice Boltzmann Model for Nonlinear Diffusion and Image Denoising

Federal Research Center “Computer Science and Control” of the Russian Academy of Sciences, Vavilova-44, 2, 119333 Moscow, Russia
Mathematics 2023, 11(22), 4601; https://doi.org/10.3390/math11224601
Submission received: 9 October 2023 / Revised: 1 November 2023 / Accepted: 7 November 2023 / Published: 10 November 2023
(This article belongs to the Special Issue Mathematical Modeling, Optimization and Machine Learning, 2nd Edition)

Abstract

:
In the present paper, a novel approach for image denoising based on the numerical solution to the nonlinear diffusion equation is proposed. The Perona–Malik-type equation is solved by employing a hybrid lattice Boltzmann model with five discrete velocities. In this method, the regions with large values of the diffusion coefficient are modeled with the lattice Boltzmann scheme for which hyper-viscous defects are reduced, while other regions are modeled with the conventional lattice Boltzmann model. The new method allows us to solve Perona–Malik-type equations with relatively large time steps and good accuracy. In numerical experiments, the removal of salt and pepper, speckle and Gaussian noise is considered. For salt and pepper noise, the novel scheme yields better peak signal-to-noise ratios in image denoising problems compared to the standard lattice Boltzmann approach. For certain non-small values of time steps, the novel model shows better results for speckle and Gaussian noise on average.

1. Introduction

The lattice Boltzmann (LB) method has been widely used for solving interdisciplinary problems [1,2]. In particular, one of several possible applications for the LB method is in the simulation of advection–diffusion and diffusion processes with applications to image processing.
The first LB models for the convection–diffusion processes have been studied by Flekkøy [3]. Stability domains and the accuracy of several single-relaxation time LB models have been explored in a very detailed way [4,5,6]. Numerical experiments show that four and five-velocity LB models are at least as accurate as nine-velocity models for small Peclet numbers; hence, for diffusion problems, the models with small velocity sets can be preferred. In addition, for pure diffusion, LB models are stable for any time step. This property distinguishes popular LB models from the explicit difference schemes for the diffusion equation, which have an upper stability limit on time steps. Although one can take large time steps while modeling a diffusion process using conventional LB models, the limitations emerge from accuracy issues. For non-small values of the time step and diffusion coefficient, LB models suffer from hyper-viscous errors (or truncation errors) [7,8,9,10,11,12]. In order to improve the accuracy of LB schemes, several approaches can be adopted. A fourth-order accurate scheme for the diffusion equation in one dimension can be obtained by modifying the proportion of the rest particles in three velocity LB models [13]. In a multiple-relaxation-time (MRT) LB method [14,15], a combination of relaxation times can be tuned in such a way that hyper-viscosity defects are diminished and stability is improved [16,17,18,19,20]. A fourth-order accurate MRT scheme has been proposed recently [21,22], but this model operates in one spatial dimension. On the other hand, in several spatial dimensions, hyper-viscosity reduction can be achieved only for special values of diffusion coefficient, and in some cases, single-relaxation-time models show better precision than two-relaxation-time models [23].
The LB method has found several applications in image processing. The Perona–Malik diffusion [24,25], which is widely used for denoising, can be modeled using nine-velocity and five-velocity LB schemes [26,27]. Appropriate LB schemes with external source terms can be constructed for the more complicated variants of the Perona–Malik model, in which the diffusion coefficient is also governed by an additional reaction–diffusion equation [28]. The Rudin–Osher–Fatemi total variation minimization equation for image restoration [29] can be efficiently modeled using a nine-velocity LB scheme [30]. Nine-velocity MRT LB models are applied for the simulation of a nonlinear reaction–advection–diffusion equation with a constant diffusion coefficient and non-constant advection velocity, with applications in sea clutter denoising in marine radar images [31]. LB models simulating anisotropic diffusion equations can be used for the smoothing and segmentation of two-dimensional and three-dimensional medical images [32,33,34,35]. A four-velocity LB model with a non-local collision step has been applied for the development of Gaussian pyramids and feature detection [23].
In this study, a five-velocity single-nonconstant-relaxation-time LB model with non-local collision term is considered. Since the similar LB scheme with constant relaxation time shows small hyper-viscous errors while modeling the diffusion equation with non-small diffusion coefficients and time step, one expects that this model can be used for the efficient simulation of the Perona–Malik diffusion. By using the Chapman–Enskog expansion, it is shown that this model simulates the nonlinear diffusion equation, which is very similar to the Perona–Malik equation. As an application of the model, a hybrid approach is proposed: boundaries and regions with small diffusion coefficients are modeled with the conventional five-velocity LB model, while other regions are modeled with the model with a non-local collision term. Numerical experiments for the removal of salt and pepper, speckle and Gaussian noise are carried out. It is demonstrated that in several cases, the proposed hybrid method shows a better peak signal-to-noise ratio (PSNR) than the conventional five-velocity LB model applied with the same time step.
This paper is organized as follows. In Section 2, the LB model with non-local collision kernel is considered. By applying multiple-scale Chapman–Enskog expansion, it is demonstrated that the model simulates a Perona–Malik-type diffusion. In Section 3, the application of the hybrid method and boundary conditions are discussed. Numerical experiments are carried out: four test images are considered, the results of the removal of salt and pepper, speckle and Gaussian noise are presented. In the Discussion section, the main results and possible directions for the future work are outlined.

2. Lattice Boltzmann Scheme for Nonlinear Diffusion

Assume that a two-dimensional grayscale image color intensity is given by a function I 0 ( x , y ) . The application of the nonlinear Perona–Malik filter is equivalent to solving the following equation [24,25]:
I ( t , x , y ) t = d i v D ( | I ( t , x , y ) | 2 ) I ( t , x , y ) ,
where I ( t , x , y ) | t = 0 = I 0 ( x , y ) and
D ( s ) = 1 1 + ( s / λ ) 2
and λ is a positive constant (contrast parameter). It is worth mentioning that the Perona–Malik evolution is a forward diffusion in the regions of a constant or small color intensity but is backward diffusion in the regions of sharp intensity variations (larger than λ ). As a result, similarly to the Gaussian filtering, the Perona–Malik filter blurs small-scale fluctuations, but in contrast, tends to preserve edges; this is obviously a very desirable property.
In order to solve Equation (1) numerically, the appropriate LB model is adopted. For modeling the advection–diffusion process in two dimensions, it is sufficient to employ LB models with four, five or nine velocities. Recent studies indicate that a five-velocity model ( D 2 Q 5 ) shows a good trade-off between accuracy, stability and computation complexity, for instance, five-velocity models with small Peclet numbers show better precision than nine-velocity LB models [6,36].
Then, in the present study, a five-velocity lattice is adopted, and the corresponding particle velocities are as follows: c 0 = ( 0 , 0 ) , c 1 = ( 0 , c ) , c 2 = ( c , 0 ) , c 3 = ( 0 , c ) , c 4 = ( c , 0 ) . Each time iteration for the LB equation is solved in two steps:
f i * ( t , x ) = f i ( t , x ) + J i ( t , x ) , i = 0 4 ,
f i ( t + Δ t , x + c i Δ t ) = f i * ( t , x ) , i = 0 4 ,
where Equation (3) is the collision step and (4) is the streaming step; f i , f i * are the distribution functions before and after collision, respectively; and J i is the collision kernel, whose explicit form is discussed further. Moreover, x = ( x , y ) , t are the spatial and time variables, whereby Δ t is the time step and the lattice spacing is Δ x = c Δ t or c = Δ x / Δ t . In this study, Δ x = 1 .
In the present study, two collision kernels are used. The first one is a standard local collision term for the D 2 Q 5 model:
f i * ( t , x ) = f i ( t , x ) + Δ t τ ( f e q ( t , x f i ( t , x ) ) , i = 1 5 ,
where f e q = ( 1 / 5 ) I . Using Chapman–Enskog expansion, it has been shown that the LB model with the collision term in the form (5) reproduces the Perona–Malik equation [27], where D = 2 c 2 5 ( τ Δ t / 2 ) .
The second considered model has a non-local collision term [23]:
f i * ( t , x ) = 1 Δ t τ 2 f i ( t , x c i Δ t ) + + Δ t τ 1 Δ t τ f e q ( t , x c i Δ t ) + Δ t τ f e q ( t , x ) , i = 1 5 .
Note that the model (6) shows significantly better accuracy for diffusion problems than the conventional model (5) for non-small time steps Δ t due to reduced hyper-viscous errors. Compared with the conventional single-relaxation-time LB models [4,5,6], the model (6) differs only in the form of their collision step. In the streaming step, the local equilibrium (with zero advecting velocity) remains unchanged. The conventional LB models employ local collisions, i.e., the after-collision values of the distribution function depend on the data at the same spatial node. In the present approach, the collision step also depends on the values of the distribution function in the adjacent nodes. In the Chapman–Enkog expansion, this affects the proportional second-order and higher-derivative terms, this affects the second-order and higher-derivative terms and allows us to reduce hyper-viscous errors in the corresponding nonlinear diffusion equation.
One needs to evaluate the equation for I in the case of non-constant τ , depending on I ( t , x , y ) . Following the Chapman–Enskog method [37], let us assume that the solution can be expanded on a small parameter ε :
f i ( t , x ) = I ( t 1 , t 2 , , x 1 ) / 5 + ε f i ( 1 ) ( t 1 , t 2 , , x 1 ) + ε f i ( 2 ) ( t 1 , t 2 , , x 1 ) + ,
where t m ε m t , x m ε m x . The non-equilibrium components f i ( n ) , n > 0 of the distribution function f i satisfy the relations: (a) the functions f i ( n ) , n > 0 depend on time and spatial variables via I and its derivatives; and (b) i f i ( n ) = 0 , n > 0 .
By applying the Taylor expansion on ε up to the second-order terms, and collecting the terms standing at the same powers of ε , one obtains a set of equations. The terms proportional ε 1 yields
Δ t ( t 1 + ( c i , 1 ) ) f ( e q ) + f i ( 1 ) = Δ t 1 Δ t τ ( c i , 1 ) f ( e q ) + 1 Δ t τ 2 f i ( 1 ) ,
where 1 x 1 . Taking the sum over the index i and remembering that i f i ( 1 ) = 0 , one concludes that t 1 f e q = 0 , and hence, t 1 f i ( n ) = 0 , n > 0 ; and then one obtains from (8):
f i ( 1 ) = τ ( c i , 1 ) f e q .
Collecting the ε 2 terms, one obtains
f i ( 2 ) + Δ t ( c i , 1 ) f i ( 1 ) + Δ t t 2 + 1 2 Δ t 2 ( c i , ) 2 f ( e q ) = 1 Δ t τ 2 × × f i ( 2 ) Δ t ( c i , ) f i ( 1 ) + 1 2 Δ t 2 ( c i , ) 2 f ( e q ) + Δ t 2 2 Δ t τ 1 Δ t τ ( c i , ) 2 f ( e q ) .
Using (9) and taking the sum over i, the following equation is derived from (10):
t I = c 2 Δ t 2 5 τ d i v ( I ) + 2 c 2 5 2 2 Δ t τ + Δ t τ 2 d i v τ I ,
where t ε t 1 + ε 2 t 2 and ε 1 .
In the case of constant τ , one obtains the diffusion equation
t I = 4 c 2 τ 5 1 Δ t 2 τ 2 d i v ( I ) ,
where the diffusion coefficient D in (12) equals
D = 4 c 2 τ 5 1 Δ t 2 τ 2 .
Using (13), one can express τ in terms of D by solving a quadratic equation (only the largest root τ is used since the other one leads to Δ t / τ > 1 , which can cause instability):
τ = Δ t 2 1 + 5 D 4 c 2 Δ t + ( 5 D 4 c 2 Δ t + 1 ) 2 1 .
Formula (14) will be also used in the case of D depending on I .
Note that (11) is not exactly the Perona–Malik Equation (1) but one can demonstrate that its properties are very similar to (1). First, one needs to introduce an orthonormal frame consisting of two vectors: one vector defined by T is parallel to I and another one N is perpendicular to I . These vectors have the components T = 1 | I | ( I x , I y ) and N = 1 | I | ( I y , I x ) . By straightforward computations, one can convince that the second-order derivatives of I along T and N are as follows:
I T T = 1 | I | 2 ( I x 2 I x x + 2 I x I y I x y + I y 2 I y y ) , I N N = 1 | I | 2 ( I y 2 I x x 2 I x I y I x y + I x 2 I y y ) .
Employing (15) and remembering that τ τ ( s ) τ ( | I | ) , one eventually obtains
t I = 4 c 2 τ 5 | I | 2 1 Δ t 2 τ 2 I N N + + 1 | I | 2 4 c 2 τ 5 1 Δ t 2 τ 2 + 2 c 2 5 | I | d τ d | I | 2 2 Δ t τ + Δ t 2 τ 2 I T T ,
or
t I = D N N I N N + D T T I T T ,
where the diffusion coefficients D N N , D T T account for the diffusion processes in the directions perpendicular and parallel to I , respectively. The relaxation time τ is expressed via D ( | I | ) using (14). The exact forms of D N N , D T T are as follows:
D N N = 4 c 2 τ 5 | I | 2 1 Δ t 2 τ 2 , D T T = 4 c 2 τ 5 | I | 2 1 Δ t 2 τ 2 + 2 c 2 5 | I | d τ d | I | 2 2 Δ t τ + Δ t 2 τ 2 .
If λ in (2) is smaller than 1, then the coefficients D N N = D = 1 / ( 1 + ( I / λ ) 2 ) , and D T T behaves similarly to the ones for the Perona–Malik equation. For the Perona–Malik equation, the corresponding coefficients are [25]:
D N N = D = 1 1 + ( I / λ ) 2 , D T T = 2 | I | d D d ( | I | ) + D .
In Figure 1, the coefficients D N N , D T T computed from (18) and (19) are presented against the values of s I , λ = 0.2 . One can see that the diffusion in the direction perpendicular to I is the same for both models. In the direction tangential to I , the diffusion is qualitatively similar: the regions with relatively small I (“ramps”) are smoothed ( D T T > 0 ), while the regions with non-small I (edges) are sharpened since D T T < 0 . Thus, one can conclude that the model (6) leads to Perona–Malik-type diffusion.

3. Applications of Image Denoising

3.1. The Hybrid Method

The previous study [23] shows that the LB model with non-local collision term (6) has significantly smaller hyper-viscous errors than the scheme with the local collision step (5). This means that the model (6) has better accuracy than (5) in the case of non-small time step and diffusion coefficient. Then, one expects that the denoising algorithm based on (6) can be applied for a larger time step, and hence, a smaller number of time iterations are required to model the diffusion process efficiently.
The simulations are performed with the following hybrid LB approach. For the regions with a small diffusion coefficient, they are defined by the condition Δ t / τ > 1 , where τ is given by (14), and for the boundary nodes, the collision step is local (5). Internal regions with a non-small diffusion coefficient ( Δ t / τ 1 ), which are prone to hyper-viscous errors are modeled using the non-local collision step (6). Note that the condition Δ t / τ 1 guarantees the stability of the LB scheme with the non-local collision step (6). For instance, for Δ t = 1 , Δ x = 1 , the condition Δ t / τ < 1 is satisfied when D > 0.2 . In practice, approximately 80–95% nodes are covered with the non-local collision term if Δ t = 1 , and 70–90% if Δ t = 0.5 .
At the moment t = 0 , the values of the distribution functions f i are taken as I 0 ( x , y ) / 5 . Now, consider the formulation of the boundary conditions. Assume that one needs to find f i ( t + 1 , x 0 ) , where index i corresponds to the lattice velocity, such that the node x 0 c i Δ t lies beyond the considered spatial domain (ghost node); therefore, the value of f i * ( x 0 c i Δ t ) is undefined. In this case, the extrapolation scheme is used
f i * ( t , x 1 ) = 2 f i * ( t , x 0 ) f i * ( t , x 1 ) ,
where 1 , 0 , 1 define the ghost node ( x 0 c i Δ t ), boundary node x , and the first node in the interior of the spatial domain ( x 0 + c i Δ t ), respectively.
In addition, the gradient I is approximated using central differences. It is worth mentioning that it is also possible to approximate first spatial derivatives using the non-equilibrium part of the distribution function f i n e q = f i f i e q [28]. The contrast parameter λ is taken to equal 0.1 max ( | I | ) .

3.2. Numerical Experiments

In the simulations, peak signal-to-noise ratio (PSNR) is measured over time variable t, and PSNR is defined as
P S N R ( t ) = 10 log 10 max ( I ) 2 1 N x N y x y ( I ( t , x , y ) I ¯ ( x , y ) ) 2 ,
where N x , N y is image resolution, I ¯ is the image gray color intensity without noise, and I | t = 0 I 0 is the considered image with noise. Salt and pepper, speckle and Gaussian noise are considered. The salt and pepper noise distorts images in such way that some pixels take maximal or minimal color intensity values (randomly); in this study, the proportion of such pixels is set to 5 % . The speckle noise modifies the image by adding the term n I , where n takes values from the uniform distribution, with zero mean and variance equal to 0.05 in this study. The applied Gaussian noise has zero mean and the variance is taken to equal 0.05 . The simulations are performed in Matlab (version 2022b), where four Matlab built-in gray-colored demo images are considered, and noise is added using the imnoise function.
The visual difference between the results of the hybrid method based on the collision terms (5), (6), and the standard model (5) is mostly pronounced for the synthetic images and salt and pepper noise. In Figure 2 and Figure 3, the applications of the hybrid method and the model D 2 Q 5 for the denoising of the synthetic images “circlesBrightDark” ( 512 × 512 ) and ”threads” ( 500 × 500 ) are presented. Clearly, the hybrid approach provides significantly better results for non-small time steps. For the other images—“Flamingos” ( 1296 × 972 ), “Parkavenue” ( 2048 × 1536 ), and speckle and Gaussian noise (Figure 4 and Figure 5)—the results of both methods are very similar. Hence, in order to assess noise reduction quantitatively, PSNR values are recorded for the different time steps Δ t . It is important to mention that the reduction in Δ t leads to larger values of Δ t / τ and therefore to a smaller proportion of nodes covered by the non-local collision term (6). In addition, for small Δ t , hyper-viscous errors are negligible. As a result, for Δ t < 0.5 , both LB methods present similar results. Then, the simulations are performed for Δ t = 0.5 , 1 . Also, PSNR data are provided for the D 2 Q 5 model, where Δ t = 0.1 ; the smaller values of Δ t present a small improvement in PSNR values.
From Figure 6, Figure 7, Figure 8 and Figure 9, one can observe that the hybrid model shows significantly better results than the D 2 Q 5 model when salt and pepper noise is considered. For speckle noise and Gaussian noise, the hybrid model yields noticeably greater PSNR values if Δ t = 1 , except for the “Parkavenue” image and speckle noise (Figure 9b). In the case of Δ t = 0.5 and the speckle and Gaussian noise, the hybrid LB method is only slightly better than the D 2 Q 5 model. Moreover, one should mention that for salt and pepper, speckle noise PSNR values for the hybrid method reaches the maximum value faster.

4. Discussion

In the present paper, the hybrid LB method for solving Perona–Malik-type diffusion with applications to image noise removal is proposed. Obviously, it is desirable to solve diffusion equations with as small a number of time steps as possible. Although there are no stability restrictions on the upper bound of the time step, the LB scheme’s accuracy still suffers from hyper-viscous errors when the time step and diffusion coefficients are not small. Then, boundary nodes and regions with small diffusion coefficients are modeled with the conventional five-velocity LB model (5), while other regions are modeled with that which contains a non-local collision term (6), for which hyper-viscous errors are reduced. By employing the multiple-scale Chapman–Enskog expansion, it is demonstrated that the model (6) simulates the nonlinear diffusion Equation (11), which can also be rewritten in the form (17), (18). Similarly to the Perona–Malik equation for Equations (17) and (18), the regions with relatively small intensity gradients are smoothed, while the regions with non-small intensity gradients are sharpened. In the numerical experiments, the removal of salt and pepper, speckle, and Gaussian noise from four test images is considered. On average, the hybrid approach yields better PSNR values for non-small time steps than the standard D 2 Q 5 model for the Perona–Malik equation—the improvement is very pronounced for salt and pepper noise, and Δ t = 1 . In addition, PSNR profiles for the hybrid method tends to achieve maximum values faster. The time is kept constant during any particular simulation, although it is possible to apply a variable time step; for example, at initial moments of time, the diffusion process can be modeled with small time steps, while in later moments, Δ t can be increased since the image is smoothed further. Compared to deep learning algorithms [38], the LB method (and other methods based on numerical solution to the diffusion equations) does not require large training samples of real noisy images and is able to filter either multiplicative or additive noise, while several types of neural networks are aimed to remove only a special type of noise like Gaussian with known a priori variance. In addition, the training process of neural networks requires significant memory resources.
As a final remark, one should mention that the Perona–Malik diffusion has some drawbacks. Its solutions are susceptible to “staircase” effects, i.e., the solutions form piecewise constant profiles. To alleviate this effect, several approaches have been proposed. For instance, in a ramp-preserving Perona–Malik model [39], the contrast parameter λ in (2) depends on the edge indicator function. In ramp regions, λ is increased, which results in a smoothing effect, while for the edge regions, λ takes smaller values and this yields a sharpening effect. Fourth-order diffusion equations for noise removal, on the other hand, are free from staircase effects [40,41]. In the directional Laplacian Perona–Malik model [42,43], smoothing is performed along the edge direction. Due to the flexibility of the LB method for anisotropic diffusion–advection [44,45,46,47], applications of LB schemes for the modifications of the Perona–Malik equation (especially for anisotropic extensions like the directional Laplacian method [42,43]) are of interest for future study.

Funding

This research is funded by the Ministry of Science and Higher Education of the Russian Federation, Project No. 075-15-2020-799.

Data Availability Statement

The simulation code that supports the findings of this study is available from the author upon reasonable request.

Conflicts of Interest

The author declares no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
LBLattice Boltzmann
MRTMultiple relaxation time
PSNRPeak signal-to-noise ratio

References

  1. Succi, S. The Lattice Boltzmann Equation: For Complex States of Flowing Matter; Oxford University Press: Oxford, UK, 2018. [Google Scholar]
  2. Succi, S. Lattice Fluid Dynamics: Thirty-five Years Down the Road. C. R. Méc. 2022, 350, 1–12. [Google Scholar] [CrossRef]
  3. Flekkøy, E. Lattice Bhatnagar–Gross–Krook models for miscible fluids. Phys. Rev. E 1993, 47, 4247–4257. [Google Scholar] [CrossRef] [PubMed]
  4. Wolf-Gladrow, D. Lattice Boltzmann Equation for Diffusion. J. Stat. Phys. 1995, 79, 1023–1032. [Google Scholar] [CrossRef]
  5. Suga, S. Numerical Schemes Obtained from Lattice Boltzmann Equations for Advection Diffusion Equations. Int. J. Mod. Phys. C 2006, 17, 1563–1577. [Google Scholar] [CrossRef]
  6. Suga, S. Stability and accuracy of lattice Boltzmann schemes for anisotropic advection-diffusion equations. Int. J. Mod. Phys. C 2009, 20, 633–650. [Google Scholar] [CrossRef]
  7. Holdych, D.; Noble, D.; Georgiadis, J.; Buckius, R. Truncation error analysis of lattice Boltzmann methods. J. Comput. Phys. 2004, 193, 595–619. [Google Scholar] [CrossRef]
  8. Zhao, F. Optimal relaxation collisions for lattice Boltzmann methods. Comput. Math. Appl. 2013, 65, 172–185. [Google Scholar] [CrossRef]
  9. Silva, G.; Semiao, V. Truncation errors and the rotational invariance of three-dimensional lattice models in the lattice Boltzmann method. J. Comput. Phys. 2014, 269, 259–279. [Google Scholar] [CrossRef]
  10. Bauer, M.; Silva, G.; Rúde, U. Letter to the Editor: Truncation errors of the D3Q19 lattice model for the lattice Boltzmann method. J. Comput. Phys. 2020, 405, 109111. [Google Scholar] [CrossRef]
  11. Wissocq, G.; Sagaut, P. Hydrodynamic limits and numerical errors of isothermal lattice Boltzmann schemes. J. Comput. Phys. 2022, 450, 110858. [Google Scholar] [CrossRef]
  12. Silva, G. Discrete effects on the source term for the lattice Boltzmann modelling of one-dimensional reaction–diffusion equations. Comput. Fluids 2023, 251, 105735. [Google Scholar] [CrossRef]
  13. Suga, S. An Accurate Multi-level Finite Difference Scheme for 1D Diffusion Equations Derived from the Lattice Boltzmann Method. J. Stat. Phys. 2010, 140, 494–503. [Google Scholar] [CrossRef]
  14. d’Humières, D. Generalized lattice-Boltzmann equations. In Proceedings of the 18th International Symposium, Rarefied Gas Dynamics, Vancouver, BC, Canada, 26–30 July 1992; Volume 159, pp. 450–458. [Google Scholar]
  15. Chai, Z.; Shi, B. Multiple-relaxation-time lattice Boltzmann method for the Navier-Stokes and nonlinear convection-diffusion equations: Modeling, analysis, and elements. Phys. Rev. E 2020, 102, 023306. [Google Scholar] [CrossRef] [PubMed]
  16. Servan-Camas, B.; Tsai, F. Lattice Boltzmann method with two relaxation times for advection–diffusion equation: Third order analysis and stability analysis. Adv. Water Res. 2008, 31, 1113–1126. [Google Scholar] [CrossRef]
  17. Servan-Camas, B.; Tsai, F.T.-C. Non-negativity and stability analyses of lattice Boltzmann method for advection–diffusion equation. J. Comput. Phys. 2009, 228, 236–256. [Google Scholar] [CrossRef]
  18. Ginzburg, I.; d’Humières, D.; Kuzmin, A. Optimal Stability of Advection-Diffusion Lattice Boltzmann Models with Two Relaxation Times for Positive/Negative Equilibrium. J. Stat. Phys. 2010, 139, 1090–1143. [Google Scholar] [CrossRef]
  19. Kuzmin, A.; Ginzburg, I.; Mohamad, A. The role of the kinetic parameter in the stability of two-relaxation-time advection–diffusion lattice Boltzmann schemes. Comput. Math. Appl. 2011, 61, 3417–3442. [Google Scholar] [CrossRef]
  20. Ginzburg, I. Truncation Errors, Exact and Heuristic Stability Analysis of Two-Relaxation-Times Lattice Boltzmann Schemes for Anisotropic Advection-Diffusion Equation. Commun. Comput. Phys. 2012, 11, 1439–1502. [Google Scholar] [CrossRef]
  21. Lin, Y.; Hong, N.; Shi, B.; Chai, Z. Multiple-relaxation-time lattice Boltzmann model-based four-level finite-difference scheme for one-dimensional diffusion equations. Phys. Rev. E 2021, 104, 015312. [Google Scholar] [CrossRef]
  22. Chen, Y.; Chai, Z.; Shi, B. Fourth-order multiple-relaxation-time lattice Boltzmann model and equivalent finite-difference scheme for one-dimensional convection-diffusion equations. Phys. Rev. E 2023, 107, 055305. [Google Scholar] [CrossRef]
  23. Ilyin, O. Lattice Boltzmann model for diffusion equation with reduced truncation errors: Applications to Gaussian filtering and image processing. Appl. Math. Comput. 2023, 456, 128123. [Google Scholar] [CrossRef]
  24. Perona, P.; Malik, J. Scale-space and edge detection using anisotropic diffusion. IEEE Trans. Pattern Anal. Mach. Intell. 1990, 12, 629–639. [Google Scholar] [CrossRef]
  25. Weickert, J. Anisotropic Diffusion in Image Processing; B.G. Teubner: Stuttgart, Germany, 1998. [Google Scholar]
  26. Jawerth, B.; Lin, P.; Sinzinger, E. Lattice Boltzmann Models for Anisotropic Diffusion of Images. J. Math. Imag. Vis. 1999, 11, 231–237. [Google Scholar] [CrossRef]
  27. Zhang, W.; Shi, B. Application of Lattice Boltzmann Method to Image Filtering. J. Math. Imag. Vis. 2012, 43, 135–142. [Google Scholar] [CrossRef]
  28. Chen, J.; Chai, Z.; Shi, B.; Zhang, W. Lattice Boltzmann method for filtering and contour detection of the natural images. Comput. Math. Appl. 2014, 68, 257–268. [Google Scholar] [CrossRef]
  29. Rudin, L.; Osher, S.; Fatemi, E. Nonlinear total variation based noise removal algorithms. Phys. D 1992, 60, 259–268. [Google Scholar] [CrossRef]
  30. Chang, Q.; Yang, T. A Lattice Boltzmann Method for Image Denoising. IEEE Trans. Image Process. 2009, 18, 2797–2802. [Google Scholar] [CrossRef]
  31. Michelet, J.; Tekitek, M.M.; Berthier, M. Multiple Relaxation Time Lattice Boltzmann Schemes for Advection-Diffusion Equations with Application to Radar Image Processing. J. Comput. Phys. 2022, 471, 111612. [Google Scholar] [CrossRef]
  32. Chen, Y.; Yan, Z.; Qian, Y. An Anisotropic Diffusion Model for Medical Image Smoothing by Using the Lattice Boltzmann Method. In Proceedings of the 7th Asian-Pacific Conference on Medical and Biological Engineering, Beijing, China, 22–25 April 2008; Volume 19, pp. 255–259. [Google Scholar]
  33. Chen, Y.; Yan, Z.; Shi, J. Application of Lattice Boltzmann Method to image segmentation. In Proceedings of the 29th Annual International Conference of the IEEE Engineering in Medicine and Biology Society, Lyon, France, 23–26 August 2007; pp. 6561–6564. [Google Scholar]
  34. Chen, Y.; Navarro, L.; Wang, Y.; Courbebaisse, G. Segmentation of the thrombus of giant intracranial aneurysms from CT angiography scans with lattice Boltzmann method. Med. Image Anal. 2014, 18, 1–8. [Google Scholar] [CrossRef]
  35. Nguyen, K.L.; Tekitek, M.M.; Delachartre, P.; Berthier, M. Multiple Relaxation Time Lattice Boltzmann Models for Multigrid Phase-Field Segmentation of Tumors in 3D Ultrasound Images. SIAM J. Image Sci. 2019, 12, 1324–1346. [Google Scholar] [CrossRef]
  36. Li, L.; Mei, R.; Klausner, J. Lattice Boltzmann models for the convection-diffusion equation: D2Q5 vs. D2Q9. Int. J. Heat Mass Transf. 2017, 108, 41–62. [Google Scholar] [CrossRef]
  37. Krüger, T.; Kusumaatmaja, H.; Kuzmin, A.; Shardt, O.; Silva, G.; Viggen, E. The Lattice Boltzmann Method. Principles and Practice; Springer: Cham, Switzerland, 2017. [Google Scholar]
  38. Tian, C.; Fei, L.; Zheng, W.; Xu, Y.; Zuo, W.; Lin, C.-W. Deep learning on image denoising: An overview. Neural Netw. 2020, 131, 251–275. [Google Scholar] [CrossRef] [PubMed]
  39. Chen, Q.; Montesinos, P.; Sun, Q.S.; Xia, D.S. Ramp preserving Perona–Malik model. Signal Process. 2010, 90, 1963–1975. [Google Scholar] [CrossRef]
  40. You, Y.-L.; Kaveh, M. Fourth-order partial differential equations for noise removal. IEEE Trans. Image Process. 2000, 9, 1723–1730. [Google Scholar] [CrossRef]
  41. Didas, S.; Weickert, J.; Burgeth, B. Properties of Higher Order Nonlinear Diffusion Filtering. J. Math. Imag. Vis. 2009, 35, 208–226. [Google Scholar] [CrossRef]
  42. Wang, Y.Q.; Guo, J.; Chen, W.; Zhang, W. Image denoising using modified Perona–Malik model based on directional Laplacian. Signal Process. 2013, 93, 2548–2558. [Google Scholar] [CrossRef]
  43. Wang, N.; Shang, Y.; Chen, Y.; Yang, M.; Zhang, Q.; Liu, Y.; Gui, Z. A Hybrid Model for Image Denoising Combining Modified Isotropic Diffusion Model and Modified Perona-Malik Model. IEEE Access 2018, 6, 33568–33582. [Google Scholar] [CrossRef]
  44. Ginzburg, I. Equilibrium-type and link-type lattice Boltzmann models for generic advection and anisotropic-dispersion equation. Adv. Water Res. 2005, 28, 1171–1195. [Google Scholar] [CrossRef]
  45. Li, Q.; Chai, Z.; Shi, B. Lattice Boltzmann model for a class of convection–diffusion equations with variable coefficients. Comput. Math. Appl. 2015, 70, 548–561. [Google Scholar] [CrossRef]
  46. Chai, Z.; Shi, B.; Guo, Z. A Multiple-Relaxation-Time Lattice Boltzmann Model for General Nonlinear Anisotropic Convection–Diffusion Equations. J. Sci. Comput. 2016, 69, 355–390. [Google Scholar] [CrossRef]
  47. Chen, Y.; Chai, Z.; Shi, B. A General Fourth-Order Mesoscopic Multiple-Relaxation-Time Lattice Boltzmann Model and Equivalent Macroscopic Finite-Difference Scheme for Two-Dimensional Diffusion Equations. arXiv 2023, arXiv:2308.05280. [Google Scholar]
Figure 1. Diffusion coefficients D T T and D N N computed from D in the form (2); s I ; λ = 0.2 ; (a) the model (6); the coefficients are given by Formula (18); Δ t = 0.5 ; (b) Perona–Malik model (1); the coefficients are given by Formula (19).
Figure 1. Diffusion coefficients D T T and D N N computed from D in the form (2); s I ; λ = 0.2 ; (a) the model (6); the coefficients are given by Formula (18); Δ t = 0.5 ; (b) Perona–Malik model (1); the coefficients are given by Formula (19).
Mathematics 11 04601 g001
Figure 2. “circlesBrightDark” image and application of salt and pepper noise: (a) initial image, (b) addition of salt and pepper noise with density 0.05 , (c) the result of the application of D 2 Q 5 model at moment t = 4 and Δ t = 1 , (d) the result of the application of hybrid D 2 Q 5 model at the moment t = 4 and Δ t = 1 .
Figure 2. “circlesBrightDark” image and application of salt and pepper noise: (a) initial image, (b) addition of salt and pepper noise with density 0.05 , (c) the result of the application of D 2 Q 5 model at moment t = 4 and Δ t = 1 , (d) the result of the application of hybrid D 2 Q 5 model at the moment t = 4 and Δ t = 1 .
Mathematics 11 04601 g002
Figure 3. “Threads” image and application of salt and pepper noise: (a) initial image, (b) addition of salt and pepper noise with density 0.05 , (c) the result of the application of D 2 Q 5 model at moment t = 4 and Δ t = 1 , (d) the result of the application of hybrid D 2 Q 5 model at the moment t = 4 and Δ t = 1 .
Figure 3. “Threads” image and application of salt and pepper noise: (a) initial image, (b) addition of salt and pepper noise with density 0.05 , (c) the result of the application of D 2 Q 5 model at moment t = 4 and Δ t = 1 , (d) the result of the application of hybrid D 2 Q 5 model at the moment t = 4 and Δ t = 1 .
Mathematics 11 04601 g003
Figure 4. “Parkavenue” image and application of speckle noise: (a) initial image, (b) addition of speckle noise with variation 0.05 , (c) the result of the application of D 2 Q 5 model at moment t = 5 and Δ t = 1 , (d) the result of the application of hybrid D 2 Q 5 model at the moment t = 5 and Δ t = 1 .
Figure 4. “Parkavenue” image and application of speckle noise: (a) initial image, (b) addition of speckle noise with variation 0.05 , (c) the result of the application of D 2 Q 5 model at moment t = 5 and Δ t = 1 , (d) the result of the application of hybrid D 2 Q 5 model at the moment t = 5 and Δ t = 1 .
Mathematics 11 04601 g004
Figure 5. “Flamingos” image and application of Gaussian noise: (a) initial image; (b) addition of Gaussian noise with variation 0.05 ; (c) result of application of D 2 Q 5 model at moment t = 5 , Δ t = 1 ; (d) result of application of hybrid D 2 Q 5 model at the moment t = 5 , Δ t = 1 .
Figure 5. “Flamingos” image and application of Gaussian noise: (a) initial image; (b) addition of Gaussian noise with variation 0.05 ; (c) result of application of D 2 Q 5 model at moment t = 5 , Δ t = 1 ; (d) result of application of hybrid D 2 Q 5 model at the moment t = 5 , Δ t = 1 .
Mathematics 11 04601 g005
Figure 6. PSNR for “circlesBrightDark” image: (a) salt and pepper noise, (b) speckle noise, (c) Gaussian noise. The solutions to the hybrid LB model are defined as: (—□—) for Δ t = 1 and (--□--) for Δ t = 0.5 . The solutions to LB D 2 Q 5 model are defined as: (—▵—) for Δ t = 1 , (--▵--) for Δ t = 0.5 , and (—) for Δ t = 0.1 .
Figure 6. PSNR for “circlesBrightDark” image: (a) salt and pepper noise, (b) speckle noise, (c) Gaussian noise. The solutions to the hybrid LB model are defined as: (—□—) for Δ t = 1 and (--□--) for Δ t = 0.5 . The solutions to LB D 2 Q 5 model are defined as: (—▵—) for Δ t = 1 , (--▵--) for Δ t = 0.5 , and (—) for Δ t = 0.1 .
Mathematics 11 04601 g006
Figure 7. PSNR for “threads” image: (a) salt and pepper noise, (b) speckle noise, (c) Gaussian noise. The solutions to the hybrid LB model are defined as: (—□—) for Δ t = 1 and (--□--) for Δ t = 0.5 . The solutions to LB D 2 Q 5 model are defined as: (—▵—) for Δ t = 1 , (--▵--) for Δ t = 0.5 , and (—) for Δ t = 0.1 .
Figure 7. PSNR for “threads” image: (a) salt and pepper noise, (b) speckle noise, (c) Gaussian noise. The solutions to the hybrid LB model are defined as: (—□—) for Δ t = 1 and (--□--) for Δ t = 0.5 . The solutions to LB D 2 Q 5 model are defined as: (—▵—) for Δ t = 1 , (--▵--) for Δ t = 0.5 , and (—) for Δ t = 0.1 .
Mathematics 11 04601 g007
Figure 8. PSNR for “parkavenue” image: (a) salt and pepper noise, (b) speckle noise, (c) Gaussian noise. The solutions to the hybrid LB model are defined as: (—□—) for Δ t = 1 and (--□--) for Δ t = 0.5 . The solutions to LB D 2 Q 5 model are defined as: (—▵—) for Δ t = 1 , (--▵--) for Δ t = 0.5 , and (—) for Δ t = 0.1 .
Figure 8. PSNR for “parkavenue” image: (a) salt and pepper noise, (b) speckle noise, (c) Gaussian noise. The solutions to the hybrid LB model are defined as: (—□—) for Δ t = 1 and (--□--) for Δ t = 0.5 . The solutions to LB D 2 Q 5 model are defined as: (—▵—) for Δ t = 1 , (--▵--) for Δ t = 0.5 , and (—) for Δ t = 0.1 .
Mathematics 11 04601 g008
Figure 9. PSNR for “flamingos” image: (a) salt and pepper noise, (b) speckle noise, (c) Gaussian noise. The solutions to the hybrid LB model are defined as: (—□—) for Δ t = 1 and (--□--) for Δ t = 0.5 . The solutions to LB D 2 Q 5 model are defined as: (—▵—) for Δ t = 1 , (--▵--) for Δ t = 0.5 , and (—) for Δ t = 0.1 .
Figure 9. PSNR for “flamingos” image: (a) salt and pepper noise, (b) speckle noise, (c) Gaussian noise. The solutions to the hybrid LB model are defined as: (—□—) for Δ t = 1 and (--□--) for Δ t = 0.5 . The solutions to LB D 2 Q 5 model are defined as: (—▵—) for Δ t = 1 , (--▵--) for Δ t = 0.5 , and (—) for Δ t = 0.1 .
Mathematics 11 04601 g009
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Ilyin, O. Hybrid Lattice Boltzmann Model for Nonlinear Diffusion and Image Denoising. Mathematics 2023, 11, 4601. https://doi.org/10.3390/math11224601

AMA Style

Ilyin O. Hybrid Lattice Boltzmann Model for Nonlinear Diffusion and Image Denoising. Mathematics. 2023; 11(22):4601. https://doi.org/10.3390/math11224601

Chicago/Turabian Style

Ilyin, Oleg. 2023. "Hybrid Lattice Boltzmann Model for Nonlinear Diffusion and Image Denoising" Mathematics 11, no. 22: 4601. https://doi.org/10.3390/math11224601

APA Style

Ilyin, O. (2023). Hybrid Lattice Boltzmann Model for Nonlinear Diffusion and Image Denoising. Mathematics, 11(22), 4601. https://doi.org/10.3390/math11224601

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