On Sat, Nov 7, 2009 at 5:38 PM, Pierre GM <pgmdevl...@gmail.com> wrote: > Linear interpolation with the delaunay package doesn't work great for > my data. I played with the radial basis functions, but I'm afraid > they're leading me down the dark, dark path of parameter fiddling. In > particular, I'm not sure how to prevent my interpolated values to be > bounded by the min and max of my actual observations. > Ralf' suggestion of smoothing the values afterwards is tempting, but > sounds a bit too ad-hoc (btw, Ralf, those were relative differences of > monthly average precipitation between El NiƱo and Neutral phases for > November).
That was me, not Ralf. :) And I agree, I the interpolated field does look a bit noisy for such data. I've been doing the smoothing on top of natural neighbor for doing some of my own meteorological analysis. Using the Gaussian kernel isn't really *that* ad hoc considering the prevalence of Barnes/Cressman weighting for spatial averaging typically used in meteorology. And if you have no idea what I'm talking about, Google them, and you'll see. :) Ryan -- Ryan May Graduate Research Assistant School of Meteorology University of Oklahoma _______________________________________________ NumPy-Discussion mailing list NumPy-Discussion@scipy.org http://mail.scipy.org/mailman/listinfo/numpy-discussion