式を実装する方法2つのラスターグリッドセルを実装する方法（各グリッドセルにuとvの風力コンポーネントを使用して平均風を計算する）

Question

風速のuとv成分を持つ2つのラスタスタックファイルurとvrがあります。 各グリッドセルの平均風速と風向はどのように計算できますか？ は、私が使用して、以下の式

windir<-calc(atan2(vwind, uwind) * 360/2/pi) + 180 

と平均風を利用して風の方向を算出することができる知っている

winav<- ((mean.u^2 + mean.v^2)^0.5) 

私の問題は私が呼ばurとvr の各グリッドセルにこれらの式を実装する方法でありますthis questionおよび他のリンクは、この計算ではまだストックされています。あなたが本当にあなたの数式を使用することができますRasterStackのvrとur

vr <- stack(system.file("external/rlogo.grd", package="raster")) 
ur <- flip(vr, 'y') 

で

Answer 1

：

windir <- atan2(vr, ur) * 180/pi + 180 

winav <- (ur^2 + vr^2)^0.5 

また、あなたは平均的な風の方向を算出することを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)