library(raster) library(zoo) # 读取tif文件 files <- list.files(path = "E:/LAI/GLASS_LAI/IWEMS_LAI/8118x", pattern = "*.tif", full.names = TRUE) rasters <- lapply(files, raster) rasters_list <- list(rasters) # 将栅格数据转换为时间序列 ts_list <- lapply(rasters_list, function(x) { ts <- as.zoo(x) index(ts) <- as.POSIXct(strptime(index(ts), format = "%Y-%m-%d")) ts }) # 生成一个每小时的时间序列 start_time <- as.POSIXct(strptime("2021-01-01 00:00:00", format = "%Y-%m-%d %H:%M:%S")) end_time <- as.POSIXct(strptime("2021-12-31 23:00:00", format = "%Y-%m-%d %H:%M:%S")) hourly_ts <- zoo(, seq(start_time, end_time, by = "hour")) # 将每8天的数据插值为每小时的数据 for (i in seq_along(ts_list)) { hourly_ts <- merge(hourly_ts, na.approx(ts_list[[i]]), all = TRUE) } # 将每小时的数据导出为tif文件 for (i in seq_along(hourly_ts)) { tif_name <- paste0("hourly_", format(index(hourly_ts)[i], "%Y%m%d%H%M"), ".tif") writeRaster(raster(hourly_ts[i]), filename = tif_name, format = "GTiff") } 请根据刚才的问题优化这个代码

好的，根据你之前的问题，我建议你在将rasters转换为rasters_list时去掉列表，直接将rasters转换为zoo对象。这样你就可以直接在栅格数据上进行时间序列插值和导出。以下是修改后的代码：

library(raster) 
library(zoo) 

# 读取tif文件
files <- list.files(path = "E:/LAI/GLASS_LAI/IWEMS_LAI/8118x", pattern = "*.tif", full.names = TRUE) 
rasters <- lapply(files, raster)

# 将栅格数据转换为时间序列
ts_list <- lapply(rasters, function(x) {
  ts <- as.zoo(x)
  index(ts) <- as.POSIXct(strptime(index(ts), format = "%Y-%m-%d"))
  ts
})

# 生成一个每小时的时间序列
start_time <- as.POSIXct(strptime("2021-01-01 00:00:00", format = "%Y-%m-%d %H:%M:%S")) 
end_time <- as.POSIXct(strptime("2021-12-31 23:00:00", format = "%Y-%m-%d %H:%M:%S")) 
hourly_ts <- zoo(, seq(start_time, end_time, by = "hour"))

# 将每8天的数据插值为每小时的数据
for (i in seq_along(ts_list)) {
  ts_hourly <- na.approx(ts_list[[i]], xout = index(hourly_ts))
  hourly_ts <- merge(hourly_ts, ts_hourly, all = TRUE)
}

# 将每小时的数据导出为tif文件
for (i in seq_along(hourly_ts)) {
  tif_name <- paste0("hourly_", format(index(hourly_ts)[i], "%Y%m%d%H%M"), ".tif")
  writeRaster(raster(hourly_ts[[i]]), filename = tif_name, format = "GTiff")
}

这个代码现在应该能够成功将每8天的栅格数据插值为每小时的数据，并将每小时的数据导出为tif文件。