计算无边界相交的多边形交集：不同年份美国县界相交问题求解

预处理修正方案（适配已尝试函数的正确用法） #

st_make_valid()、st_snap()、st_set_precision()未生效多为参数或流程顺序错误，正确执行流程如下：

• 先将所有年份的shapefile统一转换为同一米制投影（推荐美国本土Albers等面积投影），避免地理坐标系下角度单位带来的精度计算偏差
• 对所有图层执行st_make_valid()修复原始数据的几何拓扑错误
• 调用st_set_precision(geom, 1000)设置几何精度，米制投影下该参数代表所有坐标对齐到1米精度网格，可自动抹平不同版本边界的微小偏移
• 精度设置完成后再次执行st_make_valid()修复精度调整引发的几何异常，再执行相交操作

带容差的直接相交方案 #

不需要单独做吸附操作，可直接在相交步骤指定吸附容差，一步完成对齐+相交：

• 所有图层完成投影统一、几何修复后，调用st_intersection(layer_year1, layer_year2, snapping_tolerance = 10)，其中容差参数单位与投影单位一致，米制投影下10代表10米以内的边界偏移会自动吸附，不会生成无效碎面

结果过滤优化方案 #

如果相交后仍有少量残留碎面，可通过面积规则直接过滤：

• 计算所有相交结果的面积，删除面积小于预设阈值的几何（比如阈值设为1000平方米，这类极小面基本都是边界偏移带来的无效结果）
• 若过滤后出现微小缺口，可调用st_buffer(geom, 1)做1米级的微小正缓冲补全，再通过st_simplify()平滑边界

小幅度负缓冲区的方案如果缓冲距离控制在5米以内，对县域级分析的面积统计误差可忽略，若对边界精度要求不高可继续使用。

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