Local Lagrange interpolation schemes for quadratic C-1-splines on arbitrary triangulations with Powell-Sabin splits are constructed. By using the concept of weak interpolation, it is proved that the interpolation method yields optimal approximation order. We test our method by interpolating scattered data and show how the method can be applied for terrain modelling. We compare the interpolating splines on fine and coarse triangulations obtained from thinning strategies and analyze the data reduction.