[Numpy-discussion] Vector interpolation on a 2D grid

Ryan May rmay31 at gmail.com
Fri Nov 6 17:33:23 EST 2009


On Fri, Nov 6, 2009 at 4:20 PM, Pierre GM <pgmdevlist at gmail.com> wrote:
> All,
> I have a vector of observations (latitude,longitude,value) that I need
> to interpolate over a given area.
> So far, I used Jeff W's basemap to transform the lat/lon in arbitrary
> x/y and construct a regular grid, Robert K's delaunay triangulation
> routine to project the points on the grid and its nearest-neighbor
> interpolator to get a 2D array of estimates.
> Have a look for yourself at http://grab.by/gY3
> Some colleagues complained that the result looked a bit too choppy,
> meaning that too much weight was given to the actual observations.
> What are my other options ?
> Thx in advance for any idea

You could try smoothing the grid with some kind of weighted average
(like a gaussian kernel).

Out of curiosity, what values are you plotting there? I

Ryan

-- 
Ryan May
Graduate Research Assistant
School of Meteorology
University of Oklahoma



More information about the NumPy-Discussion mailing list