Vairāk

Izmantojot cilpu NDVI aprēķināšanai R?

Izmantojot cilpu NDVI aprēķināšanai R?


Es šobrīd izmantoju MOD09Q1 un vēlos aprēķināt NDVI visa gada datus. Vienu gadu man ir 46 1. joslas attēli un 46 2. joslas attēli.

  1. list.files

Es izveidoju 1. grupas un 2. grupas sarakstu.

Band1list <- list.files (path = ".", Pattern = "*.b01.tif", all.files = FALSE, full.names = FALSE, rekursīvs = FALSE, ignore.case = FALSE) Band2list <- list. faili (ceļš = ".", modelis = "*.b02.tif", all.files = FALSE, full.names = FALSE, rekursīvs = FALSE, ignore.case = FALSE)
  1. rastrs :: kaudze

Visu 46 1. joslas (2. joslas) attēlu sakraušana StackBand1 (StackBand2)

StackBand1 <- kaudze (Band1list) StackBand2 <- kaudze (Band2list)
  1. NDVI aprēķins, izmantojotpriekšcilpa

Es mēģinu sakraut viena datuma attēlus (piemēram, "2012.01.01.Smb01.tif" un "2012.01.01.Smb02.tif") un aprēķināt NDVI

priekš (i in 1: garums (Band1list)) {StackBand1 [[i]] == i StackBand2 [[i]] == i Stack <- ķieģelis (StackBand1 [[i]], StackBand2 [[i]])}

Šeit es nevarēju izpildīt cilpu visiem 46 attēliem. To var palaist visu, un rezultāts tiek parādīts tikai pēdējam pārim.


Nav nepieciešams apriekš(vai jebkura cita veida) cilpa, jo šāda informācija ir pilnībā iekļauta sadaļā rastrs iepakojums. Ja vēlaties aprēķināt NDVI, tad būtībā viss, kas jums nepieciešams, ir

## nepieciešamā pakotne (rastrs) y+x)})

Skatīties video: How to download multiple MODIS NDVI images for time series