R实现pm2.5地图数据展示

library(rvest)
library(leafletCN)
Sys.setlocale("LC_CTYPE","eng")
doc = read_html("http://www.pm25s.com/cn/rank/")
cities = doc %>% html_nodes(".cityrank a") %>%
  html_text()
AQI = doc %>% html_nodes("span[class^='lv']") %>%
  html_text() %>% .[c(F,F,T)] %>% as.numeric
dat = data.frame(city=cities,AQI=AQI)
dat
geojsonMap(dat,"city",
           popup = paste0(dat$city,":",dat$AQI),
           palette = "Reds",legendTitle = "AQI")