R语言中基于另一栅格修改指定数量栅格单元值的技术咨询

Hey there! 我帮你把你的土地利用/覆盖变化模拟相关内容整理成清晰的Markdown格式啦，方便你后续调试、分享或者进一步扩展：

我目前正在用R语言开展土地利用/覆盖变化模拟工作，手里有两幅关键栅格数据：

• 土地利用/覆盖类型栅格
• 毁林风险信息栅格

我的核心需求是：借助毁林风险栅格，从土地利用类型中识别出那些更容易发生毁林的森林像素（森林是其中一类土地利用/覆盖类型）。

目前我已经写出了可运行的基础代码，以下是完整可复现的示例：

# 创建土地利用/覆盖栅格
real <- raster(nrows=25, ncols=25, vals=round(rnorm(625, 3), 0))
real[ real > 3 ] <- 3 # 仅为确保土地利用类型值不超过3（这里假设数值3代表森林类型）

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