Advances in Mathematical Physics

Volume 2014 (2014), Article ID 494237, 6 pages

http://dx.doi.org/10.1155/2014/494237

Research Article

Adaptive Fractional Differentiation Harris Corner Detection Algorithm for Vision Measurement of Surface Roughness

Rui-Yin Tang1,2 and Zhou-Mo Zeng¹

¹State Key Laboratory of Precision Measuring Technology and Instruments, Tianjin University, 92 Weijin Road, Nankai District, Tianjin 300072, China

²College of Electrical Engineering, Hebei United University, 46 Xinhua Road, Tangshan 063009, China

Received 13 March 2014; Accepted 19 March 2014; Published 17 April 2014

Academic Editor: Xiao-Jun Yang

Copyright © 2014 Rui-Yin Tang and Zhou-Mo Zeng. 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.

The Harris algorithm via fractional order derivative (the adaptive fractional differentiation Harris corner detection algorithm), which adaptively adjusts the fractal dimension parameter, has been investigated for an analysis of image processing relevant to surface roughness by vision measurements. The comparative experiments indicate that the algorithm allows the edge information in the high frequency areas to be enhanced, thus overcoming shortcomings. The algorithm permits real-time measurements of surface roughness to be performed with high precision, superior to the conventional Harris algorithm.

1. Introduction

Fractional calculus theory [1–5] has increasingly focused on various applied areas. In this context, fractional-order viscoelasticity [6, 7], continuous-time finance [8], population dynamics [9], fractional-order feedback control [10], fractional-order fluid mechanics [11], the fractional-order supersymmetric theory in quantum mechanics [12], fractional conservation laws [13], fractional-order operators to electrocardiogram signals [14], fractional Lorenz system [15], subdiffusion phenomena [16], and granular mechanics [17] have been intensively studied.

In the context of surface analysis, fractals have been applied to model the behaviors of complex surface structures [18–23]. As a power mathematical tool the fractional calculus has been suggested to be adequate dealing with phenomena with fractals complexities (see [24–27] and the cited references therein).

The rapid development of computer technologies allows a variety of algorithms based on fractional calculus to 2D image processing to be developed [28–31]. Recently, analyzing surface roughness of a hard disk by visual measurements and assessment of the measurement precision (less 20 μm), the effects on the processing of images captured by a camera have been developed through improved Harris algorithm employing fractional calculus [30]. The aim of this paper is to propose the adaptive fractional differentiation Harris corner detection algorithm for vision measurement of surface roughness.

The paper is organized as follows. Section 2 presents the system of visual measurement for surface roughness. Section 3 gives the theoretical basis for image corner detection and edge detection, too. Section 4 develops an adaptive fractional differentiation Harris algorithm for image corner detection. Section 5 outlines the main results and the conclusions.

2. The System of Visual Measurement for Surface Roughness

The system for visual measurement of surface roughness consisting of a projector creating a thin line of light, camera, position system, and a hard disk is shown in Figure 1.

Figure 1: Structure of the system of vision measurement.

The procedure used to measure structured light is as follows.

(i)First, a thin luminous straight line is generated and projected onto the surface of the object that is to be measured using the laser line projector.(ii)Then, a stripe image modulated by the object height is formed by the camera. The camera coordinate system (see Figure 2) is defined by the world coordinate system () of the hard disk measurement system. The point is the center of the perspective projection, while is the optical axis of the camera. Additionally, is the image coordinate system, where is the point of intersection between the optical axis and the image plane. The distance between points and is the focal length.(iii)When a linear structured laser light is projected to the plane a certain point of the studied surface is projected through the center of the lens to the point in the image plane.(iv)The points such as depend on the height of the object. Thus, if the object onto which the laser line is projected differs in height, the line is not imaged as a straight line but represents a profile of the object allowing the differences in heights of the object to be determined. In order to gather many object profiles the object moves by means of a scanning system.

Figure 2: Structured light vision measurement principle.

The transformations from the pixel coordinates to the world coordinate system can be expressed mathematically as follows:

where

In order to calculate 3D world coordinates based on 2D image coordinates the light plane equation should be defined by the equation

The precision with which the 2D image coordinates are processed determines the precision of the 3D world coordinates that directly affects the precision in the surface roughness measurements. Therefore, precise extraction of the 2D coordinates in the image processing section (see Figure 1) is of primary importance.

In the hard disk image, the highlighted areas have to be measured which invokes both the corner detection and the edge detection algorithms to be developed and consequently their image coordinates to be determined.

3. Theoretical Analyses of the Image Corner Detection and Edge Detection

There exist a variety of corner and edge detection algorithms, such as the Moravec corner detection algorithm [28], the Harris corner detection algorithm [29], and the SUSAN corner detection algorithm [30]. In this context, Schmid [31] reported a comparative analysis of various corner detection algorithms and revealed that the Harris algorithm has the best detection performance. The Harris algorithm has been widely applied due to its simplicity in extraction of image corners. The detection results are not affected by factors such as image rotation and light intensity, for example. Moreover, the operator introduces a Gaussian smoothing template in processing calculations which makes the algorithm robust.

In accordance with the Harris corner detection algorithm a buttery hatch centered by pixel point moves with velocity along the direction and moves with velocity along the direction. The following analytical expression describes the gray level change:

where

Let us denote

where

is the Gaussian smoothing template, which improves the algorithm’s ability to resist noise. Usually, both and

result in the quadratic form of being

where

The partial autocorrelation function is very close.

The corner response function can be defined as

where is the value of the determinant of the matrix and denotes the matrix trace.

Let and be two eigenvalues of . If and are large and , then the local white correlation function has a peak which is a corner point. Otherwise, if and , then the point belongs to the image edge; if it is shown the point is in an area of unconspicuous gray change. Figure 3 shows results from processing results before and after application of the Harris corner detection algorithm.

(a) Original hard disk image

(b) Hard disk image processed by the Harris algorithm

Figure 3: Processing results before and after application of the Harris corner detection algorithm: (a) original hard disk image; (b) hard disk image processed by the Harris algorithm.

4. The Adaptive Fractional Differentiation Harris Corner Detection Algorithm

Developing the adaptive fractional differentiation Harris corner detection algorithm we have to introduce some definitions of fractional-order differentiation [30].

Let the signal exit the order derivative. For , the fractional differentiation of order is defined as [30]

Because the signal is digital, the first order signal is extended to 2D image signal . In this context, the fractional differentiation of signal is

Detecting the corner of high complexity texture image, it was detected that the Harris algorithm provides false corner points. Precisely, as a consequence of the first step of the Harris algorithm (see (5)) the first order differentiation can result in missing of details of many low-frequency region textures and consequently a lot of texture information cannot be retrieved. Moreover, many false corners simultaneously appear as it is shown in Figure 4. Therefore, the low-frequency extraction precision of corner points by the Harris algorithm needs to be improved. The improvement developed in this paper is based on adaptive fractional differentiation.

Figure 4: Fake corners in an image processed by the Harris algorithm.

As a first step in the improvement of the conventional Harris algorithm the Prewitt operator is replaced by fractional differentiation mask (the fractional differentiation mask in the and directions is given in Table 1).

Table 1: Fractional differentiation mask applied in the direction shown in (a) and direction shown in (b).

(a) Mask applied in the direction

0

0

0
1
0

0

0

(b) Mask applied in the direction

0
0
0

1

0
0
0

The fractal dimension represents the texture complexity and determines the order number of the fractional differentiation. As a consequent step a score box dimension algorithm [30] is used to calculate the order of fractional differentiation.

The experimental results reveal a nonuniform distribution of the gray image complexity that if the score box dimension is used [2,3], the fractal dimension of the complex texture image can be recovered [2.7,3]. In order to recover the fractional differentiation order at [0,1] taking into account the complexity of the image texture, the interval [2,2.7] is amplified and mapped to [0.3,1]; inasmuch the interval [2,2.7) is simultaneously zoomed out and mapped to [0,0.3).

Because differential operation details can be easily lost in images having rich textures there exists a negative correlation between the complexities and the image texture. The order number is defined as

The proposed results relevant to corner extraction developed by the adaptive fractional differentiation Harris algorithm are shown in Figure 5.

Figure 5: Results of corner extraction using the adaptive fractional differentiation Harris corner detection algorithm.

A comparison between analysis of the original Harris algorithm and the adaptive fractional differentiation Harris corner detection is illustrated by Figure 6. In general, details of many low-frequency region textures are missing if the original Harris algorithm is applied (see Figure 6(a)); many texture details remain (indicated by the yellow ellipse) invisible. The results of the adaptive fractional differentiation Harris corner detection algorithm are shown in Figure 6(b). The image can be subsequently processed by application of common methods, such as the morphology and gray centroid methods; the results of light stripe center extraction are shown in Figure 7. The light stripe center region marked by A in Figure 7(a) is obviously inaccurate, whereas the light stripe center marked B in Figure 7(b) is more accurate than Figure 7(a).

(a) Result of the original Harris algorithm

(b) Result of the improved Harris algorithm

Figure 6: Comparison between the results of the original Harris algorithm and the improved Harris algorithm.

(a) Light stripe center extraction result based on Figure 6(a)

(b) Light stripe center extraction result based on Figure 6(b)

Figure 7: Light stripe center extraction comparison.

The results reported in the previous section reveal that the fractional differentiation has advantages with respect to the first order differential used in the Prewitt operator. It is well-known that the differential operation promotes high-frequency signals and weakens the low-frequency signals. As the order number increases, the more intensely high-frequency signals are enhanced and low-frequency signals are weakened. The order number of fractional differentiation is lower than that of integer order counterpart and the fractional differentiation is strengthening the high-frequency signals and the low-frequency parts of the signals remain as a nonlinearity. Comparing to the integer order differentiation, the fractional differentiation allows retaining more texture information in the low-frequency areas, thus maintaining the target contour of the image. Moreover, the fractional differentiation can enhance the edge information of high-frequency areas.

The data in Table 2 present the heights measured when the original Harris algorithm (V1) and the adaptive fractional differentiation Harris corner detection algorithm (V2) are used to process the stripe image. These data indicate that the highest measurement precision error of V2 is about 0.023 mm. It sufficiently satisfies the measurement requirements of the hard disk plane and is less than that of V1, too.

Table 2: Comparison of the measurement precision.

Measurement point
Measured height V1 (mm)
Measured height V2 (mm)
Standard height (mm)
V2 error (mm)

P1
23.281
23.398
23.412
−0.014

P2
4.368
4.224
4.201
0.023

P3
23.350
23.397
23.420
−0.023

P4
21.782
21.701
21.611
0.090

P5
9.435
9.403
9.418
−0.015

P6
16.673
16.518
16.541
−0.023

P7
8.425
8.497
8.507
−0.010

P8
6.247
6.334
6.313
0.021

P9
6.269
6.345
6.327
0.018

P10
23.359
23.429
23.446
−0.017

5. Conclusion

Large number of false corner points generated by the conventional Harris algorithm when image textures are complex have been avoided by application of the adaptive fractional differentiation. The Harris corner detection algorithm with implemented adaptive fractional differentiation is for the first time suggested in this paper as a principle contribution improving the data processing. The adaptive fractional differentiation allows the improved Harris corner detection to retain more texture information in low-frequency areas, thus maintaining the target contour of images and enhances the edge information of high-frequency areas in a manner better than the conventional Harris algorithm. The proposed algorithm enhances the Harris algorithm for visual measurements relevant to surface roughness of hard disks.

Conflict of Interests

The authors declare that they have no competing interests in this paper.

Acknowledgment

This work was supported by National Natural Science Foundation of China (no. 51105273).

References

J. T. Machado, A. M. Galhano, and J. J. Trujillo, “On development of fractional calculus during the last fifty years,” Scientometrics, vol. 98, no. 1, pp. 577–582, 2014.

A. Kilbas, H. M. Srivastava, and J. J. Trujillo, Theory and Applications of Fractional Differential Equations, Elsevier, Amsterdam, The Netherlands, 2006.

J. A. Tenreiro Machado, A. C. J. Luo, and D. Baleanu, Nonlinear Dynamics of Complex Systems: Applications in Physical, Biological and Financial Systems, Springer, New York, NY, USA, 2011.

V. E. Tarasov, Fractional Dynamics: Applications of Fractional Calculus to Dynamics of Particles, Fields and Media, Springer, Berlin, Germany, 2011.

F. Mainardi, Fractional Calculus and Waves in Linear Viscoelasticity, Imperial College Press, London, UK, 2010.

R. C. Koeller, “Applications of fractional calculus to the theory of viscoelasticity,” Journal of Applied Mechanics, vol. 51, no. 2, pp. 299–307, 1984.

F. Mainardi, “An historical perspective on fractional calculus in linear viscoelasticity,” Fractional Calculus and Applied Analysis, vol. 15, no. 4, pp. 712–717, 2012.

F. Mainardi, M. Raberto, R. Gorenflo, and E. Scalas, “Fractional calculus and continuous-time finance. II: the waiting-time distribution,” Physica A: Statistical Mechanics and Its Applications, vol. 287, no. 3-4, pp. 468–481, 2000.

A. M. A. El-Sayed, S. Z. Rida, and A. A. M. Arafa, “Exact solutions of fractional-order biological population model,” Communications in Theoretical Physics, vol. 52, no. 6, pp. 992–996, 2009.

I. Petráš, “Fractional—order feedback control of a DC motor,” Journal of Electrical Engineering, vol. 60, no. 3, pp. 117–128, 2009.

V. V. Kulish and J. L. Lage, “Application of fractional calculus to fluid mechanics,” Journal of Fluids Engineering, vol. 124, no. 3, pp. 803–806, 2002.

D. Baleanu and S. I. Muslih, “About fractional supersymmetric quantum mechanics,” Czechoslovak Journal of Physics, vol. 55, no. 9, pp. 1063–1066, 2005.

G. S. F. Frederico and D. F. M. Torres, “Fractional conservation laws in optimal control theory,” Nonlinear Dynamics, vol. 53, no. 3, pp. 215–222, 2008.

M. Benmalek and A. Charef, “Digital fractional order operators for R-wave detection in electrocardiogram signal,” IET Signal Processing, vol. 3, no. 5, pp. 381–391, 2009.

I. Grigorenko and E. Grigorenko, “Chaotic dynamics of the fractional Lorenz system,” Physical Review Letters, vol. 91, article 034101, 2003.

J. Hristov, “Approximate solutions to fractional subdiffusion equations,” European Physical Journal: Special Topics, vol. 193, no. 1, pp. 229–243, 2011.

R. A. Pfaffenzeller, M. K. Lenzi, and E. K. Lenzi, “Modeling of granular material mixing using fractional calculus,” International Review of Chemical Engineering, vol. 3, no. 6, pp. 814–817, 2011.

B. B. Mandelbrot, The Fractal Geometry of Nature, Macmillan, London, UK, 1983.

T. Vicsek, Fractal Growth Phenomena, World scientific, Singapore, 1989.

E. Goldfain, “Fractional dynamics, cantorian space-time and the gauge hierarchy problem,” Chaos, Solitons and Fractals, vol. 22, no. 3, pp. 513–520, 2004.

B. J. West, M. Bologna, and P. Grigolini, Physics of Fractal Operators, Springer, New York, NY, USA, 2003.

X.-J. Yang, Local Fractional Functional Analysis and Its Applications, Asian Academic, Hong Kong, China, 2011.

X.-J. Yang, Advanced Local Fractional Calculus and Its Applications, World Science, New York, NY, USA, 2012.

C. F. Liu, S. S. Kong, and S. J. Yuan, “Reconstructive schemes for variational iteration method within Yang-Laplace transform with application to fractal heat conduction problem,” Thermal Science, vol. 17, no. 3, pp. 715–721, 2013.

A. K. Golmankhaneh, A. K. Golmankhaneh, and D. Baleanu, “Lagrangian and Hamiltonian mechanics on fractals subset of real-line,” International Journal of Theoretical Physics, vol. 52, no. 11, pp. 4210–4217, 2013.

X. -J. Yang and D. Baleanu, “Fractal heat conduction problem solved by local fractional variation iteration method,” Thermal Science, vol. 17, no. 2, pp. 625–628, 2013.

A. Carpinteri, B. Chiaia, and P. Cornetti, “The elastic problem for fractal media: basic theory and finite element formulation,” Computers and Structures, vol. 82, no. 6, pp. 499–508, 2004.

P. M. Hans, “Towards automatic visual obstacle avoidance,” in Proceedings of International Joint Conference on Artificial Intelligence, pp. 584–590, Cambridge, MA, USA, 1977.

C. Harris and M. Stephens, “A combined corner and edge detector,” in Proceedings of the 4th Alvey Vision Conference, pp. 147–151, Manchester, UK, 1998.

S. M. Smith and J. M. Brady, “SUSAN—a new approach to low level image processing,” International Journal of Computer Vision, vol. 23, no. 1, pp. 45–78, 1997.

N. Sarkar and B. B. Chauduri, “Efficient differential box-counting approach to compute fractal dimension of image,” IEEE Transactions on Systems, Man and Cybernetics, vol. 24, no. 1, pp. 115–120, 1994.