使用 R 语言计算各个城市距离北上广深的距离

根据培训班的一个小伙伴的需要，我再帮大家计算一下中国各个地级市距离北上广深的距离，首先加载相关的 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
#>    省代码 省     市代码 市       类型   距离北京市的距离  经度  纬度
#>                             
#>  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。