Início / coding / Perguntas / 79542480
Accepted
nukubiho
Asked: 2025-03-29 05:56:48 +0800 CST

Como calcular NDVI a partir de um raster quadridimensional?

  • 772

Tenho várias cenas do satélite Sentinel-2 para o mesmo local de datas diferentes. Cada cena contém bandas espectrais básicas (incluindo vermelho "B04" e infravermelho próximo "B08"). Carreguei os dados como um SpatRasterDatasetusando o terrapacote. Como posso calcular o NDVI definido como NDVI = (NIR - RED) / (NIR + RED) para cada data?

library(terra)
r <- rast(system.file("ex/logo.tif", package="terra")) / c(1,2,3)
r <- c(r, r[[1]] + r[[2]])
names(r) <- paste0("B0", c(2:4, 8)) 
dates <- c("2023-02-08", "2023-03-02", "2023-04-06")
s <- sds(r, r, r)
names(s) <- dates

1 respostas

  1. Best Answer

    Dados de exemplo

    library(terra)
r <- rast(system.file("ex/logo.tif", package="terra")) / c(1,2,3)
r <- c(r, r[[1]] + r[[2]])
names(r) <- paste0("B0", c(2:4, 8)) 
dates <- c("2023-02-08", "2023-03-02", "2023-04-06")
s <- sds(r, r, r)
names(s) <- dates

    Solução (aqui assumindo que "BO2" é "NIR" e "B04" é "RED")

    ndvi <- function(NIR, RED) (NIR - RED) / (NIR + RED)
v <- lapply(s, \(x) lapp(x[[c("B02", "B04")]], ndvi))
v <- rast(v)
panel(v, nr=1)

    insira a descrição da imagem aqui

    Você pode usar uma lista em vez de um SpatRasterDataset

    ss <- list(r, r, r)
names(ss) <- dates
vv <- lapply(ss, \(x) lapp(x[[c("B02", "B04")]], ndvi)) |> rast()
    • 1

relate perguntas