Interpoler les données des degrés à l'aide de QGIS?

Je veux créer des cartes de vent. Il existe une introduction à ArcGIS et à QGIS . Comme je n'ai pas accès à Spatial Analyst, j'ai essayé la version QGIS. Cependant, je ne sais pas comment interpoler les données de degré.

Par exemple, la direction du vent pourrait être de 359 ° à un point et de 1 ° au deuxième point. Les valeurs entre ce point devraient probablement être de 359 °, 0 ° ou 1 °.

Comme l'interpolation ne sait pas qu'il s'agit de données de degré, il semble que les valeurs interpolées soient d'environ 180 °.

Comment puis-je interpoler des données de degrés à l'aide de QGIS?

Interpoler le sinus et le cosinus de l'angle, puis reconvertir en un angle avec la fonction atan. Ces fonctions sont disponibles dans le moteur d'expression de QGIS. Il y a une atan2(dy,dx)fonction comme celle utilisée dans RI ci-dessous ...

Voici une fonction R pour illustrer. J'ai utilisé meanici pour donner l'interpolation:

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

qui donne:

> 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

Dans QGIS, en supposant que vous créez un raster interpolé à partir d'un ensemble de données ponctuelles, votre flux de travail sera le suivant:

• Créez deux nouveaux attributs à chaque point sinwindet coswind, à partir de votre direction du vent en degrés. Assurez-vous de convertir en radians si nécessaire.
• Créez deux rasters interpolés, un de sinwindet un de coswind.
• Utilisez la calculatrice raster pour revenir à la direction en calculant atan2les rasters sinus et cosinus interpolés.
• Convertissez en degrés si nécessaire.

Vous pouvez créer un workflow de traitement pour cela - il peut également y avoir déjà un plugin. Notez que bien que cela vous donne une interpolation de la direction du vent, toute évaluation de l' incertitude de cette interpolation est un peu délicate et nécessite des méthodes statistiques circulaires.