根據(jù)培訓(xùn)班的一個小伙伴的需要，我再幫大家計(jì)算一下中國各個地級市距離北上廣深的距離，首先加載相關(guān)的 R 包：

library(tidyverse)
library(sf)

讀取各個城市的地理數(shù)據(jù)：

read_sf("2019行政區(qū)劃/市.shp") -> city
計(jì)算質(zhì)心：

city %>% 
  st_centroid() -> city_centroid

先以北京市為例，計(jì)算所有的地級市距離北京的距離：

# 計(jì)算所有地級市距離北京的距離
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("經(jīng)度", "緯度"))
) %>% 
  rename(距離北京市的距離 = dist)
#> # A tibble: 371 x 8
#>    省代碼 省     市代碼 市       類型   距離北京市的距離  經(jīng)度  緯度
#>     <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() 計(jì)算距離北上廣深四個城市的距離并分別保存為四個文件：

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("經(jīng)度", "緯度"))
  ) %>% 
    writexl::write_xlsx(paste0("各個地級市質(zhì)心距離", x, "質(zhì)心的距離.xlsx"))
})
#> [[1]]
#> [1] "/Users/ac/Desktop/各個地級市質(zhì)心距離北京市質(zhì)心的距離.xlsx"
#> 
#> [[2]]
#> [1] "/Users/ac/Desktop/各個地級市質(zhì)心距離上海市質(zhì)心的距離.xlsx"
#> 
#> [[3]]
#> [1] "/Users/ac/Desktop/各個地級市質(zhì)心距離廣州市質(zhì)心的距離.xlsx"
#> 
#> [[4]]
#> [1] "/Users/ac/Desktop/各個地級市質(zhì)心距離深圳市質(zhì)心的距離.xlsx"

這樣就完成了，所有文件里面的距離單位都是 km。