在R中迭代合并人口未达阈值的普查区至邻近单元

• 核心需求：以迭代方式将人口低于设定阈值（如≥200）的普查区与人口最少的邻近单元合并，合并时对单元的属性值（如总人口、西班牙裔人口）进行求和；需处理无邻近单元的"孤立单元"（岛屿型），最终用于空间自相关分析。
• 工具偏好：优先使用R语言编写可重复运行的代码；也熟悉QGIS、Geoda，可手动调整邻接矩阵，但更希望通过代码自动查找最近邻。
• 示例数据：纽约州里士满县（斯塔滕岛）2020年普查区数据，其中7个普查区人口低于200阈值，数据获取代码如下：

library(tidycensus)
richmondtracts <- get_decennial(
  geography = "tract",
  state = "NY",
  county = "Richmond",
  variables = c(
    Hispanic = "P2_002N" #Hispanic population per tract
  ),
  summary_var = "P2_001N", #total population per tract
  year = 2020,
  geometry = TRUE
)

• 当前困境：作为编程新手，尝试ChatGPT生成的while循环方案无法正常运行；曾找到一个基于面积合并小型多边形至最大邻近单元的论坛方案，但与当前按人口阈值合并至人口最少邻近单元的需求不匹配。

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