基于DBF表匹配栅格值并结合嵌套流域实现条件过滤的技术咨询

我手头有一张覆盖澳大利亚的大型栅格数据，已经用多面要素的sf对象（嵌套流域）对它做了crop()和mask()处理。我用foreign包读取了DBF表，用terra::rast()读取了TIFF栅格，并且已经确保两者的CRS一致。

但当我尝试把DBF表和栅格的像元匹配时，结果和预期不符——因为我的栅格本身没有数值，只包含“森林分类（原生林和非森林）”的信息。我需要把DBF表的数据关联到栅格上，这样才能对像元进行筛选。

举个具体的需求例子：我想筛选出满足以下条件的像元：

• 来自Forest_CATEGO列的“原生林”
• HT_CODE和COVER_CODE列的值属于c(1,2,3)

之后，我还想进一步筛选那些**满足上述条件的像元占比超过60%**的流域（这里要注意流域是嵌套结构的）。我不确定这种筛选是否可行？

备注：内容来源于stack exchange，提问作者yataya