Abstract
From the observed datasets, we should be able to produce curve surfaces that have the same characteristics as the original datasets. For instance, if the given data are positive, then the resulting curve or surface must be positive on entire given intervals, i.e., everywhere. In this study, a new partial blended rational bi-quartic spline with C1 continuity is constructed through the partially blended scheme. This rational spline is defined on four corners of the rectangular meshes. The sufficient condition for the positivity of rational bi-quartic spline is derived on four boundary curve networks. There are eight free parameters that can be used for shape modification. The first-order partial derivatives are estimated by using numerical techniques. We also show that the proposed scheme is local quadratic reproducing such that it can exactly reproduce the quadratic surface. We test the proposed scheme to interpolate various types of positive surface data. Based on statistical indicators such as the root mean square error (RMSE) and coefficient of determination (R2), we found that the proposed scheme is on par with some established schemes. In fact, it requires less CPU times (in seconds) to generate the interpolating surface on rectangular meshes. Furthermore, by combining the statistical indicators’ result and graphically visualizing the test functions, the proposed method has the capability to reconstruct very comparable smoothing interpolating positive surfaces compared to some existing schemes. This finding is significant in producing a better interpolating surface for computer graphics applications since the proposed scheme has a smaller error compared with existing schemes.
Introduction
Computer-aided geometric design (CAGD) provides the mathematical basis when dealing with geometric datasets. The term “CAGD” was coined by Barnhill and Riesenfeld in 1974 [1–5]. One of the main tools in CAGD is a spline function. The spline function was introduced by Schoenberg in 1967 for application in statistics. Since then, the spline has been used extensively in shape designing especially in car industries, in the design of airplane fuselages and wings at Boeing Company, as well as in the development of cartoon movies. General Motors developed its first computer-aided design/computer-aided manufacturing (CAD/CAM) system called DAC-I by utilizing the spline curve and surface techniques initiated by de Boor and Gordon [5]. Furthermore, in some applications, such as data visualization, there is a need to reconstruct the curve or surface for the data collected from the measurement or experimental laboratory. This task can be achieved by using interpolation or approximation schemes.
Shape-preserving interpolation is an important topic in CAGD that can be used to preserve the geometric shape of the dataset such as positivity, convexity, and monotonicity. For instance, given positive data, then the final interpolating curve or surface must be positive everywhere. Any negativity is meaningless and unacceptable. Some examples of positive data are the amount of gas discharge during the experiment (Brodlie and Butt [6], Butt and Brodlie [7], Brodlie et al. [8], Sarfraz et al. [9]). Meanwhile, monotonicity-preserving concerns about preserving the monotone datasets, i.e., either monotonically increasing or monotonically decreasing. For instance, the approximation of couples and quasi-couples in statistics and the empirical option model in finance are always involving monotone datasets [10]. The blood uric acid in patients suffering from gout is also an example of monotone data. Meanwhile, convexity always arises in finance and engineering-based problems such as the optimal control, non-linear programming, and parameter estimation problems.
In general, the standard cubic Hermite polynomial and cubic spline interpolation functions are unable to produce the interpolating curve or surface with shape-preserving properties. Therefore, to cater this weakness, many researchers have proposed several methods for shape-preserving approximation and interpolation. One of the simplest ways is to use rational splines in many different forms such as cubic/cubic, cubic/quadratic, cubic/linear, and quartic/quadratic. Those schemes provided more design parameters which can be used to modify the interpolating curves or surfaces. In the next paragraphs, we provide some literature review on positivity-preserving interpolation schemes.
The modified quadratic Shepard (MQS) method has been used widely for visualization of scattered data. For example, Asim et al. [11, 12] improved the original MQS method for positivity-preserving on scattered datasets. Their strategy is to reduce the deviation from the original shape while at the same time preserving positivity in the sense of least square. Brodlie et al. [8] constructed the MQS method to interpolate scattered data of any dimensionality. Their scheme preserves the positivity of the data for curves or surfaces by forcing the quadratic basis functions to be positive. They also extended the method to handle modeling other types of constraints such as the lower bound of 0 and upper bound of 1 and generalized the constraint to any arbitrary functions as lower and upper bounds. But from the numerical results, they indicated that their schemes have 10% error bound associated with the data values. Wu et al. [13] discussed about positivity-preserving for curve and surface approximation and interpolation by using compactly supported radial basis functions (CSRBFs). To preserve the positivity of data, the optimization problem needs to be solved. In [14, 15], the multi-quadric (MQ) function has been used as the quasi-interpolation operator and the CSRBFs are used to construct the interpolation function that preserves the monotonicity and convexity of planar datasets.
Most shape-preserving surface interpolation methods are constructed by extending the work of Casciola and Romani [16]. Nowadays, their schemes are known as partially blended rational bi-cubic functions. Brodlie and Butt [6] constructed the piecewise cubic Hermite interpolant with continuity to preserve positivity, monotonicity, and convexity, respectively. Their schemes required one or two extra knots to be inserted in the interval where shape violation exists without needs for derivative modification. However, this technique will increase the total number of data points, and indeed, the computation times will be increased.
Brodlie et al. [8] extended the Butt and Brodlie [7] idea to construct the positivity-preserving scheme of positive surface data over a rectangular grid. The sufficient conditions on first partial derivatives and twist values are derived where the obtained values will be projected to the valid interval through the efficient knot insertion algorithm. However, their scheme has increased unnecessary computation time.
Schmidt and Hess [17, 18] derived the necessary and sufficient conditions for the positivity of rational quadratic splines and for the positivity of cubic spline interpolation on intervals. Schmidt [19] described the positive, monotone, and -convex surface interpolation for data arranged on rectangular grids. The rational quadratic spline has been used, and the necessary conditions for positivity, monotonicity, and -convex surface are derived. However, no numerical results were presented. Hussain et al. [20] constructed shape-preserving surface schemes for 3D positive and convex surface data. The rational bi-quadratic function with eight parameters was used to construct the positive and monotone rational interpolant. But their scheme was unable to produce smooth surfaces as well as no free parameters to alter the final resulting interpolating surfaces. Chan and Ong [21] described a local scheme for range-restricted scattered data interpolation by using cubic triangular Bézier patches. The piecewise interpolating surface is obtained through a convex combination of three cubic Bézier triangular patches. The sufficient conditions for the non-negativity of a cubic Bézier triangle on the Bézier ordinates were derived with respect to a lower bound. The gradients (first-order partial derivatives) at the data sites are modified if necessary, to ensure that the non-negativity conditions are fulfilled. Luo and Peng [22] described the rational spline as a piecewise rational convex combination of three cubic Bézier triangular patches sharing the same boundary Bézier ordinates. The sufficient conditions for non-negativity were derived on the boundary Bézier ordinates of the adjacent triangle and the normal derivatives at the data points. If in any triangular patch non-negativity was lost, then the gradients at the data points and normal derivatives at the edge knots will be modified to ensure positivity is preserved. Indeed, their main scheme also requires the modification of the first-order partial derivative and the normal derivatives if shape violation is found.
Piah et al. [23] improved the lower bound of Bézier ordinates obtained in the work of Chan and Ong [21] by proposing an alternative scheme which is simpler and has more relaxed conditions of positivity. The new lower bound for Bézier ordinates can become infinitely smaller compared to the lower bound of Bézier ordinates in [21]. Some numerical results were presented to show the capability of the alternative scheme. The main drawback of their scheme is that the actual first-order partial derivative needs to be adjusted if positivity on each triangular patch is violated. Hussain and Hussain [24] constructed positive scattered data interpolation by using a rational cubic spline defined by the side-vertex scheme.
Saaban et al. [25] constructed the positivity-preserving scheme by using quintic triangular Bézier patches. The surface is constructed by using a convex combination comprising of three local quintic triangular Bézier patches. They implemented their proposed method by using rainfall data collected at various meteorology stations in West (Peninsular) Malaysia.
Hussain et al. [26, 27] extended the application of rational quartic spline given by Wang and Tan [28] to preserve positivity, constrained data, and convexity. The data-dependent sufficient conditions for the rational interpolant to satisfy the shape-preserving properties are derived. It was proved by Harim et al. [29] that their sufficient conditions may not successfully produce positive and convex interpolating curves on entire given intervals.
Liu et al. [30] described a new bivariate rational quartic spline (quartic/quadratic) for positivity- and monotonicity-preserving interpolation. The rational bi-quartic spline has four parameters, and the sufficient conditions for positivity and monotonicity are derived on all parameters. Thus, there is no free parameter for shape modification. They have compared the performance of their schemes against that of Hussain and Sarfraz’s [31] for positivity-preserving interpolation. No error measurement is given to show the superiority of the proposed scheme. Han [32, 33] discussed the rational quartic spline with a quadratic denominator. Harim et al. [29, 34, 35] constructed a new rational quartic spline (quartic/quadratic) with three parameters for shape modification. Karim et al. [31, 36, 37] discussed the positivity-preserving interpolation and constrained surface modeling by using rational bi-cubic spline interpolation with 12 parameters. Qin et al. [38] proposed a new method to derive the sufficient condition for the positivity of rational bi-cubic spline interpolation. However, their method requires extra computation time.
In this section, we have given a comprehensive review on the existing schemes of positivity-preserving interpolation. Table 1 summarizes the comparison of various shape-preserving interpolation schemes related to our studies.
TABLE 1
| Bil | Scheme | Features |
|---|---|---|
| 1 | Gregory [39], Hussain and Ali [37], Hussain and Hussain [40], and Tian et al. [41] | No free parameters for shape refinement. Therefore, we cannot control the final interpolating curves and surfaces. |
| 2 | Hussain et al. [20, 42] | Cannot produce positivity-preserving interpolation on entire given intervals as shown by Harim et al. [29]. |
| 3 | Asim and Brodlie [43], Asim et al. [11], Brodlie and Butt [6], Butt and Brodlie [7], Fiorot and Tabka [44], and Peng et al. [45] | Additional knots are inserted. This will increase the CPU times. |
| 4 | Fristch and Butland [46], Fristch and Carlson [47], and Wang and Tan [28, 48, 49] | Shape-preserving with modification of derivatives. Cannot interpolate first derivatives. |
| 5 | Hussain et al. [20] | The interpolating curve and surface are C0 continuous. Thus, the surface is not smooth and not suitable for shape-preserving interpolation. |
| 6 | Abbas et al. [50, 51] and Karim et al. [36] | RMSE is a bit higher and smaller R2 value. They did not prove whether their schemes are local quadratic polynomial reproducing or not. |
| 7 | Wu et al. [13, 14] | Requires solving an optimization problem that will increase unnecessary CPU times (in seconds) to produce the curves. |
| 8 | Ibraheem et al. [52] | The interpolating surface is not smooth as well as not visually pleasing as seen in their study. |
| 9 | Ashraf et al. [43] | Use the subdivision scheme to preserve the shape of the data. |
Comparison of shape-preserving schemes.
From
Table 1, we have come out with our main objectives of the study. There are three main objectives given as follows:
1) We extend the univariate rational quartic spline with three parameters from Harim et al. [34] to the bivariate rational quartic spline interpolation. We employ the partially blended scheme initiated by Casciola and Romani [16]. There are three free parameters on each boundary of the rectangular meshes. This will give 12 parameters that can be used to control the shape of the interpolating surface.
2) To derive the sufficient condition for the positivity of the rational bi-quartic spline on each boundary, we also derive the quadratic polynomial reproducing.
3) Finally, to compare the performance between the proposed scheme for positivity-preserving interpolation and some established schemes such as Abbas et al. [50] and Karim et al. [36], the root mean square error (RMSE) and coefficient of determination (R2) are used as statistical goodness-of-fit measurements [53].
Furthermore, we identified several advantages of the proposed partially blended rational bi-quartic spline interpolation for positivity-preserving:
1) There are 12 parameters in the description of the rational bi-quartic spline, and eight are free parameters for shape modification. Meanwhile, Hussain and Ali [37], Hussain and Hussain [40], Tian et al. [41], and Sarfraz et al. [9] schemes have no free parameters for shape modification.
2) Unlike existing schemes such as Karim et al. [36], Abbas et al. [50, 51], Hussain and Ali [37], Hussain and Hussain [40], Tian et al. [41], and Sarfraz et al. [9], the proposed scheme is local quadratic reproducing. This is one of the main novelties of the proposed scheme.
3) Based on statistical indicators, i.e., root mean square error (RMSE) and coefficient of determination (R2), the proposed partially blended rational bi-quartic spline is the best compared with the works of Abbas et al. [50] and Karim et al. [36]. In fact, the proposed scheme has a smaller RMSE error and higher R2 compared with Brodlie et al.’s scheme [7].
4) Furthermore, the proposed scheme is direct and not involving any iterative scheme via the subdivision scheme as discussed by Ashraf et al. [43].
The remainder of this paper is organized as follows. Materials and Methods is devoted to materials and methods used in this study. This includes the method to estimate first derivatives using the arithmetic mean method (AMM) and geometric mean method (GMM), a review on the rational quartic spline, the derivation of the rational bi-quartic spline, and local control and quadratic reproducing properties. Meanwhile, Positivity-Preserving Interpolation discusses the construction of positive rational bi-quartic spline interpolation on each boundary curve network. Furthermore, an efficient algorithm is presented in detail. Results and Discussion is devoted to the results obtained and discussion. Conclusion gives the concluding remarks.
Materials and Methods
Derivative Estimation
In most applications, the first derivative needs to be estimated by using some numerical methods. Three common methods are the arithmetic mean method (AMM), geometric mean method (GMM), and harmonic mean method (HMM) [54].
Arithmetic Mean Method
Given scalar datasets where. Let the step size and the gradient , . The first derivatives at the data point , can be estimated as follows.
At both end points and , the first derivatives are [54]Meanwhile, at the interior points, , the values of are estimated using the following formula:Formulas in Equations (1)–(3) can be extended to 3D data as follows.
Given 3D datasets . Let , and . The first partial derivatives are given as follows.
The first partial derivative in the direction is given asfor .
The first partial derivative in the direction is given as for . Meanwhile, the twists (mixed derivative) at the interior points can be estimated by
Geometric Mean Method
The first derivatives also can be calculated by using the GMM. Given 2D data , at the end points and , the first derivatives are given aswith andwith
Meanwhile, at interior points, , the values of are estimated by
The advantages of using this choice are the values of are always positive. This method is suitable for monotonically increasing data.
Rational Quartic Spline Interpolation
Given scalar datasets where and the first derivative, at the respective point. Then, the RQS scheme with three shape control parameters, andon the interval is given by [29, 34, 35]wherewith ,, and .
The RQS scheme defined in Eq. 14 gives the C1 continuity at the knots and satisfies the following conditions:From Eq. 15 and after some derivations, the unknown , and are given as [34]
Partially Blended Rational Bi-Quartic Spline Interpolation
We construct a new rational bi-quartic spline based on the partially blended scheme of Casciola and Romani [16]. There are four curve networks on each edge of the rectangular domain. Figure 1 shows the arrangement of the partially blended rational bi-quartic spline interpolation on a rectangular domain.
FIGURE 1
Let the function be arranged over a rectangular domain . The partition of the intervals and is defined as and Let , , and indicate the data values and the first partial derivatives with respect to and at the node . The partially blended rational bi-quartic function over each rectangle , is defined as follows:withwhereand, and are rational quartic functions defined in Equations (18)–(21) arranged on the edges of the rectangular domain (see Figure 1). They are defined aswithwithwithwithThe partial derivatives and are estimated by using the AMM given in Equations (4)–(10).
Local Control Shape Property Analysis
The proposed rational bi-quartic spline defined by
Eq. 17has local control properties. Some observations can be made as follows:
1) The proposed scheme is symmetrical, i.e., the surface can be constructed either in x and y directions or in y and x directions. This is one of the main properties in CAGD.
2) When the parameters , and , the partially blended rational bi-quartic spline in Eq. 17 reduced to a standard bi-quartic polynomial spline. Usually, in this case, the resulting interpolating surface is not shape-preserving.
3) When or the partially blended rational bi-quartic spline defined by Eq. 17 is reduced to a bilinear surface, i.e., flat surface, since all four curve networks are reduced to a straight line (see [34]):
• Similarly, we can obtain the following observations:
Based on
Equations (22)–
(25), we conclude that the final interpolating surface is a bilinear surface.
Theorem 1. The partially blended rational bi-quartic spline interpolation defined in Eq. 17 is C1 continuous everywhere provided that all free parameters are positive.
Theorem 2. The proposed rational bi-quartic spline interpolation given in Eq. 17 is shape-preserving provided that all four curve networks defined in Eqs. 18–21 are shape-preserving.
Theorem 3. The partially blended rational bi-quartic spline defined by Eq. 17 is a continuous degree-nine piecewise rational surface (degree seven in the numerator and degree two in the denominator).The proof of Theorems 1–3 can be obtained by generalizing Proposition 3 in the study of Casciola and Romani [16].The behavior of local control of the proposed partially blended rational bi-quartic spline based on the observations given in Equations (21)–(25) is given in Example 1.
Example 1. Data from the following function are truncated to five decimal places [50]:Figure 2 shows the interpolating surface by using the proposed scheme for data listed in Table 2.Figure 2A shows the interpolating surface with , and Meanwhile, Figure 2B shows the surface with , and Clearly, the surfaces are bilinear or flat. A similar effect can be seen in Figure 2C, i.e., when , and Meanwhile, Figure 2D shows the behavior of the interpolating surface when , become a larger number, i.e., equal to 50. Finally, Figure 2E shows the surface when become a smaller number, i.e., equal to 0.005.Furthermore, the parameters cannot be negative because the interpolating surface will become a rational surface, i.e., having poles on the surface. This will result in the discontinuities of the surface. This fact can be seen in Figure 3. Thus, the parameters are restricted to having positive values only. This is in line with the existing works such as Karim et al. [36], Abbas et al. [50, 51], Hussain and Ali [37], Hussain and Hussain [40], Tian et al. [41], and Sarfraz et al. [9].Finally, as the parameters , are increased and are fixed, we found that the rational interpolant is stable, i.e., the produced surfaces are almost identical to each other. Meanwhile, as the parameters keep increasing while parameters , we also found that the rational interpolant is stable. These results are shown in Figure 4.
FIGURE 2
TABLE 2
| 0 | 2 | 4 | 6 | |
|---|---|---|---|---|
| 0 | 1.33000 | 0.011261 | 0.10505 | 0.41710 |
| 2 | 1.79240 | 0.619300 | 0.39739 | 0.45990 |
| 4 | 0.41370 | 0.020814 | 0.16294 | 0.33635 |
| 6 | 0.39537 | 0.281670 | 0.30087 | 0.33560 |
Data from function .
FIGURE 3
FIGURE 4
Local Quadratic Reproducing
The proposed partially blended rational bi-quartic spline defined in Eq. 17 is local quadratic reproducing. To prove this property, we use numerical simulations by using the following quadratic test function:
Figure 5A shows the true function. Meanwhile, Figures 5B–G show the interpolating surfaces with various parameters and their values. We can see that all figures are almost identical to each other. Table 3 shows the values of RMSE, R2, and maximum absolute error for each figure. Based on Table 3, we found that as or , the proposed rational bi-quartic spline is local quadratic reproducing. This property is important since we can use our scheme to reconstruct any conic sections. Meanwhile, the works by Karim et al. [36] and Abbas et al. [50] are not local quadratic reproducing.
FIGURE 5
Positivity-Preserving Interpolation
The rational bi-quartic spline defined in Eq. 17 is not positivity-preserving for positive datasets. As seen in Figure 2A, the partially blended rational bi-quartic spline is not positivity-preserving interpolation. We can obtain a positive surface by manipulating all free parameters. However, this approach is not practical and really time-consuming. Therefore, we want to find an automatic method to calculate the parameter values that will produce a positive rational interpolant. To do this, we derive the sufficient condition for the positivity of the rational quartic splines on all four boundary curve networks (see Figure 1). On each boundary, there will be two free parameters. This results in eight free parameters for shape modification. The derivation of the sufficient condition is described as follows.
Let be strictly positive data (to avoid dividing by zero at later stages) defined over a rectangular grid , such thatThe partially blended rational bi-quartic spline in Eq. 17 is positive if all boundary curves defined in Eqs. 18–21 are positive, i.e., , and respectively. The following theorem summarizes the construction of positivity-preserving interpolation using the proposed rational bi-quartic spline.
Theorem 4. Given strictly positive surface data the piecewise partially blended rational bi-quartic function in Eq. 17 defined over the rectangular mesh , preserves the positivity of the data if the parameters satisfy the following sufficient conditions:
Proof.To prove the theorem, we adopted the idea from Karim et al. [22]. As conjectured by Casciola and Romani [16], the surface is shape-preserving if all four curve networks are shape-preserving. This can be summarized as that the rational bi-quartic surface is positive if and only if all four curves are positive.From Equations (18)–(21), if Therefore, if and only if and Both conditions can be simplified as follows:Combining conditions in Eqs. 28, 29 results in the following sufficient conditions for the positivity of in inequality form:Likewise, the remaining three boundary curves , and are positive if the following sufficient conditions are satisfied. if and only if and if and only if and : if and only if and :Combining the inequality conditions stated in Eqs. 31–36 leads us to the following conditions:Then from Eqs. 30, 37, the sufficient conditions for the positivity of all the four curve networks , and are obtained. This completes the proof.The sufficient conditions in Eq. 27 can be rewritten aswhere the parameters
Theorem 5. The positive partially blended rational bi-quartic spline constructed by using shape parameters calculated using Eq. 38 is a continuous degree-nine piecewise rational surface with positivity-preserving properties with respect to the positive datasets.
Proof. This is a consequence from Theorem 3.The following algorithm can be used to implement the proposed positivity-preserving interpolation.
Algorithm 1
Input: Strictly positive data points,
Output: Positive interpolating surfaces and RMSE and R2 values.
Step 1: Calculate the first partial derivative values by using numerical techniques, i.e., AMM.
Step 2: Choose suitable positive values for free parameters, and and compute the values of parameters by using Theorem 4.
Step 3: Construct the positive surfaces by substituting all the required parameters into the partially blended rational bi-quartic spline function in Eq. 17.
Repeat Steps 2 and 3 to generate a different positive surface through different positive datasets.
Results and Discussion
We test the proposed scheme by using two positive datasets. Example 2 uses the same data as in Example 1. All computational, simulation, and graphical results are obtained by using MATLAB Version 2019a installed on Windows 10, AMD Ryzen 3 2200G with Radeon Vega Graphics 3.50 GHz. As validation, we calculate RMSE and R2 values and compare the results with those of Karim et al. [36] and Abbas et al. [50].
Example 2. Positive data from the following function are truncated to five decimal places [50]:Figure 6A shows the default bi-quartic Hermite spline for the positive data given in Table 4. Figures 6B,C show the xz-view and yz-view for Figure 6A, respectively. Figure 6D shows positivity-preserving by using the proposed rational bi-quartic spline with . Meanwhile, Figures 6E,F show the xz-view and yz-view for Figure 6D, respectively. Figure 6G shows the other view for Figure 6D. Clearly, the interpolating surfaces are positive on entire given intervals. Meanwhile, Figures 6H,I show the positive interpolating surfaces using Karim et al. [36] and Abbas et al. [50] schemes.
FIGURE 6
TABLE 4
| 0 | 2 | 4 | 6 | |
|---|---|---|---|---|
| 0 | 1.33000 | 0.011261 | 0.10505 | 0.41710 |
| 2 | 1.79240 | 0.619300 | 0.39739 | 0.45990 |
| 4 | 0.41370 | 0.020814 | 0.16294 | 0.33635 |
| 6 | 0.39537 | 0.281670 | 0.30087 | 0.33560 |
Positive surface data from function .
Example 3. Positive data from the following function are truncated to four decimal places [51]:Figure 7A shows the default bi-quartic Hermite spline for the positive data given in Table 5. Figures 7B,C show the xz-view and yz-view for Figure 7A, respectively. Figure 7D shows positivity-preserving by using the proposed rational bi-quartic spline with . The parameters are obtained by applying Theorem 4, i.e., Eq. 27. Meanwhile, Figures 8E,F show the xz-view and yz-view for Figure 8D, respectively. Figure 7G shows the other view for Figure 7A. Meanwhile, Figure 7H shows the other view for the surface in Figure 8D. We found that the interpolating surfaces are positive everywhere. Figures 7I,J show the interpolating surfaces using Karim et al. [36] and Abbas et al. [50] schemes.
FIGURE 7
TABLE 5
| −3 | −2 | −1 | 0 | 1 | 2 | 3 | |
|---|---|---|---|---|---|---|---|
| −3 | 0.0401 | 0.0404 | 0.1755 | 1.0401 | 0.1755 | 0.0404 | 0.0401 |
| −2 | 0.0583 | 0.0586 | 0.1936 | 1.0583 | 0.1936 | 0.0586 | 0.0583 |
| −1 | 0.4078 | 0.4082 | 0.5432 | 1.4079 | 0.5432 | 0.4082 | 0.4078 |
| 0 | 1.0400 | 1.0403 | 1.1753 | 2.0400 | 1.1753 | 1.0403 | 1.0400 |
| 1 | 0.4078 | 0.4082 | 0.5432 | 1.4079 | 0.5432 | 0.4082 | 0.4078 |
| 2 | 0.0583 | 0.0586 | 0.1936 | 1.0583 | 0.1936 | 0.0586 | 0.0583 |
| 3 | 0.0401 | 0.0404 | 0.1755 | 1.0401 | 0.1755 | 0.0404 | 0.0401 |
Positive surface data from function .
FIGURE 8
We compare the performance between the proposed partially blended rational bi-quartic spline interpolation and the works of Abbas et al. [50] and Karim et al. [36] by calculating the root mean square error (RMSE) and coefficient of determination (R2). Tables 6, 7 summarize the main results. It can clearly be seen that the proposed rational bi-quartic spline interpolation is on par with the works of Abbas et al. [50] and Karim et al. [36]. Based on R2, for both examples, the proposed scheme is the best. Furthermore, according to Renka and Brown [53] criteria, the proposed scheme is excellent since the R2 value is in between 0.91 and 0.997. Meanwhile, based on the RMSE value, we found that the proposed scheme is the best when the parameters are From all numerical simulations, we found that the proposed scheme requires less computation times (in seconds) compared with Abbas et al. [50] and Karim et al. [36] schemes. The speedup is around 10% faster than the other schemes. Furthermore, based on graphical results, all schemes are looking competent to each other. Most probably, the proposed scheme is the best for both examples. Based on this, we conclude that the rank from the best scheme to the lowest is as follows:
TABLE 6
| Method | Example 2 | Example 3 | |
|---|---|---|---|
| Abbas et al. [50] | 0.0252 | 0.0024 | |
| Karim et al. [36] | 0.0233 | 0.0012 | |
| The proposed scheme | A | 0.0230 | 0.0012 |
| B | 0.0231 | 0.0012 | |
| C | 0.0229 | 0.0013 | |
| D | 0.0205 | 0.0007 | |
Root mean square error (RMSE) values.
TABLE 7
| Method | Example 2 | Example 3 | |
|---|---|---|---|
| Abbas et al. [50] | 0.8991 | 0.9846 | |
| Karim et al. [36] | 0.9001 | 0.9952 | |
| The proposed scheme | A | 0.9013 | 0.9946 |
| B | 0.9007 | 0.9947 | |
| C | 0.9016 | 0.9945 | |
| D | 0.9129 | 0.9965 | |
Coefficient of determination (R2) values.
Parameters’ arrangement: A- , and
B- , and
C- , and
D- , and
The Proposed Scheme > Karim et al. [36] > Abbas et al. [50]
Our final example shows that the proposed rational bi-quartic spline for positivity-preserving has the capability to reduce the CPU times (in seconds) to construct the interpolating surface compared with interpolation without positivity-preserving. Table 8 shows the positive data from the following function:Figure 8 shows the examples of various interpolating surfaces for data given in Table 8. Figure 8A shows the interpolating surface when the free parameters are , and Figures 8B,C show the xz-view and yz-view for Figure 8A. Clearly, the surfaces are not positive everywhere. Figure 8D shows the positive interpolating surface after applying Theorem 4. The free parameters are . Meanwhile, are obtained through Eq. 27. Figures 8E,F show the xz-view and yz-view for Figure 8D. From both graphs, the surfaces are positive and lie above both x and y axes. Figure 8G shows another interpolating surface with parameters . However, the surface is not positive as seen in Figure 8H. Finally, Figure 8I shows the interpolating surface using the proposed scheme. From Figure 8J, it is seen that the interpolating surface shown in Figure 8I is positive everywhere. Therefore, Theorem 4 is confirmed to produce a positive interpolating surface on entire given intervals. This is in contrast with what was claimed by Qin et al. [38]. Therefore, the proposed scheme indeed has reduced the number of mathematical derivations needed to produce a positive rational interpolant. This is very significant in numerical analysis aspects. We only require lesser computation time compared with the work of Qin et al. [38]. CPU times (in seconds) are measured by using MATLAB tic and toc commands.
TABLE 8
| 0 | 2 | 4 | 6 | |
|---|---|---|---|---|
| 0 | 0.2500000 | 0.0147059 | 0.000972763 | 0.00019275 |
| 2 | 0.0147059 | 0.003846154 | 0.000623441 | 0.00015615 |
| 4 | 0.000972763 | 0.0006234414 | 0.000243902 | 9.2421e−05 |
| 6 | 0.000192752 | 0.0001561524 | 9.242144e−05 | 4.8216e−05 |
Positive surface data from function .
As seen in Table 9, by applying the positivity-preserving interpolation scheme to the given data, the CPU times (in seconds) needed to construct the interpolating surface are reduced. This is true for all set of parameters. We also noticed that the proposed sufficient condition for the positivity of the rational quartic interpolant is guaranteed to produce a positive interpolating surface everywhere, i.e., on entire given intervals. However, we also found that, by applying positivity-preserving interpolation, R2 values will be reduced, while at the same time, the RMSE value will be increased. However, the values are not much different.
TABLE 9
| The proposed scheme | R2 | RMSE | Number of points below x and y axes | CPU times (in seconds) | |
|---|---|---|---|---|---|
| I | Non-PP | 0.9227 | 8.8280e−05 | 110 | 2.435965 |
| PP | 0.9067 | 1.0117e−04 | 0 | 2.008921 | |
| II | Non-PP | 0.9280 | 8.2205e−05 | 56 | 2.379029 |
| PP | 0.9062 | 1.0173e−04 | 0 | 1.971622 | |
| III | Non-PP | 0.9126 | 9.9883e−05 | 160 | 2.968662 |
| PP | 0.9153 | 9.5954e−05 | 0 | 1.795824 | |
| IV | Non-PP | 0.9333 | 7.6204e−05 | 28 | 2.28106 |
| PP | 0.9335 | 7.5327e−05 | 0 | 1.567297 | |
Error estimation.
PP: positivity-preserving.
Note: sets of parameters.
I- and
II- and
III- and
IV- and
Final Remark: Qin et al. [
38] have proposed a new approach to derive the sufficient condition for the positivity of the partially blended rational spline. They claimed that their condition is guaranteed to produce a positive interpolating surface everywhere. However, their method is mathematically more complicated to apply since it will increase the unnecessary computation time to produce the positive surface. If we apply the Qin et al. [
38] method, the sufficient condition for the positivity of the rational quartic spline will take the following form:
It is very complicated to solve
Eq. 43until
Eq. 46. This is because we need to simplify all equations to form a new rational function, i.e., quintic numerator and quadratic denominator. Compared with our derivation for the positivity of the rational quartic spline which is direct and easy as shown in
Eq. 28until
Eq. 36. Furthermore, based on all numerical and graphical results presented in this study, we found that the sufficient condition on all four curve networks is sufficient to produce positive interpolating surfaces on the entire given domain and not just on the four boundary curves only. This fact can be seen clearly in the final example. When we applied the proposed scheme, we found that, without positivity-preserving interpolation, the z values are negative. Meanwhile, after we apply the sufficient condition stated in Theorem 4, we obtain the interpolating surface as positive everywhere as seen in
Figures 6–
8as well as in
Table 9. Therefore, we encouraged the user to use our proposed rational bi-quartic spline for positivity-preserving interpolation with the following merits:
1) It will guarantee to produce a positive interpolating surface on entire given intervals and not just on the four boundary curve networks as claimed by Qin et al. [38]. From all numerical examples, we found that, after the positivity conditions are applied, all points on the surface are positive.
2) It can reproduce an exactly quadratic polynomial unlike Abbas et al. [50], Karim et al. [36], and Qin et al. [38].
3) Overall, based on R2 and RMSE values, the proposed scheme is the best compared with some existing schemes.
Conclusion
In this study, a new partially blended rational bi-quartic spline interpolation with C1 continuity is constructed for positivity-preserving application. The sufficient condition for the positivity of the rational quartic spline is derived on four boundary curves defined on rectangular meshes. This condition results in eight free parameters for shape modification. The proposed scheme is tested for many positive surface datasets. Based on numerical analysis and graphical comparison, the proposed scheme is the best, i.e., higher R2 and smaller RMSE. We also found that the proposed scheme has the capability to reduce the CPU times (in seconds) when we want to construct a positive interpolating surface. Furthermore, the proposed positivity-preserving interpolation scheme is guaranteed to deliver a positive surface on the entire given interval. The derivation for the positivity sufficient condition in this study is easier than that in the study of Qin et al. [38]. By using our proposed scheme, we have reduced the required numbers of mathematical derivations, and our scheme is easy to use. On top of that, the proposed scheme is local quadratic polynomial reproducing. This is significant when the user wants to reconstruct any quadric surfaces. Future studies will be focusing on the construction of general shape-preserving for scattered data arranged on triangular meshes as well as its applications in image processing and image denoising such as in the studies of Karim and Saaban [55], Karim et al. [56], Walther and Schmidt [57], and Zulkifli et al. [58].
Statements
Data availability statement
The original contributions presented in the study are included in the article/supplementary material, and further inquiries can be directed to the corresponding author.
Author contributions
SA, AS, and VN curated the data, involved in formal analysis, obtained the resources, and wrote the original draft. SA and VN performed the methodology. SA acquired the funding and investigated the data. VN involved in project administration. SA and AS ran the software, validated and visualized the results, and reviewed and edited the paper.
Funding
The authors would like to acknowledge Universiti Teknologi PETRONAS (UTP) and Ministry of Education (MOE), Malaysia, for the financial support received in the form of a research grant: FRGS/1/2018/STG06/UTP/03/1/015MA0-020 and YUTP: 015LC0-315 (Uncertainty estimation based on quasi-Newton methods for Full Waveform Inversion (FWI)).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
1.
DuanQDjidjeliKPriceWGTwizellEH. Constrained Control and Approximation Properties of a Rational Interpolating Curve. Inf Sci (2003) 152:181–94. 10.1016/s0020-0255(02)00409-7
2.
DuanQWangLTwizellEH. A New Bivariate Rational Interpolation Based on Function Values. Inf Sci (2004) 166:181–91. 10.1016/j.ins.2003.12.001
3.
DuanQZhangHZhangYTwizellEH. Bounded Property and Point Control of a Bivariate Rational Interpolating Surface. Comput Math Appl (2006) 52:975–84. 10.1016/j.camwa.2006.04.021
4.
DuanQBaoFZhangY. Shape Control of a Bivariate Interpolating Spline Surface. Int J Comp Math (2008) 85(No. 5):813–25. 10.1080/00207160701463245
5.
FarinGE. Curves and Surfaces for CAGD: A Practical Guide. 5th ed.San Francisco, CA: Morgan Kaufmann Publishers (2002).
6.
BrodlieKWButtS. Preserving Convexity Using Piecewise Cubic Interpolation. Comput Graphics (1991) 15(1):15–23. 10.1016/0097-8493(91)90026-e
7.
BrodlieKMashwamaPButtS. Visualization of Surface Data to Preserve Positivity and Other Simple Constraints. Comput Graphics (1995) 19(4):585–94. 10.1016/0097-8493(95)00036-c
8.
ButtSBrodlieKW. Preserving Positivity Using Piecewise Cubic Interpolation. Comput Graphics (1993) 17(1):55–64. 10.1016/0097-8493(93)90051-a
9.
SarfrazMHussainMZNisarA. Positive Data Modeling Using Spline Function. Appl Math Comput (2010) 216:2036–49. 10.1016/j.amc.2010.03.034
10.
BeliakovG. Monotonicity Preserving Approximation of Multivariate Scattered Data. Bit Numer Math (2005) 45(4):653–77. 10.1007/s10543-005-0028-x
11.
AsimMRMustafaGBrodlieKW. Constrained Visulaization of 2D Positive Data Using Modified Quadratic Shepard Method. MSCG POSTERS Proceedings, WSCG’2004. Plzen, Czech Republic. Curve drawing subject to positivity and more general constraints. Comput Graphics (2004) 27:469–85.
12.
BaoFSunQDuanQ. Point Control of the Interpolating Curve with a Rational Cubic Spline. J Vis Commun Image Representation (2009) 20:275–80. 10.1016/j.jvcir.2009.03.003
13.
WuJZhangXPengL. Positive Approximation and Interpolation Using Compactly Supported Radial Basis Functions. Math Probl Eng (2010) 2010. 10.1155/2010/964528
14.
WuJLaiYZhangX. Radial Basis Functions for Shape Preserving Planar Interpolating Curves. J Inf Comput Sci (2010) 7(7):1453–8.
15.
ZhuY. C2 Positivity-Preserving Rational Interpolation Splines in One and Two Dimensions. Appl Math Comput (2018) 316:186–204. 10.1016/j.amc.2017.08.026
16.
CasciolaGRomaniL. Rational Interpolants with Tension Parameters, In: TLycheM-LMazureLLSchumaker, editors Curve and Surface Design. Nashboro Press (2003). pp. 41–50.
17.
SchmidtJWHessW. Positive Interpolation with Rational Quadratic Splines. Computing (1987) 38:261–7. 10.1007/bf02240100
18.
SchmidtJWHessW. Positivity of Cubic Polynomials on Intervals and Positive Spline Interpolation. BIT (1988) 28:340–52. 10.1007/bf01934097
19.
SchmidtJW. Positivity, Monotone and -convex Interpolation on Rectangular Grids. Computing (1992) 48(3-4):363–71. 10.1007/bf02238643
20.
HussainMZSarfrazMShakeelA. Shape Preserving Surfaces for the Visualization of Positive and Convex Data Using Rational Bi-quadratic Splines. Int J Compu Appl (2011) 27(10):12–20.
21.
ChanESOngBH. Range Restricted Scattered Data Interpolation Using Convex Combination of Cubic Bézier Triangles. J Comput Appl Math (2001) 136:135–47. 10.1016/s0377-0427(00)00580-x
22.
LuoZPengX. A C1-Rational Spline in Range Restricted Interpolation of Scattered Data. J Comput Appl Math (2006) 194(2):255–66. 10.1016/j.cam.2005.07.010
23.
PiahARMGoodmanTNTUnsworthK. Positivity-preserving Scattered Data Interpolation. Lecture Notes Comp Sci (LNCS) (2005) 3604:336–49. 10.1007/11537908_20
24.
HussainMZHussainM. C1 Positive Scattered Data Interpolation. Comput Math Appl (2010) 59:457–67. 10.1016/j.camwa.2009.06.019
25.
SaabanAMajidAAPiahARM. Visualization of Rainfall Data Distribution Using Quintic Triangular Bezier Patches. Bull Malays Math Sci Soc (2009) 32(2):137–50.
26.
HussainMZHussainMSamreenS. Visualization of Shaped Data by Rational Quartic Spline Interpolation. Int J Appl Math Stat (2008) 13(J08):34–46.
27.
HussainMZHussainM. Visualization of 3D Data Preserving Convexity. J Appl Math Comput (2007) 23(2):397–410. 10.1007/bf02831986
28.
WangQTanJ. Rational Quartic Spline Involving Shape Parameters. J Inf Comput Sci (2004) 1(1):127–30.
29.
HarimNAKarimAAAOthmanOSaabanAGhaffarANisarKSet alPositivity Preserving Interpolation by Using Rational Quartic Spline. AIMS Math (2020) 5(4):3762–82. 10.3934/math.2020244
30.
LiuXZhuYLiuS. Positivity and Monotonicity Preserving Biquartic Rational Interpolation Spline Surface. J Appl Math (2014) 2014:11Article ID 987076. 10.1155/2014/987076
31.
HussainMZSarfrazM. Positivity-Preserving Interpolation of Positive Data by Rational Cubics. J Comput Appl Math (2008) 218:446–58. 10.1016/j.cam.2007.05.023
32.
HanX. Convexity-preserving Piecewise Rational Quartic Interpolation. SIAM J Numer Anal (2008) 46(2):920–9. 10.1137/060671577
33.
HanX. Shape-preserving Piecewise Rational Interpolant with Quartic Numerator and Quadratic Denominator. Appl Math Comput (2015) 251:258–74. 10.1016/j.amc.2014.11.067
34.
HarimANKarimSAAOthmanMSaabanAGhaffarANisarKSet alPositivity Preserving Interpolation by Using Rational Quartic Spline. AIMS Math (2020) 5(4):3762–82. 10.3934/math.2020244
35.
HarimAKarimSAAOthmanMSaabanA. High Accuracy Data Interpolation Using Rational Quartic Spline with Three Parameters. Int J Scientific Tech Res (2019) 1219–27432.
36.
KarimSAAKongVPSaabanA. Positivity Preserving Interpolation Using Rational Bi-cubic Spline. J Appl Math (2015) 2015. 10.1155/2015/572768
37.
HussainMZAliJM. Positivity Preserving Piecewise Rational Cubic Interpolation. Matematika (2006) 22(2):147–53.
38.
QinXQinLXuQ. C1 Positivity-Preserving Interpolation Schemes with Local Free Parameters. IAENG Int J Comp Sci (2016) 43(2).
39.
GregoryJA. Shape Preserving Spline Interpolation. Computer-Aided Des (1986) 18(1):53–7. 10.1016/s0010-4485(86)80012-4
40.
HussainMZHussainM. Visualization of Data Subject to Positive Constraint. J Inf Comput Sci (2006) 1(3):149–60.
41.
TianMZhangYZhuJDuanQ. Convexity-Preserving Piecewise Rational Cubic Interpolation. J Inf Comput Sci (2005) 2(4):799–803.
42.
HussainMZSarfrazMShaikhTS. Shape Preserving Rational Cubic Spline for Positive and Convex Data. Egypt Inform J (2011) 12:231–6. 10.1016/j.eij.2011.10.002
43.
AshrafPSabirMGhaffarANisarKSKhanI. Shape-Preservation of the Four-Point Ternary Interpolating Non-stationary Subdivision Scheme. Front Phys (2020) 7:241. 10.3389/fphy.2019.00241
44.
FiorotJCTabkaJ. Shape-Preserving C 2 Cubic Polynomial Interpolating Splines. Math Comput (1991) 57:291–8. 10.2307/2938674
45.
PengXLiZSunQ. Nonnegativity Preserving Interpolation by Bivariate Rational Spline Surface. J Apllied Math (2012) 2012:Article ID 624978.
46.
FritschFNButlandJ. A Method for Constructing Local Monotone Piecewise Cubic Interpolants. SIAM J Sci Stat Comput (1984) 5:300–4. 10.1137/0905021
47.
FritschFNCarlsonRE. Monotone Piecewise Cubic Interpolation. SIAM J Numer Anal (1980) 17:238–46. 10.1137/0717021
48.
WangQTanJ. Shape Preserving Piecewise Rational Biquartic Surface. J Inf Comput Sci (2006) 3:295–302.
49.
WangYTanJLiZBaiT. Weighted Rational Quartic Spline Interpolation. J Inf Comput Sci (2013) 10(9):2651–8. 10.12733/jics20101820
50.
AbbasMMajidAAAliJM. Shape Preserving Rational Bi-cubic Function for Positive Data. World Appl Sci J (2012) 18(11):1671–9. 10.5829/idosi.wasj.2012.18.11.3527
51.
AbbasMMajidAAAliJM. Positivity-preserving Rational Bi-cubic Spline Interpolation for 3D Positive Data. Appl Math Comput (2014) 234:460–76. 10.1016/j.amc.2014.02.031
52.
IbraheemFHussainMHussainMZBhattiAA. Positive Data Visualization Using Trigonometric Polynomials. J Appl Math (2012) 2012:19Article ID 247120. 10.1155/2012/247120
53.
RenkaRGBrownR. Algorithm 792: Accuracy Tests of ACM Algorithms for Interpolation of Scattered Data in the Plane. ACM Trans Math Softw (1999) 25:1. 10.1145/305658.305745
54.
DelbourgoRGregoryJA. The Determination of Derivative Parameters for a Monotonic Rational Quadratic Interpolant. IMA J Numer Anal (1985) 5:397–406. 10.1093/imanum/5.4.397
55.
KarimSAASaabanA. Shape Preserving Interpolation and its Application in Image Processing. Math Problem Eng (2017) 12:7459218. 10.1155/2017/7459218
56.
KarimSAASaabanASkalaV. Range-Restricted Surface Interpolation Using Rational Bi-cubic Spline Functions with 12 Parameters. IEEE Access (2019) 7:104992–5007. 10.1109/ACCESS.2019.2931454
57.
WaltherMBSchmidtJW. Range Restricted Interpolation Using Gregory’s Rational Cubic Splines. J Comput Appl Math (1999) 103:221–37.
58.
ZulkifliNAKarimSAAShafieASarfrazMGhaffarANisarKS. Image Interpolation Using a Rational Bi-cubic Ball. Mathematics (2019) 7:1045. 10.3390/math7111045
Summary
Keywords
rational bi-quartic spline, surface interpolation, quadratic reproducing, positivity-preserving, partially blended, numerical analysis
Citation
Abdul Karim SA, Saaban A and Nguyen VT (2021) Construction of C1 Rational Bi-Quartic Spline With Positivity-Preserving Interpolation: Numerical Results and Analysis. Front. Phys. 9:555517. doi: 10.3389/fphy.2021.555517
Received
05 May 2020
Accepted
30 June 2021
Published
19 August 2021
Volume
9 - 2021
Edited by
Alex Hansen, Norwegian University of Science and Technology, Norway
Reviewed by
Muhammad Abbas, Ton Duc Thang University, Vietnam
Abdul Ghafar, Ton Duc Thang University, Vietnam
Updates
Copyright
© 2021 Abdul Karim, Saaban and Nguyen.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Samsul Ariffin Abdul Karim, samsul_ariffin@utp.edu.my
†ORCID:Samsul Ariffin Abdul Karimorcid.org/0000-0001-6518-6705
This article was submitted to Computational Physics, a section of the journal Frontiers in Physics
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.