Jak mogę interpolować dane georeferencyjne w pythonie?

Question

Mam wiele wierszy georeferencyjnych danych hydrologicznych z tygodniowym rozdzielczość:

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

Niestety, ja też mam pewne dane z jedynego miesięcznego uchwały:

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

zamiast „wyważania otwartych drzwi , "czy ktoś może polecić ulubiony moduł, pakiet lub technikę, która zapewni rozsądną interpolację cotygodniowych wartości z wartości miesięcznych? Liniowy będzie dobrze, ale byłoby miło, gdybyśmy mogli użyć współrzędnych, aby poprawić interpolację na podstawie pobliskich stacji.

Tagowałem ten post w python, ponieważ jest to język, którego ostatnio używałem (choć nie jest to jego funkcja statystyczna). Jeśli odpowiedź brzmi "użyj programu statystycznego, takiego jak r", niech to będzie, ale jestem ciekawy, co tam jest dla Pythona. Dzięki!

Answer 1

Nie miałem okazję zagłębić się nim, ale hpgl (High Performance Geostatystyka Library) zawiera szereg kriging (geoprzestrzennych interpolacji) metod:

Algorytmy

• Proste kriging (SK)
• Zwyczajne kriging (OK)
• Wskaźnik kriging (IK)
• lokalny Różnej Mean kriging (LVM kriging)
• Proste CoKriging (modele Markowa 1 & 2)
• Sequential Wskaźnik Simulation (SIS)
• Corellogram Local Różnie Mean SIS (CLVM SIS)
• Local Różnie Mean SIS (LVM SIS)
• Sequential Gaussa Simulation (SGS)

Answer 2

Jeśli są Państwo zainteresowani w rozwija swoje doświadczenia do badań, istnieje wiele dobrych, dobrze wykorzystywane i udokumentowanych pakietów tam. Zacznę od obejrzenia Spatial Taskview, który wyszczególnia, jakie pakiety można wykorzystać do danych przestrzennych. Jeden z akapitów dotyczy interpolacji. Jestem zaznajomiony z automap/gstat (napisałem automap), gdzie szczególnie gstat to potężny pakiet geostatystyczny, który obsługuje szeroki zakres metod.

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

Integracja Pyton i R mogą być wykonane na różne sposoby, na przykład Używanie wywołań systemowych lub łącza do pamięci za pomocą Rpy. Zobacz także:

Python interface for R Programming Language

Answer 3

szukam do robi to samo, i znalazłem this kriging module napisany przez Sat Kumar Tomer w AMBHAS.

Wydaje się, że istnieją metody wytwarzania variogramów i wykonywania krigingu zwyczajnego.

Zaktualizuję tę odpowiedź, jeśli użyję tego i dokonam dalszych odkryć.

Answer 4

Ponieważ pierwotnie opublikowany na to pytanie aktywnie rozwinięty moduł Pythona kriging został wydany https://github.com/bsmurphy/PyKrige

Jest też to starsza opcja: https://github.com/capaulson/pyKriging