Range restricted positivity-preserving scattered data interpolation


  • Abd. Rahni Mt. Piah
  • Azizan Saaban
  • Ahmad Abd. Majid




Scattered data, Interpolation, Positivity, Range-restricted,


The construction of a range restricted bivariate C1 ( or G1 ) interpolant to scattered data is considered in which the interpolant is positive everywhere if the original data are positive. This study is motivated by earlier work in which sufficient conditions are derived on Bézier points in order to ensure that surfaces comprising cubic Bézier triangular patches are always positive and satisfy C1 continuity conditions. In the current work, simpler and more relaxed conditions are derived on the Bézier points. The gradients at the data sites are then calculated (and modified if necessary) to ensure that these conditions are satisfied. Each triangular patch of the interpolating surface is formed as a convex combination of three quartic Bézier triangular patches. Its construction is local and easily extended to include as upper and lower constraints to the interpolant surfaces of the form z = P(x,y) where P is a polynomial of degree less or equal to 4. Moreover, C1 ( or G1 ) piecewise polynomial surfaces consisting of polynomial pieces of the form z = P(x,y) on the triangulation of the data sites are also admissible constraints. A number of examples are presented.


Brodlie, K. W., Butt, S. and Mashwama, P., Computers and Graphics 19 (1995) 585-594.

Chan, E.S. and Ong, B.H., J. Comp. Appl. Math. 136 (2001) 135-147.

Conte, S. D., de Boor, C., Elementary Numerical Analysis, McGraw-Hill, Tokyo, 1992.

Fang, T.P., Piegl, L.A., Computer Aided Design, 24 (1992) 425-436.

Farin, G., Curves and Surfaces for Computer Aided Geometric Design-A Practical Guide 4th Edition,

Academic Press, San Diego, 1996.

Goodman, T.N.T., Said, H.B. and Chang, L.H.T., Appl. Math. Comp. 80 (1994) 1-10.

Halstead, M., Kass, M., DeRose,T.D., SIGGRAPH’93 Conference Proceedings, New York (1993), 35– 44.

Herrmann, M., Mulansky, B. & Schmidt, J. W., Journal of Computational and Applied Mathematics 73(1996) 209 – 223.

Kobbelt, L., Proceedings of Seventh IMA Conference on the Mathematics of Surfaces (1996), 101–131.

Lancaster, P. & Salkauskas, K., Curve and Surface Fitting: An Introduction, Academic Press, San Diego,1986.

Mulansky, B. and Schmidt, J. W., Computing 53 (1994) 137-154.

Ong, B. H. and Wong, H. C., A C1 positivity preserving scattered data interpolation scheme, In Advanced Topics in Multivariate Approximation, Fontanella, F. et. al. (eds) (1996) 259-274.

Piah, A.R.M., Goodman, T.N.T. and Unsworth, K., Positivity preserving scattered data interpolation. In Mathematics of Surfaces, LNCS 3604, R. Martin et. al. (eds), Springer-Verlag Berlin, Heidelberg (2005) 336-349.

Piper. B., Visually smooth interpolation with triangular Bézier patches, In Geomertic Modelling: Algorithm

and New Trends, Farin. G (Ed.), SIAM Philadelphia (1987) 221-233.

Saaban, A., Piah, A.R.M., Majid, A.A. and Chang, L.H.T., G1 scattered data interpolation with minimized sum of squares of principal curvatures, In Computer Graphics, Imaging and Visualization: New Trends, Sarfraz, M. et. al. (eds), IEEE Computer Society, Los Alamitos (2005) 385-390.

Walton, D.J. and Meek, D.S., Computer Aided Design 28(2) (1996) 113–123.