Texture Enhancement for Medical Images Based on Fractional Differential Masks

Texture enhancement for medical images is the most important technique in medical image diagnosis. This paper introduces a texture enhancement technique for medical images by using fractional differential (FD) masks based on Srivastava-Owa fractional operators. We also construct a 2D isotropic gradient mask based on generalized fractional operators. Texture enhancement performance is measured by applying experiments according to visual perception and by using Sobel/Canny edge filters and gray-level co-occurrence matrix. We discuss the capability of the FD mask for texture enhancement. The experiments and analysis show that the operator can extract subtle information and make the edges prominent.


Introduction
Fractional calculus and its applications (e.g., the theory of derivatives and integrals of any arbitrary, real, or complex order) are significant for different areas of mathematics, physics, and engineering.Fractional calculus generalizes the ideas of integer-order differentiation and -fold integration.Fractional derivatives introduce an excellent instrument for describing the general properties of various materials and processes.This feature is the main advantage of fractional derivatives compared with classical integer-order models, in which such effects are neglected.The advantages of fractional derivatives become apparent in modeling, describing the mechanical, electrical, and physical properties of real materials [1][2][3][4], and removing noise and irrelevant information from the absorption spectra [5].
Studies on fractional calculus and fractional differential (FD) equations that involve different operators, such as Riemann-Liouville operators, Erd' elyi-Kober operators, Weyl-Riesz operators, Caputo operators, and Grünwald-Letnikov operators, have evolved during the past four decades and have been extended in other fields [6].
Fractional calculus (integral and differential operators) has been employed recently in signal and image possessing.
Fractional calculations are capable of enhancing image quality for revealing faint objects in astronomical image analysis [7,8] and signal analysis [9]; this approach provides interesting possibilities for edge detection and image restoration.FD has received considerable attention in various image enhancement applications such as texture enhancement [10] and image denoising [6].
Medical images are a key investigation tool in medical diagnosis.The quality of a medical image is determined by the imaging method, equipment characteristics, and imaging variables selected by the operator.A large number of medical image acquisition devices exist such as computed tomography scanners, magnetic resonance imagers, ultrasound probes, and nuclear imagers [11,12].
Medical images usually have low resolutions, high noise, and high textures.Moreover, medical images represent an enormous amount of data.The noise in medical images renders the interpretation of images difficult.Thus, denoising is necessary before the analysis of medical images.Denoising techniques have been exploited in the field of medicine [13,14].
Texture is an important feature of natural images.Thus, various image texture applications have been subjected to intense studies by researchers [15].Image texture is defined as the function of the spatial variation in pixel intensities (gray values).Thus, image texture is one of the fundamental factors required in medical imaging analysis.Castellano et al. [16] clarified the principles of texture analysis in medical images and concluded that texture parameters are simply mathematical representations of image features, which can be categorized as smooth, rough, grainy, and so on.This result implies that texture analysis can be applied to any set of image regions that can be differentiated by such descriptions.Smith and Chang [17] defined texture as a visual pattern with homogenous properties that cannot result from the presence of only a single color or intensity.Therefore, texture parameter analysis is a useful method of increasing obtainable information from medical images.
Texture enhancement is one of the most important techniques in medical image processing.The most common method for enhancing texture is by using masks (also referred to as operators, filters, or templates), which is basically a small 2D array.Each pixel is modified in the mask according to the value of the neighborhood around the pixel of interest.The texture enhancement technique based on mask operations has become increasingly important in image processing and analysis [18].
The frequencies of an image determine the changes in gray values with changing distance.High-frequency components are characterized by large changes in gray values over small distances; edges and noise are examples of highfrequency components.By contrast low-frequency components are characterized by small changes in gray values; backgrounds and textures are examples of low-frequency components [19].
The FD mask preserves low-frequency contour features in smooth areas, nonlinearly maintains high-frequency marginal features in areas wherein the changes in gray level are considerable, and enhances texture details in areas wherein the changes in gray level are insignificant [17].
Several works have used mask-based texture enhancement by using FD.Qing-li et al. [20] constructed an 2D isotropic nonlinear filter algorithm based on a FD algorithm for image contrast enhancement.The structure and parameters of the algorithm were presented in eight directions, and fractional-order differential operators were used to control the image contrast enhancement.This nonlinear filter algorithm was implemented on various low-contrast digital images.
Yang et al. [21] showed that image edge detection, which is based on FD operators and not traditional first and second differential operators, has a higher signal-to-noise ratio.This work was a preliminary attempt in exploring fractional-order differential operators in image processing applications.Garg and Singh [22] proposed a 2D FD operator, which is an improved version of the Grümwald-Letnikov FD operator.The improved Grümwald-Letnikov FD operator enhances both the texture and brightness of the image.Pu et al. [23] proposed six FD masks.They demonstrated the most efficient fractional mask by theoretic and experimental analysis.Their experiments show that the FD-based approach is better than traditional integral differential-based algorithms when nonlinearly enhancing comprehensive texture details in the smooth areas of texture-rich digital images.
In our previous work, we proposed a generalized FD based on the generalized Savitzky-Golay filtering of Srivastava-Owa fractional operators for image texture enhancement.This algorithm was used to control the degree of texture enhancement with fractional power values.
We aim to develop a differential mask based on Srivastava-Owa fractional operators for the texture enhancement of grayscale medical images and demonstrate the capability of the FD mask by varying the FD power parameters.
This paper is organized as follows.Section 2 explains the discrete FD mask.Section 3 presents the simulation results.Finally, Section 4 concludes.

Discreet FD Mask
This section describes the mathematical background of the differential operator mask in the proposed algorithm.We constructed a generalized fractional mask by using the following generalized FD operator [24]: where the function () is analytic in simply connected region of the complex -plane C that contains the origin.The multiplicity of ( +1 −  +1 ) − is removed by requiring log( +1 −  +1 ) to be real when ( +1 −  +1 ) > 0.
We can compute the value of the FD operator in (1) by using the numerical calculation that references the discrete form.If the interval is [0, ],  is real,  is divided into equal parts, and the differential of each part is calculated and combined.Thus, we observe the following:

× (𝜇 + 1)
−1 ( − 1) However, an image is a function of two variables; thus, we can generalize these definitions to include both the  and  values.If the image is required to complete the nonlinear filter, we must convolute the pixel values of the image with a ( × ) mask: where  = ,  is the pixel value, and  is the value of the mask.The mask is designed into  ×  size matrix , which  has  layers ( is odd). has eight directions: 0 ∘ , 45 ∘ , 90 ∘ , 135 ∘ , 180 ∘ , 225 ∘ , 270 ∘ , and 315 ∘ (Figure 1).The distinct th layer of  is 2 × .From the above construction, we conclude the following: If we assume that  = 0,  ∈ (0, 1), and  ≥ 0, we obtain the following:  Specifically, when  = 0, the case of the Riemann-Liouville fractional operator is obtained: For texture enhancement, the sharpening edges must add to the value of the ordinal pixel.Thus, we impose that   =   , where  > 0 is the intensity factor:  For digital images, the coefficients of 2D FD masks can be obtained in eight directions: 180 ∘ , 90 ∘ , 0 ∘ , 270 ∘ , 45 ∘ , 135 ∘ , 315 ∘ , and 225 ∘ (Figure 1).These masks are located on the negative -coordinate, negative -coordinate, positive -coordinate, positive -coordinate, left upward diagonal, right upward diagonal, left downward diagonal, and right downward diagonal, respectively.These masks are rotation invariant and are used to describe the textures and edges in images.
The values of   are provided in (7), and  0 is calculated as follows [25]: The nine values of each FD mask are obtained by sliding the mask window over the image from the top left corner of the image block through all pixels until the FD mask fits entirely within the image boundaries.The output of each image block is eight values, which represents the texture information of each image.The dimension of the FD mask is an odd natural number.To obtain a good result, the size of the mask should be small.

Experimental Results and Discussion
This section demonstrates the texture-enhancing capability of the 2D FD discrete gradient operator based on Srivastava-Owa fractional operators.We analyze the texture-enhancing Figure 5 shows the summary of the results of the proposed enhancement for different values of  and Riemann-Liouville cases where  = 0.By using qualitative analysis, we observe that the first image with  = 0.9 is clearer than all images, even for the Riemann-Liouville case where  = 0.The proposed algorithm not only enhances the texture quality of the image but also maintains the edges present in the image.The proposed algorithm for image texture enhancement provides satisfactory results.The qualitative analysis of the proposed algorithm serves as one of the most important parameters in judging performance.
For the quantity analysis, we use GLCM to describe the comprehensive information of the texture details. Figure 6 compares the texture enhancement of the present study with that of Pu et al. [23] for  = 0.1, 0.5, and 0.9.Pu et al. [23] developed six FD masks based on Grümwald-Letnikov operators and presented the structures and parameters of each mask in eight directions.The main differences between the present study and the study of Pu et al. [23] are as follows.
(1) The FD masks of Pu et al. [23] are based on Grümwald-Letnikov operators.In the present study, texture enhancement is introduced by employing a generalized 2D FD discrete gradient operator based on Srivastava-Owa fractional operators.The discrete approximations derived from the Grünwald-Letnikov fractional derivatives present several limitations [28] as follows.
(i) Unstable numerical methods are frequently generated.(ii) The order of accuracy is never higher than one.
(2) The algorithm of Pu et al. [23] uses only one FD parameter.By contrast, a key contribution of our work is the multi-FD model.Changing the value of the fractional powers allows the FD coefficients to be adjusted for each image according to various characteristics.
The present study and the study of Pu et al. [23] have the following characteristics.
(i) The contrast of GLCM increases with increasing .
An increasing  corresponds to more pixels, greater contrast, and improved image clarity.
(ii) The correlation of GLCM decreases with increasing , and texture details are less observable.
(iii) The energy of GLCM decreases with increasing , and texture details become less noticeable.
(iv) The homogeneity of GLCM decreases with increasing .Thus, intensive textures are increased.

Conclusion
A texture enhancement technique for medical images is introduced by using FD masks based on Srivastava-Owa fractional operators.In the FD spectra, the texture and lightness of the image are enhanced by the FD mask-based algorithm.The proposed algorithm controls the degree of texture enhancement by using the fractional power parameters of the multi-FD model.Our proposed technique is not only limited to medical images but can also be applied to different image applications by considering the limitations of each imaging method.We aim to prevent the effect of noise and control the degree of texture enhancement by using filter mask parameters.The experiment results demonstrate the efficacy of the proposed algorithm.
image Sobel for FD enhanced image Canny for original image Canny for FD enhanced image
image Sobel for FD enhanced image Canny for original image Canny for FD enhanced image