Interpolar datos de grado usando QGIS?

9

Quiero crear mapas de viento. Hay una introducción para ArcGIS y para QGIS . Como no tengo acceso al Spatial Analyst, probé la versión QGIS. Sin embargo, no sé cómo interpolar datos de grado.

Por ejemplo, la dirección del viento podría ser 359 ° en un punto y 1 ° en el segundo punto. Los valores entre este punto probablemente deberían ser 359 °, 0 ° o 1 °.

Como la interpolación no sabe que se trata de datos de grados, parece que los valores interpolados están alrededor de 180 °.

¿Cómo puedo interpolar datos de grado usando QGIS?

Alex
fuente

Respuestas:

7

Interpolar el seno y el coseno del ángulo, y luego convertir de nuevo a un ángulo con la función atan. Estas funciones están disponibles en el motor de expresión de QGIS. Hay una atan2(dy,dx)función como la que utilizo en RI a continuación ...

Aquí hay una función R para ilustrar. He usado meanaquí para dar la interpolación:

 dinterp = function(d){
             r=d*pi/180
             (180/pi)*atan2(mean(sin(r)),mean(cos(r)))
           }

lo que da:

> dinterp(c(10,20))
[1] 15                *as expected*
> dinterp(c(340,20))
[1] 3.384676e-15      *approximately zero*
> dinterp(c(0,180))   
[1] 90                *could be +/- 90
> dinterp(c(359,180))
[1] -90.5             *small adjustment makes it go to ~-90

En QGIS, suponiendo que está creando un ráster interpolado a partir de un conjunto de datos de puntos, su flujo de trabajo será el siguiente:

  • Cree dos nuevos atributos en cada punto sinwindy coswind, a partir de la dirección del viento en grados. Asegúrese de convertir a radianes si es necesario.
  • Cree dos rásteres interpolados, uno desde sinwindy otro desde coswind.
  • Use la calculadora ráster para volver a convertir en dirección calculando atan2los rásteres seno y coseno interpolados.
  • Convierte a grados si es necesario.

Puede crear un flujo de trabajo de procesamiento para esto; también puede haber un complemento. Tenga en cuenta que aunque esto le proporciona una interpolación de la dirección del viento, cualquier evaluación de la incertidumbre de esa interpolación es un poco complicada y requiere algunos métodos de estadísticas circulares.

Hombre espacial
fuente
Gracias por esa respuesta Spacedman. Funciona bien. ¿No tienes algunas referencias para eso? Sería genial poder decir, ... se ajustó de acuerdo con alguien et al. (2000) ...
Alex
Este documento menciona la técnica, así como otros: pphmj.com/abstract/7133.htm y cita referencias: sospecho que los libros de Mardia y Cressie serán útiles, pero no los tengo a mano en este momento.
Spacedman