Complex Uncertainty of Surface Data Modeling via the Type-2 Fuzzy B-Spline Model

This paper discusses the construction of a type-2 fuzzy B-spline model to model complex uncertainty of surface data. To construct this model, the type-2 fuzzy set theory, which includes type-2 fuzzy number concepts and type-2 fuzzy relation, is used to define the complex uncertainty of surface data in type-2 fuzzy data/control points. These type-2 fuzzy data/control points are blended with the B-spline surface function to produce the proposed model, which can be visualized and analyzed further. Various processes, namely fuzzification, type-reduction and defuzzification are defined to achieve a crisp, type-2 fuzzy B-spline surface, representing uncertainty complex surface data. This paper ends with a numerical example of terrain modeling, which shows the effectiveness of handling the uncertainty complex data.


Introduction
Data points are the representation of the visible object in a digital system. The shape modification of items is carried out by tweaking the data points as the designer desires. However, the data points representing an object are not necessarily precise due to their characteristics, which are uncertain, leading to complexity. The complexities that are found in the data are due to various possible reasons, e.g., inappropriate analysis, human perception and logical assumption [1,2], limitations of tool accuracy, and the nature of data collection itself [3,4] with different types of uncertainty being based on BIPM (Bureau International des Poids et Mesures) concepts. The specific meaning of 'complex' in complex uncertainty data is the stack of the uncertainty of two arguments of the collected data points. It is impossible to accurately model the complex uncertainty data using an appropriate standard curve or surface function, such as the B-spline function, unless we formulate a new definition of the B-spline function with the complex uncertainty meaning.
To make the complex uncertainty data useful in modeling, we define complex uncertainty data in this paper. The complex uncertainty data can be explained by using type-2 fuzzy number (T2FN) concepts [5][6][7], which are defined by type-2 fuzzy set theory (T2FST) [8][9][10], especially interval type-2 fuzzy number (IT2FN) concepts [5,11]. The IT2FN concept is implemented to define the complex uncertainty data of real numbers subtly transformed into type-2 fuzzy data points (T2FDP), in line with the definition of type-2 fuzzy relation (T2FR). This is then followed by representing T2FDP as a type-2 fuzzy control point (T2FCP), which can be used for modeling the B-spline surface function [12][13][14]. Hence, we end up with a standard definition denoted as the type-2 fuzzy B-spline surface.
The interval type-2 fuzzy set (IT2FS) concept has been widely used to model higherorder uncertainty, which has been proven to be more suitable compared to the interval type-

Previous Work
The T2FST is useful in defining the uncertainty complex data to create a type-2 fuzzy curve and surface using the B-spline curve and surface functions. Regarding the designers' requirements in modeling the complex uncertainty data, the designers cannot choose and decide which data points are essential among the collectives of complex uncertainty data in modeling curves and surfaces. Therefore, we need the T2FST to define the complex uncertainty data and then model them through the B-spline curve and surface functions. This approach will make the complex uncertainty data become T2FDP/T2FCP, which can be modeled after integrating with the B-spline curve and surface functions.
Many types of research are carried out using the type-1 fuzzy set theory to model surface in dealing with uncertainty issues. Examples include the surface model proposed by Gallo et al. [21,22] of Mount Etna, Zakaria et al. [4] of lLakebed's modeling of Kenyir Lake, and Sarwar and Akram [23] proposed the fuzzy tensor product of Bezier surface. Note that the uncertainty level of these models is at level one. Suppose that the level of uncertainty increases due to specific errors, in particular the uncertainty at level two (complex uncertainty data). In that case, the current proposed approach for level one

Definition 3 ([25]).
The following definition is T2FDP which is defined using the previous definitions as stated before [25]. Definition 4. Let P = {x|x type − 2 fuzzy point} and ↔ ↔ P = {P i |P i data point} be the set of the type-2 fuzzy data point with P i ∈ P ⊂ X, where X is a universal set and µ P (P i ) : P → [0, 1] is the membership function defined as µ P (P i ) = 1 and formulated as , µ P ( f i ) and µ P (c i ) are right-right, right, left-right membership grade values, which can be written as follows: For every i, After T2FDP has been defined, the next procedure is the fuzzification process that applies the alpha-cut of IT2FN [25]. This definition is determined based on the fuzzification process of type-1 fuzzy data points discussed in [4]. Therefore, the description of the After T2FDP has been defined, the next procedure is the fuzzification process that applies the alpha-cut of IT2FN [25]. This definition is determined based on the fuzzification process of type-1 fuzzy data points discussed in [4]. Therefore, the description of the fuzzification process against T2FDP can be given through Definition 5, as follows. P where i = 0, 1, . . . , n − 1. Then ↔ ↔ P i α is the alpha-cut operation of T2FDPs with i = 0, 1, 2, . . . , n which is given as follows.

Definition 5. Let
This definition is illustrated in Figure 3. After T2FDP has been defined, the next procedure is the fuzzification process that applies the alpha-cut of IT2FN [25]. This definition is determined based on the fuzzification process of type-1 fuzzy data points discussed in [4]. Therefore, the description of the fuzzification process against T2FDP can be given through Definition 5, as follows.
; ; This definition is illustrated in Figure 3.   Figure 3 shows the implication of the alpha-cut operation against T2FDP, which is the fuzzification process with the specific value of alpha (membership value). The alpha value of this operation is 0.5. If the alpha value increases to one, then the crisp data point is obtained. This can be illustrated by Figure 4 as follows.
After performing the fuzzification process, then the following process is type-reduction. Type-reduction is defined, then used against T2FDP to allow the defuzzification of the type-1 fuzzy set. The type-reduction is an approach to simplify type-2 defuzzification. The first type-reduction had been proposed by Nie and Tan [27]. On the other hand, type-reduction on discretized interval type-2 fuzzy sets have been discussed in [28]. The consistent linear and quadratic type-reduction methods have been introduced by [29].
The type-2 defuzzification has been already discussed, such as the Karnik-Mendel algorithm [30]. Subsequently, type-2 defuzzification with explicit models of the uncertainty was proposed in [31]. Moreover, the various methods of type-2 defuzzification have different mathematical properties discussed by [31].
In this paper, the proposed methods for type-reduction and type-2 defuzzification against T2FDP are based on the centroid method based on Definition 6 and Definition 7, respectively. Figure 3 shows the implication of the alpha-cut operation against T2FDP, which is the fuzzification process with the specific value of alpha (membership value). The alpha value of this operation is 0.5. If the alpha value increases to one, then the crisp data point is obtained. This can be illustrated by Figure 4 as follows. After performing the fuzzification process, then the following process is type-reduction. Type-reduction is defined, then used against T2FDP to allow the defuzzification of the type-1 fuzzy set. The type-reduction is an approach to simplify type-2 defuzzification. The first type-reduction had been proposed by Nie and Tan [27]. On the other hand, typereduction on discretized interval type-2 fuzzy sets have been discussed in [28]. The consistent linear and quadratic type-reduction methods have been introduced by [29].
The type-2 defuzzification has been already discussed, such as the Karnik-Mendel algorithm [30]. Subsequently, type-2 defuzzification with explicit models of the uncertainty was proposed in [31]. Moreover, the various methods of type-2 defuzzification have different mathematical properties discussed by [31].
In this paper, the proposed methods for type-reduction and type-2 defuzzification against T2FDP are based on the centroid method based on Definition 6 and Definition 7, respectively.
; 0,1, 2,..., P i α which is represented as ↔ P α can be defined as follows: where P i is crisp data points and ↔ P ← i α and ↔ P → i α are left and right fuzzified type-reduction T2FDP respectively with their formulation given by the following: Through the implementation of Definition 6, the fuzzified type-reduction T2FDP is obtained. Then, the next procedure to get the crisp T2FDP is the defuzzification process of type-1. The defuzzification process has been defined in [4,32,33] and described as follows.

Definition 7. Let
↔ P i α be the fuzzified type-reduction T2FDP with i = 0, 1, 2, . . . , n. Then, P i α is the defuzzification process of where each P i α can be formalized as: The process of obtaining T2FDP defuzzified from defining T2FDP, including the fuzzification, type-reduction and defuzzification processes can be summarized and illustrated in the following Figure 5.
The process of obtaining T2FDP defuzzified from defining T2FDP, including the fuzzification, type-reduction and defuzzification processes can be summarized and illustrated in the following Figure 5.

Results: Type-2 Fuzzy B-Spline Model
This section will discussing type-2 fuzzy B-spline model which specifically used the interpolation method. When creating a type-2 fuzzy curve and surface, the T2FDP are integrated into the B-spline curve and surface function where the end result is known as a type-2 fuzzy B-spline curve and surface. This type-2 fuzzy model meets the processes of fuzzification, type-reduction and defuzzification to obtain the crisp type-2 fuzzy curve and surface solution (single curve and surface solution).
The construction of the type-2 fuzzy B-spline model is based on the studies carried out by Zakaria et al. [4], Zakaria and Wahab [32], Wahab et al. [33], and Karim et al. [34]. These studies discussed the construction of the type-1 fuzzy B-spline model. The type-2 fuzzy interpolation B-spline model [25] is defined in Definition 8 and illustrated by Figure 6.

Definition 8. Let
↔ ↔ D i ∈ R be a list of T2FDP with i = 0, 1, 2, . . . , n, then the type-2 fuzzy interpolation B-spline curve (T2FIBsC) can be defined as follows: where ↔ ↔ P i are T2FCP, N i,h (t) is a basic function of B-spline and t is crisp knot sequences t 1 , t 2 , . . . , t m=d+n+1 in which d represents the degree of B-spline function and n represents the numbers of control points.

Definition 8. Let
∈   i D R be a list of T2FDP with = 0,1, 2,..., i n, then the type-2 fuzzy interpolation B-spline curve (T2FIBsC) can be defined as follows: where   i P are T2FCP,

Definition 9. Let
↔ ↔ D a,b be given as a set of T2FDP with a = 0, 1, . . . , m and b = 0, 1, . . . , n which has degree p and q. Then, the type-2 fuzzy interpolation B-spline surface (T2FIBsS) model with degree p and q can be given as follows: where ↔ ↔ P k,l is the set of T2FCP which is the unknown value.
Equation (10) can be rewritten as with ↔ ↔ Mathematics 2021, 9, 1054 9 of 18 Therefore, Equation (12) is resolved first which used ↔ ↔ D a,b before N k,p (s a ). It is followed by solving ↔ ↔ P k,l using Equation (10). Therefore, we obtain T2FCP values that allow the surface to interpolate the T2FDP. Thus, the illustration of the T2FIBsS model can be illustrated in Figure 7 as follows.
Therefore, Equation (12) (10). Therefore, we obtain T2FCP values that allow the surface to interpolate the T2FDP. Thus, the illustration of the T2FIBsS model can be illustrated in Figure 7 as follows.  Figure 7 shows that the T2FIBsS model in bicubic shapes was constructed based on the T2FCP net obtained using Equation (10). The T2FDP net of T2FIBsS of Figure 7a illustrates the first impression of the against the type-2 fuzzy surface.
The next process is the fuzzification process. This fuzzification is performed against the T2FDPs before being integrated with the basic B-spline surface function. Therefore, the fuzzification of T2FIBsS can be defined as the following equation  Figure 7 shows that the T2FIBsS model in bicubic shapes was constructed based on the T2FCP net obtained using Equation (10). The T2FDP net of T2FIBsS of Figure 7a illustrates the first impression of the against the type-2 fuzzy surface.
The next process is the fuzzification process. This fuzzification is performed against the T2FDPs before being integrated with the basic B-spline surface function. Therefore, the fuzzification of T2FIBsS can be defined as the following equation where α j is the value of alpha-cut operation of the type-2 triangular fuzzy number represents the fuzzification process with α ∈ (0, 1], j = 1, 2, . . . , y. Therefore, based on Figure 7, the illustration after fuzzification process can be shown in Figure 8. After the fuzzification process has been applied, the next step is the type-reduction process. Therefore, the reduction of fuzzified T2FIBsS can be formalized as follows, based on Definition 6.
The defuzzification of type-1 fuzzy then follows it to obtain the final result as the crisp T2FIBsS. Therefore, the defuzzification of type-reduced fuzzified T2FIBsS is given by the following equation.
where α j is the value of alpha-cut operation of the type-2 triangular fuzzy number represents the fuzzification process with α ∈ (0,1] , = 1, 2,..., j y. Therefore, based on Figure  7, the illustration after fuzzification process can be shown in Figure 8. After the fuzzification process has been applied, the next step is the type-reduction process. Therefore, the reduction of fuzzified T2FIBsS can be formalized as follows, based on Definition 6. Then, the illustration of this type-reduction process can be illustrated in Figure 9 as follows. The defuzzification of type-1 fuzzy then follows it to obtain the final result as the crisp T2FIBsS. Therefore, the defuzzification of type-reduced fuzzified T2FIBsS is given by the following equation. BsSi s t N s N t P D (15) The illustration of the defuzzification process against the type-reduced fuzzified T2FIBsS is shown in Figure 10. Then, the illustration of this type-reduction process can be illustrated in Figure 9 as follows. After the fuzzification process has been applied, the next step is the type-reduction process. Therefore, the reduction of fuzzified T2FIBsS can be formalized as follows, based on Definition 6. BsSi s t N s N t P D (14) Then, the illustration of this type-reduction process can be illustrated in Figure 9 as follows. The defuzzification of type-1 fuzzy then follows it to obtain the final result as the crisp T2FIBsS. Therefore, the defuzzification of type-reduced fuzzified T2FIBsS is given by the following equation.
The illustration of the defuzzification process against the type-reduced fuzzified T2FIBsS is shown in Figure 10. The illustration of the defuzzification process against the type-reduced fuzzified T2FIBsS is shown in Figure 10.  Figure 10 shows that the defuzzification-reduced T2FIBsS model along with crisp interpolation B-spline surface. Both surfaces were built by finding the control points that interpolates data. The crisp data points were marked by red and the defuzzfication-reduced T2FDPs were marked by cyan.

Application: Seabed Modeling
This section discusses the practical application of the T2FIBsS model of seabed modeling. In the seabed modeling, multiple uncertainties or errors occur due to the nature of collecting data points, i.e., the wavy water surface and the collector's uncertain perception and truth level. Thus, the errors of the data collected are bounded by two uncertainties. This scenario is illustrated in Figure 11, which is also the extension of uncertainty data for  Figure 10 shows that the defuzzification-reduced T2FIBsS model along with crisp interpolation B-spline surface. Both surfaces were built by finding the control points that interpolates data. The crisp data points were marked by red and the defuzzfication-reduced T2FDPs were marked by cyan.

Application: Seabed Modeling
This section discusses the practical application of the T2FIBsS model of seabed modeling. In the seabed modeling, multiple uncertainties or errors occur due to the nature of collecting data points, i.e., the wavy water surface and the collector's uncertain perception and truth level. Thus, the errors of the data collected are bounded by two uncertainties. This scenario is illustrated in Figure 11, which is also the extension of uncertainty data for lakebed modeling [4].  Figure 10 shows that the defuzzification-reduced T2FIBsS model along with crisp interpolation B-spline surface. Both surfaces were built by finding the control points that interpolates data. The crisp data points were marked by red and the defuzzfication-reduced T2FDPs were marked by cyan.

Application: Seabed Modeling
This section discusses the practical application of the T2FIBsS model of seabed modeling. In the seabed modeling, multiple uncertainties or errors occur due to the nature of collecting data points, i.e., the wavy water surface and the collector's uncertain perception and truth level. Thus, the errors of the data collected are bounded by two uncertainties. This scenario is illustrated in Figure 11, which is also the extension of uncertainty data for lakebed modeling [4]. Figure 11. Illustration of procedure in taking depth data point (in meter) of the seabed, which consists of uncertainty complex data. Figure 11 shows that obtaining the depth sea by means of an echo sounder where the data points obtained have the complex uncertainty properties. Therefore, it necessary to use the T2FIBsS method to model those complex uncertainty data after the data points are defined through the T2FST.
The following algorithm shows a step-by-step process of defining uncertainty complex data until the final stage of defuzzification-reduced T2FDP B-spline surface interpolation function.
The result of Algorithm 1 is illustrated in Figures 12-14 as follows. Figure 11. Illustration of procedure in taking depth data point (in meter) of the seabed, which consists of uncertainty complex data. Figure 11 shows that obtaining the depth sea by means of an echo sounder where the data points obtained have the complex uncertainty properties. Therefore, it necessary to use the T2FIBsS method to model those complex uncertainty data after the data points are defined through the T2FST.
The following algorithm shows a step-by-step process of defining uncertainty complex data until the final stage of defuzzification-reduced T2FDP B-spline surface interpolation function.
The result of Algorithm 1 is illustrated in Figures 12-14 as follows. Figure 12 shows the implementation of T2FIBsS against the seabed data modeling. In getting the final output named as crisp T2FIBsS, the fuzzifying process involves the output of T2FIBsS for seabed surface from the complex uncertainty data. Then, the fuzzification process is applied to shorten the T2FIBsS interval with the alpha value of 0.2 by using the alpha-cut of a triangular type-2 fuzzy number. After the fuzzification process, the following process applied is the type-reduction process, which reduces the type-2 fuzzy to become type-1 fuzzy to allow the defuzzfication process. Later, the defuzzification process is applied to obtain the crisp type-2 fuzzy solution that shows the distinct surface.
For Figures 13 and 14, getting the crisp T2FIBsS is the same as the processes in    Figure 12 shows the implementation of T2FIBsS against the seabed data modeling. In getting the final output named as crisp T2FIBsS, the fuzzifying process involves the output of T2FIBsS for seabed surface from the complex uncertainty data. Then, the fuzzification process is applied to shorten the T2FIBsS interval with the alpha value of 0.2 by using the alpha-cut of a triangular type-2 fuzzy number. After the fuzzification process, the following process applied is the type-reduction process, which reduces the type-2 fuzzy to Step 1: Define the uncertainty complex data of seabed by using Definition 4.
Step 3: Find the values of N k,p (s a ) after Equation (12) has been solved.
Step 4: Find the points of ↔ ↔ P k,l through Equation (10) by using Equation (12).
Step 9: Find and plot the error between defuzzification-reduced T2FDPs and crisp data points of seabed depth data using the following equation:

Discussion and Conclusions
This paper proposed a new method for defining complex uncertainty data, and modeled it on the hybrid model constructed called T2FIBsS. The complex uncertainty data was defined through IT2FN concepts obtaining T2FDPs and integrated with the B-spline surface function to produce the T2FIBsS model. This model has an advantage in modeling complex uncertainty data, as shown in seabed modeling. Meanwhile, the error and the percentage error have been calculated between defuzzification-reduced T2FDPs and crisp data points demonstrating the feasibility of the proposed model. This developed model as T2FIBsS also deals with uncertainty modeling, which type-1 fuzzy modeling can do when the T2FIBsS first is reduced to the type-1 fuzzy model. However, the type-1 fuzzy model, which has been discussed in the literature, cannot be used in defining and modeling the complex uncertainty data. Therefore, T2FIBsS can be used to define and model the uncertainty and complex uncertainty compared to the type-1 fuzzy model, which only can be used to define and model the uncertainty data, but not applicable for complex uncertainty data.
The limitation of this study is the properties of the data that we want to model. The data that we obtained is secondary data, which had been filtered from raw data. Next is the surface function, which generates a desired surface based on the properties of the data. This function generates the surface entirety, rather than as a suitable surface patch. This surface that is generated entirety will make the surface sometimes unreasonable due to the data point position.
This research can be improved in the future by taking into account the complex uncertainty data, which has more uncertainty based on decision making or perceptions. This research can be expanded by using a complicated surface function, such as the Non-Uniform Rational B-spline (NURBS) surface, for modeling parts. The NURBS surface function has the advantage of local control on the surface other than the control points and knots, but the added value as weight in the NURBS formulation can modify the surface locally.
Author Contributions: R.Z., A.F.W. and M.I.E.Z. conceived of the presented the idea, conceptualization and methodology. I.I., gives the seabed data based on GIS from ArcGIS. All authors have read and agreed to the published version of the manuscript.