R 空间约束的多变量聚类

我有一个以多边形表示的区域空间数据集。每个多边形都有一个人口属性，表示居住在那里的人数

我必须将这些地区划分为更大的空间区域，同时考虑到由此产生的区域的人口不应超过阈值

看起来ArcGIS有一个解决此问题的算法：

因为我没有访问ArcGIS的权限，所以我正在寻找使用R Stats的解决方案。是否有人知道具有此功能的软件包