Как я могу найти максимальное значение ячейки из сложенного растра.
Rmax <- maxValue(RAD1998.all[[1]])
работает нормально но
Rmax <- maxValue(RAD1998.all[[2]])
дает Н.А.
Конечно, не в сложенном растре.
Вот мой код:
RAD1998 <- raster(paste(getwd(), "/1998bil/1998ASC5min_ppt_spas1214_0001_19980202_0810_UTC.asc.bil", sep = ""))
list.ras <- mixedsort(list.files(paste(getwd(), "/1998bil/", sep = ""), full.names = T, pattern = ".asc.bil"))
RAD1998.all <- stack(list.ras)
maxValue
правильный путь. Согласно странице справки, вам лучше использовать дополнительный аргумент... Additional argument: layer number (for RasterStack or RasterBrick objects)
Ответы:
В следующем примере показаны два способа получить максимальное значение растра в стеке. Первый использует,
max()
который также дает вам множество другой полезной информации. Второй метод используетmaxValue()
, который дает только максимальное значение обоих растров в стекеисточник
max(MaxValue(raster))
обеспечит максимальное значение всех растров в стеке.источник