此代码用于计算流苏帽,但我不知道为什么它不起作用。
错误是:
lsat8 [[c(1:6)]]中的错误:尝试在其中选择多个元素 vectorIndex
这是代码
library("RStoolbox")
library(raster)
lsat8 <- list.files(path = "E:/ndvi_1/", pattern = "*.tif", full.names=FALSE)
lsat8_tc <- tasseledCap(lsat8[[c(1:6)]], sat = "Landsat8OLI")
lsat8_tc
plot(lsat8_tc)
outpath <- "E:\ndvi_1\tc output test/"
outfiles <- paste0(outpath,"NDVI_",inputfiles)
for (f in 1:length(files)){
r <- raster(files[f])
rc <- crop(r, e)
rm <- mask(rc, ppoly)
rw <- writeRaster(rm,outfiles[f],overwrite=TRUE)