R数据可视化－动态、交互式地图神器（一）概述与实现

这几天在网上看到一个了不起的R包：Leaflet/LeafletCN，忍不住跟着大神学习、动手实践一番，先看个效果图（交互式界面，可实现拖拽）：

下面是一些简单的用法与介绍

一、Leaflet/LeafletCN实例

Leaflet是一个国外动态交互地图，超赞；而LeafletCN适应了“中国国情”，划分更能细到县市级，虽然不一定完全准，但是够用就行啦~

回到前面的效果图，先看如何实现：

【1】载入包，读取网页数据（这里使用网上提供的各城市AQI指数）

library(XML)
library(leafletCN)
# Sorry for 爬了你家网站
table = readHTMLTable("http://www.pm25.in/rank",  encoding = "UTF-8", stringsAsFactors = F)[[1]]

【2】整理数据

# 整理数据并命名
dat = table[ , 2:3]
names(dat) = c("city","AQI")
dat$AQI = as.numeric(dat$AQI)

【3】使用geojsonMap绘制城市的污染情况

# 调用geojsonMap进行绘制
geojsonMap(dat, "city",
           popup =  paste0(dat$city, ":", dat$AQI),
           palette = "Reds", legendTitle = "AQI")

我去，丧心病狂，7行代码搞定 ~

二、函数使用

简单的COPY谁都会，各种函数、参数选择多少也要会一点，下面介绍常用函数

【regionNames函数】

如果传入省份级，会输出该省份下的城市

> regionNames("江苏")
 [1] "盐城市"   "徐州市"   "南通市"   "淮安市"   "苏州市"   "宿迁市"   "连云港市" "扬州市"  
 [9] "南京市"   "泰州市"   "无锡市"   "常州市"   "镇江市" 

若传入市级，则对应输出区县级

> regionNames("无锡")
[1] "北塘区" "滨湖区" "崇安区" "惠山区" "江阴市" "南长区" "锡山区" "宜兴市"

【geojsonmap函数】

dat = data.frame(name = regionNames("china"), value = runif(34))
geojsonMap(dat,"china")

# 可以看到有很多参数，默认颜色是蓝色，地图类型是高德地图，可自己调参
> geojsonMap
function (dat, mapName, namevar = NULL, valuevar = NULL, palette = "Blues", 
    colorMethod = "numeric", na.color = "#808080", popup = NULL, 
    stroke = T, smoothFactor = 1, weight = 1, fillOpacity = 0.7, 
    legendTitle = "Legend", tileType = amap, ...) 

三、其他实例

【1】单点标注：地图+标点

leaflet() %>%amap() %>%addMarkers(lng=121.587687, lat=31.201832,popup="张江老司机")

# %>%：管道操作符，它的主要作用就是把前面的语句（变量）传递给下一个语句，并且作为第一个参数使用

【2】底图地图多样选择

# 使用addProviderTiles("XXX")来调用这些源, 部分源需要注册
卫星地图:
Esri.WorldImagery
带海拔的地图:
OpenTopoMap
黑白地图:
Stamen.Toner
NASA气象地图:
NASAGIBS.ModisTerraChlorophyll
NASA夜景:
NASAGIBS.ViirsEarthAtNight2012

这么炫，简直没朋友啊

待续未完...