> --- Ursprüngliche Nachricht --- > Hope this helps answer your question. Thanks, this is basically what I expected - sum of weighted *squared* errors + smoothness term. Btw, my understanding of thin plate splines (TPS) is, that the scattered data are basically fitted by a sum of N radial basis functions (r^2*ln(r)), each of the RBFs centered at one of the N scattered data points (plus optional linear/polynomial terms and an optional smoothness term). Does this apply to RSPLs as well, or are there some differences? Does the RSPL grid simply represent a resampling of the TPS at regular intervals? And is my assumption correct, that the spline_interp() method does not use the TPS for interpolation, but only the scattered data fitting involves TPS? Thank you, Gerhard -- Gerhard Fuernkranz nospam456@xxxxxx 5 GB Mailbox, 50 FreeSMS http://www.gmx.net/de/go/promail +++ GMX - die erste Adresse für Mail, Message, More +++