Positivity-Preserving Scattered Data Interpolating Surface using C^1 Piecewise Cubic Triangular Patches

The construction of a bivariate C¹ 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 Bezier points in order to ensure that surfaces comprising cubic Bezier triangular patches are always positive and satisfy C¹ continuity conditions. Initial gradients at the data sites are estimated and then modified if necessary to ensure that these conditions are satisfied. The construction is local and easy to implement. Graphical examples are presented using two test functions