同 RasterStack 小号 vr 和 ur
RasterStack
vr
ur
vr <- stack(system.file("external/rlogo.grd", package="raster")) ur <- flip(vr, 'y')
你确实可以使用你的公式:
windir <- atan2(vr, ur) * 180/pi + 180 winav <- (ur^2 + vr^2)^0.5
或者,您可以使用 overlay
overlay
windir2 <- overlay(vr, ur, fun=function(x,y) atan2(x,y) *180/pi + 180) winav2 <- overlay(vr, ur, fun=function(x, y) (x^2 + y^2)^0.5 )