leaflet的小搭档leaflet.minicharts来了，从此动态地图又多了一些乐趣~~~

本文内容取材自leaflet.minicharts包官方主页的案例介绍，本篇案例虽然是关于leaflet在线地图的辅助包，但是该包的出现对于leaflet生态系统来说，确是有着划时代的意义。

该包大大扩充了leaflet包所能呈现的图表形式，打破了散点图、路径图、热力图三类图表对于传统地图数据呈现形式的垄断地位。

该包为此提供了气泡饼图、玫瑰图、mini柱形图以及时间维度控制器等多种数据可视化元素。

以下便是作者对该包的简要介绍及案例演示：

leaflet.minicharts包简介：

一直以来，借助于Rstudio团队开发的交互式地图工具包——leaflet，在R语言中创建交互式地图已非难事。但是该包仅仅提供了很少的集中图表呈现类型，因而在一张单独的地图上呈现的数据信息量被大大限制了。

倘若你要在地图上呈现一些点数据，你只能需要使用两个变量定位出这些点在地图上的位置，然后更改点的半径和点的填充颜色来呈现数据。
然后leaflet.minicharts包的出现大大改变了这一格局。你可以在利用其提供的两个附加函数，在leaflet包的交互地图上增加更多的mini图表。

这些mini图表可以呈现比之前多得多的变量数据（理论上不受变量数量的限制），目前为止，仅支持三种图表类型：

barcharts (the default), pie charts and polar area charts.

数据预览：

该包内提供了一个法国自2010~2017年包含电力产品、电力消费及结构变动趋势的数据集。

library("ggplot2")
library("leaflet.minicharts")
library(geojsonio)
library(rgdal)
data("eco2mix")
head(eco2mix)

地图数据及作图数据整理：

geojson <-readOGR("D:/R/mapdata/State/France.json","OGRGeoJSON",stringsAsFactors=FALSE)
library(dplyr)
prod2016 <- eco2mix %>%
  mutate(
    renewable = bioenergy + solar + wind + hydraulic,
    non_renewable = total - bioenergy - solar - wind - hydraulic
  ) %>%
  filter(grepl("2016",month) & area!= "France") %>%
  select(-month) %>%
  group_by(area,lat,lng) %>%
  summarise_all(sum) %>%
  ungroup()
``` 
  
### 制作背景地图：
抽取了法国2016年的能源消费结构数据信息：
```r
library(leaflet)
library(leafletCN)
basemap <- leaflet(geojson,width ="100%",height="400px") %>% amap() %>%
addPolygons(weight=1,color="#444444",opacity=1,fillColor ="white",fillOpacity =0.7,smoothFactor=0.5)

默认图表（柱形图）

renewable2016 <- prod2016 %>% select(hydraulic,solar,wind)
colors <- c("#3093e5","#fcba50","#a0d9e8")
basemap %>%
  addMinicharts(
    prod2016$lng, prod2016$lat,
    chartdata = renewable2016,
    colorPalette = colors,
    width = 45, height = 45
  )

气泡饼图

colors <- c("#4fc13c", "#cccccc")
basemap %>% 
  addMinicharts(
    prod2016$lng, prod2016$lat,
    type = "pie",
    chartdata = prod2016[, c("renewable", "non_renewable")], 
    colorPalette = colors, 
    width = 60 * sqrt(prod2016$total)/sqrt(max(prod2016$total)),transitionTime=0)

玫瑰图1：

basemap %>%
addMinicharts(
prod2016$lng, prod2016$lat,
type = "polar-radius",
chartdata = prod2016[, c("renewable", "non_renewable")],
colorPalette = colors,
width = 60 * sqrt(prod2016$total)/sqrt(max(prod2016$total)),transitionTime=0)

玫瑰图2：

basemap %>%
addMinicharts(
prod2016$lng, prod2016$lat,
type = "polar-area",
chartdata = prod2016[, c("renewable", "non_renewable")],
colorPalette = colors,
width = 60 * sqrt(prod2016$total)/sqrt(max(prod2016$total)),transitionTime=0)

添加时间趋势：

柱形图（带时间趋势）

prodRegions <- eco2mix %>% filter(area != "France")
basemap %>% 
  addMinicharts(
    prodRegions$lng, prodRegions$lat, 
    chartdata = prodRegions[, c("hydraulic", "solar", "wind")],
    time = prodRegions$month,
    colorPalette = colors,
    width = 45, height = 45
  )

饼图（带时间趋势）

prodRegions <- eco2mix %>% filter(area != "France")
basemap %>% 
  addMinicharts(
    prodRegions$lng, prodRegions$lat, 
    chartdata = prodRegions[, c("hydraulic", "solar", "wind")],
    time = prodRegions$month,
    colorPalette = colors,
    type="pie",
    width = 45, height = 45
  )

玫瑰图（带时间趋势）

prodRegions <- eco2mix %>% filter(area != "France")
basemap %>% 
  addMinicharts(
    prodRegions$lng, prodRegions$lat, 
    chartdata = prodRegions[, c("hydraulic", "solar", "wind")],
    time = prodRegions$month,
    colorPalette = colors,
    type="pie",
    width = 45, height = 45
  )

联系方式：
wechat：ljty1991
Mail:578708965@qq.com
个人公众号：数据小魔方（datamofang）

qq交流群：[魔方学院]298236508

个人简介：
杜雨
财经专业研究僧；
伪数据可视化达人；
文科背景的编程小白；
喜欢研究商务图表与地理信息数据可视化，爱倒腾PowerBI、SAP DashBoard、Tableau、R ggplot2、Think-cell chart等诸如此类的数据可视化软件，创建并运营微信公众号“数据小魔方”。
Mail:578708965@qq.com