I have some GeoTiff rasters with over 100 bands in total, but they are not huge: resolution is about 1200x600 and total size is less than 7 mb. I'm trying to do a simple NDVI calculation:

raster <- stack('path_to_raster_1', 'path_to_raster_2')
red <- 60 # red band number
NIR <- 75 # NIR band number
f_NDVI <- function(x){
NDVI <- calc(raster, f_NDVI)

The calculation runs extremely slow and takes about 10 minuets.

Is there a way to optimise it and why is the calculation so slow in my case?


1 Answer 1


I expect it will be much quicker this way (using a vectorized approach in the NDVI function)

f_NDVI <- function(x) {

# subset the two layers you need
s <- raster[[c(60, 75)]]

NDVI <- calc(s, f_NDVI)

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service and acknowledge that you have read and understand our privacy policy and code of conduct.

Not the answer you're looking for? Browse other questions tagged or ask your own question.