Shape Preserving Interpolation Using Rational Cubic Ball Function and Its Application in Image Interpolation

New rational cubic Ball interpolation with one parameter is proposed for shape preserving interpolation such as positivity, monotonicity, and convexity preservations and constrained data lie on the same side of the given straight line. To produce shape preserving interpolant, the data dependent sufficient condition is derived on the parameter. The rational bicubic Ball function is constructed by using tensor product approach and it will be used for application in image upscaling. Numerical and graphical results are presented by using Mathematica and MATLAB including comparison with some existing scheme.


Introduction
Shape preserving interpolation and approximation is important in many sciences and engineering applications.Common strategy in shape preserving interpolation is the construction of rational or nonrational interpolant that satisfies some shape preserving criteria.When the interpolant satisfies the conditions, then the resulting interpolating curve or approximating curve preserves the shape of the given data.For instance, if the given data is positive, then the interpolating curve is positive on the whole interval.There is many good research in this area.For instance, Abbas et al. [1] discussed the application of rational cubic spline for image interpolation based on genetic algorithm (GA).Abbas et al. [2] use similar idea as in [1] but they replace the rational cubic spline with rational cubic Ball interpolation.
Majeed et al. [3] discussed the craniofacial reconstruction by using rational cubic Ball with cubic denominator without any free parameters.They utilized the original Ball function from Ball [4].The main difference between their scheme and the proposed scheme in this study is that we construct rational cubic Ball spline with quadratic denominator with one parameter.Thus the computation time will be less compared with the work of Majeed et al. [3] (Table 7).Besides that, there are different types of rational cubic Ball interpolant.For instance, Karim [5,6] uses different form from what we proposed in this study.The main contributions of this study can be summarized as follos: (i) The proposed scheme does not require the modification of the first derivative if shape preserving is violated as shown in the work of Brodlie and Butt [7].(ii) The proposed scheme is guaranteed to produce positive and convex curves for positive and convex data sets while the scheme by Hussain et al. [8,9] might be producing nonpositivity interpolating curve on some interval.
(iii) The proposed scheme requires less computation time compared to the bicubic spline interpolation for image upscaling.
This paper is organized as follows.After the Introduction to the subject matter, Methodology will be discussed in Section 2. This includes the construction of the new rational cubic Ball interpolation and the derivation of the sufficient condition of the rational cubic Ball for constrained data interpolation, monotonicity, and convexity preservation, respectively.Numerical and graphical results are given in Section 3. Section 4 is devoted for application of rational bicubic Ball function for image interpolation.Conclusions are given in the final section.
In this study we adopted AMM since it is simple to use and suitable for all types of data set.The derivation is given as follows [5,6].

Arithmetic Mean Method.
The values of   are given as The derivative value of  0 and   is given as ) .
Shape control analysis can further be elucidated by using a graphical approach.To see the effect of the shape control, by varying the value of shape parameter   ≥ 0,  = 0, 1, . . ., −1, the data given in Table 1 is used.

Application in Shape Preserving
Interpolation.In this Section, we discuss the application of the proposed rational cubic Ball function for shape preserving interpolation.The data under consideration are positive, monotone, convex, and on the same side as the constraint line.
The following theorem gives the sufficient conditions for the positivity of the rational cubic Ball interpolant.It is data dependent and has one free parameter to alter the final positive interpolating curves.
Theorem 1.For a strictly positive data defined in (11), the rational cubic Ball interpolant on [ 0 ,   ] is positive if, in each subinterval [  ,  +1 ],  = 0, 1, . . .,  − 1, if the following sufficient conditions are satisfied: Proof.The sufficient condition of the positivity is derived from condition (13) as follows: Both conditions can be further simplified as follows: The sufficient condition the positivity of the rational cubic Ball defined by ( 1) can be obtained by combining conditions (17) and ( 18) into one condition given as follows: For the purpose of computer implementation, the above condition can be rewritten as follows: such that   > 0 satisfies 0 <   ≤ 0.25.
Algorithm for computer implementation is as follows.
2.6.Constrained Data Interpolation.This section discusses the construction of the sufficient condition for constrained data interpolation for data that lies on the same side as the given constraint line.There are three types for constraint line such as above straight line, below straight line, and in between two lines.In this paper, we will only study the constrained interpolation for the data lie above arbitrary constraint line.
The problem statement can be read as follows.
where a positive value   satisfies the inequality 0 <   ≤ 0.25.

Monotonicity Preserving.
There are two cases for strictly monotonic preserving interpolation.
Theorem 4. The rational cubic Ball interpolation defined in (1) preserves the monotonic increasing data on each subinterval [  ,  +1 ],  = 0, 1, . . .,  − 1 if the parameter satisfies the following condition: Proof.Shape preserving interpolation for strictly monotonic data can be stated as follows.
This completes the proof for Theorem 4.
Thus  (2) () > 0 if and only if  0 > 0 and  3 > 0. Both conditions provide the following inequalities: Simplifying conditions (50) lead to These two conditions give the sufficient condition for the convexity of the rational cubic Ball interpolant The result is summarized as Theorem 1 below.
Theorem 1.The rational cubic Ball defined by (1) preserves the convexity of the data if the parameter satisfies the following condition: with 0 <   ≤ 0.25.

Results and Discussion
This section is devoted to the numerical results for the constrained data interpolation subject to constraint line above straight line, monotonicity preservation and convexity preservation.

Numerical Example 2.
The positive data sets from Hussain et al. [8] are listed in Table 2.

Numerical Example 3. The data in Table
Figure 3(a) shows that the cubic Ball cannot produce the interpolating curve that lies above the given straight line.3.
Applying Theorem 3, we obtain the interpolating curve that lies above a straight line with   = 0.1 for Figure 3(b) and   = 0.1 for Figure 3(c).

Numerical Example 5.
Table 5 shows the convex data from Brodlie and Butt [7] Figure 5(a) shows the cubic Ball cannot produce the convex interpolating curve.Condition (52) guarantees to produce convex interpolating curve on the whole interval.This can be verified by plotting the second derivative graph for Figures 5(a) and 5(b), respectively.From Figure 6(a) the second derivative for cubic Ball polynomial is a straight line and has some negative value which contradicts with the sufficient condition given in (49) and (54).Figure 6(b) shows that the proposed rational cubic Ball interpolant produces the convex interpolating curve since its second derivative is greater than zero.

Application in Image
Processing.An application of the proposed rational cubic Ball interpolation in image processing is discussed in this section.Firstly the univariate spline given in (1) is extended to the bivariate cases by using tensor product approach.The rational bicubic function over each rectangular patch [  ,  +1 ] × [  ,  +1 ],  = 0, 1, . . .,  − 1;  = 0, 1, . . .,  − 1 is defined as follows: where where where   , and   , are partial derivatives on  and  directions, respectively.Meanwhile   , is mixed partial derivatives (twists) at the interior points [12].All can be estimated by using the method discussed in [12].
We apply the rational bicubic Ball function to image interpolation application.The main spotlight here is the image upscaling.Given the original image 256 × 256, then we can upscale the image with factor two or four and so on.To achieve this, we employ the following algorithm.Step 1.Let (, ),  = 1, 2, . . ., ,  = 1, 2, . . .,  represent input pixel indexes and   (0-255) be their corresponding grayscale intensity.Our aim is to find the function; that is,   = (, ) which interpolate the given input pixels.
Step 2. Construct rectangular mesh for given input pixels using our rational bicubic Ball surface defined by (55).
Output Image. by  pixels will be obtained by the following transformation: ) . (58) Step 4. Identify the rectangles with the missing grayscale intensity values.
The original pixels and their intensity values of input image are at vertex of rectangular mesh.
Step 5. Estimate the derivative at each vertex of rectangular (input pixels) by method discussed in [12].
Step 6. Estimate missing grayscale intensity value in Step 4 by using the proposed rational cubic Ball scheme defined in (55).
Step 7. Display the result and compare the performance of proposed method against nearest neighbor, bilinear interpolation, and bicubic spline interpolation by using peak signal to noise ratio (PSNR).
We test the proposed rational bicubic Ball functions for image upscaling by using five different images shown in Figure 7.We compare the performance with nearest neighbor and bilinear and bicubic spline interpolation.Table 6 summarized the results including the value of PSNR.The value of PSNR will reflect the quality of image that is the larger PSNR value the higher quality of images [13].It is measured in decibels (dB).
Peak signal to noise ratio (PSNR) is defined as where MSE is mean square error given as From all PSNR value, we conclude that the proposed scheme is at par with the standard scheme for image interpolation, that is, bicubic spline.In fact, for some images, the proposed scheme gives higher PSNR compared with bicubic spline interpolation.To obtain the result we use parameter value as   = α = 0.1.Finally, for all tested images, the proposed rational bicubic Ball surface gives less computation Problems in Engineering

Conclusion
This study has proposed new rational bicubic Ball function with one parameter.This rational cubic has been used for shape preserving interpolation, that is, positivity, monotonicity, convexity, and constrained data interpolation.The univariate rational cubic Ball is extended to the bivariate cases.The rational bicubic Ball function is used for image upscaling problem in image processing.From all results, we conclude that the proposed scheme is at par with standard scheme, that is, bicubic spline interpolation.Further works can be explored such as to determine the optimum value of the parameters   and α for image processing application.This can be achieved by using genetic algorithm (GA) or neural network (NN) algorithm such as in Hussain et al. [14].Furthermore, we intend to apply the proposed scheme for image resizing and image zooming in medical imaging.
The proposed rational cubic Ball can also is extended to  2 continuous interpolating curve in line with the work of Karim and Kong [15].

Figure 2 :
Figure 2: Positivity preserving using the proposed rational cubic Ball.

Figure 4 (
Figure 4(a) shows the cubic Ball cannot preserve the monotonicity of the data set.This flaw has been recovered nicely by using condition (45), Figure 4(b) with   = 0.1 and Figure 4(c) with   = 0.25.Clearly the resulting monotonic interpolating curve is visually pleasing.

Figure 5 :
Figure 5: Interpolating curve for convex data set.

Algorithm 3 .
Input.Image  by  pixel and parameter values  , and α, .Output.Upscaling image, computation time, and PSNR value.
−        2 (60) with   and   defined as in Step 2 in Algorithm 3.

Table 1 :
Data set for shape control analysis.

Table 6 :
PSNR value for image upscaling.

Table 7 :
Computation time.Mathematical Problems in time comparing with bicubic spline interpolation.Computation time is calculated based on CPU time to generate the upscaling images with factor 4.