之前转载过徐老师的一篇：三维人口密度分布图的制作和数据分享，感觉很有意思：
这种 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 柱形图
是不是非常炫酷！
