We describe a numerical method to construct volume matching C-1 splines sat
isfying the frequently occurring requirement of nonnegativity. The approxim
ated histogram is assumed to be defined on a partition of a polygonal domai
n into polygonal subdomains. Our idea consists of applying cubic C-1 spline
s on Clough Tocher refinements of a triangulation that is compatible with t
he given partition. The proposed algorithm always works successfully for no
nnegative histograms.