- About this Journal ·
- Abstracting and Indexing ·
- Aims and Scope ·
- Annual Issues ·
- Article Processing Charges ·
- Articles in Press ·
- Author Guidelines ·
- Bibliographic Information ·
- Citations to this Journal ·
- Contact Information ·
- Editorial Board ·
- Editorial Workflow ·
- Free eTOC Alerts ·
- Publication Ethics ·
- Reviewers Acknowledgment ·
- Submit a Manuscript ·
- Subscription Information ·
- Table of Contents

International Journal of Distributed Sensor Networks

Volume 2013 (2013), Article ID 203719, 6 pages

http://dx.doi.org/10.1155/2013/203719

## Circle Fitting Using a Virtual Source Localization Algorithm in Wireless Sensor Networks

^{1}Xi'an University of Technology, Xi'an 710048, China^{2}University of Florida, Gainesville, FL 32611, USA^{3}College of Management, Shenzhen University, Shenzhen 518060, China

Received 16 August 2012; Accepted 14 January 2013

Academic Editor: Arunita Jaekel

Copyright © 2013 Junli Liang et al. 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

A novel circle fitting algorithm is proposed in this paper. The key points of this paper are given as follows: (i) it formulates the circle fitting problem into the special source localization one in wireless sensor networks (WSN); (ii) the multidimensional scaling (MDS) analysis is applied to the data points, and thus the propagator-like method is proposed to represent the circle center parameters as the functions of the circle radius; (iii) the virtual source localization model can be rerepresented as special nonlinear equations of a unique variable (the circle radius) rather than the original three ones (the circle center and radius), and thus the classical fixed-point iteration algorithm is applied to determine the radius and the circle center parameters. The effectiveness of the proposed circle fitting approach is demonstrated using the simulation and experimental results.

#### 1. Introduction

Circle fitting receives considerable attention because it plays an important role in computer vision, observational astronomy, structural geology, industry inspection, medical diagnosis, Iris recognition, military, security, and so forth [1–8]. For instance, to meet the increasing demand for manufacturing automation, the circle fitting technique is often applied to measure the diameter of the processing product in the manufacturing systems.

The fitting problem can be viewed as follows: estimate the parameters of a circle from a set of coplanar points. Several classical approaches [1–8], including the Hough transform (HT) methods [4, 5] and the least square (LS) approaches [6–8], have been developed to solve this problem. The former are actually to carry out a voting procedure in a three-dimensional (3D) Hough accumulator space, where every point represents a circle of a certain size. The corresponding coordinate of the local maxima is obtained as the estimated parameters of the circle. In comparison, the latter attempt to find the parameters of a circle by minimizing an error metric between the primitive and the data points.

In this paper, we develop a novel circle fitting approach by borrowing the idea from source localization in wireless sensor networks (WSN) [9, 10]. It is worthwhile to highlight the main contributions of this paper here.(i)It formulates the circle fitting problem into special source localization one in WSN, where each data point should be understood as an abstract sensor node in sensor networks, and the circle center represents the localized target. However, the propagation delays are unknown, and thus the existing source localization algorithms in WSN cannot be applied to solve the special source localization problem.(ii)The multidimensional scaling (MDS) analysis [11] is applied to the data points, and a special covariance-like matrix is constructed. Thus, we propose the propagator-like method to represent the circle center as the functions of the circle radius.(iii)The virtual source localization model can be rerepresented as special nonlinear equations, where the radius is the unique variable rather than the original three ones (the circle center and radius), and thus the classical fixed-point iteration algorithm [12] is applied to determine the radius and circle center.

The rest of this paper is organized as follows. The circle fitting problem is described in Section 2. A novel circle fitting approach is developed in Section 3. Simulated and experimental results are presented in Section 4. The paper is concluded in Section 5.

#### 2. Problem Formulation

The equation for a circle centered at with a radius in coordinates has the following form:

The circle fitting problem (CFP) [1–8] can be described in Figure 1, that is, given data points , , the objective of circle fitting is to estimate circle parameters that best fit to these data points. However, since in the actual application noise is introduced by some operations (e.g., the segmentation and edge detection operations in the image processing application), these obtained points are not completely precise, that is, where is the introduced noise. The objective of this paper is to estimate the circle parameters from the given data points , .

#### 3. Proposed Algorithm

In this section, we first reformulate CFP into a virtual source localization problem in wireless sensor networks (WSN) [9, 10] and then develop a novel circle fitting algorithm in this framework.

Let us review the source localization model in WSN [9, 10]: where is the propagation speed of light; is the unknown location of the emitter; is the location of the sensor node in wireless sensor networks; is the (ideal) signal propagation delay from the target to the sensor node, but is the available delay measurement, containing noise . In addition, denotes the Euclidean norm of vector . The objective of source localization in WSN is to estimate from the given measurement delays , .

To reformulate CFP into the source localization problem in WSN, we rewrite (2) in another form as where , and .

The source localization model in (3) is quite similar to the circle fitting model in (4), especially when : (i) each data point should be understood as a virtual “sensor node” in sensor networks; (ii) the circle center represents the virtual “emitter" or localized “target"; and (iii) the circle radius is the virtual “propagation delay”, which are clearly described in Figure 2.

By comparing (3) with (4), we can easily observe their differences; that is, all “propagation delays” from the “target” to “sensor nodes” or are unknown. Therefore, the existing source localization algorithms cannot solve since they require the knowledge of “propagation delays” , , or .

In the rest of this section, we will develop a novel algorithm for estimating the “target position” and the “propagation delays” of all sensor nodes.

Let And then define an similarity matrix [11]: the element of which can be represented as

Under the ideal (without noise) case, . Note that is unknown, and thus is actually the function of , that is, .

Since the rank of equals 2, the rank of is also 2. From , we introduce the following partition: where and are the first two and last columns of , respectively.

Similar to the conventional propagator method [13], we define the propagator which satisfies

Let and stand for the first and second column of . Based on the propagator method [13], we have where denotes the vector with all elements 1.

From (11), we can solve and as follows:

Note that depends on and thus and are not determined from (12) directly and can only be represented as the functions of the unknown “propagation delay” .

Plugging (12) (i.e., and ) into (4), we can obtain which implies that is the root of the equations previously mentioned in absence of noise.

According to the fixed-point iteration theory [12], is the fixed point of the function , that is, the value that the function “locks onto” in the iterative process. Thus, we can compute “propagation delay” using the following iteration procedure: which are combined to yield the following iteration process:

Once circle radius is obtained from the previously mentioned iterative procedure, the circle center can be solved from (12).

#### 4. Simulation Results

In this section, some experiments are conducted to evaluate the performance of the proposed method. For comparison, we simultaneously implement the HT method [4, 5] and the LS approach [6–8].

##### 4.1. Experiment 1

The first experiment is implemented on data points of a circle shown in Figure 3, where the noise is stochastic additive white Gaussian noise with zero mean and variance 0.01, and the true circle parameters are , and . The initial value of is 0.1. Figure 4 shows the realized , and using the proposed algorithm with 50 iterations. We can see from Figure 4 that the realization of , and generally stays around 10, 8, and 15, respectively, after 3 iterations in this problem, which shows that the proposed algorithm converges rapidly. Table 1 gives the estimated results using the proposed algorithm, the LS method, and the Hough transform approach. We can see from Table 1 that the proposed algorithm has the approximate estimation accuracy as those of HT and LS and is slightly more accurate in estimating than the LS method. Figure 5 shows the fitting results using different algorithms, which further shows that the fitting result obtained by the proposed algorithm approaches the true circle.

##### 4.2. Experiment 2

In this experiment, the proposed algorithm is applied to the real data. Figure 6 (resolution ) shows an example which computes the diameter of the bright hole after orientating the clouds of points to get the hole plane parallel to the projection plane (Available from http://www.aqsense.com/docs/docu/Compatibility.html). We can obtain the inner edge points of the hole via threshold segmentation, edge detection, and spectral clustering and then fit these points. Figure 7 shows the realized , and using the proposed algorithm with 50 iterations. Table 2 lists the estimation results using the proposed algorithm, the HT method, and the LS approach.

##### 4.3. Experiment 3

Iris recognition is a biometric identification technique based on images of the irides of an individual’s eyes. Since the Iris area lies between the pupil region (a dark ellipse with the lowest intensity) and limbus region, determining the pupil region is an important preprocessing step of Iris localization. In the third experiment, we implement the proposed algorithm on the Iris image, as shown in Figure 8. Via the thresholding segmentation and Sobel edge detection, edge points are given in Figure 9. Table 3 gives the estimation results from three different methods. Figure 10 displays the fitting results by the proposed algorithm that are marked by red points, and Figure 11 shows the realized , and using the proposed algorithm with 50 iterations, which shows that the proposed algorithm can fit the Pupil’s boundary well.

Although the HT method is of the highest estimation accuracy, it needs to be pointed out that the HT method requires (i) quantizing the three-dimensional space finely enough; otherwise the peaks in the transform plane will be broadened and (ii) the overwhelming burden of the three-dimensional search in the , and plane.

#### 5. Conclusion

In this paper, we propose a novel circle fitting algorithm by borrowing the idea from source localization in wireless sensor networks. Since the virtual propagation delays of all sensor nodes are unknown, the existing source localization algorithms cannot be applied. This paper formulates the virtual source localization model of three unknown parameters into special nonlinear equations of a unique parameter, that is, the circle radius , using the MDS analysis and propagator method, and then it employs the classical fixed-point iteration theory to determine the circle radius and circle center.

#### Acknowledgments

This work was supported in part by the National Natural Science Foundation of China under Grants 61172123 and 60901059, by Excellent Youth Research Star (2012KJXX-35), and educational Department Foundation of Shaanxi Province (12JK0526).

#### References

- U. M. Landau, “Estimation of a circular arc center and its radius,”
*Computer Vision, Graphics and Image Processing*, vol. 38, no. 3, pp. 317–326, 1987. View at Publisher · View at Google Scholar · View at Scopus - J. F. Crawford, “A non-iterative method for fitting circular arcs to measured points,”
*Nuclear Instruments and Methods in Physics Research*, vol. 211, no. 1, pp. 223–225, 1983. View at Scopus - G. Coath and P. Musumeci, “Adaptive arc fitting for ball detection in robocup,” in
*Proceedings of APRS Workshop on Digital Image Computing*, pp. 63–68, Brisbane, Australia, February 2003. - R. O. Duda and P. E. Hart, “Use of the hough transformation to detect lines and curves in pictures,”
*Communications of the ACM*, vol. 15, no. 1, pp. 11–15, 1972. View at Publisher · View at Google Scholar · View at Scopus - D. J. Kerbyson and T. J. Atherton, “Circle detection using hough transform filters,” in
*Proceedings of the 5th International Conference on Image Processing and its Applications*, pp. 370–374, July 1995. View at Publisher · View at Google Scholar · View at Scopus - H. Spath, “Least-square fitting by cirles,”
*Computing*, vol. 57, no. 2, pp. 179–185, 1996. View at Publisher · View at Google Scholar - I. D. Coope, “Circle fitting by linear and nonlinear least squares,”
*Journal of Optimization Theory and Applications*, vol. 76, no. 2, pp. 381–388, 1993. View at Publisher · View at Google Scholar · View at Scopus - L. Moura and R. Kitney, “A direct method for least-squares circle fitting,”
*Computer Physics Communications*, vol. 64, no. 1, pp. 57–63, 1991. View at Scopus - H. C. So and F. K. W. Chan, “A generalized subspace approach for mobile positioning with time-of-arrival measurements,”
*IEEE Transactions on Signal Processing*, vol. 55, no. 10, pp. 5103–5107, 2007. View at Publisher · View at Google Scholar · View at Scopus - A. Beck, P. Stoica, and J. Li, “Exact and approximate solutions of source localization problems,”
*IEEE Transactions on Signal Processing*, vol. 56, no. 5, pp. 1770–1778, 2008. View at Publisher · View at Google Scholar · View at Scopus - T. F. Cox and M. A. A. Cox,
*Multidimensional Scaling*, Chapman & Hall/CRC, Boca Raton, Fla, USA, 2001. - R. L. Burden and J. D. Faires,
*Numerical Analysis: Fixed-Point Iteration*, PWS Publishers, 2010. - S. Marcos, A. Marsal, and M. Benidir, “The propagator method for source bearing estimation,”
*Signal Processing*, vol. 42, no. 2, pp. 121–138, 1995. View at Scopus