Interpolasi sinus dan kosinus sudut, lalu konversikan kembali ke sudut dengan fungsi atan. Fungsi-fungsi ini tersedia di mesin ekspresi QGIS. Ada atan2(dy,dx)
fungsi seperti yang digunakan di RI di bawah ini ...
Inilah fungsi R untuk diilustrasikan. Saya telah menggunakan di mean
sini untuk memberikan interpolasi:
dinterp = function(d){
r=d*pi/180
(180/pi)*atan2(mean(sin(r)),mean(cos(r)))
}
pemberian yang mana:
> 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
Di QGIS, dengan asumsi Anda membuat raster yang diinterpolasi dari kumpulan data titik, alur kerja Anda adalah sebagai berikut:
- Buat dua atribut baru di setiap titik,
sinwind
dan coswind
, dari arah angin Anda dalam derajat. Pastikan Anda mengonversi ke radian jika perlu.
- Buat dua raster interpolasi, satu dari
sinwind
dan satu dari coswind
.
- Gunakan kalkulator raster untuk mengkonversi kembali ke arah dengan menghitung
atan2
raster sinus dan kosinus yang diinterpolasi.
- Konversikan ke derajat jika diperlukan.
Anda bisa membuat alur kerja pemrosesan untuk ini - mungkin juga sudah ada plugin. Perhatikan bahwa meskipun ini memberi Anda interpolasi dari arah angin, penilaian apa pun tentang ketidakpastian interpolasi itu agak rumit, dan memerlukan beberapa metode statistik melingkar.