2020 年中国人口密度 3D 蜂窝地图

2020 年中国人口密度 3D 蜂窝地图

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

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

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

首先打开 RStudio,加载下面的 R 包:

1
2
3
4
library(mapdeck)
library(tidyverse)
library(raster)
library(sf)

设置 mapbox token,我使用的是默认的 “pk.eyJ1IjoiY3p4YSIsImEiOiJjazMwYTZndGUwNGI5M2RsZ20ydTM2azhoIn0.gvP_XXYhLV7hkD0e3-9AKw”

1
2
3
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

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

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

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

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

通过聚合降低分辨率:

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

栅格数据转换成点数据(经度、纬度、人口密度):

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
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 底图上没有九段线,所以我们自己添加一下:

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

然后就可以绘图了:

1
2
3
4
5
6
7
8
9
10
11
12
13
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 柱形图

是不是非常炫酷!

阅读更多:https://mp.weixin.qq.com/s/N06zVKdVIB3VxW1ySjQrBQ
下载代码和数据:https://rstata.duanshu.com/#/brief/course/8f921af271a04b6db79c52ce48a04d49

2020 年中国人口密度 3D 蜂窝地图

https://tidyfriday.cn/posts/57187/

作者

Painter

发布于

2021-01-23

更新于

2021-02-09

许可协议

评论