Table of Contents Author Guidelines Submit a Manuscript
Journal of Applied Mathematics
Volume 2014 (2014), Article ID 285045, 8 pages
Research Article

Fuzzy B-Spline Surface Modeling

School of Informatics & Applied Mathematics, Universiti Malaysia Terengganu (UMT), 21030 Kuala Terengganu, Terengganu, Malaysia

Received 23 October 2013; Revised 10 January 2014; Accepted 18 February 2014; Published 6 July 2014

Academic Editor: Jong Hae Kim

Copyright © 2014 Rozaimi Zakaria 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.


This paper discusses the construction of a fuzzy B-spline surface model. The construction of this model is based on fuzzy set theory which is based on fuzzy number and fuzzy relation concepts. The proposed theories and concepts define the uncertainty data sets which represent fuzzy data/control points allowing the uncertainties data points modeling which can be visualized and analyzed. The fuzzification and defuzzification processes were also defined in detail in order to obtain the fuzzy B-spline surface crisp model. Final section shows an application of fuzzy B-spline surface modeling for terrain modeling which shows its usability in handling uncertain data.

1. Introduction

Data points are collected from physical objects to capture its geometric entity and representation in a digital system, that is, CAD systems. This data is collected by using specific devices such as scanning tools. However, the recorded data do not necessarily represent error-free data. It is due to the fact that the errors are produced by the limitations of the tools, environmental factors, human errors, and so forth. Usually, these kinds of data which have uncertainty characteristics cannot be used directly to produce digitized models. Hence, designers use certain type of digital filter to remove and amend the errors involved which is a painstaking process [1]. One may not truly capture the digital model of a scanned model due to the reasons stated above regardless of executing the time consuming digital filter.

In order to make the uncertain data useable for analysis and model building, these data have to be defined in a different approach which will incorporate uncertainties of the measurements. In this paper, we propose fuzzy set theory which was introduced by Zadeh in 1965 [2]. It has been widely used in dealing with uncertain matters for wise decision making processes. Readers are referred to [35] for detailed explanation regarding the subject matter.

Natural spline, Bezier, and B-spline functions are example of functions which can be used to create CAD models with data points [68]. These curves and surfaces created with the stated functions are the standard approach to represent a set of collected data points. These curves and surfaces are used to visualize and analyze the CAD models. B-spline functions can be used to design curves and surface using either approximation or interpolation methods to model the real data points [915]. The reason why B-spline functions with its weight known as nonuniform rational B-spline (NURBS) being used in many applications is that the designer can easily tweak the control points to obtain a desired shape easily.

This paper discusses modeling of interpolating B-spline surface using fuzzy set theory and it is organized as follows. Section 2 reviews the previous work on modeling uncertain data via B-spline function in the form of curves and surface. Section 3 discusses the representation of data points using fuzzy set theory, fuzzy number, and fuzzy relation concepts. Section 4 discusses blending of fuzzy data/control point with B-spline curve and surface function where the end results in fuzzy B-spline curve and surface model. This section also defines fuzzification and defuzzification processes. To show the application of fuzzy B-spline surface model, we apply the proposed method to model lakebed from a set of uncertain data in Section 5. This section also compares the result of statistical analysis to show the effectiveness of proposed fuzzy B-spline model.

2. Previous Works

The requirement of fuzzy set theory is essential in handling ambiguous data in order to create a model using B-spline curve and surface function. The designers are unable to choose the appropriate control points which are exposed to errors and uncertainties due to the reasons stated above. Therefore, this is the reason why we define the uncertain data by using fuzzy set theory and model them through B-spline curve and surface function.

There are a number of methods that have been developed in dealing with uncertain data. Examples include modeling surface of Mount Etna which was proposed by Gallo et al. [16, 17]. In this paper, they developed a fuzzy B-spline model which was conceptualized from fuzzy numbers in the form of intervals. They also proposed alpha level within . The data of Mount Etna was represented in the form of interval fuzzy surface based on alpha value. Anile et al. [18] further enhanced this method for modeling uncertain and sparse data.

Both of the methods discussed above do not have defuzzification process to obtain a crisp fuzzy surface (defuzzified surface). Although both methods employed fuzzification phase, it is still in interval surface form which is not in single surface form. Furthermore, the fuzzification and defuzzification process of B-spline fuzzy system [19] were not elaborated in detail. Therefore, this paper elucidates defuzzification process upon the application of fuzzification process.

3. The Process of Defining Fuzzy Data

This section defines uncertain data based on the concept of fuzzy number [35, 20] based on the interval of fuzzy number.

Definition 1. Let be a universal set in which is a real number and is subset of . A Fuzzy set in is called a fuzzy number and expressed using the -level with various -cut; that is, if for every , there exist set in where and [21].
Definition 1 provides the basis to define uncertain data in which these data are in the form of real numbers. We use triangular fuzzy number for defining the uncertain data in the form of interval. Therefore, the triangular fuzzy number can be defined in Definition 2 as follows [21].

Definition 2. If triangular fuzzy number is represented as and is a -cut operation of triangular fuzzy number, then crisp interval by -cut operation is obtained as with where the membership function, , is given by where and are left and right fuzzy number which form the interval fuzzy number and is the crisp point in the interval. The symbol means the -level values of triangular fuzzy -cut (see Figure 1).

Figure 1: Triangular fuzzy number, .

In Figure 1, the value is the crisp fuzzy number which has full membership function; that is, it is equal to 1.

Before we define the uncertain data points, we must define the fuzzy relation which is used as a converter from the definition of fuzzy number to the definition of fuzzy data points in real numbers. Then, the definition of fuzzy relation and also the definition of the relation between two fuzzy points can be given in Definition 3 till Definition 5.

Definition 3. Let be universal sets; then is called a fuzzy relation on [4, 21].

Definition 4. Let and and are two fuzzy sets. Then is a fuzzy relation on and if , and , [4].

Definition 5. Let with and represent two fuzzy data. Then, the fuzzy relation between both fuzzy data is given by .
The uncertain data points can be defined after fuzzy number and fuzzy relation had been defined. This uncertain data point becomes fuzzy data point (FDP) as shown in Definition 6.

Definition 6. Let and is the set of FDPs which is with is the universal set and is membership function defined as in which . Therefore, with where and are left-grade and right-grade membership values, respectively. This can be written as for all , with , , and being left FDP, crisp data point, and right FDP, respectively [20]. The procedure in defining FDP is illustrated in Figure 2.

Figure 2: The process of defining FDPs.

Figure 2 shows the approach to transform ordinary data point to Fuzzy data point (FDP). The membership grades of FDPs in the form of are illustrated in Figure 3.

Figure 3: FDP form after being defined by fuzzy relation.

Figure 3 shows the formation of FDP by using the definitions of fuzzy relation and fuzzy number. The construction of FDPs is in - and -axis.

Definition 6 gives us the definition of FDP in 2D form and for FDP in 3D form similar concept is applicable given by (5) (based on (4)). Consider

4. The Proposed Method

This section discusses blending of FDPs into B-spline function to produce fuzzy B-spline curves and surfaces. The fuzzification and defuzzification process towards fuzzy B-spline model in the form of either curves or surface are also discussed.

Based on B-spline function [68], we can define the fuzzy B-spline model [22, 23] as follows.

Definition 7. A fuzzy B-spline curve is a function which represents a curve to the set of real fuzzy numbers and it is defined as where are fuzzy control points which are also known as fuzzy data point and are B-spline basic function with crisp knot sequences where represents the degree of B-spline function and represents the numbers of control points.

Definition 8. A fuzzy B-spline surface is defined by the following equation: where (i) and are B-spline basic function of degrees and with crisp parameters of and in ; (ii) each vector knot must satisfy the conditions and ; (iii) are fuzzy control point in th row and th column.
Therefore, both of Definitions 7 and 8 can be illustrated in the form of numerical examples with Figure 4.
For fuzzy B-spline curve model, the resultant curve interpolates the first and last fuzzy control points. This fuzzy curve was designed based on five fuzzy control/data points (). The same concept is applied to surface which uses diagonal fuzzy control/data points having 16 fuzzy control/data points.
For fuzzy B-spline curve model, the fuzzy control/data points are defined for -element and for fuzzy B-spline surface model, and the fuzzy control/data points are defined at -element. Therefore, defining uncertain data can be done by combining either the tuple axis or one of the axes while maintaining other axes as crisp values [21].
Upon defining the fuzzy B-spline model, we use the alpha-cut operation of triangular fuzzy number to do the fuzzification process based on Definition 2. Therefore, the fuzzification process of fuzzy control points is given by the following definition which is intact with the B-spline surface function.

Figure 4: The fuzzy B-spline models in the form of (a) curves and (b) surface.

Definition 9. Let be the set of fuzzy control points where and . Then, is the alpha-cut operation of fuzzy control point which is given as the following equation where with :
Upon fuzzification, the next procedure is the defuzzification process. Defuzzification process is applied to obtain a fuzzy solution in a single value. The result of defuzzification process is also known as fuzzy crisp solution. Therefore, the defuzzification process is defined as in Definition 10.

Definition 10. The defuzzification of and can be given as
The illustration of fuzzification and defuzzification processes based on Definitions 9 and 10, respectively, is illustrated in Figure 5 with the alpha value being 0.5.
Figure 5 shows the fuzzy B-spline surface after fuzzification (Figure 5(a)) and defuzzification of B-spline surface (Figure 5(b)). From Figure 5(b), the fuzzification process was applied by means of alpha-cut operation with the value of alpha as 0.5. Finally, the defuzzification of B-spline surface is equal to crisp B-spline surface because the left and right interval of fuzzy control points is equal.
The constructed fuzzy B-spline surface along with defuzzification process has its merits. The advantage includes the effectiveness model which can be used to model either the crisp data (exact data) or fuzzy data (with various -cuts) compared to the crisp model which can be used to model the crisp data only but not for fuzzy data.

Figure 5: Fuzzy B-spline surface after fuzzification and defuzzification processes.

5. Lakebed Modeling

This section illustrates an application example of fuzzy B-spline surface model proposed in this paper. The proposed model is used to generate the lakebed by using collected data points which is exposed to various kinds of errors.

These errors which occurred during data point retrieval include the wavy water surface condition which gives the uncertain data reading. Therefore, every set of data points in modeling underwater ground surface has the error of accuracy to a certain extent. Figure 6 shows the scenario.

Figure 6: The illustration of collecting the depth of lake.

Figure 6 shows the process of getting uncertain depth lake data which has been taken by echo sounder. We can clearly comprehend the reason for modeling the uncertain data indicating the depth of lake with fuzzy B-spline surface model. This uncertainty in the data exists for -elements (depth).

Fuzzy number concept and fuzzy relation definition are utilized to define this uncertain data. We represent the following algorithm to illustrate the steps to be modeled with fuzzy B-spline from fuzzification to defuzzification processes.

Algorithm 11.
Step  1. Define the uncertain data of lakebed by using Definition 6.
Step  2. Blend the FDPs of lakebed together with B-spline surface function which is given as which is in the form of bicubic surface.
Note that the set of FDPs of lakebed of Kenyir Lake has 64 data points where with .
Step  3. Apply alpha-cut operation as fuzzification process (Definition 9) towards (10) with the alpha value being 0.5. Consider where are the fuzzified data points of lakebed after alpha-cut operation was applied with
Step  4. Use Definition 10 to defuzzify B-spline surface model of lakebed which is where are defuzzified data points of lakebed with
The result of Algorithm 11 can be illustrated in Figure 7.

Figure 7: The process of modeling fuzzy data of Kenyir Lake lakebed.

Figure 7 shows the processes of defining the uncertain data of lakebed which is then modelled by using fuzzy B-spline surface. The fuzzy data points of lakebed can be defined as the fuzzy control points because we used the approximation method to create the fuzzy B-spline surface which used fuzzy data points as fuzzy control points. After achieving fuzzy B-spline surface modeling, we then apply the fuzzification process which utilized the alpha-cut operation of triangular fuzzy number with the alpha value as 0.5. Then, we defuzzify the fuzzy lake surface based on the definition of defuzzification. By setting the new alpha values that are 0.2 and 0.9, we may obtain the result of the defuzzified lakebed model as illustrated in Figures 8 and 9, respectively.

Figure 8: The defuzzification process of fuzzy data of a lakebed for Kenyir Lake with .
Figure 9: The defuzzification process of fuzzy data of a lakebed for Kenyir Lake with .

In order to investigate the effectiveness of the output of lakebed, we find the errors between the defuzzification data points and crisp data points of lakebed which can be given through (15). Therefore, the error between those data can be illustrated by Figures 10, 11, and 12 for different alpha values as follows:

Figure 10: The error between defuzzified and crisp data points of lakebed with .
Figure 11: The error between defuzzified and crisp data points of lakebed with .
Figure 12: The error between defuzzified and crisp data points of lakebed with .

Figure 10 until Figure 12 shows the errors between defuzzified and crisp data points of lakebed in order to see the effectiveness of proposed model in modeling uncertain data of lakebed. The average percentage of the errors is 0.0413623 m, 0.0661796 m, and 0.00827245 m, respectively. These errors are acceptable for terrain modeling.

6. Conclusion

In this paper, we proposed a new paradigm in modeling the uncertain data by using the hybrid method between fuzzy set theory and B-spline function surface function. This model has an upper hand in dealing with uncertain problem compared to the existing model which only can be used in dealing with and modeling crisp data.

The fuzzification and defuzzification processes were also elucidated exclusively for fuzzy B-spline surface model. For fuzzification process, the alpha-cut operation was applied which is in the form of triangular functions. This fuzzification process was applied to obtain the fuzzy interval of fuzzy data points where the crisp fuzzy solution is in this interval. It is then followed by the defuzzification process to find crisp B-spline surface which focused on the defuzzification of fuzzy data points.

Finally, to identify the effectiveness of fuzzy B-spline surface model, this model was applied to modeling of the uncertain data of lakebed. The process of defining, fuzzification, and defuzzification of uncertain data of lakebed can be represented by an algorithm which is applicable in dealing with various fuzzy data. The errors produced in the case of lakebed modeling indicate that it can be used for terrain modeling.

This work can expand further to solve Hermite data problems which occur in designing aesthetic splines [24]. One may create a system to optimize and propose a suitable alpha-cut which satisfies given Hermite data which may facilitate the designers in creating aesthetic shapes efficiently.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.


The authors would like to acknowledge Research Management and Innovation Centre (RMIC) of Universiti Malaysia Terengganu and Ministry of Higher Education (MOHE), Malaysia, for their funding (FRGS: 59244) and providing the facilities to conduct this research. They further acknowledge anonymous referees for their constructive comments which improved the readability of this paper.


  1. K. T. Miura, R. Shirahata, S. Agari, S. Usuki, and R. U. Gobithaasan, “Variational formulation of the log aesthetic surface and development of discrete surface filters,” Computer Aided Design and Application, vol. 9, no. 6, pp. 901–914, 2012. View at Publisher · View at Google Scholar
  2. L. A. Zadeh, “Fuzzy sets,” Information and Computation, vol. 8, pp. 338–353, 1965. View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet
  3. D. Dubois and H. Prade, Fuzzy Sets and Systems: Theory and Applications, vol. 144 of Mathematics in Science and Engineering, Academic Press, New York, NY, USA, 1980. View at Zentralblatt MATH · View at MathSciNet
  4. H. J. Zimmermann, Fuzzy Set Theory and Its Applications, Kluwer Academic, New York, NY, USA, 1985.
  5. G. J. Klir and B. Yuan, Fuzzy Sets and Fuzzy Logic: Theory and Application, Prentice Hall, New York, NY, USA, 1995. View at Zentralblatt MATH · View at MathSciNet
  6. D. F. Rogers, An Introduction to NURBS: With Historical Perspective, Academic Press, New York, NY, USA, 2001.
  7. G. Farin, Curves and Surfaces for CAGD: A Practical Guide, Academic Press, New York, NY, USA, 5th edition, 2002.
  8. F. Yamaguchi, Curves and Surfaces in Computer Aided Geometric Design, Springer, Berlin, Germany, 1988. View at MathSciNet
  9. B. Aggarwal, B-spline finite elements for plane elasticity problems [M.S. thesis], Texas A&M University, 2006.
  10. C.-C. Leung, P. C.-K. Kwok, K.-Y. Zee, and F. H.-Y. Chan, “B-spline interpolation for bend intra-oral radiographs,” Computers in Biology and Medicine, vol. 37, no. 11, pp. 1565–1571, 2007. View at Publisher · View at Google Scholar · View at Scopus
  11. H.-B. Jung, “An interpolation method of b-spline surface for hull form design,” International Journal of Naval Architecture and Ocean Engineering, vol. 2, no. 4, pp. 195–199, 2010. View at Publisher · View at Google Scholar · View at Scopus
  12. T. Maekawa, T. Noda, S. Tamura, T. Ozaki, and K.-I. Machida, “Curvature continuous path generation for autonomous vehicle using B-spline curves,” CAD Computer Aided Design, vol. 42, no. 4, pp. 350–359, 2010. View at Publisher · View at Google Scholar · View at Scopus
  13. S. Okaniwa, A. Nasri, L. Hongwei, A. Abbas, Y. Kineri, and T. Maekawa, “Uniform B-spline curve interpolation with prescribed tangent and curvature vectors,” IEEE Transactions on Visualization and Computer Graphics, vol. 18, no. 9, pp. 1474–1487, 2012. View at Publisher · View at Google Scholar
  14. H. Pungotra, G. K. Knopf, and R. Canas, “An alternative methodology to represent B-spline surface for applications in virtual reality environment,” Computer-Aided Design and Applications, vol. 10, no. 4, pp. 711–726, 2013. View at Publisher · View at Google Scholar
  15. A. A. Nada, “Use of B-spline surface to model large-deformation continuum plates: procedure and applications,” Nonlinear Dynamics, vol. 72, no. 1-2, pp. 243–263, 2013. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet
  16. G. Gallo, M. Spagnuolo, and S. Spinello, “Fuzzy B-splines: a surface model encapsulating uncertainty,” Graphical Models, vol. 62, no. 1, pp. 40–55, 2000. View at Publisher · View at Google Scholar · View at Scopus
  17. G. Gallo, M. Spagnuolo, and S. Spinello, “Rainfall estimation from sparse data with fuzzy B-splines,” Journal of Geographic Information and Decision Analysis, vol. 2, pp. 194–203, 1998. View at Google Scholar
  18. A. M. Anile, B. Falcidieno, G. Gallo, M. Spagnuolo, and S. Spinello, “Modeling uncertain data with fuzzy B-splines,” Fuzzy Sets and Systems, vol. 113, no. 3, pp. 397–410, 2000. View at Publisher · View at Google Scholar · View at MathSciNet
  19. T. Yanhua and L. Hongxing, “Faired MISO B-spline fuzzy systems and its applications,” Mathematical Problems in Engineering, vol. 2013, Article ID 870595, 9 pages, 2013. View at Publisher · View at Google Scholar · View at MathSciNet
  20. G. J. Klir, U. S. Clair, and B. Yuan, Fuzzy Set Theory: Foundation and Application, Prentice Hall, Upper Saddle River, NJ, USA, 1997. View at MathSciNet
  21. R. Zakaria and A. B. Wahab, “Fuzzy B-spline modeling of uncertainty data,” Applied Mathematical Sciences, vol. 6, pp. 6971–6991, 2012. View at Google Scholar
  22. A. F. Wahab, J. M. Ali, A. A. Majid, and A. O. M. Tap, “Fuzzy set in geometric modeling,” in Proceedings of the International Conference on Computer Graphics, Imaging and Visualization (CGIV '04), pp. 227–232, July 2004. View at Scopus
  23. N. A. A. Karim, A. F. Wahab, R. U. Gobithaasan, and R. Zakaria, “Model of fuzzy B-spline interpolation for fuzzy data,” Far East Journal of Mathematical Sciences, vol. 72, pp. 269–280, 2013. View at Google Scholar
  24. K. T. Miura, D. Shibuya, R. U. Gobithaasan, and S. Usuki, “Designing log-aesthetic splines with G2 continuity,” Computer-Aided Design and Applications, vol. 10, no. 6, pp. 1021–1032, 2013. View at Publisher · View at Google Scholar