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 mean
aquí 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
sinwind
y 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
sinwind
y otro desde coswind
.
- Use la calculadora ráster para volver a convertir en dirección calculando
atan2
los 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.