MPE Mathematical Problems in Engineering 1563-5147 1024-123X Hindawi Publishing Corporation 149289 10.1155/2013/149289 149289 Research Article Texture Enhancement Based on the Savitzky-Golay Fractional Differential Operator Jalab Hamid A. 1 Ibrahim Rabha W. 2 Tenreiro Machado José 1 Faculty of Computer Science & Information Technology University of Malaya 50603 Kuala Lumpur Malaysia um.edu.my 2 Institute of Mathematical Sciences University of Malaya 50603 Kuala Lumpur Malaysia um.edu.my 2013 24 2 2013 2013 13 08 2012 05 01 2013 05 01 2013 2013 Copyright © 2013 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.

Texture enhancement for digital images is the most important technique in image processing. The purpose of this paper is to design a texture enhancement technique using fractional order Savitzky-Golay differentiator, which leads to generalizing the Savitzky-Golay filter in the sense of the Srivastava-Owa fractional operators. By employing this generalized fractional filter, texture enhancement is introduced. Consequently, it calculates the generalized fractional order derivative of the given image using the sliding weight window over the image. Experimental results show that the operator can extract more subtle information and make the edges more prominent. In general, the capability of the generalized fractional differential will be high because it is sensitive to the subtle fluctuations of values of pixels.

1. Introduction

Texture is an important feature of natural images; hence, a variety of image texture applications has been intensively studied by many researchers . Image texture is defined as a function of the spatial variation in pixel intensities (gray values). Smith and Chang  have defined texture as visual patterns, which have properties of homogeneity and not resulting from only a single color or intensity.

In the image, texture features capture information about repeating patterns. Texture analysis can be classified into three models: structural, statistical, and signal theoretic methods . Therefore, the analysis of texture parameters is a useful approach for increasing the information accessible from images. In texture enhancement technique, which is based on mask operation, each pixel is modified according to value of the neighbourhood around the pixel of interest. One important aspect of an image, which enables us to perform this, is the notion of frequencies. Fundamentally, the frequencies of an image are the amount, by which the gray values change with distance. High-frequency components are characterized by huge changes in gray values over small distances; examples of high frequency components are edges and noise. On the other hand, low-frequency components are parts of the image, which are characterized by little change in the gray values . Fractional differential mask can further preserve the low-frequency contour feature in those smooth areas, and nonlinearly keep the high-frequency marginal feature in those areas, where the gray-level changes heavily, and also enhances texture details in those areas, where the gray-level does not change evidently.

Fractional integration and fractional differentiation are generalizations of notions of integer-order integration and differentiation and include nth derivatives and n-fold integrals as particular cases [5, 6]. Many applications of fractional calculus in physics have replaced 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. Currently, fractional calculus (integral and differential operators) is heavily used in control design [12, 13], Furthermore, in image processing , all results that are based on the fractional calculus operators (differential and integral) show that this method not only is effective, but also has good immunity.

The digital fractional order differentiator is an important topic in fractional calculus that can estimate the fractional order derivative of any given digital signals, without known function. The Savitzky-Golay filter is a simplified digital differentiator that is implemented by a local polynomial regression technique [20, 21]. Now, Savitzky-Golay digital differentiator has been one of the most popular numerical differentiation methods, due to its high computing speed and strong antinoise ability.

Recently, the interest in using texture enhancement technique based on mask operation has grown in the field of image processing. Pu and Zhou  have implemented multiscale texture segmentation by fractional differential. They have proposed two fractional differential masks and presented the structures and parameters of each mask, respectively. Then they have discussed the multi-scale texture segmentation based on the fractional mask. Pu  has proposed a fractional calculus approach to enhance the texture of digital image. He has found that the textural detail enhancing capability of fractional derivative-based texture operator is much better than integer derivative. Zhang et al.  have used the fractional differential masks based on the classical Riemann-Liouville definition. They have concluded that the fractional order between 1 and 2 can enhance the texture and edges in multiscale, by controlling the fractional order.

In this paper, we have used a generalized fractional differential based on the generalized Savitzky-Golay filter in sense of Srivastava-Owa fractional operators for image texture enhancement. The Savitzky-Golay filter has become a powerful signal and image processing tool, which has found application in many scientific areas. Moreover, the Savitzky-Golay filter method is considered to be a good approach in image texture enhancement, which is used as an alternative to classical techniques. The rest of the paper is organized as follows: Sections 2 and 3 explain the fractional calculus and the generalized fractional integral operator, respectively, Section 4 describes the construction of fractional differential Savitzky-Golay filter, Section 5 elucidates the experimental results, and Section 6 concludes the paper.

2. Fractional Calculus

The idea of the fractional calculus (i.e., calculus of integrals and derivatives of any arbitrary real or complex order) was found over 300 years ago. Abel in 1823 scrutinized the generalized tautochrone problem and for the first time applied fractional calculus techniques in a physical problem.

2.1. The Riemann-Liouville Operators

The Riemann-Liouville fractional derivative strongly poses the physical interpretation of the initial conditions required for the initial value problems involving fractional differential equations. Moreover, this operator possesses advantages of fast convergence, high stability, and higher accuracy to derive different types of numerical algorithms .

The fractional (arbitrary) order integral of the function f of order α>0 is defined by (1)Iaαf(t)=at(t-τ)α-1Γ(α)f(τ)dτ. When a=0, we write Iaαf(t)=f(t)*ϕα(t), where (*) denoted the convolution product, ϕα(t)=tα-1/Γ(α),   t>0 and ϕα(t)=0,  t0, and ϕαδ(t) as α0, where δ(t) is the delta function.

The fractional (arbitrary) order derivative of the function f of order 0α<1 is defined by (2)Daαf(t)=ddtat(t-τ)-αΓ(1-α)f(τ)dτ=ddtIa1-αf(t). When a=0, we have (3)Dαtμ=Γ(μ+1)Γ(μ-α+1)tμ-α,μ>-1,0<α<1,Iαtμ=Γ(μ+1)Γ(μ+α+1)tμ+α,μ>-1,α>0. The Caputo fractional derivative of order μ>0 is defined, for a smooth function f, by (4)  cDaμf(t)=at(t-τ)n-μ-1Γ(n-μ)f(n)(τ)dτ, where n=[μ]+1 (the notation [μ] stands for the largest integer not greater than μ). Note that there is a relationship between the Riemann-Liouville differential operator and the Caputo operator: (5)Daμf(t)=cDaμf(t)+1Γ(1-μ)f(a)(t-a)μ, and they are equivalent in a physical problem (i.e., a problem which specifies the initial conditions).

2.2. The Srivastava-Owa Operators

In , Srivastava and Owa defined and studied fractional operators (derivative and integral) in the complex z-plane for analytic functions.

The fractional derivative of order β is defined, for a function f(z) by (6)Dzβf(z):=1Γ(1-β)ddz0zf(ζ)(z-ζ)βdζ,0β<1, where the function f(z) is analytic in simply-connected region of the complex z-plane containing the origin and the multiplicity of (z-ζ)-β is removed by requiring log(z-ζ) to be real when (z-ζ)>0. Furthermore, for nβ<n+1, the fractional differential operator is defined as (7)Dzβf(z)=dndznDzβ-nf(z),n.

The fractional integral of order β is defined, for a function f(z), by (8)Izβf(z):=1Γ(β)0zf(ζ)(z-ζ)β-1dζ,β>0, where the function f(z) is analytic in simply connected region of the complex z-plane () containing the origin and the multiplicity of (z-ζ)β-1 is removed by requiring log(z-ζ) to be real when (z-ζ)>0: (9)Dzβzμ=Γ(μ+1)Γ(μ-β+1)zμ-β,μ>-1,0β<1Izβzμ=Γ(μ+1)Γ(μ+β+1)zμ+β,μ>-1,β>0. Note that the real case of the Srivastava-Owa operators is equivalence to the Riemann-Liouville operators.

3. Generalized Fractional Integral Operator

This section briefly describes the mathematical background for the fractional integral operator that has been used by the proposed algorithm. The usual way of representing the fractional derivatives is by the Riemann-Liouville formula Dtα. Another way to represent the fractional derivatives is by the Grünwald-Letnikov formula . The discrete approximations derived from the Grünwald-Letnikov fractional derivatives present some limitations, such as the following:

they frequently originate unstable numerical methods;

the order of accuracy of such approaches is never higher than one.

To implement the generalized fractional integral method, Ibrahim in  has imposed a formula for the generalized fractional integral. Consider, for natural n={1,2,} and real μ, the n-fold integral of the form (10)Izα,μf(z)=0zζ1μdζ10ζ1ζ2μdζ20ζn-1ζnμf(ζn)dζn. Applying the Cauchy formula for iterated integrals implies (11)0zζ1μdζ10ζ1ζμf(ζ)dζ=0zζμf(ζ)dζζzζ1μdζ1=1μ+10z(zμ+1-ζμ+1)ζμf(ζ)dζ. Repeating the above step n-1 times we obtain (12)0zζ1μdζ10ζ1ζ2μdζ20ζn-1ζnμf(ζn)dζn=(μ+1)1-n(n-1)!0z(zμ+1-ζμ+1)n-1ζμf(ζ)dζ, which yields the fractional operator type (13)Izα,μf(z)=(μ+1)1-αΓ(α)0z(zμ+1-ζμ+1)α-1ζμf(ζ)dζ, where α and μ-1 are real numbers and the function f(z) is analytic in the simply connected region of the complex z-plane containing the origin, and the multiplicity of (zμ+1-ζμ+1)-α is removed by requiring log(zμ+1-ζμ+1) to be real when (zμ+1-ζμ+1)>0. When μ=0, 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 (13), we define the generalized differential operator of order α by (14)Dzα,μf(z):=(μ+1)αΓ(1-α)ddz0zζμf(ζ)(zμ+1-ζμ+1)αdζ,0<α1, where the function f(z) is analytic in the simply connected region of the complex z-plane containing the origin and the multiplicity of (zμ+1-ζμ+1)-α is removed by requiring log(zμ+1-ζμ+1) to be real when (zμ+1-ζμ+1)>0.

Proposition 1 (see [<xref ref-type="bibr" rid="B27">27</xref>]).

The generalized derivative of the function f(z)=zν,  ν is given by the following: (15)Dzα,μf(z)=(μ+1)α-1Γ(ν/(μ+1)+1)Γ(ν/(μ+1)+1-α)z(1-α)(μ+1)+ν-1, which is later used to compute the coefficient matrix W(α,μ).

4. Construction of the Fractional Differential Savitzky-Golay Filter

The Savitzky-Golay filter has been introduced for computing the numerical derivatives and is also called a digital smoothing polynomial filter. The Savitzky-Golay method is often used to preserve higher moments in the data, thus reducing the distortion of essential features of the data.

In this section, we willgeneralize this filter for calculating the fraction derivatives which will be utilized by the proposed algorithm.

Assume a uniformly sampled signal, our aim is to estimate its nth order derivative using I-point filtering window and an n-degree polynomial : (16)fn(i)=k=0nakik, which is used to fit the given signal i=1,2,,I. In matrix notation, (16) is reduced to the system (17)Y=XA+ε, where ε is the estimate error, A is the n+1×1 coefficient matrix and X is the I×(n+1) Vandermonde matrix defined by (18)X=(1111n1212n1I1In). The coefficients of the best-fit polynomial can be obtained by minimizing the sum of the squared errors between the actual data and fitting points. Thus, (19)B=(XTX)-1XTY implies (20)Y^=WY, where W denotes the moving window’s coefficients matrix. Consequently, the nth order derivative can be estimate by (21)Y^(dn)=W(dn)Y.

Now, in view of Proposition 1, we have (22)Y^(α,μ)=W(α,μ)Y=A(XTX)-1XTY, where (23)A=[(μ+1)α-1Γ(1-α)i(1-α)(μ+1)-1,(μ+1)α-1Γ(1/(μ+1)+1)Γ(1/(μ+1)+1-α)i(1-α)(μ+1),,(μ+1)α-1Γ(n/(μ+1)+1)Γ(n/(μ+1)+1-α)i(1-α)(μ+1)+n-1].

Note that when µ=0, we have the Riemann-Liouville differential operator. Moreover, when I=n+1, the Vandermonde matrix X is a square matrix. The purpose of the Savitzky-Golay filter is to estimate W(α,μ), which can be used to calculate the nth order derivative of any given signal . The coefficient matrix W(α,μ) can be computed by (24)W(α,μ)=AX-1, where (25)X=(1111n1(1+ϕ)1(1+ϕ)n1(1+2ϕ)1(1+2ϕ)n1I1In), where ϕ is the sampling interval.

The matrix X can be assumed as the formula of computation time; therefore, the generalized Savitzky-Golay filter can be viewed as the generalization of the differentiator. However, digital image is a function of two variables, so, we can generalize these definitions to include both the x and y values.

The mask is designed into r×r size matrix M which has r layers (r is odd). The window’s size can be an arbitrary odd number, and a larger window can improve the accuracy of fractional differential, but increases the computational time. Therefore, we proceed to use moving window W(α,μ) with 3×3 size as shown in Table 1.

Moving window of fractional differential weights with 3 × 3 dimension.

 1 11 12 1 ( 1 + ϕ ) 1 ( 1 + ϕ ) 2 1 I 1 I 2

The fractional differential operator can enhance edges and contours as well as reserve the texture details. The nine values output of each fractional differential window h(x,y) is performed by sliding the mask window w(s,t) over the image f. Generally one can start at the top left corner of the image block through all the pixels, where the fractional differential mask fits entirely within the boundaries of the image. The output of each image block is nine values, which represent the texture information in each image block, that takes the following formula: (26)h(x,y)=w(s,t)f(x+s,y+t), where f is the value of an image pixel and w is the value of filter mask.

5. Experimental Results and Discussion

The reason of this experiment is to validate the correctness of the proposed algorithm.

Performance tests for the algorithm proposed by this paper were implemented using Matlab 2010a on Intel (R) Core i7 at 2.2 GHz, 4 GB DDR3 Memory, and system type 64-bit, Window 7. The computation time per image differs for each image and depends mainly on the window’s size as well as image size.

The proposed texture features enhancement algorithm includes the following steps:

set the value of n and I (≥1);

set the value of the image sampling interval ϕ;

set the values of the fractional power parameters (α,µ(0,1));

compute Savitzky-Golay moving window W(α,μ) as in  (24);

compute the Vandermonde matrix X as in (25);

apply the Savitzky-Golay fractional differential mask with the corresponding image pixels by sliding the window over the image.

By varying both the fractional powers α and µ, keeping I, and n fixed (I=9, n=2), the elements of the Savitzky-Golay moving window W(α,μ) have been computed as shown in Tables 2 and 3.

Moving window of fractional differential weights with 3 × 3 dimension for φ=1, α=0.4, and μ=0.5.

 0.2268 0.3057 −0.0060 −0.8687 1.2983 0.0420 −1.0277 1.2742 0.2016

Moving window of fractional differential weights with 3 × 3 dimension for φ=1, α=0.5, and μ=0.5.

 0 . 0732 0.3951 −0.0076 −0.8192 1.126 0.0416 −0.7999 0.9307 0.1772

Tables 2 and 3. show the coefficients of the fractional differential moving window for different values of α. All coefficient values are not equal to zero, which implies that the magnitude response of Savitzky-Golay filter is not also zero in the image region. This will likelyimprove the texture detail. However, the qualities of texture is defined by the spatial distribution of gray values for this reason, we have used gray-scale images for testing, which are shown in Figures 1(a) and 1(b).

Original images.

In order to illustrate the efficiency of the proposed algorithm in Figure 2, we have presented an illustration of the obtained results for the texture enhancement of the original images of Figure 1.

The results of the proposed algorithm for µ=0.5 and; (a) α=0.4, (b) α=0.5.

The proposed enhancement algorithm shows good enhancement performance for both, testing images by different degrees of fractional power values α and µ which are experimentally fixed at α=0.4, 0.5 and µ=0.5, and the value of the image sampling interval φ=1. It is seen that, the proposed enhancement algorithm using fractional differential masks, can extract more texture information and sharpen edges more efficiently. The eye’s qualitative analysis of the proposed algorithm acts as one of the important parameters to judge its performance.

Other metrics used to judge the algorithm performance are the statistical measures. In this paper, among the statistical features, the following second-order statistics are used as texture features in representing images. The gray-level co-occurrence matrix (GLCM) is a statistical method used to describe textures in an image, by modeling texture as a two-dimensional gray level variation . Four statistical measures are extracted to evaluate the images texture enhancement; these are entropy, homogeneity, contrast, and energy.

Entropy measures the amount of information, and the larger value of entropy is the greater amount of information carried by image, but inversely correlated to energy. Entropy feature of gray-scale cooccurrence matrix is one of the features having the best discriminatory power, which is given in following equation: (27)Entropy=i,jp(i,j)logp(i,j),

where p(i,j) is the probability for gray-scale i and j and occurs at two pixels.

Homogeneity measures the closeness of the distribution of elements in the GLCM to the GLCM diagonal: (28)Homogeneity=i,jp(i,j)1+|i-j|.

Contrast measures the intensity contrast between a pixel and its neighbour over the whole image: (29)Contrast=i,jp(i-j)2.

Energy measures the sum of squared elements in the gray-level cooccurrence matrix (GLCM): (30)Energy=i,jp(i,j)2.

Figure 2 shows the results of the proposed enhancement algorithm for µ=0.5 with (a) α=0.4, and (b) α=0.5. The variation of the image texture is observed when α is increased from 0.4 to 0.5. So, the selection of differential order is important.

Tables 4 and 5 and Figures 3 and 4 show the performance evaluation of the proposed algorithm for image (a) and image (b) according to those four statistical measures of gray-level cooccurrence matrix (GLCM). It can be clearly seen that there has been a large increase in the value of entropy, which means the greater amount of information is carried by image due to texture enhancement. The entropy values for image (a) climbed to approximately 0.224 for α=0.4 and to 0.26 for α=0.5 for all texture enhancement cases. Moreover, for image (a), it inversely correlated to energy, which decreases to approximately 0.881 for α=0.4 and to 0.860 for α=0.5 and from its value of the original testing image. The homogeneity steadiness is reduced with the increase of texture enhancement, which means more divergence of the distribution of elements of information carried by image due to texture enhancement process. While the contrast showed diverse tendencies for all texture enhancement cases, it is conclude that the intensity contrast between a pixel and its neighbour over the whole image are changed too. This variation in the statistical measures makes the proposed algorithm capable to control the degree of texture enhancement of the image by controlling the fractional order parameters α,µ and φ.

Performance evaluation of the proposed algorithm compared to four statistical measures for image (a).

Images Entropy Homogeneity Contrast Energy
Original image (a) 0.036 0.999 0.025 0.998
Proposed algorithm with (µ=0.5, α=0.4) 0.224 0.949 2.767 0.881
Proposed algorithm with (µ=0.5, α=0.5) 0.260 0.941 3.219 0.860

Performance evaluation of the proposed algorithm compared to four statistical measures for image (b).

Images Entropy Homogeneity Contrast Energy
Original image (b) 0 1 0 1
Proposed algorithm with (µ=0.5, α=0.4) 0.017 0.997 0.112 0.994
Proposed algorithm with (µ=0.5, α=0.5) 0.032 0.995 0.214 0.989

The evaluation performance of the proposed algorithm compared to four statistical measures for image (a).

The evaluation performance of the proposed algorithm compared to four statistical measures for image (b).

6. Conclusion

In this paper, a texture enhancement technique using fractional order Savitzky-Golay differentiator, which leads to generalize Savitzky-Golay filter in sense of Srivastava-Owa fractional operators, have been introduced. The new algorithm presented in this paper can control the degree of texture enhancement of the image with the fractional power values. The new approach can control the degree of texture enhancement of the image with fractional order of the parameters α,µ, and φ. However, the technique is by no means limited only to images, instead, it can be applied in the setting of different image applications, taking into consideration the limitations of each imaging method. Furthermore, our goal is to keep away from the effect of the noise that caused in the texture enhancement of the image and to control the degree of texture enhancement of the image with the filter mask parameters. The experiment results had demonstrated the efficacy of this algorithm according to the metrics used to judge the algorithm performance.

Acknowledgments

The authors would like to thank the reviewers for their comments on earlier versions of this paper. This research has been funded by university of Malaya, under Grant no. UMRG 104-12ICT.

Jalab H. A. Image retrieval system based on color layout descriptor and Gabor filters Proceedings of the IEEE Conference on Open Systems (ICOS '11) 2011 32 36 Smith J. R. Chang S. F. Automated binary texture feature sets for image retrieval Proceedings of the IEEE International Conference on Acoustics, Speech, and Signal Processing (ICASSP '96) May 1996 2239 2242 2-s2.0-0029725713 Bouguila N. Hammoud R. I. Color texture classification by a discrete statistical model and feature selection Proceedings of the IEEE International Conference on Image Processing (ICIP '09) November 2009 1381 1384 2-s2.0-77951970553 10.1109/ICIP.2009.5413376 McAndrew A. An introduction to digital image processing with matlab notes for SCM2511 image processing School of Computer Science and Mathematics 2004 Victoria University of Technology 1 264 Samko S. G. Kilbas A. A. Marichev O. I. Fractional Integrals and Derivatives 1993 Yverdon, Switzerland Gordon and Breach Science Publishers xxxvi+976 MR1347689 Miller K. S. Ross B. An Introduction to the Fractional Calculus and Fractional Differential Equations 1993 New York, NY, USA John Wiley & Sons xvi+366 A Wiley-Interscience Publication MR1219954 Oustaloup A. Levron F. Mathieu B. Nanot F. M. Frequency-band complex noninteger differentiator: characterization and synthesis IEEE Transactions on Circuits and Systems I 2000 47 1 25 39 2-s2.0-0033888365 10.1109/81.817385 Malti R. Aoun M. Levron F. Oustaloup A. Analytical computation of the H2-norm of fractional commensurate transfer functions Automatica 2011 47 11 2425 2432 10.1016/j.automatica.2011.08.021 MR2886871 Malti R. Moreau X. Khemane F. Oustaloup A. Stability and resonance conditions of elementary fractional transfer functions Automatica 2011 47 11 2462 2467 10.1016/j.automatica.2011.08.029 MR2886877 ZBL1228.93111 Victor S. Melchior P. Oustaloup A. Robust path tracking using flatness for fractional linear MIMO systems: a thermal application Computers & Mathematics with Applications 2010 59 5 1667 1678 10.1016/j.camwa.2009.08.008 MR2595939 ZBL1189.93049 Sommacal L. Melchior P. Dossat A. Petit J. Cabelguen J. M. Oustaloup A. Ijspeert A. J. Improvement of the muscle fractional multimodel for low-rate stimulation Biomedical Signal Processing and Control 2007 2 3 226 233 2-s2.0-35348866111 10.1016/j.bspc.2007.07.013 Gruel D. N. Lanusse P. Oustaloup A. Robust control design for multivariable plants with time-delays Chemical Engineering Journal 2009 146 3 414 427 2-s2.0-58349094549 10.1016/j.cej.2008.09.041 Lanusse P. Oustaloup A. Sabatier J. Robust design of an anti-windup compensated 3rd-generation CRONE controller Advances in Fractional Calculus 2007 Dordrecht, The Netherlands Springer 527 542 10.1007/978-1-4020-6042-7_36 MR2432161 ZBL1125.93418 Hu J. Pu Y. Zhou J. A novel image denoising algorithm based on riemann-liouville definition Journal of Computers 2011 6 7 1332 1338 2-s2.0-79960730033 10.4304/jcp.6.7.1332-1338 Jalab H. A. Ibrahim R. W. Denoising algorithm based on generalized fractional integral operator with two parameters Discrete Dynamics in Nature and Society 2012 2012 14 529849 MR2923390 ZBL1244.94007 Jalab H. Ibrahim R. Texture feature extraction based on fractional mask convolution with cesáro means for content-based image retrieval Proceedings of the 12th Pacific Rim Conference on Artificial Intelligence 2012 170 179 Liu Y. Pu Y. Zhou J. Design of image denoising filter based on fractional integral Journal of Computational Information Systems 2010 6 9 2839 2847 2-s2.0-77958106219 Hu J. Pu Y. Zhou J. Fractional integral denoising algorithm and implementation of fractional integral filter Journal of Computational Information Systems 2011 7 3 729 736 2-s2.0-79953748688 Mathieu B. Melchior P. Oustaloup A. Ceyral C. Fractional differentiation for edge detection Signal Processing 2003 83 11 2421 2432 2-s2.0-0141927230 10.1016/S0165-1684(03)00194-4 ZBL1145.94309 Schafer R. W. What is a savitzky-golay filter? IEEE Signal Processing Magazine 2011 28 4 111 117 2-s2.0-79959394502 10.1109/MSP.2011.941097 Chen D. Xue D. Pan F. Digital fractional order Savitzky-Golay differentiator 58 IEEE Transactions on Circuits and Systems II 2011 758 762 Pu Y.-F. Zhou J.-L. A novel approach for multi-scale texture segmentation based on fractional differential International Journal of Computer Mathematics 2011 88 1 58 78 10.1080/00207160903134263 MR2754676 ZBL1209.94013 Pu Y. Fractional calculus approach to texture of digital image Proceedings of the 8th International Conference on Signal Processing (ICSP '06) November 2006 2-s2.0-34249335716 10.1109/ICOSP.2006.345713 Zhang Y. Pu Y. Zhou J. Construction of fractional differential masks based on Riemann-Liouville definition Journal of Computational Information Systems 2010 6 10 3191 3200 2-s2.0-78149381405 Srivastava H. M. Owa S. Univalent Functions, Fractional Calculus, and Their Applications 1989 Chichester, UK Ellis Horwood 404 MR1199135 Sousa E. How to approximate the fractional derivative of order 1<α2 International Journal of Bifurcation and Chaos in Applied Sciences and Engineering 2012 22 4 13 10.1142/S0218127412500757 MR2926051 Ibrahim R. W. On generalized Srivastava-Owa fractional operators in the unit disk Advances in Difference Equations 2011 2011 1 10 MR2855254 Selvarajah S. Kodituwakku S. Analysis and comparison of texture features for content based image retrieval International Journal of Latest Trends in Computing 2011 108