Mathematics Stack Exchange is a question and answer site for people studying math at any level and professionals in related fields. It's 100% free, no registration required.

Sign up
Here's how it works:
  1. Anybody can ask a question
  2. Anybody can answer
  3. The best answers are voted up and rise to the top

I am looking for a way to fit a surface given a set of measured data $(x, y) \mapsto z$. A typical example would consist of anywhere between $10$ and $30$ measurements spread evenly over a disc. Until now it was sufficient to fit a plane through these points with a least square error fit. This model is not sufficient anymore since the accuracy has to be increased. The first thing we tried was to move to more general polynomial fitting. However, low order polynomials ($2$, $3$, $4$) cannot model the surface well enough in general and with higher orders we quickly seem to run into the "Runge phenomenon" that the global fitting error is hard to control.

Therefore I'm currently thinking in the direction of triangulation and piecewise fitting. I can imagine that even piecewise linear fitting could be an option. I'm not sure how to proceed from here however. Are there ways to construct a simple (low count) triangulation and piecewise fit? (The triangulation vertices may be located anywhere and need not coincide with the locations of the measurements.)

Any (other) suggestions are welcome!

share|cite|improve this question
up vote 1 down vote accepted

A cubic surface, which maintains first derivative continuity is a Clough Tocher finite element approach. This does require an estimate of the gradients though. I'm not sure if you have this information easily available or not.

share|cite|improve this answer
Thanks for that reference! – WimC Sep 10 '12 at 10:45

Your Answer


By posting your answer, you agree to the privacy policy and terms of service.

Not the answer you're looking for? Browse other questions tagged or ask your own question.