Açının sinüsü ve kosinüsü enterpolasyonu yapın ve sonra atanan fonksiyon ile bir açıya geri dönün. Bu işlevler QGIS'in ifade motorunda mevcuttur. atan2(dy,dx)
Aşağıda RI kullanımındaki gibi bir fonksiyon var ...
İşte açıklamak için bir R işlevi. Ben mean
burada enterpolasyonu verdim :
dinterp = function(d){
r=d*pi/180
(180/pi)*atan2(mean(sin(r)),mean(cos(r)))
}
hangi verir:
> 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
QGIS'de, bir nokta veri kümesinden enterpolasyonlu bir raster oluşturduğunuzu varsayarsak, iş akışınız aşağıdaki gibi olacaktır:
- Her noktada
sinwind
ve coswind
rüzgar yönünüzden derece olarak iki yeni özellik oluşturun . Gerekirse radyana dönüştürdüğünüzden emin olun.
- İki interpolasyon rasterları, bir tane oluşturun
sinwind
ve bir tane coswind
.
atan2
Enterpolasyonlu sinüs ve kosinüs rasterlerini hesaplayarak yönüne geri dönüştürmek için raster hesap makinesini kullanın .
- Gerekirse dereceye dönüştürün.
Bunun için bir işleme iş akışı oluşturabilirsiniz - zaten bir eklenti de olabilir. Bu size rüzgar yönünün enterpolasyonunu vermesine rağmen, bu enterpolasyonun belirsizliğinin herhangi bir değerlendirmesinin biraz zor olduğunu ve bazı dairesel istatistik yöntemleri gerektirdiğini unutmayın.