Jak znaleźć maksymalną wartość komórki ze stosu rastrowego.
Rmax <- maxValue(RAD1998.all[[1]])
działa dobrze, ale
Rmax <- maxValue(RAD1998.all[[2]])
daje NA.
Oczywiście nie w ustawionym w stos rastrze.
Oto mój kod:
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
prawidłowego sposobu. Według strony pomocy lepiej użyć dodatkowego argumentu... Additional argument: layer number (for RasterStack or RasterBrick objects)
Odpowiedzi:
Poniższy przykład pokazuje dwa sposoby uzyskania maksymalnej wartości rastra na stosie. Pierwszy wykorzystuje,
max()
który daje również wiele innych przydatnych informacji. Druga metoda wykorzystujemaxValue()
, która daje tylko maksymalną wartość obu rastrów w stosieźródło
max(MaxValue(raster))
zapewni maksymalną wartość wszystkich rastrów w stosie.źródło