在openair中不能用"timeAverage“函数计算风向

Question

我正在做一个使用openair项目的项目。具体来说，我想使用"timeAverage“函数在每天的时间步骤中平均数据集，但是在执行这样的任务时，将不会提取可变的风向。下面是我使用的R代码：

library("pmetar")
library("openair")
library("dplyr")
library("stringr")
NFTL<-metar_get_historical("NFTL", start_date = "2022-01-14", end_date = "2022-01-18", from = "iastate")
decoded_NFTL <- metar_decode(NFTL, metric = TRUE, altimeter = FALSE)

NFTL_obs<-select(decoded_NFTL, METAR_Date, Wind_direction)
NFTL_obs1 <- NFTL_obs 
NFTL_obs1$Wind_direction <- str_replace_all(NFTL_obs1$Wind_direction, 'Variable', 'NA')
NFTL_obs1$Wind_direction<-gsub(",.*","",NFTL_obs1$Wind_direction)
names(NFTL_obs1)[names(NFTL_obs1) == "METAR_Date"] <- "date"
names(NFTL_obs1)[names(NFTL_obs1) == "Wind_direction"] <- "wd"
daily <- timeAverage(NFTL_obs1, avg.time = "day")

在本例中，您可以检查在执行来自openair的最后一个命令时，没有提取风向(wd)变量，我如何解决这个问题？

Answer 1

将Wind_direction保留为字符的原因是，有时值是文本“变量”或“从n1到n2的变量”。您可以使用函数metar_dir作为numeric提取风向。

NFTL_obs$Wind_direction <- metar_dir(NFTL, numeric_only = TRUE)