如何筛选sf对象中的行：采样点及已知点位筛选方法问询

问题1：通过几何匹配筛选sf对象中的单个点 #

你的代码无效是因为st_equals()返回的是列表格式的索引匹配结果，直接用于子集化sf对象会出错。以下是两种可行的修正方案：

方案1：利用st_equals结合逻辑向量转换 #

library(sf)
data(meuse)
meuse_sf <- st_as_sf(meuse, coords = c("x", "y"), crs = 28992)

# 提取单个点的几何并转换为sf对象
temp <- meuse_sf[1, ]
# 用sapply将st_equals的列表结果转为逻辑向量
meuse_sf[sapply(st_equals(meuse_sf, temp), function(x) length(x) > 0), ]

方案2：更简洁的st_filter方法 #

st_filter可以直接按空间关系筛选要素，代码更直观：

temp <- meuse_sf[1, ]
st_filter(meuse_sf, temp, .predicate = st_equals)

问题2：通过已知坐标匹配sf要素 #

首先要把单独的st_point转换为带对应CRS的sf对象（必须和meuse_sf的CRS一致），再用空间匹配筛选：

library(sf)
data(meuse)
meuse_sf <- st_as_sf(meuse, coords = c("x", "y"), crs = 28992)

# 定义带CRS的点sf对象
location_point <- st_sfc(st_point(c(181027, 333363)), crs = 28992)
location_sf <- st_as_sf(location_point)

# 筛选匹配的要素
st_filter(meuse_sf, location_sf, .predicate = st_equals)

sf与ArcGIS对比学习资源推荐 #

• 专栏类文章：聚焦sf包实现ArcGIS核心工具的对应方法，对比两种工具在缓冲区分析、叠加运算、空间连接等操作上的流程差异与代码实现
• 分步教程：按空间数据处理全流程（导入、清洗、分析、输出），逐一对比sf与ArcGIS的操作逻辑和实现步骤
• 案例解析：针对点模式分析、邻域分析、空间插值等具体任务，对比两种工具的算法实现、参数设置与结果输出

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