Interpolazione dei dati di laurea utilizzando QGIS?

Voglio creare mappe del vento. C'è un'introduzione per ArcGIS e per QGIS . Dato che non ho accesso all'analista spaziale, ho provato la versione di QGIS. Tuttavia, non so come interpolare i dati di laurea.

Ad esempio, la direzione del vento potrebbe essere 359 ° in un punto e 1 ° nel secondo punto. I valori tra questo punto dovrebbero probabilmente essere 359 °, 0 ° o 1 °.

Poiché l'interpolazione non sa che si tratta di dati di gradi, sembra che i valori interpolati siano di circa 180 °.

Come posso interpolare i dati di laurea usando QGIS?

Interpolare il seno e il coseno dell'angolo, quindi riconvertire in un angolo con la funzione atan. Queste funzioni sono disponibili nel motore di espressione di QGIS. C'è una atan2(dy,dx)funzione come quella in uso RI sotto ...

Ecco una funzione R per illustrare. Ho usato meanqui per dare l'interpolazione:

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

che dà:

> 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

In QGIS, supponendo che si stia creando un raster interpolato da un set di dati puntuali, il flusso di lavoro sarà il seguente:

• Crea due nuovi attributi in ciascun punto sinwinde coswind, dalla direzione del vento in gradi. Assicurati di convertire in radianti, se necessario.
• Crea due raster interpolati, uno da sinwinde uno da coswind.
• Usa la calcolatrice raster per tornare alla direzione calcolando atan2i raster seno e coseno interpolati.
• Converti in gradi se necessario.

È possibile creare un flusso di lavoro di elaborazione per questo - potrebbe anche esserci già un plug-in. Si noti che, sebbene ciò fornisca un'interpolazione della direzione del vento, qualsiasi valutazione dell'incertezza di tale interpolazione è un po 'complicata e richiede alcuni metodi statistici circolari.