if (FALSE) {
# Get station table
(stations <- isd_stations())
## plot stations
### remove incomplete cases, those at 0,0
df <- stations[complete.cases(stations$lat, stations$lon), ]
df <- df[df$lat != 0, ]
### make plot
library("leaflet")
leaflet(data = df) %>%
addTiles() %>%
addCircles()
}
Run the code above in your browser using DataLab