根据培训班的一个小伙伴的需要，我再帮大家计算一下中国各个地级市距离北上广深的距离，首先加载相关的 R 包：
library(tidyverse)
library(sf)
读取各个城市的地理数据：
read_sf("2019行政区划/市.shp") -> city
计算质心：
city %>%
st_centroid() -> city_centroid
先以北京市为例，计算所有的地级市距离北京的距离：
# 计算所有地级市距离北京的距离
st_distance(
city_centroid,
city_centroid %>%
dplyr::filter(市 == "北京市")
) -> mat
mat[,1]
#> Units: [m]
#> [1] 0.0 126702.0 284720.8 172266.5 236542.1 435058.8 358031.5
#> [8] 167576.3 139156.7 160799.7 215046.3 104499.9 273788.2 432058.3
#> [15] 231742.6 344875.2 512764.2 603197.9 328554.4 436454.0 727999.2
...
#> [358] 2368608.1 2533359.2 2960674.1 3168140.8 2412249.6 2405613.5 2582123.9
#> [365] 2825578.0 2963079.5 3211553.9 2632945.0 1872284.5 1985789.5 2016356.2
city_centroid %>%
mutate(dist = units::set_units(mat[,1], km),
dist = as.numeric(dist)) -> df
bind_cols(
df %>%
st_drop_geometry(),
df %>%
st_coordinates() %>%
as_tibble() %>%
set_names(c("经度", "纬度"))
) %>%
rename(距离北京市的距离 = dist)
#> # A tibble: 371 x 8
#> 省代码 省 市代码 市 类型 距离北京市的距离 经度 纬度
#> <dbl> <chr> <dbl> <chr> <chr> <dbl> <dbl> <dbl>
#> 1 110000 北京市 110000 北京市 直辖市 0 116. 40.2
#> 2 120000 天津市 120000 天津市 直辖市 127. 117. 39.3
#> 3 130000 河北省 130100 石家庄市 地级市 285. 114. 38.1
#> 4 130000 河北省 130200 唐山市 地级市 172. 118. 39.7
#> 5 130000 河北省 130300 秦皇岛市 地级市 237. 119. 40.1
#> 6 130000 河北省 130400 邯郸市 地级市 435. 115. 36.6
#> 7 130000 河北省 130500 邢台市 地级市 358. 115. 37.2
#> 8 130000 河北省 130600 保定市 地级市 168. 115. 39.0
#> 9 130000 河北省 130700 张家口市 地级市 139. 115. 40.9
#> 10 130000 河北省 130800 承德市 地级市 161. 118. 41.3
#> # … with 361 more rows
然后我们可以使用 lapply() 计算距离北上广深四个城市的距离并分别保存为四个文件：
lapply(c("北京市", "上海市", "广州市", "深圳市"), function(x){
st_distance(
city_centroid,
city_centroid %>%
dplyr::filter(市 == x)
) -> mat
mat[,1]
city_centroid %>%
mutate(dist = units::set_units(mat[,1], km),
dist = as.numeric(dist)) -> df
bind_cols(
df %>%
st_drop_geometry(),
df %>%
st_coordinates() %>%
as_tibble() %>%
set_names(c("经度", "纬度"))
) %>%
writexl::write_xlsx(paste0("各个地级市质心距离", x, "质心的距离.xlsx"))
})
#> [[1]]
#> [1] "/Users/ac/Desktop/各个地级市质心距离北京市质心的距离.xlsx"
#>
#> [[2]]
#> [1] "/Users/ac/Desktop/各个地级市质心距离上海市质心的距离.xlsx"
#>
#> [[3]]
#> [1] "/Users/ac/Desktop/各个地级市质心距离广州市质心的距离.xlsx"
#>
#> [[4]]
#> [1] "/Users/ac/Desktop/各个地级市质心距离深圳市质心的距离.xlsx"
这样就完成了，所有文件里面的距离单位都是 km。
