QGIS kullanarak derece verilerini enterpolasyon mu?


9

Rüzgar haritaları oluşturmak istiyorum. ArcGIS ve QGIS için bir giriş var . Uzamsal Analiste erişimim olmadığından, QGIS sürümünü denedim. Ancak, derece verilerinin enterpolasyonunu nasıl yapacağımı bilmiyorum.

Örneğin rüzgar yönü bir noktada 359 ° ve ikinci noktada 1 ° olabilir. Bu nokta arasındaki değerler muhtemelen 359 °, 0 ° veya 1 ° olmalıdır.

İnterpolasyon, derece verisi olduğunu bilmediğinden enterpolasyonlu değerlerin 180 ° civarında olduğu görülmektedir.

QGIS kullanarak derece verilerini nasıl enterpolasyon yapabilirim?

Yanıtlar:


7

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 meanburada 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 sinwindve coswindrü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 sinwindve bir tane coswind.
  • atan2Enterpolasyonlu 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.


Bu cevap için teşekkürler Spacedman. İyi çalışıyor. Bunun için referansınız yok mu? Söyleyebilmek harika olurdu, ... birisine ve ark. (2000) ...
Alex

Bu makalede, diğer tekniklerin yanı sıra teknikten de bahsedilmektedir: pphmj.com/abstract/7133.htm ve atıflar - Mardia ve Cressie kitaplarının yararlı olacağından şüpheleniyorum, ancak şu anda teslim etmeleri gerekmiyor.
Spacedman
Sitemizi kullandığınızda şunları okuyup anladığınızı kabul etmiş olursunuz: Çerez Politikası ve Gizlilik Politikası.
Licensed under cc by-sa 3.0 with attribution required.