Displaying 1 result from an estimated 1 matches for "navalue".
Did you mean:
avalue
2010 Oct 06
1
R getting slower until it breaks...
...#########################
## Fonction 1 - NDVI a partir de 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(ce...