누적 래스터의 최대 셀 값

누적 래스터에서 최대 셀 값을 어떻게 찾을 수 있습니까?

Rmax <- maxValue(RAD1998.all[[1]]) 

잘 작동하지만

Rmax <- maxValue(RAD1998.all[[2]]) 

NA를 제공합니다.

물론 쌓아 놓은 래스터에는 없습니다.

내 코드는 다음과 같습니다.

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)

다음 예제는 스택에서 최대 래스터 값을 얻는 두 가지 방법을 보여줍니다. 첫 번째는 max()다른 유용한 정보를 제공합니다. 두 번째 방법 maxValue()은를 사용하여 스택에있는 두 래스터의 최대 값 만 제공합니다.

library(raster)  

# Generate some georeferenced raster data
x = matrix(rnorm(400),20,20)
rast = raster(x)
extent(rast) = c(36,37,-3,-2)
projection(rast) = CRS("+proj=longlat +datum=WGS84")

y = matrix(rnorm(400),20,20)
rast2 = raster(y)
extent(rast2) = c(36,37,-3,-2)
projection(rast2) = CRS("+proj=longlat +datum=WGS84")

raster = stack(rast, rast2)

# Now run the statistics
max(raster) # Provides min, max and additional details  # Example 1

maxValue(raster)  # Gives both values                   # Example 2...
maxValue(raster)[[1]] # Gives first in stack max value
maxValue(raster)[[2]] # Gives second in stack max value

> maxValue(raster)  # Gives both values
[1] 2.688376 2.971443
> maxValue(raster)[[1]] # Gives first in stack max value
[1] 2.688376
> maxValue(raster)[[2]] # Gives second in stack max value
[1] 2.971443
> 
> max(raster) # Provides min, max and additional details
class       : RasterLayer 
dimensions  : 20, 20, 400  (nrow, ncol, ncell)
resolution  : 0.05, 0.05  (x, y)
extent      : 36, 37, -3, -2  (xmin, xmax, ymin, ymax)
coord. ref. : +proj=longlat +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0 
data source : in memory
names       : layer 
values      : -1.457908, 2.971443  (min, max)

max(MaxValue(raster)) 스택에있는 모든 래스터의 최대 값을 제공합니다.