Estoy trabajando con múltiples imágenes MODIS(una sola banda, extraída previamente) en R. Me gustaría encontrar la manera de, en este caso, extraer el píxel x = 1 e y = 1 de todas ellas e introducirlos en un raster stack para después aplicarles otras operaciones. He intentado diferentes cosas pero no doy con la tecla. Por ejemplo, he estado probando cosas de este estilo.
library(raster)
library(rgdal)
setwd("/Users/guillermorodriguezlopez/Desktop/ndsi_modis") #Directorio
# Carga de imagenes
images <- list.files(path=getwd(),
pattern="*.tif$|*.TIF$")
imagestack <- stack(images) #Stack con los archivos raster.
for (i in seq_along(imagestack)) {
r = stack(crop(imagestack, extent(imagestack, 1,1,1,1)))
}
Si a alguien tiene alguna idea, será bienvenida! Gracias de antemano y un cordial saludo!!