Research Article  Open Access
Implementation of Steiner Point of Fuzzy Set
Abstract
This paper deals with the implementation of Steiner point of fuzzy set. Some definitions and properties of Steiner point are investigated and extended to fuzzy set. This paper focuses on establishing efficient methods to compute Steiner point of fuzzy set. Two strategies of computing Steiner point of fuzzy set are proposed. One is called linear combination of Steiner points computed by a series of crisp αcut sets of the fuzzy set. The other is an approximate method, which is trying to find the optimal αcut set approaching the fuzzy set. Stability analysis of Steiner point of fuzzy set is also studied. Some experiments on image processing are given, in which the two methods are applied for implementing Steiner point of fuzzy image, and both strategies show their own advantages in computing Steiner point of fuzzy set.
1. Introduction
Associated with every closed bounded convex set in is a point known as Steiner point or curvature centroid [1]. Though the Steiner point of smooth convex curves was defined and investigated by Seiner already in 1840, some properties, like additivity, were discovered almost in 1960s by Shephard [2]. From then on, Steiner point was known as one of the important geometry points and was studied more and more. Nowadays, the Steiner point is regarded as the only pointvalued function defined for all convex bodies, which is additive and uniformly continuous and commutes with similarity transformations [1].
The utility of the Steiner point extends beyond its definition as a robust center of a set of static points. Locating the Steiner point of an object is helpful for many tasks, because Steiner point is an invariant point of an object, while a transform is used on it in certain ways, such as growing uniformly in all directions, moving in a line, and rotating around an axis [3]. By finding the Steiner point of an object, one can analyze some properties of an image [4]. To detect or recognize an object in an image, Steiner point can help us in some cases [5]. If two objects are similar but have different Steiner points, one can distinguish them in this way [6]. Tracking moving objects is now a popular approach for research workers [7–10]; if Steiner points of objects are referenced, they could save a large amount of computation.
In the early years, much work has been done on some algebraic and analytic structure and behavior of Steiner points, such as linear translation, continuity, and even affine translation of an object. Three important properties were studied and known as basic properties of Steiner points, which are shortly denoted by commutation, addition, and continuity [2, 11, 12]. Furthermore, the definition of Steiner point was generalized from a polytope to a nonempty compact subset of [13]. Furthermore, the stability and eccentricity of Steiner point were researched and applied to mobile facility location [14]. To implement calculation of Steiner point, there are several alternatives [4]; one of the efficient ways refers to [1], which is based on the exterior angle of convex points in a polytope. In recent years, Steiner point has been extended to fuzzy set and provides an alternative strategy of defuzzication, which is regarded as the center of the fuzzy set [3, 15].
This approach focuses on implementing Steiner point of fuzzy set. The motivation is trying to find an efficient method to calculate the Steiner point of a fuzzy set. This paper is arranged as follows. In the second part, some definitions and properties of Steiner point of fuzzy set are investigated referring to the literatures. The third part discusses strategies to compute the Steiner point of a fuzzy set. Two main methods are proposed for calculating the Steiner point of a fuzzy set. In the fourth part, stability analysis of Steiner point of fuzzy set is proposed. In the fifth part, some experiments on image processing are presented. The last part of the paper contains the conclusions.
2. Definitions of Steiner Point and Properties
In the following let us suppose that is an integer. We denote by the set of nonempty compact convex subsets of . The set is endowed with a linear structure in which the addition of two subsets and the multiplication of a subset by a positive real number are defined pointwise. We furthermore endow with the Hausdorff metric . Let denote the unit sphere in , and let denote the space of continuous functions from to , endowed with the supremum norm. Now, for , we define the support function of ; see, for example, [16], by where denotes the usual inner product of . The following definition is due to [3, 11].
Definition 1. The Steiner point of in is defined as where varies over the unit vectors of , is the Lebesgue measure on , and is the volume of the unit ball of . Notice that .
We denote by the set of all functions from to which are (i) decreasing and (ii) left continuous on and continuous at 0. For a fuzzy set and a rigid motion we set , .
Definition 2. Let be a measure function. For , let where is the Steiner point of crisp sets. Then is a Steiner point of fuzzy set .
Based on Definition 1, the following properties refer to [2, 12].
Theorem 3. Let have the following properties:(S1)for any , ;(S2)for and any rigid motion , one has ;(S3) is continuous.
Then is the Steiner point. These three properties are described in [11] as addition, commutation, and continuity of Steiner point.
The following theorem is due to [3], which is an extension of Theorem 3 to the case of fuzzy set.
Theorem 4. A function is called a Steiner point if it has the following properties:(SF0)for any , ;(SF1)for any , ;(SF2)for and any rigid motion , one has ;(SF3) is continuous.
As mentioned in [17, 18], a Steiner point of fuzzy set is not defined unambiguously by the properties (SF0)–(SF3). It is amazingly difficult to impose further properties on to obtain uniqueness; it is an open question if this is possible in some reasonable way. For the purpose of calculating Steiner point of fuzzy set, we introduce the following definitions and lemmas [3].
Definition 5. Let be a division of , which is . Then one calls a fuzzy set a step fuzzy set if it is constant on , respectively. One denotes by the set of all step fuzzy sets.
Definition 6. Let be a division of . Let be a function fulfilling the properties (SF0)–(SF3) of Theorem 4. Then one calls a Steiner point a step Steiner point if it is constant on , respectively. One denotes by the set of all step Steiner points.
Lemma 7. Let be a division of . Let be a function fulfilling the properties (SF0)–(SF3) of Theorem 4. Then there are unique real numbers such that and, for all ,
Lemma 8. Let be a division of . Let be a function fulfilling the properties (SF0)–(SF3) of Theorem 4. Let be the unique real numbers fulfilling (4). Then .
For the general case, the Steiner point of a fuzzy set can be calculated as follows.
Lemma 9. Let be a Steiner point. Then there is a measure function such that, for all ,
3. Calculation of Steiner Point of Fuzzy Set
From Lemma 9, we know that computing Steiner point of a fuzzy set can be transformed to computing Steiner point of cut sets of a fuzzy set, which means that defuzzification is necessary in fuzzy set Steiner point computing. Defuzzification methods may be divided into two classes, considering either the horizontal or the vertical representation of the fuzzy set. In the former, one assigns to each possibility value a set of elements of the universe in agreement with that possibility value, that is, cut representation, while in the latter one assigns to each element of the universe a possibility value. In this paper we focus on the former and discuss a little about the later.
Definition 10. Let be a scale. Let be a fuzzy set. Then one calls an cut and , which is a crisp set.
Now there are two strategies to compute Steiner point of a fuzzy set. One is to find a series of cut sets of the fuzzy set, compute Steiner point for each cut set, and combine them in linear form as in (4). The other strategy is to try to find a crisp set in the cut sets, which has the same Steiner point with the fuzzy set. We will discuss these two methods in detail in the following.
In the case of step fuzzy set, we have fixed the number of cut sets, so it is easy to transform a fuzzy set into a series of crisp sets. But there may be a large number of cut sets, which are not necessary in calculating Seiner point of a fuzzy set. So, in this case, we prefer the step Steiner point, which is as defined in Definition 6.
From Lemma 7, it is known that once cut sets of fuzzy set are confirmed, which can be denoted by , computing Steiner point of fuzzy set is equivalent to computing Steiner points of the series of cut sets, namely, , and combining them by a series of weights which satisfy and . One of the choices is simply taking where . In the view of numerical computing, calculating Steiner point of a fuzzy set in this way can be much more complex, while a lot of level sets are taken for the fuzzy set.
For the sake of convenience in computing Steiner point of a fuzzy set, we introduce the second strategy (which is similarly approximate). Firstly, consider the following definition.
Definition 11. Let be a cut of fuzzy set . Then one calls the optimal approximate crisp set to if where is the distance between crisp set and fuzzy set .
Theorem 12. Let , one of cut sets of fuzzy set , be the optimal approximate crisp set to fuzzy set . Then , as defined in Definition 1, namely, is a Steiner point of fuzzy set .
Proof. We prove that satisfies (S1)–(S3) in Theorem 3. Denote by a real number and denote by and the cut sets of fuzzy sets , respectively. According to [17], and ; then
So (S1) and (S2) are satisfied. Rather, more is continuous, which is (S3). This completes the proof.
Another motivation is from [15, 17] and Lemma 9. If we rewrite (5) as it is clear that Steiner point of a fuzzy set is the average of the Steiner points of all the level sets . Note that is continuous with respect to ; by the Mean Value Theorem, there exists an such that
Theorem 13. If is a convex fuzzy set and is continuous with respect to , then in (12) satisfies
It follows from [15] that the Steiner point is a characteristic point of a fuzzy set in the sense of where is the metric on fuzzy space . This implies the following.
Theorem 14. For any Steiner point of fuzzy set there is , where is an cut set of fuzzy set and satisfies (8).
Theorem 15. A Steiner point of , which is an cut set of fuzzy set and satisfies (8), is also a Steiner point of fuzzy set ; namely,
If it is reasonable, we can introduce the following definitions.
Definition 16. For , one defines the support function of by where is a membership function of the fuzzy set , and
Definition 17. For , the Steiner point of fuzzy set is given by where varies over the unit vectors of , is the Lebesgue measure on , and is the volume of the unit ball of .
Now what we need is computing Steiner point based on crisp set, for a 2dimensional case, based on polygon or convex polygon. For a 2D set, there are two steps to compute its Steiner point in numerical sense. The first step is to find all the convex points of the set and to form a convex polygon , which is proved to be linear computational complex in [19]. The second step is referring to [2, 20]. Consider where , for , is the proportion to of the external angle of convex polygon at .
4. Stability Analysis of Steiner Point
From [15], we know the following fact. For any cut of a is a convex compact subset of and can be uniquely characterized by its socalled support function. Therefore, the fuzzy set itself is uniquely characterized by the function Now, the Steiner point of a fuzzy set with is given by where is the Lebesgue measure on with . So if we denote the Steiner point of a fuzzy set takes the following form:
Supposing are any two values, the corresponding cut sets are and . Now let us discuss their differences in computing the Steiner points: Considering that is bounded, namely, there exists an upper boundary for each such that Therefore where is a Lipschitz constant.
Let us consider an example in the following. Suppose that is a crisp set and is the corresponding fuzzy set with the following membership: That is, . The cut of fuzzy set is given by Notice the following fact: and (25); then This leads to
5. Experimental Examples
In this section, we investigate four fuzzy images as an example, which were proposed in [21]; see Figure 1. Figure 1(A) shows a synthetic fuzzy set, described by a membership function radially nonincreasing from the centroid. Figure 1(B) presents a fuzzy segmented slice of a threedimensional magnetic resonance angiography (MRA) image of a human aorta at the position where it splits into the two iliac arteries. Figure 1(C) shows a part of a histological light microscope image of a bone implant (inserted in a leg of a rabbit). The selected part of the fuzzy segmented image contains a bone area, surrounded by a nonbone area. Figure 1(D) shows a synthetic fuzzy set caused by motion, as described in [22], with and . We illustrate the implementation of the Steiner point of the fuzzy object in Figure 1 with the method presented in Section 3 and the comparison of the linear combination Steiner point with the proposed approximate Steiner point.
There are some technical issues that should be interpreted here. First, commonly not always all grays appear in an image, so we can find the minimum gray (denoted by ) and the maximum one (denoted by ) of the image by computing histogram of the image. Then, we divide the interval into several levels with equal metric; that is, , and each gray corresponds to an cut of the fuzzy image. Second, in order to define the distance between the fuzzy image (the original image) and certain level image, we unify the grays of the fuzzy image by , where is the gray matrix of the original fuzzy image, unify all the image grays of level , , which is denoted by , and compute the distance by . Third, considering that the Steiner point of an object in an image is controlled by the shape of the object and the shape of the object depends greatly on boundary detecting, in most cases, boundaries of an object can be detected well in certain interval of gray, for example, in the middle part of interval . So we can choose part of the interval and divide it into several levels to compute the Steiner points. Here, we choose half of the interval as the considered domain; namely, or .
In this paper, every degree of cut image is given according to the gray level of the image. Steiner points of the synthetic images, Figures 1(A) and 1(D), are shown in Figures 1(a) and 1(d), and Steiner points of the two real fuzzy segmented images, Figures 1(B) and 1(C), are given in Figures 1(b) and 1(c). In each image, the linear combination of Steiner point with equal weight of for all levels is given and marked by “o.” Also, given in Figure 1, the approximate Steiner point is marked by “+,” which has the minimum distance between the fuzzy image and those gray images in all cut sets.
In Figure 2, we plot the distance between the Steiner point of each gray level and the linear combined Steiner point in Figures 1(A), 1(B), 1(C), and 1(D). Also, we present their corresponding cut sets in Figure 1. The minimal value of distance between the approximate Steiner point and the combination of Steiner point is indicated by “*” in the plot. The corresponding value of provides the optimal approximate cut. The three Steiner points calculated by cut set with the minimum distance, with the second smallest distance, and at , respectively, are shown the last three columns in Figure 1. Table 1 presents the Steiner point of the object obtained by different methods and the distance between two Steiner points. var1 in Table 1 represents the variance of distance between the Steiner point of each gray level and the linear combination of Steiner point, while var2 presents the variance of distance between Steiner points of two adjacent gray levels. Results, expressed in var1 and var2, show that the distance between Steiner point of each gray level and the linear combination of Steiner point is more stable than the distance between Steiner points of two adjacent gray levels.
Experiment Results Analysis. Both two strategies of implementing Steiner point of fuzzy set have their own advantages and shortcuts. With the growing of the distance of image gray, the distance between the linear combined Steiner point and the Steiner points of each gray level do not enlarge rapidly. This means that the linear combined Steiner point shows more stability than the approximate Steiner point. Unfortunately, the former needs more computational time than the later, especially when more levels of gray are chosen to compute Steiner point. In our opinion, a suitable weight for combining the Steiner points of all the levels of fuzzy image is not reasonable to perform without having a particular application in mind. Also, from Figure 2, we see that the cut image is obtained by cut at optimal value of in second column, which is chosen as the approximate Steiner point.
6. Conclusion
This approach focuses on implementing Steiner point of fuzzy set and some properties of Steiner point on fuzzy set. We try to find some efficient methods to compute Steiner point of fuzzy set. Two strategies of computing Steiner point of fuzzy set are proposed, namely, the linear combination of Steiner point, which calculates the Steiner point based on the approximate cut set. We also discuss some stable properties of Steiner point of fuzzy set and give some experiments on image processing. However, there are still some open problems which need more investigation while implementing Steiner point of a fuzzy set, such as how to choose each level of a fuzzy set, which is suitable for defuzzificating a fuzzy set? Does there exist unique choice of levels for fuzzy set in computing Steiner point? How to choose suitable weights for combining Steiner point of different levels? Those problems have high potential value in image processing.
Appendix
Implementations of Algorithm
See Algorithm 1.

Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
This work was supported by CMP Laboratory, Department of Cybernetics, Faculty of Electrical Engineering, Czech Technical University. One of the authors was supported by the Agreement between Czech Ministry of Education and Chinese Ministry of Education. This work was also supported by Blue Project of Universities in Jiangsu Province Training Young Academic Leaders Object and the National Natural Science Foundation of China (no. 61170121).
References
 B. Grunbaum, Convex Polytopes, Springer, New York, NY, USA, 2nd edition, 2003.
 G. C. Shephard, “A uniqueness theorem for the Steiner point of a convex region,” Journal of the London Mathematical Society, vol. 43, no. 1, pp. 439–444, 1968. View at: Publisher Site  Google Scholar
 T. Vetterlein and M. Navara, “Defuzzification using Steiner points,” Fuzzy Sets and Systems, vol. 157, no. 11, pp. 1455–1462, 2006. View at: Publisher Site  Google Scholar
 S. R. Sternberg, “Grayscale morhpology,” Computer Vision, Graphics, and Image Processing, vol. 35, no. 3, pp. 333–355, 1986. View at: Google Scholar
 J. Z. Liang and W. Song, “Clustering based on Steiner points,” Machine Learning, vol. 3, no. 2, pp. 141–148, 2011. View at: Google Scholar
 H. Alt, O. Aichholzer, and G. Rote, “Matching shapes with a reference point,” International Journal of Computational Geometry and Applications, vol. 7, no. 4, pp. 349–363, 1997. View at: Google Scholar
 A. J. Lipton, H. Fujiyoshi, and R. S. Patil, “Moving target classification and tracking from real time video,” in Proceedings of the 4th IEEE Workshop on Applications of Computer Vision (WACV'98), pp. 8–14, 1998. View at: Google Scholar
 S. S. Intille, J. W. Davis, and A. F. Bobick, “Realtime closedworld tracking,” in Proceedings of the IEEE Computer Society Conference on Computer Vision and Pattern Recognition ( CVPR '97), pp. 697–703, 1997. View at: Google Scholar
 I. Haritaoglu, D. Harwood, and L. S. Davis, “W^{4}S: a real time system for detecting and tracking people in $2.5\text{D}$,” in Proceedings of the 5th European Conference on Computer Vision (ECCV '98), 1998. View at: Google Scholar
 C. Isaac and M. Gerard, “Detecting and tracking moving objects for video surveillance,” Computer Vision and Pattern Recognition, vol. 2, pp. 23–25, 1999. View at: Google Scholar
 W. J. Meyer, “Characterization of the Steiner point,” Pacific Journal of Mathematics, vol. 35, no. 3, pp. 717–725, 1970. View at: Publisher Site  Google Scholar
 R. Schneider, “On steiner points of convex bodies,” Israel Journal of Mathematics, vol. 9, no. 2, pp. 241–249, 1971. View at: Publisher Site  Google Scholar
 J. Mattioli, “Minkowski operations and vector spaces,” SetValued Analysis, vol. 3, no. 1, pp. 33–50, 1995. View at: Publisher Site  Google Scholar
 D. Stephane and D. Kirkpatrick, “The Steiner centre of a set of points: Stability, eccentricity, and applications to mobile facility location,” International Journal of Computational Geometry and Applications, vol. 16, no. 4, pp. 345–371, 2006. View at: Publisher Site  Google Scholar
 R. Korner and W. Nather, “Linear regression with random fuzzy variables: extended classical estimates, best linear estimates, least squares estimates,” Information Sciences, vol. 109, no. 14, pp. 95–118, 1998. View at: Google Scholar
 R. Schneider, Convex Bodies: The BrunnMinkowski Theory, Cambridge University Press, Cambridge, UK, 1993.
 D. Butnariu, M. Navara, and T. Vetterlein, “Linear space of fuzzy vectors,” in Fuzzy Logics and Related Structed, Liuz Austria, pp. 23–26, 2005. View at: Google Scholar
 R. T. Rockaellar, Convex Analysis, Princeton University Press, Princeton, NJ, USA, 1970.
 J. Luis and R. Muniz, “Revisiting Bobylevąŕs definition,” in Proceedings of the 5th Conference of the European Society for Fuzzy Logic and Technology, pp. 11–14, Stept, 2007. View at: Google Scholar
 C.L. Chen, “Computing the convex hull of a simple polygon,” Pattern Recognition, vol. 50, no. 22, pp. 561–565, 1989. View at: Google Scholar
 N. Sladoje, http://www.inf.uszeged.hu/ssip/2006/lectures/fuzzy_application.pdf.
 N. Sladoje, http://www.inf.uszeged.hu/ssip/2011/PDF/LectureN1hour_Sladoje.pdf.
Copyright
Copyright © 2014 Jiuzhen Liang and Dejiang Wang. 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.