使用随机森林制作最终预测图

Question

对于随机森林模型，我有22个植被指数栅格作为我的预测变量。运行模型时，我最终确定只能使用五个植被指数（栅格）进行模型预测。因此，仅想使用这五个栅格图层，我想制作一个预测图，对于该预测图，我对每个图都拥有具有落叶百分比的现场数据以训练模型。在这里，我不知道如何将植被指数栅格值与野外数据合并在一起，以便我可以使用自己创建的模型制作最终的预测图。

我已经搜索了几个代码，并试图合并我的问题，但是有些行根本没有运行，而有些却让我感到困惑。

###########################################################################
library(rgdal)
library(raster)
setwd('C:\\Users\\rajeev.bhattarai\\Documents\\Research\\3_date_change_stack')
image='prediction_stack.tif'######it is a 5 layered raster file(predictor variables for RF model).
Stack=new("GDALReadOnlyDataset", image)#####when run shows up in data as a Formal class GDALReadonlyDataset
head(stack)
dim(stack) # it shows null
width<-dim(stack)[2]####I don't know what this code is doing, shows null(empty) when run
height<-dim(stack)[1]####I don't know what this code is doing, shows null(empty) when run
imagedata<-data.frame(getRasterTable(stack))#####when run, shows error saying object is not a member of class GDALReadonlyDataset

##############我现在被卡住了

我希望获得使用随机森林的预测图