`Discrete Dynamics in Nature and SocietyVolume 2012, Article ID 529849, 14 pageshttp://dx.doi.org/10.1155/2012/529849`
Research Article

## Denoising Algorithm Based on Generalized Fractional Integral Operator with Two Parameters

1Faculty of Computer Science and Information Technology, University of Malaya, 50603 Kuala Lumpur, Malaysia
2Institute of Mathematical Sciences, University of Malaya, 50603 Kuala Lumpur, Malaysia

Received 23 January 2012; Accepted 21 February 2012

Academic Editor: Garyfalos Papaschinopoulos

Copyright © 2012 Hamid A. Jalab and Rabha W. Ibrahim. 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.

#### Abstract

In this paper, a novel digital image denoising algorithm called generalized fractional integral filter is introduced based on the generalized Srivastava-Owa fractional integral operator. The structures of fractional masks of this algorithm are constructed. The denoising performance is measured by employing experiments according to visual perception and PSNR values. The results demonstrate that apart from enhancing the quality of filtered image, the proposed algorithm also reserves the textures and edges present in the image. Experiments also prove that the improvements achieved are competent with the Gaussian smoothing filter.

#### 1. Introduction

Fractional integration and fractional differentiation are generalizations of notions of integer-order integration and differentiation and include th derivatives and -fold integrals as particular cases. Many applications of fractional calculus in physics amount to replace the time derivative in an evolution equation with a derivative of fractional order. Fractional calculus has been applied to a variety of physical phenomena, including anomalous diffusion, transmission line theory, problems involving oscillations, nanoplasmonics, solid mechanics, astrophysics, and viscoelasticity [16].

Nowadays, fractional calculus (integral and differential operators) is utilized in signal processing and image possessing. The fractional calculation enhances the quality of images, with interesting possibilities in edge detection and image restoration, to reveal faint objects in astronomical images and is devoted to astronomical images analysis [7, 8]. Furthermore, fractional calculus is employed in design problems of variables [9] and in different applications in engineering [10]. Finally, the fractional calculus (integral operators) is used in image denoising [11]. All results based on the fractional calculus operators (differential and integral) show that this method is not only effective, but also has good immunity. Therefore, the fractional calculus in the field of image processing and signal prosecuting has broad application prospect.

Many studies on fractional calculus and fractional differential equations, involving different operators, such as the Riemann-Liouville operators, the Erdélyi-Kober operators, the Weyl-Riesz operators, the Caputo operators, and the Grünwald-Letnikov operators, have evolved during the past three decades with its applications in other field. Moreover, the existence and uniqueness of holomorphic solutions for nonlinear fractional differential equations such as Cauchy problems and diffusion problems in complex domain are established and posed [1221].

Denoising is one of the most fundamental image restoration problems in computer vision and image processing. In this paper, we have introduced an image denoising algorithm called generalized fractional integral image denoising algorithm based on the Srivastava-Owa fractional integral operator. The structures of fractional masks of this algorithm are constructed. The denoising performance is measured by employing experiments according to standard of visual perception and PSNR values. This paper is organized as follows. In Section 2, we introduce the generalized integral operator. In Section 3 construction of fractional integral mask, which is the novelty of this work, is presented. The experimental results are shown in Section 4. Finally, conclusion is presented in Section 5.

#### 2. Generalized Integral Operator

In [22], Srivastava and Owa have defined fractional operators (derivative and integral as follows.

The fractional derivative of order is defined, for a function , by where the function is analytic in simply connected region of the complex -plane containing the origin and the multiplicity of is removed by requiring to be real when .

The fractional integral of order is defined, for a function , by where the function is analytic in simply connected region of the complex -plane containing the origin and the multiplicity of is removed by requiring to be real when .

In [23], Ibrahim derived a formula for the generalized fractional integral, considering for natural and real the -fold integral of the form

By employing the Cauchy formula for iterated integrals yields Repeating the above step times, we have which implies the fractional operator type where and are real numbers and the function is analytic in simply connected region of the complex -plane containing the origin and the multiplicity of is removed by requiring to be real when . When , we arrive at the standard Srivastava-Owa fractional integral, which is used to define the Srivastava-Owa fractional derivatives.

Corresponding to the generalized fractional integrals (2.6), we define the generalized differential operator of order by where the function is analytic in simply connected region of the complex -plane containing the origin and the multiplicity of is removed by requiring to be real when .

#### 3. Construction of Fractional Integral Mask

Using the generalized fractional integral operator defined in (2.6), we have proceeded to construct the generalized fractional integral mask which is the main contribution of this work. Since the fractional differential operator is performed if the order is positive and the fractional integral operator is performed if the order is negative; for analytic function , we assume that , and then

For converting continuous integral into discrete sum. Integral section [] is equally transformed into the sum of the integrals of parts, and when is greater enough, an approximate formula can be obtained: Moreover, we can derive that

Substituting (3.4) into (3.2), we obtain However, in the context of image processing, (3.5) is applied uniformly in the whole digital image and therefore should be in two directions of and . Thus for two variables on the negative direction of and coordinates, we have

The next two formulae show the construction of the numerical implementation algorithm for generalized fractional integral operation in sense of the Srivastava-Owa operators

The nonzero values of corresponding terms in formula (3.7) are which are all fractional coefficients according to the generalized Srivastava-Owa fractional integral operator. Note that when , we have the Riemann-Liouville integral operator while for two variables on the positive direction of and coordinates, we have obtained

The nonzero values of corresponding terms in formula (3.10) are For digital images, 2-dimensional fractional integral filter coefficients can be obtained in eight directions of 180°, 0°, 90°, 270°, 45°, 135°, 315°, 225° as shown in Figure 1. These filters are rotation invariant and are used to describe the edges present and for removing noise. They are, respectively, on the directions of negative -coordinate, positive -coordinate, negative -coordinate, positive -coordinate, right upward diagonal, left upward diagonal and right downward diagonal, and left downward diagonal, which are all implemented. In this paper, we have applied fractional mask convolution on eight directions with the gray value of corresponding digital grayscale image pixels, adding all product terms to obtain weighting sum on eight directions. For digital color images, the same algorithm which is used for gray image can be applied but it performs separately for each of the R, G, B color components (RGB).

Figure 1: Fractional integral masks on directions of 180°, 0°, 90°, 270°, 45°, 135°, 315°, and 225°.

The fractional integral filter coefficients are labeled as , , , , , , , and , respectively, where represents the location of pixel inside each mask.

The magnitude for each filter can be obtained as follows: The proposed image denoising algorithm for grayscale images includes the following steps:(i)set the mask window size and the values of the fractional powers and ;(ii)apply fractional mask convolution on eight directions with the gray value of corresponding image pixels, adding all product terms to obtain weighting sum on eight directions;(iii)find the arithmetic mean of the weighting sum value on the eight directions as approximate value of fractional integral for image pixel;(iv)repeat steps (ii) and (iii) for all image pixels;(v)measure the PSNR for the result image.

The fractional mask convolution is performed by sliding the mask window over the image, generally starting at the top left corner of the image through all the pixels where the fractional fits entirely within the boundaries of the image. The size of the mask window and the values of the fractional powers and are chosen to achieve the requirements of image denoising. For testing we have added a Gaussian noise to the original image, and the noisy image is then used for image smoothing; therefore the PSNR can be calculated for the restored images.

#### 4. Experimental Results

This section aims at demonstrating that the proposed image denoising algorithm using fractional integral masks has better capability than the traditional approaches for image denoising.

The test image employed here is the grayscale images “Lena,” “Cameraman,” “Boat” and “Peppers” with 512 × 512 pixels. The default Gaussian noise is added into the image with different noise variances. All filters considered operate using 3 × 3 processing window mask. The values of the fractional powers are taken with the range and . The performance of filters was evaluated by computing the peak signal to noise ratio (PSNR) which has been wildly used in the literature [11, 24]. The value of PSNR depends totally on the size of the mask window and the values of the fractional powers and . PSNR is defined via the mean squared error (MSE) for two images and , where one of the images is considered the original noisy image and the other is the filtered image: where max is the maximum possible pixel value of the image. In grayscale image, this is equal to 255. Table 1 shows the results of PSNR values of our proposed method as compared with Gaussian smoothing filter with different Gaussian variance noise values. As it can be seen, the maximum PSNR value was obtained by our proposed approach. From the human visual system effect, Figures 2, 3, 4, and 5, illustrate that the proposed denoising algorithm using fractional integral masks has good denoising performance for both testing images by different degrees of noise. The proposed algorithm not only enhances the quality of filtered image but also reserves the textures and edges present in the image. Table 2 shows the comparison of experimental results of the proposed algorithm with other denoising algorithm with the variance of noise () of 10, 15, 20 and 25. The proposed algorithm for image denoising algorithm provides satisfactory results. The good PSNR of the proposed algorithm acts as one of the important parameters to judge its performance.

Table 1: Results of denoising by Gaussian smoothing filter and fractional integral filter with different noise variance values.
Table 2: Comparison of the experimental results with other standard methods.
Figure 2: Results of denoising by fractional integral and Gaussian smoothing filter with noise variance = 0.01.
Figure 3: Results of denoising by fractional integral and Gaussian smoothing filter with noise variance = 0.02.
Figure 4: Results of denoising by fractional integral and Gaussian smoothing filter with noise variance = 0.03.
Figure 5: Results of denoising by fractional integral and Gaussian smoothing filter with noise variance = 0.05.

#### 5. Conclusion

In this paper, a novel digital image denoising algorithm called generalized fractional integral filter based on generalized Srivastava-Owa fractional integral operator was used. Fractional mask convolution on eight directions with the gray value had been applied on eight directions. The results proved that the proposed algorithm not only enhances the quality of filtered image but also reserves the textures and edges present in the image. Changing the size of the mask window or any of the values of the fractional powers and will allow adjusting the fractional integral filter coefficients to each image according to it characteristics. Several experiments proved that the improvements achieved were comparable to Gaussian smoothing filter. Besides, the proposed algorithm exhibited better PSNR than the other two standard denoising methods.

#### Acknowledgment

This research has been funded by the University of Malaya, under the Grant no. RG066-11ICT.

#### References

1. K. S. Miller and B. Ross, An Introduction to the Fractional Integrals and Derivatives-Theory and Application, John Wiley & Sons, New York, NY, USA, 1993.
2. I. Podlubny, Fractional Differential Equations, vol. 198 of Mathematics in Science and Engineering, Academic Press, San Diego, Calif, USA, 1999.
3. R. Hilfer, Applications of Fractional Calculus in Physics, World Scientific, Singapore, 2000.
4. B. J. West, M. Bologna, and P. Grigolini, Physics of Fractal Operators, Institute for Nonlinear Science, Springer, New York, NY, USA, 2003.
5. A. A. Kilbas, H. M. Srivastava, and J. J. Trujillo, Theory and Applications of Fractional Differential Equations, vol. 204 of North-Holland Mathematics Studies, Elsevier Science, Amsterdam, The Netherlands, 2006.
6. J. Sabatier, O. P. Agrawal, and J. A. T. Machado, Advance in Fractional Calculus: Theoretical Developments and Applications in Physics and Engineering, Springer, Dordrecht, The Netherlands, 2007.
7. A. C. Sparavigna, “Using fractional differentiation in astronomy,” Computer Vision and Pattern Recognition (2010), http://arxiv.org/abs/0910.2381v3.
8. R. Marazzato and A. C. Sparavigna, “Astronomical image processing based on fractional calculus: the AstroFracTool,” Instrumentation and Methods for Astrophysics (2009), http://arxiv.org/abs/0910.4637v2.
9. C. C. Tseng, “Design of variable and adaptive fractional order fir differentiators,” Signal Processing, vol. 86, no. 10, pp. 2554–2566, 2006.
10. J. A. T. Machado, M. F. Silva, R. S. Barbosa et al., “Some applications of fractional calculus in engineering,” Mathematical Problems in Engineering, vol. 2010, Article ID 639801, 34 pages, 2010.
11. J. Hu, Y. Pu, and J. Zhou, “A novel image denoising algorithm based on riemann-liouville definition,” Journal of Computers, vol. 6, no. 7, pp. 1332–1338, 2011.
12. R. W. Ibrahim and M. Darus, “Subordination and superordination for analytic functions involving fractional integral operator,” Complex Variables and Elliptic Equations. An International Journal of Elliptic Equations and Complex Analysis, vol. 53, no. 11, pp. 1021–1031, 2008.
13. R. W. Ibrahim and M. Darus, “Subordination and superordination for univalent solutions for fractional differential equations,” Journal of Mathematical Analysis and Applications, vol. 345, no. 2, pp. 871–879, 2008.
14. S. M. Momani and R. W. Ibrahim, “On a fractional integral equation of periodic functions involving Weyl-Riesz operator in Banach algebras,” Journal of Mathematical Analysis and Applications, vol. 339, no. 2, pp. 1210–1219, 2008.
15. R. W. Ibrahim, “Solutions of fractional diffusion problems,” Electronic Journal of Differential Equations, vol. 2010, no. 147, pp. 1–11, 2010.
16. R. W. Ibrahim, “On the existence for diffeo-integral inclusion of Sobolev-type of fractional order with applications,” ANZIAM Journal, vol. 52, no. (E), pp. E1–E21, 2010.
17. D. Băleanu, O. G. Mustafa, and R. P. Agarwal, “On the solution set for a class of sequential fractional differential equations,” Journal of Physics A, vol. 43, no. 38, Article ID 385209, 2010.
18. R. W. Ibrahim and M. Darus, “On analytic functions associated with the Dziok-Srivastava linear operator and Srivastava-Owa fractional integral operator,” Arabian Journal for Science and Engineering, vol. 36, no. 3, pp. 441–450, 2011.
19. R. W. Ibrahim, “Existence and uniqueness of holomorphic solutions for fractional cauchy problem,” Journal of Mathematical Analysis and Applications, vol. 380, no. 1, pp. 232–240, 2011.
20. S. G. Samko, A. A. Kilbas, and O. I. Marichev, Fractional Integrals and Derivatives, Gordon and Breach Science, Yverdon, Switzerland, 1993.
21. J. T. Machado, V. Kiryakova, and F. Mainardi, “Recent history of fractional calculus,” Communications in Nonlinear Science and Numerical Simulation, vol. 16, no. 3, pp. 1140–1153, 2011.
22. H. M. Srivastava and S. Owa, Univalent Functions, Fractional Calculus, and Their Applications, Ellis Horwood Series: Mathematics and Its Applications, John Wiley & Sons, New York, NY, USA, 1989.
23. R. W. Ibrahim, “On generalized Srivastava-Owa fractional operators in the unit disk,” Advances in Difference Equations, vol. 2011, article no. 55, 2011.
24. E. Cuesta, M. Kirane, and S. Malik, “Image structure preserving denoising using generalized fractional time integrals,” Signal Processing, vol. 92, no. 2, pp. 553–563, 2012.