Abstract
Fractional differential models are playing a vital role in many applications such as diffusion, probability potential theory, and scattering theory. In this study, the variableorder space and time fractional diffusion model is employed for denoising the medical images. The finite difference approach is implemented to find the numerical solution of the proposed model. Convergence and stability of the numerical method are presented. The experimental outcomes of the variableorder model are analyzed with those of the fractional and integerorder diffusion models. It was noticed that the peak signaltonoise ratio (PSNR) value is increased considerably for the proposed model.
1. Introduction
The field of fractional differential equations has drawn immense consideration towards theoretical [1–3] and applied research studies [4–6]. Fractional differential equations are found to be an effective tool used in certain mathematical models such as hydrology [7], finance [8], physics [9], and signal and image processing [10–14]. Various theories of fractional integrals and derivatives were developed by many authors, for instance, RiemannLiouville [2], Grnwald [15], Caputo [16], and Riesz [17].
Many researchers are focusing on image denoising based on total variation, wavelet transform, bilateral filter, histogram of gradient, fractional , anisotropic diffusion filter, and primaldual algorithm [18–25]. However, the drawback of using a total variation, secondorder, and fourthorderbased image denoising models suffered from the staircase effect, too much of smoothness and preserving in discontinues. The integerorder fractional derivatives are not efficient in specifying some complex diffusion processes. So, our work aims to further development of the variableorder fractional model in the context of image denoising to overcome the above issues.
Nowadays, variableorder fractional calculus is especially acknowledged as a helpful and hopeful approach in the modeling of a diffusion process. The pioneering work of variableorder operators can be traced. Chen et al. [26] introduced Caputofractional derivatives for removing noisy signal using wavelet transform. In recent years, a finite difference technique for solving multivariableorder fractional integrodifferential equations is introduced using Bernstein basis functions together with the NewtonCotes collocation points in [27]. Furthermore, they have found a numerical scheme to solve variableorder diffusionwave and differential equation. Also, in [28], Guo et al. proposed threedimensional fractional total variation under tensor algebrabased model for image denoising. Recently, in [29], Garrappa et al. investigated the main properties of the emerging variableorder operators and discussed some practical applications of the variableorder Scarpi integral and derivative. Gu et al. [30] proposed an unconditionally stable implicit difference scheme for solving generalized timespace fractional diffusion equations with variable coefficients with numerical scheme that utilizes the type formula for the generalized Caputofractional derivative in time discretization and the secondorder weighted and shifted Grnwald difference formula in spatial discretization. Also, Fang et al. [31] developed a fast finite difference method for solving a class of variableorder time fractional diffusion equations. In [32], Gu and Wu considered a class of Volterra partial integrodifferential problems and proposed a novel iterative algorithm for parallelintime pattern (PinT) computation.
In this research, we study the variableorder fractional diffusion model for medical image denoising using the Caputo finite difference scheme for the proposed problem. The experiments demonstrate the advantage of the variableorder fractional model that achieves the highest PSNR value, and it indicates the quality enhancement of the medical images. This research is organized as follows. In Section 2, few types of variableorder fractional derivatives are explored. In Section 3, the discrete approximation of the proposed model is developed. In Section 4, the stability and convergence of this model are discussed. In Section 5, the numerical results are compared with fractional and integerorder diffusion models. This study concludes in Section 6 with directions for further research.
2. VariableOrder Fractional Operators
The research on fractional calculus is pursued over a long time in different disciplines such as biomedical, computational biology, economics, and control engineering. The classical calculus extends the definition of fractional calculus, where the orders need not be positive integers. On the other hand, the variableorder calculus is a hereditary extension of the integer order calculus. In this regard, the order may use in any variable such as time and space variables or a system of other variables. The basic definitions for variableorder fractional derivatives are provided in the following sequel.
Definition 1 (See [16]). The variableorder Caputo derivative with respect to time is defined as follows:where .
Definition 2 (See [16]). The variableorder Caputo derivative with respect to space is defined as follows:where .
3. VariableOrder Fractional Diffusion EquationBased Image Denoising
Consider is a closed domain in . The format of noisy image would be mathematically modeled aswhere and are the observed and clean images, respectively, is the Gaussian noise, and takes the location with a rectangular image domain . The proposed work approximates the required clean image by solving the following variableorder fractional diffusion equation:where , , and , , where , , , and . Define is in the case , as the grid size in direction, as the grid size in direction, and , where is an image domain with initial and Neumann boundary conditions:
Now, the above Neumann boundary conditions are discretized aswhere is the dimension of the image . We use these Neumann boundary conditions to preserve continuity of the image boundary. The boundary values of the image can be obtained by reflecting the closest samples inside the image region. The solution to the proposed problem can be obtained by discretization of space and time variables. The discretization of the Caputotype variableorder space fractional derivative can be stated as follows [33].where and .where and .
Similarly, the discretization of the Caputotype variableorder time fractional derivative can be stated as follows [2].where and .
Discretize equation (4) in the following form:
Now, equation (8) can be written as
Equation (9) can be entered at time in the matrix form:where , , and is a matrix with the following coefficients:
The computational solutions of and are described from the identified initial and symmetric boundary conditions (3) and (4). The solution of (8) is meant as a denoised image.
The implementation of the proposed method is expressed in the following Algorithm 1:

4. Stability and Convergence Analysis
4.1. Stability Analysis
Consider and be the numerical solutions of (6) with initial values given by and , respectively.
Theorem 1. The explicit approximation method defined by (9) to variableorder spacetime diffusion equation (2) is stable, that is,where is a positive number independent of , , and .
Proof. Define . From (6),whereWe note that for any .
Let . From equation (15), . Then,The stability via mathematical induction is analyzed in [16]. From (1), , where C is a constant.
From (11), , which proves that the explicit scheme is stable. If there is a perturbation in , then the small change would not cause large error in numerical solution.
4.2. Convergence Analysis
Theorem 2. Let be the numerical solution of equations (4)–(6) at mesh point and computed by equations(5)(6), and (8). Then, there is a positive constant C, such that , where , and .
Proof. Define and . Using and , substitution into equation (16) leads toand , where is a positive constant . Using mathematical induction, we provide the convergence analysis as follows. For ,We have ,Suppose that , and letbecauseHence, there is a constant C, . If is finite.
Thus, we see that for any , and as , and approach 0 in such a way that tends to . This proves that converges to as , and tends to zero [34]. Hence, the conclusion follows.
5. Experimental Results and Discussion
The effect of this variableorder model is estimated by calculating the numerical indicator PSNR value, which is frequently used to find the quality of the processed image and also compute high quality of restoration in image firmness. The PSNR is described via the mean square error (MSE) for two images, that is to say, (corrupted image) and (denoised image), respectively.
The PSNR and MSE values of the denoised image are determined which act as a computable quality for comparison of the proposed model, fractionalorder, and integerorder diffusion models (Table 1). Figure 1 shows the relationship between , and PSNR values. Here, = 250, that is, after 250 iterations, we obtained the resulting image which is very close to the input image. The computational time (CPU) is given in Table 1. To this end, we performed on the benchmark image databases where a noisy image is evolved by using the fractionalorder and integerorder diffusion models and the proposed model with () fractionalorder derivative would attain the spotting of the inflection point without any relocation. The parameter shows Gaussian noise, and each pixel in the image will be changed from its original value by a small amount. The terminating formula for the iteration is . We studied in the numerical experiment. In experimental and numerical analyses, we have used the version of MATLAB R2012a. The proposed model has been studied with four types of medical images, namely, OCT (optical coherence tomography) image, MRI (magnetic resonance image), CT (computed tomography) scan, and Xray image. From Figures 2(e), 3(e), 4(e), and 5(e), the proposed model provides affirmative evidence for better performance in maintaining the details and edge information than the fractional and integerorder models.
(a)
(b)
(c)
(d)
(e)
(a)
(b)
(c)
(d)
(e)
(a)
(b)
(c)
(d)
(e)
(a)
(b)
(c)
(d)
(e)
Figure 2 (OCT image) shows us beneath the surface of the retina. The OCT is the most valuable advancement in the retinal diagnostic image, and we can greater understand the very fine changes which can be indicated abnormally. It constructs a crosssectional view of ocular structures accurate to less than 10 microns. So, denoising these types of images is still challenging. Hence, the retinal diagnostic image (Figure 2(a)) is used for the first experiment. Figure 2(a) is a corrupted image by an additive Gaussian noise at four various intensities such as 10, 15, 20, and 25, as shown in Figure 2(b). The denoised image 1(e) by the presented model is added clear, and it has a larger PSNR value than Figure 2(c) and Figure 2(d). Also, Figure 2(d) shows the information of retina that is critical in picking up the most subtle changes same from Figure 2(a). Therefore, the proposed model can be used for enhanced eye examination and in addition to all the normal tests.
Figure 3 (MRI image) shows a brain tumor image. Recently, biomedical images are taking an important role in identifying anatomy and the physiological process of the body in both health and disease. This MRI image has proven successful for the diagnosis of all parts of the body including cancer, stroke, heart and vascular disease, breast disease, and musculoskeletal disorders. Detection of tumors from MRI data is tedious for physicians and challenging for computers in all medical disciplines. MRI scans detect a large number of lesions and define the location more readily, and they are also the best at detecting spread to the meninges, the lining around the brain (and spinal cord). Figure 3(a) is the brain image affected by the tumor. Figure 3(b) is affected by additive Gaussian noise at four different values such as 10, 15, 20, and 25. Figures 3(c) and 3(d) have lost the information as well as the appearance of the tumor is not clear. Figure 3(e) shows original information without any loss and can be used for diagnosis and surgical purposes.
Figure 4 (CT image) shows a lung cancer image. This CT image is a diagnostic image that is used to test and construct detailed images of internal organs, soft tissue, bones, and blood vessels. A chest CT scan can assist determine the cause of lung symptoms such as shortness of breathing or chest pain or check lung problems such as a tumor and excess fluid around the lung. As shown in Figure 4(a), lung image affected by cancer and CT lung screening is capable of detecting lung nodules as small as 2 or 3 millimeters. Figure 4(b) is corrupted by additive Gaussian noise at four different levels such as 10, 15, 20, and 25. In Figure 4(c) and Figure 4(d), small nodules and affected areas are not visible. Figure 4(e) shows malignant tumors when they are still small and can be removed before the disease spreads to other areas of the body.
Figure 5 (Xray image) shows a spider image. Xray image has enhanced an important method for visualizing cellular and histological structures in a broad range of biological and medical studies. Figure 5(a) is a highresolution Xray image of a spider. Figure 5(b) has been perturbed by additive Gaussian noise at 10, 15, 20, and 25. Figures 5(c) and 5(d) are not making possible to see smaller details. But, Figure 5(e) is making it possible to see similar details of Figure 5(a).
6. Conclusion
In this work, the convergence and stability of the numerical method are presented. The remarkable difference between integer, fractional, and variableorder fractional models of diffusion equation can be seen easily. The variableorder fractional model describes the characteristics of denoising with more accuracy compared to integer and fractionalorder. Thus, the fundamental goal of this work to construct an image denoising algorithm 1 for the variableorder fractional diffusion equation in space and time by using finite difference approximation is achieved. The present work shows the validity and great potential of variableorder fractional diffusion equation for denoising the noised images. The experimental results show that the quality of denoising images and the highest PSNR and least MSE values are obtained by a fractional variableorder diffusion model.
Data Availability
The data used to support the findings of this study are included within the article.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
The work was supported by the National Research Foundation of Korea (NRF) grant funded by the Korea government (MSIT) (2021R1F1A1048937).