我试图在一系列连续时间点的太阳辐射强度图下找到该区域。 基本上我想要在7天内每小时过去24小时太阳辐射的积分 - 过去24小时的移动总和 - （我怀疑土壤温度是过去24小时太阳辐射的结果）< / p>
这是我正在使用的代码，它可以工作，但我想自动化它，这样我就可以轻松地更改集成的时间窗口（尝试12,18,24,36小时）并获得每小时打印/保存的表格综合的太阳辐射值，我可以根据我的每小时温度数据绘制，以确定是否存在关系）
这里：Rg-solar辐射在10分钟内测量 num - 数据帧中的条目号 AUC_xxx - 过去24小时内的太阳总辐射量
y<-as.numeric(xx$Rg[xx$num["2015-09-13 14:10"]:xx$num["2015-09-14 14:00"]])
x<-c(1:length(y))
id <- order(x)
AUC_s14_14 <- sum(diff(x[id])*rollmean(y[id],2))
我尝试了rollapply，但我再次陷入困境：
rollapply(xx$Rg[xx$num["2015-09-13 00:10"]:xx$num["2015-09-14 00:00"]], width = 144, by = 6, FUN = **"INTEGRAL"**, na.rm = TRUE, align = "left")
感谢您的帮助！