我有一个问题与R中的rasterVis中的vectorplot函数如何根据u和v风分量计算风向有关。我在下面的代码中看到它使用sqrt(u ^ 2 + v ^ 2)计算风速,但是我也想知道如何计算方向。
更具体地说,当风向向东(东风)时,我的纬向分量(U)为正;当风向北(北风)时,子午分量(V)为正。所以我的问题是vectorplot如何将U和V的方向转换为风来自何处的复合方向?这里提供了vectorplot的完整代码:https://rdrr.io/cran/rasterVis/src/R/vectorplot.R
else if (isTRUE(isField)) {
if (isTRUE(dXY)) {
u <- subset(object, 1)
v <- subset(object, 2)
object <- sqrt(u^2 + v^2)
}
else {
object <- subset(object, 1)
}
}