r - 緯度/経度に基づいて気象データをプルする R の関数 - RNOAA パッケージ

Question

気象観測所に似た緯度と経度に基づいて、R で気象観測の種類のリストを作成しています。

## List of airports you want to include in your weather extract

airport_list <- c("KABE" , "KBWI", "KRAL")

## Drilldown of your airport locations ( I have a separate table that pulls in this info

airport_list_dd <- airport_locs[airport_locs$icao %in% airport_list,]

## Mutate the data frame to make lat/lon compatible for searching the    NOAA GHCND 

airport_list_similar <- airport_list_dd %>%
  mutate(lon_similar = str_extract(longitude, "([-0-9]+)\\.."),
         lat_similar = str_extract(latitude, "([-0-9]+)\\.."),
         lon_exact = str_extract(longitude, "([-0-9]+)\\....."),
         lat_exact = str_extract(latitude, "([-0-9]+)\\....."))

## Define your date range

 date_min <- Sys.Date() - (10 * 365)
 date_max <- Sys.Date()
 filter_year <- year(Sys.Date()) - 1

# THIS IS WHERE I AM HAVING THE ISSUE
 ## Build your weather extracts
 ghcnd_near_airport <- list()
 build_lon_table <- function(x){
  i <- 1
  for (i in 1:length(x)) {
    lon_similar <- x$lon_similar[i]
    lat_similar <- x$lat_similar[i]
    ghcnd_near_airport <- c(ghcnd_stations %>%
                                filter(str_detect(longitude, lon_similar), str_detect(latitude, lat_similar)), list(i))
  }
  return(ghcnd_near_airport)
}

ただし、これはカウントが 11 の空のリストを返します。これは、関数を適切な時間反復処理しているが、リスト内のデータを返していないことを意味します。