2012-07-23 18 views
9

tengo muchas líneas de datos hidrológicos georreferenciados con una resolución semanal:¿Cómo puedo interpolar datos georeferenciados en python?

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

Por desgracia, también tengo algunos datos con resolución única mensual:

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

En lugar de "reinventar la rueda , "¿alguien puede recomendar un módulo, paquete o técnica favorita que proporcione una interpolación razonable de los valores mensuales de los valores mensuales? Lineal estaría bien, pero estaría bien si pudiéramos usar las coordenadas para mejorar la interpolación basada en estaciones cercanas.

He etiquetado esta publicación con python porque es el idioma que he usado recientemente (aunque no son sus funciones estadísticas). Si la respuesta es "use un programa de estadísticas como r", que así sea, pero tengo curiosidad sobre qué hay para python. ¡Gracias!

Respuesta

4

no he tenido la oportunidad de profundizar en ella, pero la hpgl (High Performance Geoestadística Library) proporciona una serie de kriging (interpolación) geoespacial métodos:

Algoritmos

  • simple kriging (SK)
  • kriging ordinario (OK)
  • Indicador de kriging (IK)
  • local La variación media de Kriging (LVM Kriging)
  • Cokriging simple (Markov modelos 1 & 2)
  • Secuencial Indicador de simulación (SIS)
  • Corellogram local Variar Mean SIS (CLVM SIS)
  • Variar Local Mean SIS (LVM SIS)
  • secuencial gaussiana Simulación (SGS)
+2

Esto se ve un poco peluda - escasa en (en 2012!) documentación y ejemplos, pero aprecio su esfuerzo. Sería bueno saber de alguien que lo ha usado, por supuesto. –

+1

Probablemente encontré ese paquete de esta pregunta (no específica de Python): http://gis.stackexchange.com/questions/1041/open-source-methods-for-kriging – monkut

3

Si usted está interesado en ampliar su experiencia en I, hay un buen número de paquetes, así utilizadas y documentadas por ahí. Comenzaría por mirar Spatial Taskview, que enumera qué paquetes se pueden usar para datos espaciales. Uno de los párrafos trata de la interpolación. Estoy muy familiarizado con automap/gstat (escribí automap), donde especialmente gstat es un poderoso paquete de geoestadística que admite una amplia gama de métodos.

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

Integración Python y R se puede hacer de varias maneras, por ejemplo, Usar llamadas al sistema o un enlace en memoria usando Rpy. Ver también:

Python interface for R Programming Language

1

Estoy buscando en hacer lo mismo, y me encontré con this kriging module escrito por Sat Kumar Tomer en AMBHAS.

Parece haber métodos para producir variogramas y realizar kriging ordinarios.

Actualizaré esta respuesta si la utilizo y realizo más descubrimientos.

Cuestiones relacionadas