Re: 2D interpolation

• To: mathgroup at smc.vnet.net
• Subject: [mg72916] Re: [mg72901] 2D interpolation
• From: "Josef Otta" <josef.otta at gmail.com>
• Date: Thu, 25 Jan 2007 06:39:24 -0500 (EST)

```Hi Fabrice,
if you use Interpolation your sampling points should lie on regular grid.
That is one of possible reasons why it does not not work. In general case
you can use for example Delaunay triangulation and resample your original
data(included in DiscreteMath`ComputationalGeometry`.) There is nothing
wrong with real numbers, but problem is in choosing sampling points.

Regards,
Josef Otta
http://home.zcu.cz/~jotta

2007/1/24, Jouvenot, Fabrice <F.Jouvenot at liverpool.ac.uk>:
>
> Hello,
>
> I try to have a 2 dimensional interpolation of points {x, y, f(x,y)}.
> After some try, I wrote these lines (as an exemple of what I want to do)
> that works :
> ________________________________________
>
> data = {{1, 1, -10}, {1, 2, 2}, {1, 3, 3}, {1, 4, 4}, {1, 5, 5}, {2, 1,
> 2}, {2, 2, 4}, {2, 3, 6}, {2, 4, 8}, {2, 5, 10}, {3., 1, -9}, {3, 2, 6},
> {3, 3, 90}, {3, 4, 12}, {3, 5, 15}, {4, 1, 4}, {4, 2, 8}, {4, 3, 12},
> {4, 4, 16}, {4, 5, 20}, {5, 1, 5}, {5, 2, 10}, {5, 3, 15}, {5, 4, 20},
> {5, 5, 40.5}};
>
> MatrixForm[data]
>
> test[x_, y_] := Interpolation[data][x, y];
>
> test[1.5, 1.5]
>
> Plot[test[x, 1], {x, 1, 5}];
>
> Plot[test[x, x], {x, 1, 5}];
>
> Plot3D[test[x, y], {x, 1, 5}, {y, 1, 5}];
>
> ________________________________________
>
>
> The problem is that when I want to use a point with a real x or y like
> {1.23, 1.1, -10}, nothing works anymore.
> I think that I miss something essential and fundamental, but I don'y
> know why...