如何在 R 中为 sf 对象创建几何列

编程问答 2022-05-22

问题描述

我的数据如下：

example <- tibble(geometry = "POINT (-3.8258 43.4109)")

我想将示例从 tibble 转换为 sf 对象 我有以下解决方案，但我想知道是否有更简单的方法（因为我也会有多边形和线条）

left.border <- "POINT ("
right.border <- ")"
example %>% 
  mutate(coord =str_match(geometry,paste(left.border,'(.+)',right.border,sep=''))[,2]) %>% 
  separate( coord,c("X","Y"),sep = " ") %>% 
  mutate(X = str_sub(X,2)) %>% 
  mutate(Y = str_sub(Y,1,nchar(Y)-2)) %>%
  mutate(X= as.numeric(X),Y= as.numeric(Y) ) %>%
  sf::st_as_sf(coords = c("X","Y"))

解决方法

如果您所有的几何图形都是众所周知的文本格式，您可以使用：

example = sf::st_as_sf(example,wkt = "geometry")
geospatialrrsf