2020 年中国人口密度 3D 蜂窝地图（使用 R 语言绘制 3

之前转载过徐老师的一篇：三维人口密度分布图的制作和数据分享，感觉很有意思：

这种 3D 柱形地图也可以使用 R 语言绘制的，今天就给大家介绍一下这种 3D 柱形地图的绘制方法：

2020 年中国人口地理分布 3D 柱形图

首先打开 RStudio，加载下面的 R 包：

library(mapdeck)
library(tidyverse)
library(raster)
library(sf)

设置 mapbox token，我使用的是默认的 "pk.eyJ1IjoiY3p4YSIsImEiOiJjazMwYTZndGUwNGI5M2RsZ20ydTM2azhoIn0.gvP_XXYhLV7hkD0e3-9AKw"

set_token(Sys.getenv("MAPBOX"))
# 或者，如果你没有把 token 添加到环境变量，可以：
set_token("pk.eyJ1IjoiY3p4YSIsImEiOiJjazMwYTZndGUwNGI5M2RsZ20ydTM2azhoIn0.gvP_XXYhLV7hkD0e3-9AKw")

Mapbox 中国地图有三种：

样式	style
街景地图	mapbox://styles/mapbox/streets-zh-v1
暗黑地图	mapbox://styles/mapbox/dark-zh-v1
浅色地图	mapbox://styles/mapbox/light-zh-v1

使用暗黑模式会比较炫酷：

ms = "mapbox://styles/mapbox/dark-zh-v1"

读取 2020 年中国人口密度栅格数据：

raster("每平方公里人口密度Density/Density2020.tif") -> pop

通过聚合降低分辨率：

aggregate(pop, fact = 4, fun = mean) -> pop2

栅格数据转换成点数据（经度、纬度、人口密度）：

pop2 %>% 
  rasterToPoints() %>% 
  as_tibble() -> popdf
  
popdf
#> # A tibble: 871,421 x 3
#>        x     y Density2020
#>    <dbl> <dbl>       <dbl>
#>  1  123.  53.5      1.07  
#>  2  123.  53.5      0.854 
#>  3  123.  53.5      0.887 
#>  4  123.  53.5      1.00  
#>  5  123.  53.5      1.18  
#>  6  123.  53.5      0.350 
#>  7  123.  53.5      0.680 
#>  8  123.  53.5      0.309 
#>  9  123.  53.5      0.838 
#> 10  123.  53.5      0.0453
#> # … with 871,411 more rows

由于 Mapbox 底图上没有九段线，所以我们自己添加一下：

sf::read_sf("九段线.geojson") -> jdx

然后就可以绘图了：

mapdeck(style = ms) %>%
  add_path(data = jdx, 
           stroke_colour = "#fef0d9") %>% 
  add_hexagon(data = popdf, 
              lat = "y", lon = "x", 
              elevation = "Density2020",
              colour = "Density2020",
              elevation_scale = 500, 
              layer_id = "grid_layer",
              colour_range = c("#fef0d9", "#fdd49e", "#fdbb84", "#fc8d59", "#e34a33", "#b30000")) %>% 
  mapdeck_view(pitch = 45, 
               location = c(104, 36), 
               zoom = 3)

效果如下：

2020 年中国人口地理分布 3D 柱形图

是不是非常炫酷！