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

Question

与我之前关于calculate atan2 from two raster object in R?的问题相关。可以创建一个矢量图，其中基于以下等式用一组箭头（例如Matlab中的箭头）显示速度（斜率）和方向（方位）：

表示速度：ws＆lt; - sqrt（u ^ 2 + v ^ 2）

方向：wd＆lt; - （180 / pi）*（atan2（u，v））

我的数据集示例可以是downloaded here

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

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

library(rasterVis)
vectorplot(uv)

Answer 1

来自vectorplot的帮助页面：

  

如果'isField ='dXY'''对象'必须是带有两层的光栅*   代表水平和垂直分量，   分别

因此，解决方案是：

vectorplot(uv, isField = 'dXY')