Alex Rivera | Logout

How can I interpolate georeferenced data in python?

Asked 2012-07-23T23:15:02.157
10

I have many lines of georeferenced hydrological data with weekly resolution:

Station name, Lat, Long, Week 1 average, Week 2 average ... Week 52 average

Unfortunately, I also have some data with only monthly resolution:

Station name, Lat, Long, January average, February average ... December average

Rather than "reinventing the wheel," can anyone recommend a favorite module, package, or technique that would provide a reasonable interpolation of weekly values from monthly values? Linear would be fine, but it would be nice if we could use the coordinates to improve the interpolation based on nearby stations.

I've tagged this post with python because it's the language I've been using recently (although not its statistical functions). If the answer is "use a stats program like r" so be it, but I'm curious as to what's out there for python. Thanks!

python statistics spatial-interpolation
Edit
Report

1 Answer

4

If you are interested into expanding your experience into R, there are a number of good, well used and documented packages out there. I would start by looking at the Spatial Taskview, which lists what packages can be used for spatial data. One of the paragraphs deals with interpolation. I am most familiar with automap/gstat (I wrote automap), where especially gstat is a powerfull geostatistics package which supports a wide range of methods.

http://cran.r-project.org/web/views/Spatial.html

Integrating Python and R can be done in multiple ways, e.g. Using system calls or an in memory link using Rpy. See also:

Python interface for R Programming Language

answered 2012-07-30T15:37:35.897

Your Answer