Displaying 1 result from an estimated 1 matches for "cells1".
Did you mean:
cells
2010 Oct 06
1
R getting slower until it breaks...
...coordonnee Pixel et largeur #####
calculate_NDVI<- function(Type, object, VALUE) {
redorgreen <- ifelse(Type=="red",2,3)
list1 <- unstack(object)
rast1 <- list1[[1]]
rast2 <- list1[[redorgreen]]
NAvalue(rast1)<- -99999
NAvalue(rast2)<- -99999
cells1 <- getValuesBlock(rast1,row=VALUE[[2]],nrow=VALUE[[3]],col=VALUE[[1]],ncol=VALUE[[3]])
cells2 <- getValuesBlock(rast2,row=VALUE[[2]],nrow=VALUE[[3]],col=VALUE[[1]],ncol=VALUE[[3]])
cells1[is.na(cells1)]<-0;
cells2[is.na(cells2)]<-0;
calculNDVI <-(cells1 - cells2) / (cells...