[SciPy-user] [Fwd: 3D interpolation over irregular data]

Stéfan van der Walt stefan@sun.ac...
Wed Jul 30 16:11:36 CDT 2008


Hi Mark

2008/7/30 mark starnes <m.starnes05@imperial.ac.uk>:
> Hi everyone,
>
> I've looked through the list here and in Numpy-users, and checked the
> 'net but can't find an answer to this problem (with luck, I've missed
> something obvious!).
>
> I've an array of velocities at 80,000 points, irregularly spaced (from a
> CFD analysis).  I'd like to generate the interpolated velocity at any
> position in the domain, to map the data to an acoustics analysis on a
> different mesh.
>
> I tried a least squares approach but the errors are too large using
> polynomials and trigonometric functions.  My conclusion is that I need a
> nearest-neighbour type interpolation routine.

Also take a look at

http://www.scipy.org/Cookbook/Matplotlib/Gridding_irregularly_spaced_data

Robert Kern's delaunay package does natural neighbour interpolation, IIRC.

Regards
Stéfan


More information about the SciPy-user mailing list