R中将EPSG25832字符格式坐标转为空间点时遇解析错误求助

错误原因 #

你的代码存在两处核心逻辑错误：

1. 用st_as_sfc将字符串转为空间对象（sfc类型）后，错误地把该空间对象当成字符串进行坐标提取，导致后续处理的内容并非原始坐标值。
2. 最后用str2lang解析"x,y"格式字符串时，R会将其视为非法表达式（单独的数字, 数字不符合R语法规则），因此触发解析报错。

方案1：修正原有代码逻辑 #

直接从原始POINT字符串中拆分坐标，跳过不必要的空间对象转换步骤：

# 提取括号内的坐标文本
Total_df$coord_str <- stringr::str_extract(Total_df$Geometri_EPSG_25832, "(?<=\\().*(?=\\))")
# 拆分x、y坐标为数值列
coord_matrix <- stringr::str_split_fixed(Total_df$coord_str, "\\s+", 2)
Total_df$x_coord <- as.numeric(coord_matrix[,1])
Total_df$y_coord <- as.numeric(coord_matrix[,2])
# 生成空间多点对象
Total_df$geos_points <- lapply(1:nrow(Total_df), function(i) {
  sf::st_multipoint(matrix(c(Total_df$x_coord[i], Total_df$y_coord[i]), ncol = 2))
})

方案2：用sf原生函数处理（更简洁高效） #

sf包支持直接解析WKT格式字符串，无需手动替换字符：

# 直接将POINT字符串转为空间点对象（带坐标系）
Total_df$geometry <- sf::st_as_sfc(Total_df$Geometri_EPSG_25832, crs = 25832)
# 提取单独的x、y坐标列
coords <- sf::st_coordinates(Total_df$geometry)
Total_df$x_coord <- coords[,1]
Total_df$y_coord <- coords[,2]
# 若需转为多点对象（单条记录转多点意义有限，批量场景可调整）
Total_df$geos_points <- sf::st_cast(Total_df$geometry, "MULTIPOINT")

内容的提问来源于stack exchange，提问作者user23675742