Monotone Data Visualization Using Rational Trigonometric Spline Interpolation
Rational cubic and bicubic trigonometric schemes are developed to conserve monotonicity of curve and surface data, respectively. The rational cubic function has four parameters in each subinterval, while the rational bicubic partially blended function has eight parameters in each rectangular patch. The monotonicity of curve and surface data is retained by developing constraints on some of these parameters in description of rational cubic and bicubic trigonometric functions. The remaining parameters are kept free to modify the shape of curve and surface if required. The developed algorithm is verified mathematically and demonstrated graphically.
The technique or algorithm employed in creating images, diagrams, or animations for imparting a piece of information is termed as visualization. It has a key role to play in different fields like science, engineering, education, and medicine as it can aid experts in identifying and interpreting different patterns and artifacts in their data and provide a three-dimensional display of data for the solution of a wide range of problems.
The methods used to obtain visual representations from abstract data have been in practice for a long time. However, physical quantities often emanate distinctive features (such as positivity, convexity, and monotonicity) and it becomes imperative that the visual model must contain the shape feature to fathom the physical phenomenon, the scientific experiment, and the idea of the designer. Spline interpolating functions play elemental role in visualizing shaped data. This paper specifically addresses the problem of visualizing monotone curve and surface data.
Monotonicity is an indispensable characteristic of data stemming from many physical and scientific experiments. The relationship between the partial pressure of oxygen and percentage dissociation of hemoglobin, consumption function in economics, concentration of atrazine and nitrate in shallow ground waters, and approximation of couples and quasi couples are few phenomena which exhibit monotone trend.
Efforts have been put in by many researchers and a variety of approaches has been proposed to solve this eminent issue [1–17]. Cripps and Hussain  visualized the 2D monotone data by Bernstein-Bézier rational cubic function. The authors in  converted the Bernstein-Bézier rational cubic function to cubic Hermite by applying the continuity conditions at the end points of interval. The lower bounds of weights functions were determined to visualize monotone data as monotone curve. Hussain and Sarfraz  have conserved monotonicity of curve data by rational cubic function with four shape parameters, two of which were set free and two were shape parameters. Data dependent constraints on shape parameters were developed which assure the monotonicity but one shape parameter is dependent on the other which makes it economically very expensive. Rational cubic function with two shape parameters suggested by Sarfraz  sustained monotonicity of curves but lacked the liberty to amend the curve which makes it inappropriate for interactive design. Piecewise rational cubic function was used by M. Z. Hussain and M. Hussain  to visualize 2D monotone data by developing constraints on the free parameters in the specification of rational cubic function. The authors also extended rational cubic function to rational bicubic partially blended function. Simple constraints were derived on the free parameters in the description of rational bicubic partially blended patches to visualize the 3D monotone data. Three kinds of monotonicity preservation of systems of bivariate functions on triangle were defined and studied by Floater and Peña . Sarfraz et al.  developed constraints in the specification of a bicubic function to visualize the shape of 3D monotone data.
This paper is a noteworthy addition in the field of shape preservation when the data under consideration admits monotone trend. The suggested algorithm offers numerous advantages over the prevailing ones. Orthogonality of sine and cosine function compels much smoother visual results as compared to algebraic spline. Derivative of the trigonometric spline is much lower than that of algebraic spline. Moreover, trigonometric splines play an instrumental role in robotic manipulator path planning.
The remainder of the paper is structured as follows. Section 2 is devoted to reviewing the rational trigonometric cubic function developed in . In Section 3, rational trigonometric cubic function is extended to rational trigonometric bicubic function. Section 4 aims to develop monotonicity preserving constraints for 2D data. Section 5 submits a solution to shape preservation of 3D monotone data. In Section 6, numerical examples have been demonstrated. Section 7 draws the conclusion and significance of this research.
2. Rational Trigonometric Cubic Function
In this section, rational trigonometric cubic function  is reviewed.
Let be the given set of data points defined over the interval , where . Piecewise rational trigonometric cubic function is defined over each subinterval as where , .
The rational trigonometric cubic function (1) is ; that is, it satisfies the following properties: Here and are derivatives at the end points of the interval . The parameters and are real numbers used to modify the shape of the curve.
3. Rational Trigonometric Bicubic Partially Blended Function
Let be the 3D regular data set defined over the rectangular mesh , let be a partition of , and let be a partition of . Rational trigonometric bicubic function which is an extension of rational trigonometric cubic function (1) is defined over each rectangular patch , where , as where , , , and are rational trigonometric bicubic functions defined on the boundary of rectangular patch as
4. Monotone Curve Interpolation
Monotonicity is a crucial shape property of data and it emanates from many physical phenomenon, engineering problems, scientific applications, and so forth, for instance, dose response curve in biochemistry and pharmacology, approximation of couples and quasi couples in statistics, empirical option pricing model in finance, consumption function in economics, and so forth. Therefore, it is customary that the resulting interpolating curve must retain the monotone shape of data.
In this section, constraints on shape parameters in the description of rational trigonometric cubic function (1) have been developed to preserve 2D monotone data.
Let be the monotone data defined over the interval ; that is, The curve will be monotone if the rational trigonometric cubic function (1) satisfies the condition Now, we have where The denominator in (16) is a squared quantity, thus, positive. Hence, monotonicity of rational trigonometric cubic spline depends upon the positivity of numerator which can be attained if the coefficients of the trigonometric basis functions are all positive. This yields the following result:
The above discussion can be summarized as follows.
Theorem 1. The piecewise trigonometric rational cubic function (1) preserves the monotonicity of monotone data if in each subinterval , the parameters and satisfy the following sufficient conditions: The above constraints can be rearranged as
Step 1. Take a monotone data set .
Step 2. Use the Arithmetic Mean Method  to estimate the derivatives ’s at knots ’s (note: Step 2 is only applicable if data is not provided with derivatives).
Step 3. Compute the values of parameters ’s and ’s using Theorem 1.
Step 4. Substitute the values of variables from Steps 1–3 in rational trigonometric cubic function (1) to visualize monotone curve through monotone data.
5. Monotone Surface Interpolation
Now, is monotone if , where with Now the positivity of entirely depends on . The denominator in (22) is always positive. Since the parameter lies in first quadrant therefore the trigonometric basis functions will be positive also. This yields the following constraints on the free parameters:
is monotone if
where with The denominator in (26) is always positive. Moreover, the trigonometric basis functions are also positive for . It follows that the positivity of entirely depends upon . This yields the following constraints on the free parameters:
is monotone if . We have where Since the denomoinator of (29) is always positive and trigonometric basis functions are positive for so the positivity of . It follows that the positivity of entirely depends upon . This yields the following constraints on the free parameters:
is monotone if . We have where Finally, is positive if are positive. This yields the following constraints on the free parameters: The above discussion can be put forward as the following theorem.
Theorem 3. The bicubic partially blended rational trigonometric function defined in (4) visualizes monotone data in view of the monotone surface if in each rectangular grid , free parameters satisfy the following constraints: The above constraints are rearranged as
Step 1. Take a 3D monotone data set .
Step 2. Use the Arithmetic Mean Method to estimate the derivatives at knots (note: Step 2 is only applicable if data is not provided with derivatives).
Step 3. Compute the values of parameters using Theorem 3.
Step 4. Substitute the values of variables from Steps 1–3 in rational trigonometric cubic function (4) to visualize monotone surface through monotone data.
6. Numerical Example
This section illustrates the monotonicity preserving schemes developed in Sections 4 and 5 with the help of examples. The data in Table 1 is observed by exposing identical samples of hemoglobin to different partial pressures of oxygen which results in varying degree of saturation of hemoglobin with oxygen. The sample obtaining the highest amount is said to be saturated. The amount of oxygen combined with the remaining samples is taken as percentage of this maximum value. At a low partial pressure of oxygen, the percentage saturation of hemoglobin is very low; that is, hemoglobin is combined with only a very little oxygen. At high partial pressure of oxygen, the percentage saturation of hemoglobin is very high; that is, hemoglobin is combined with large amounts of oxygen, that is, a monotone relation, so the resulting curve must exhibit the same behavior. Figure 1 represents the curve created by assigning random values to free parameters in description of rational trigonometric cubic function (1) which does not retain the monotone nature of the data. This impediment is removed by applying monotonicity preserving schemes developed in Section 4 and is shown in Figure 2. It is evident from the figure that this curve preserves the monotone shape of hemoglobin dissociation curve. Similar investigation in Table 2 displays a series of results for percentage saturation of myoglobin and partial pressure of oxygen. Figure 3 is produced by assigning random values to free parameters in description of rational trigonometric cubic function (1) which fails to conserve the monotone trend of data. Algorithm 2 developed in Section 4 is applied to remove this drawback and Figure 4 displays the required result. Numerical results corresponding to Figures 2 and 4 are shown in Tables 3 and 4.
Figures 5 and 7 are produced by interpolating the monotone data sets in Tables 5 and 6, respectively, by rational trigonometric bicubic function for arbitrary values of free parameter. Monotone surfaces in Figures 6 and 8 are produced by interpolating the same data by the monotonicity preserving scheme developed in Section 5. Tables 7 and 8 enclose numerical results against Figures 6 and 8.
In this paper, monotonicity of data is retained by developing constraints on free parameters in the specification of rational trigonometric function and bicubic blended function. Authors in [7, 8] used algebraic function while the proposed algorithm applies trigonometric function which gives much smoother result due to orthogonality of sine and cosine function. Shape preserving techniques of Butt and Brodlie  required insertion of additional knots. In , developed scheme failed to maintain smoothness. The proposed technique is local, affirms smoothness, works well for data with derivatives, and does not require insertion of extra knots. Derivative of trigonometric spline is much lower than that of polynomial spline.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
P. Costantini and F. Fontanella, “Shape-preserving bivariate interpolation,” SIAM Journal on Numerical Analysis, vol. 27, no. 2, pp. 488–506, 1990.View at: Google Scholar
R. D. Fuhr and M. Kallay, “Monotone linear rational spline interpolation,” Computer Aided Geometric Design, vol. 9, no. 4, pp. 313–319, 1992.View at: Google Scholar
M. S. Floater and J. M. Peña, “Monotonicity preservation on triangles,” Mathematics of Computation, vol. 69, no. 232, pp. 1505–1519, 2000.View at: Google Scholar
T. N. T. Goodman, B. H. Ong, and K. Unsworth, “Constrained interpolation using rational cubic splines,” in NURBS for Curve and Surface Design, G. Farin, Ed., pp. 59–74, 1991.View at: Google Scholar
M. Sarfraz, S. Butt, and M. Z. Hussain, “Surfaces for the visualization of scientific data preserving monotonicity,” in Proceedings of the 7th IMA Mathematics for Surfaces Conference, T. N. T. Goodman and R. Martin, Eds., pp. 479–495, Dundee, UK, September 1997.View at: Google Scholar