Looking for triangulator/interpolator

Grant Edwards grante at visi.com
Sat May 27 09:11:23 EDT 2006


On 2006-05-27, Travis E. Oliphant <oliphant.travis at ieee.org> wrote:

> Not that you made a bad choice.  I do wonder, how much of your
> difficulty was with the interface to the underlying fitpack
> routines.

I've no idea.  I had never done anything with splines before,
so it's quite possible I just wasn't doing things right.  I
never got it to work at all for non-gridded data (which is what
I needed).  Since it wasn't stable even for gridded data, I
more or less gave up.

> The interface is pretty involved as there are lots of
> parameter choices to the underlying routine.  It's very
> possible the interface is broken in some strange way.

I took a look at underlying Fortran code, but that made my
dizzy.

> Given how many people want to do interpolation.  I'm surprised
> nobody has looked into fixing up the fitpack routines that do
> it.  I suppose it's possible that the underlying fitpack
> routines are faulty in 2-d.
>
> Your examples will help in that study.
>
> Pearu made some nice class-based interfaces to fitpack in
> 2003.  They are in SciPy, but, it appears that people aren't
> making much use of them.  Did you try them?

Not that I remember.  

For that particular application, the interpolation had to be
done in real time, so a triangulation scheme turned out to be a
lot faster [with some constraints on the input data to make the
triangle search O(sqrt(N))].

-- 
Grant Edwards                   grante             Yow!  Are the STEWED PRUNES
                                  at               still in the HAIR DRYER?
                               visi.com            



More information about the Python-list mailing list