r - 如何从速度数据集创建矢量图？

Question

与我之前关于从 R 中的两个栅格对象计算 atan2 的问题相关吗？. 可以创建一个矢量图，其中速度（斜率）和方向（纵横比）用一组箭头（例如 Matlab 中的 quiver）显示，基于这个方程：

速度：ws <- sqrt(u^2+v^2)

方向：wd <- (180/pi)*(atan2(u,v))

我的数据集样本可以在这里下载

u <- raster(flname, varname = 'U')
v <- raster(flname, varname = 'V')
uv <- stack(u,v)

我使用 rasterVis 包来绘制绘图，但结果不合适。

library(rasterVis)
vectorplot(uv)

Answer 1

从帮助页面vectorplot：

如果 'isField='dXY'' 'object' 必须是一个 Raster*，有两个层分别代表水平和垂直分量。

因此，解决方案是：

vectorplot(uv, isField = 'dXY')