使用R语言计算坐标点到海岸线的距离|利用sf包

去年一个项目中的一段代码，分享出来给需要的同学。

library(tidyverse)
library(sf)
library(data.table)

#Read in registry agency list and convert the coordinates
region_list <- na.omit(read.csv(file = "./data/CI5-XId/converted/registry_info.csv", 
                                header =T)[,c(1,2,5,6)])
region_list_sf <- region_list %>% st_as_sf(coords = c('longitude','latitude')) %>% 
  st_set_crs(4326)

#Read in the global coastline info
#Downloaded from https://www.naturalearthdata.com/
global_coastline <- read_sf("./data/MapShp/GlobalCoastline/ne_10m_coastline.shp", crs = 4326)

#Calculate the distances
distance <- array()
for(i in 1:length(region_list[,1])){
  result <- array()
  result <- st_distance(region_list_sf[i,], global_coastline[,3])
  distance[i] <- min(result)
  print(paste(i,": Done!"))
}

#Combine list with distance
region_info_with_distance <- cbind(region_list, distance)

#Write to file
write.csv(region_info_with_distance, 
          file = "./data/CI5-XId/region_info_with_distance.csv", quote = T, 
          row.names = F)