無法在露天使用“timeAverage”function 計算風向

Question

我正在做一個使用 openair 項目的項目。 具體來說，我想使用“timeAverage” function 在每日時間步長對數據集進行平均，但在執行此類任務時，不會提取變量風向。 這是我正在使用的 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 保留為字符的原因是有時值是文本“Variable”或“variable from n1 to n2”。 您可以使用 function metar_dir將風向提取為numeric 。

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