R中的快速地理空间采样

我有一大组多边形（大约20k），我想从中采样点。我在R中使用了

sf

包中的

st_sample

函数，但速度非常慢。从所有多边形中采样大约需要5分钟，我需要重复此任务很多次（N>=1000），因此不实用

---

有没有办法加快采样速度？

还没有检查，但可能更快？我也有同样的问题（耗时>5小时）。我认为加快速度的一个选择是从一个简单的多边形采样，然后使用接受/拒绝步骤，使用st_交点获得所需多边形内的精确点。（我还没有测试过这个）。