我试图了解我的研究区域在过去32年中NDVI值如何变化。更具体地说，我需要评估每个像素的NDVI值在32年中如何变化。为此，我需要32个Landsat图像，并为其计算了NDVI，从而获得了32个NDVI栅格。我的问题出现在尝试对R中的栅格堆栈的每个像素应用Mann-Kendall测试以生成2个新栅格，一个栅格传递每个像素的tau值，另一个传递每个像素的p值。
我试图将栅格堆栈转换为数据帧并运行mann-kendall，但是没有找到写入新栅格的方法。我也尝试利用lapply函数，但我认为这不适用于栅格。
对此表示感谢。