The glcm
function outputs a stack
, with each band containing values for each of the statistics calculated. By default, all available statistics are computed.
You won't be able to put this stack
into a single band within the brick
as you are trying with band[i]<-
. This is probably what is giving you the error message.
One solution would be to output each stack
to a list
and then look at stacking them together once the loop has completed.
Something like this should work:
out <- list()
for(i in 1:dim(LS)[3]) {
out[i] <- glcm::glcm(LS[[i]], window = c(3, 3), na_opt = 'center')
}