如何在 R 中找到最匹配多边形整体土地覆盖比例的区域？

Question

在 R 中，我有一个多边形形状文件、一个来自 NLCD 的土地覆盖数据栅格，以及一个总结每个多边形内土地覆盖类型比例的数据框。现在我想在每个多边形内找到一个特定大小的区域（例如 115 km^2），该区域与该多边形的土地覆盖类型的总体比例最匹配。

我认为包

optmatch

可能很适合这个问题，使用

fullmatch

功能，也许

grouped_by

多边形，但我不知道如何去做。

我的多边形数据是一个带有标识符列和几何列的 sf 对象。我的栅格土地覆盖数据以 TIFF 格式加载。我的数据框类似于下面创建的数据框，但具有更多的类和 ID。

id <- c(1:9)
class <- c(11,21,35, 11,21,35, 11,21,35)
proportion <- c(20, 35, 45, 20, 20, 60, 5, 90, 5)
df <- data.frame(cbind(class, id, proportion))

我将如何使用

optmatch

来实现这一点，或者我可以做更好的选择吗？