0

我正在尝试从 R 中的 rnoaa 获取一些天气数据。因为 rnoaa 仅支持一年的提取,我试图将一个循环放在一起以获得几年。使用地图功能更好吗?

它返回一个空白列表..

library(rnoaa)
options(noaakey= "somekey") 

washington_weather <- getweather("GHCND:USW00024234")
getweather <- function(stid) {
wtr<-0
for (i in 2009:2017) {
start_date <- paste0(i, "-01-01")
end_date <- paste0(i, "-12-31")
j<- i -2008
wtr[j]$tbl <- ncdc(datasetid='GHCND', stationid=stid, startdate = start_date, enddate = end_date)
}
return(wtr)
}

fahrenheit_to_celsius <- function(temp_F) {
  temp_C <- (temp_F - 32) * 5 / 9
  return(temp_C)
}
4

2 回答 2

1

rnoaa 包允许您组合使用该包获得的多个 ncdc 对象。

如果您使用ncdc_combine()函数,则可以组合您创建的多个对象。

例如:

x <- ncdc(datasetid= "GHNCD", stationid=stid, startdate = start_date, enddate = end_date)
y <- ncdc(datasetid= "GHNCD", stationid=stid, startdate = start_date1, enddate = end_date1)

z <- ncdc_combine(x,y)

这将结合您的两个 ncdc 对象,只要您将每个对象分解为不到一年。

于 2020-08-07T18:17:24.923 回答
0

ncdc函数的返回值是一个列表。理想情况下,您只想返回列表中的数据部分。

在这个脚本中,我下载了每年的数据并将数据部分信息保存在一个列表中。然后可以使用 data.frames 列表进行额外分析或将所有数据帧绑定到一个大数据帧中。

getweather <- function(stid) {
   wtr<-list()  # create an empty list
   for (i in 2009:2011) {
      start_date <- paste0(i, "-01-01")
      end_date <- paste0(i, "-12-31")
      
      #save data portion to the list (elements named for the year
      wtr[[as.character(i)]] <- ncdc(datasetid='GHCND', stationid=stid, startdate = start_date, enddate = end_date)$data
   }
   #return the full list of data frames
   return(wtr)
}


washington_weather <- getweather("GHCND:USW00024234")

#bind the dataframes in the list together into one large dataframe
dplyr::bind_rows(washington_weather)
于 2020-06-26T18:39:01.553 回答