Interpolación de 17 puntos de datos - datos de temperatura

Question

Tengo 17 puntos de localización aquí en Islandia: puntos de estaciones meteorológicas.

Es un compuesto de 10 años, mes a mes, de los valores de temperatura de cada estación, lo que significa 120 valores para cada estación. 2040 valores en total. Como sólo tengo 17 estaciones sólo tengo 17 valores visibles en mi mapa.

Lo que quiero hacer es interpolar estos puntos para crear un mapa de temperatura para la zona y quiero poder compararlos, por ejemplo, comparar junio para cada año y hacer un mapa de eso - comparar junio-sep para cada año - comparar diciembre para 5 de esos 10 años, etc. Básicamente quiero ser capaz de comparar lo que quiero.

La distancia entre cada punto varía, pero la zona es bastante pequeña, por lo que sólo hay unos pocos kilómetros entre estaciones.

Cuando la interpolación esté hecha espero poder hacer una animación temporal a partir de los datos.

Se trata de un proyecto escolar. Asisto a la Universidad de Islandia.

Entonces, ¿qué método debo utilizar y cómo?

Gracias de antemano.

Einar Ingi

Answer 1

No soy un exportador cuantitativo, pero parece que Interpolación de splines podría ser una buena opción.

Las funciones polinómicas se consideran un buen método para la interpolación de elementos climáticos mensuales y anuales, pero son menos adecuadas en resoluciones temporales más altas, como días y horas

Este documento repasa los distintos métodos de interpolación y su relación con los datos de temperatura.

http://www.knmi.nl/bibliotheek/knmipubIR/IR2009-04.pdf

Answer 2

Respuesta similar a la de la obra de arte21.

Respondí a una pregunta hace tiempo y la misma respuesta parece aplicable aquí.

Cómo crear contornos para la zonificación

No he utilizado ArcGis, sin embargo he conseguido un resultado similar con SAGA GIS, que es una aplicación SIG de código abierto.

http://www.saga-gis.org/en/index.html

Utilizar Gridding para interpolar los datos de los puntos (B-Spline, Cubic Spline, Multi Level B-Spline...) y luego usar "Contour Lines from Grid" para convertir la malla en contornos.

Sé que no es una respuesta directa a tu pregunta pero te da opciones.

Básicamente, sigue las instrucciones anteriores pero ignora las líneas de contorno de la parte de Gird y tendrás un conjunto de datos de cuadrícula bien interpolados para cada mes, y se hace usando SAGA que es de código abierto.

Answer 3

Puede utilizar GDAL y su potente comando gdal_grid para crear cuadrículas regulares de datos dispersos.

Por ejemplo, este comando creará un ráster de interpolación en GTiff a partir de un archivo shape llamado "points.shp" almacenando los valores en un campo llamado "value" utilizando el algoritmo de la distancia inversa a una potencia:

$ gdal_grid points.shp output.tif -zfield value -ot Float64 -l points