使用QGIS插值度数据？

我想创建风图。这里有关于ArcGIS和QGIS的介绍。由于我无权访问Spatial Analyst，因此尝试了QGIS版本。但是，我不知道如何对度数据进行插值。

例如，风向可以在一个点处为359°，在第二点处为1°。此点之间的值可能应为359°，0°或1°。

由于插值不知道它是度数据，因此看起来插值在180°左右。

如何使用QGIS插值学位数据？

插入角度的正弦和余弦，然后使用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插值正弦和余弦栅格，将其转换回方向。
• 如果需要，转换为度。

您可以为此创建处理工作流-可能已经有一个插件。请注意，尽管这为您提供了风向插值，但是对插值的不确定性进行任何评估都有些棘手，并且需要一些循环统计方法。