Интерполяция данных степени с использованием QGIS?

9

Я хочу создавать карты ветра. Существует введение для ArcGIS и QGIS . Поскольку у меня нет доступа к Spatial Analyst, я попробовал QGIS-версию. Тем не менее, я не знаю, как интерполировать данные степени.

Например, направление ветра может быть 359 ° в одной точке и 1 ° во второй точке. Значения между этой точкой, вероятно, должны быть 359 °, 0 ° или 1 °.

Поскольку интерполяция не знает, что это данные в градусах, похоже, что интерполированные значения составляют около 180 °.

Как я могу интерполировать данные степени, используя QGIS?

Alex
источник

Ответы:

7

Интерполировать синус и косинус угла, а затем преобразовать обратно в угол с помощью функции atan. Эти функции доступны в механизме выражений QGIS. Существует atan2(dy,dx)функция , как в использовании RI ниже ...

Вот функция R для иллюстрации. Я использовал meanздесь, чтобы дать интерполяцию:

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

который дает:

> 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, если вы создаете растр, интерполированный из набора точечных данных, ваш рабочий процесс будет следующим:

  • Создайте два новых атрибута в каждой точке sinwindи coswindот вашего направления ветра в градусах. Убедитесь, что вы переводите в радианы, если это необходимо.
  • Создайте два интерполированных растра, один из sinwindи один из coswind.
  • Используйте калькулятор растра, чтобы преобразовать обратно в направление путем вычисления atan2интерполированных растров синуса и косинуса.
  • Преобразовать в градусы, если требуется.

Для этого вы можете создать рабочий процесс обработки - там также может быть плагин. Обратите внимание, что хотя это дает вам интерполяцию направления ветра, любая оценка неопределенности этой интерполяции немного сложна и требует некоторых методов круговой статистики.

Spacedman
источник
Спасибо за этот ответ Spacedman. Работает отлично. У вас случайно нет ссылок на это? Было бы замечательно, если бы можно было сказать, что ... был настроен в соответствии с кем-то и соавт. (2000) ...
Алекс
В этой статье упоминается техника, а также другие: pphmj.com/abstract/7133.htm и приводятся ссылки - я подозреваю, что книги Mardia и Cressie будут полезны, но сейчас я не могу их передать.
Spacedman