R sf包st_intersection计算Polygon与Multipolygon相交报错求助

这个replacement has 2 rows, data has 1报错和数据体量无关，核心触发逻辑是：

• st_intersection计算两个几何的相交结果时，你的poly1是多部件MULTIPOLYGON，和poly2相交后生成了2个独立的相交几何片段，返回结果有2行；
• 但输入的sf对象本身是1行的data.frame结构，sf包尝试把2行的几何列绑定回1行的属性表时，行数不匹配直接抛出错误。

你之前把poly2转成MULTIPOLYGON没解决问题，是因为没有触碰到核心诱因；新建简单多边形求交能正常运行，只是因为简单多边形相交后大概率只返回1个几何片段，行数刚好匹配没触发报错，属于巧合。

按顺序尝试，哪步跑通就停：

1. 统一输入对象类型，避免sf类和sfc类混算
   你贴的poly2构建代码最终生成的是sfc纯几何列对象，但你打印的poly2结构显示它是sf类型的要素集，大概率是中间转换时对象类不统一，先全部提取纯几何列计算排除属性绑定干扰：

   # 提取两个对象的纯几何列做计算
   inter_geom <- st_intersection(st_geometry(poly1), st_geometry(poly2))
   

   如果这步能正常返回结果，后续要带属性计算的话，给两个sf对象补完唯一ID再算就行。

2. 修复无效几何
   大型公开空间数据集大概率存在自相交、重叠、环方向错误等无效几何问题，是空间计算报错的高频诱因，先做有效性检查和修复：

   # 检查几何有效性
   st_is_valid(poly1)
   st_is_valid(poly2)
   
   # 一键修复无效几何
   poly1_fix <- st_make_valid(poly1)
   poly2_fix <- st_make_valid(poly2)
   
   # 修复后重新计算
   inter_res <- st_intersection(poly1_fix, poly2_fix)
   

   注意：你手动构建poly2的代码要确认b_xmin < b_xmax、b_ymin < b_ymax，如果坐标顺序写反，会生成自相交的无效矩形，也会触发奇怪报错。

3. 给输入要素加唯一行标识
   如果上面两步还没解决，给两个sf对象手动加唯一ID字段，让sf在求交时能正确匹配几何和属性的对应关系，避免行绑定错误：

   poly1$poly_id <- seq_len(nrow(poly1))
   poly2$clip_id <- seq_len(nrow(poly2))
   inter_res <- st_intersection(poly1, poly2)
   

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