Menginterpolasi data derajat menggunakan QGIS?


9

Saya ingin membuat peta angin. Ada pengantar untuk ArcGIS dan untuk QGIS . Karena saya tidak memiliki akses ke Analis Tata Ruang, saya mencoba versi QGIS. Namun, saya tidak tahu bagaimana cara menginterpolasi data gelar.

Misalnya arah angin bisa 359 ° pada satu titik dan 1 ° pada titik kedua. Nilai antara titik ini mungkin harus 359 °, 0 °, atau 1 °.

Karena interpolasi tidak tahu bahwa ini adalah data derajat, sepertinya nilai yang diinterpolasi sekitar 180 °.

Bagaimana saya bisa menginterpolasi data tingkat menggunakan QGIS?

Jawaban:


7

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 meansini 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, sinwinddan coswind, dari arah angin Anda dalam derajat. Pastikan Anda mengonversi ke radian jika perlu.
  • Buat dua raster interpolasi, satu dari sinwinddan satu dari coswind.
  • Gunakan kalkulator raster untuk mengkonversi kembali ke arah dengan menghitung atan2raster 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.


Terima kasih atas jawaban itu, Spacedman. Bekerja dengan baik. Anda tidak punya referensi untuk itu? Akan lebih bagus untuk bisa mengatakan, ... disesuaikan menurut seseorang et al. (2000) ...
Alex

Makalah ini menyebutkan teknik, serta yang lain: pphmj.com/abstract/7133.htm dan mengutip referensi - Saya menduga buku-buku Mardia dan Cressie akan berguna, tetapi saya tidak punya mereka untuk menyerahkannya sekarang.
Spacedman
Dengan menggunakan situs kami, Anda mengakui telah membaca dan memahami Kebijakan Cookie dan Kebijakan Privasi kami.
Licensed under cc by-sa 3.0 with attribution required.