Come posso interpolare dati georeferenziati in python?

Question

Ho molte linee di dati idrologici georeferenziati con risoluzione settimanale:

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

Purtroppo, ho anche alcuni dati con una risoluzione di soli mensile:

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

Invece di "reinventare la ruota , "qualcuno può raccomandare un modulo, un pacchetto o una tecnica preferita che fornisca un'interpolazione ragionevole dei valori settimanali dai valori mensili? La linearità andrebbe bene, ma sarebbe bello se potessimo usare le coordinate per migliorare l'interpolazione basata sulle stazioni vicine.

Ho taggato questo post con python perché è la lingua che ho usato di recente (anche se non le sue funzioni statistiche). Se la risposta è "usa un programma di statistiche come r" così sia, ma sono curioso di sapere cosa c'è là fuori per Python. Grazie!

Answer 1

non ho avuto la possibilità di scavare in esso, ma il hpgl (High Geostatistica Biblioteca Performance) fornisce una serie di metodi kriging (geospaziale interpolazione):

Algoritmi

• Semplice Kriging (SK)
• Kriging ordinario (OK)
• Indicator Kriging (IK)
• locale Variando medio Kriging (LVM Kriging)
• semplice cokriging (Markov Models 1 & 2)
• sequenziale Indicatore di simulazione (SIS)
• Corellogram locale Variabile Media SIS (CLVM SIS)
• Variabile locale medio SIS (LVM SIS)
• sequenziale gaussiana Simulation (SGS)

Answer 2

Se siete interessati nella espandere la vostra esperienza in R, ci sono una serie di buoni pacchetti, ben utilizzati e documentati là fuori. Vorrei iniziare guardando la Spatial Taskview, che elenca quali pacchetti possono essere usati per i dati spaziali. Uno dei paragrafi riguarda l'interpolazione. Ho più familiarità con automap/gstat (ho scritto automap), dove in particolare gstat è un potente pacchetto di geostatistica che supporta una vasta gamma di metodi.

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

integrazione Python e R può essere fatto in diversi modi, ad esempio Usando le chiamate di sistema o un link in memoria usando Rpy. Vedi anche:

Python interface for R Programming Language

Answer 3

Sto cercando di fare la stessa cosa, e ho trovato this kriging module scritto da Sat Kumar Tomer a AMBHAS.

Sembrano esserci metodi per produrre variogrammi e eseguire normali kriging.

Aggiornerò questa risposta se uso questo e faccio ulteriori scoperte.

Answer 4

Da quando ho originariamente postato questa domanda un modulo Python Kriging attivamente sviluppato è stato rilasciato https://github.com/bsmurphy/PyKrige

C'è anche l'opzione più vecchio: https://github.com/capaulson/pyKriging