Comonotone parametric C1 interpolation of nongridded data

We present a local scheme for constructing a C1 interpolating function comonotone with a set of data having a tensor product topology. The constructed interpolant is a parametric surface with piecewise bicubic components which locally maintain the monotonicity of the data along curves “parallel” to the related coordinate lines. Error estimates and graphical examples are provided.