由于罗盘方向范围从0到360度,因此对标量风向数据进行平均会产生不准确的值,所以我已经将我的列表从大小和风向角转换为u和v分量。
为了获得正确的风向,以便进行平均,我需要为以下3种情况开发某种应用、ifelse函数:
V > 0...((180 / pi) * atan((Ucomp/Vcomp)) + 180)
U and V < 0...((180 / pi) * atan((Ucomp/Vcomp)) + 0)
U > 0 and V < 0...((180 / pi) * atan((Ucomp/Vcomp)) + 360)
在我要分析的数据集中,Ucomp大于0且Vcomp小于零,但无疑会有所有3个情况都会出现的情况,因此我需要一个函数来逐步解析和计算,并在每个时间步骤应用正确的公式。我以前没有使用过lapply或函数,因此我尝试使用它们没有成功。
我提供以下数据样本...
DateTime Wind.Spd Wind.Direction Air.Density Temp.C GEP.GE16XLE GCF.GE16XLE Ucomp Vcomp
1 1981 7.662370 248.3395 0.9148207 11.28967 597.7513 37.35946 5.253453 -0.7404972
2 1982 8.199412 251.6763 0.9172176 10.12751 678.8595 42.42872 5.867979 -0.6191475
3 1983 8.188782 251.7889 0.9162767 10.30619 667.9461 41.74663 5.777208 -1.0473982
4 1984 7.942632 246.7908 0.9174074 10.05093 642.6374 40.16484 5.415773 -0.6796723
5 1985 8.016558 252.7305 0.9171721 10.38414 654.2588 40.89117 5.649406 -0.9999082
6 1986 7.739984 249.6431 0.9158740 10.99859 607.0542 37.94089 5.305971 -0.9118965