我正在尝试使用R中的lulcc package来进行一些土地利用模拟。据我了解,我需要在土地使用栅格/堆栈中有一个“时间”槽,以符合包装的数据格式。例如,根据包装中提供的示例,我们有:
require(lulcc)
## load data
data(pie)
## observed maps
obs <- ObsLulcRasterStack(x=pie,
pattern="lu",
categories=c(1,2,3),
labels=c("Forest","Built","Other"),
t=c(0,6,14))
该示例运行良好,但是,在加载自己的栅格时,我在ObsLulcRasterStack()
代码中使用“ t”遇到麻烦。首先,通过查看数据,我无法找到pie
数据中“ t”的来源。其次,我无法为ObsLulcRasterStack()
使用栅格在栅格上添加槽“ t”(希望能够解决该问题)。有任何想法吗?具体来说,这是我要修复的错误消息:
Error in validityMethod(object) : timesteps do not correspond with maps
我不知道如何将这些“时间步长”添加到栅格/堆栈中。根据该软件包,t
是“具有与每个观察到的地图相对应的时间步长的数值向量”。我尝试使用rts
包中的rts
函数,但这创建了ObsLulcRasterStack()
无法处理的新型对象:
Error in (function (classes, fdef, mtable) :
unable to find an inherited method for function ‘ObsLulcRasterStack’ for signature ‘"RasterStackTS", "character"’