这几天在网上看到一个了不起的R包:Leaflet/LeafletCN,忍不住跟着大神学习、动手实践一番,先看个效果图(<small>交互式界面,可实现拖拽</small>):
下面是一些简单的用法与介绍
一、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
这么炫,简直没朋友啊
待续未完...
网友评论