Voronoi图与同半径圆的重叠面积计算问题咨询

我已经基于Excel文件中的点集生成了Voronoi图，现在想要以另一组坐标为圆心生成同半径的圆，计算每个圆与对应Voronoi单元的重叠面积。

遇到的具体问题 #

• 我通过函数生成了SpatialPolygonsDataFrame格式的Voronoi划分，但生成第二组坐标的圆时遇到了障碍：
  • 使用st_buffer并将Voronoi结果转为sf格式后没有效果；
  • 调用st_intersection时出现st_crs(x) == st_crs(y)的错误，推测是Voronoi集未定义坐标参考系（CRS）。
• 我不确定适用的最佳函数，也想了解sf包是否有工具可以直接输出各Voronoi单元与圆的重叠面积。
• 之前参考过一些相关方法，但无法和现有的Voronoi函数兼容。

测试数据集 #

Hubs <- cbind('X Coord' = c(52.37639999,52.36989975,52.86299896,52.01011658,51.67409897,50.84980011,51.88669968,52.1048943,52.0746994), 
              'Y Coord' = c(4.894589901,4.876679897,6.493730068,4.703330517,4.548630238,5.687580109,4.491449833,5.0528496,4.310130119), 
              'Info 1' = c(13,15,62,24,9,46,73,97,69))
References <- cbind('X Coord' = c(51.88529968,52.3360136,52.37440109,51.92269897,51.9192276,51.43019867,52.0780896,51.90299988,52.04100037,50.98810196), 
                    'Y Coord' = c(4.493549824,4.8750552,4.895800114,4.47382021,4.481627464,5.451980114,5.1470653,4.458670139,4.318089962,5.78059721))

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