Giá trị ô tối đa từ một raster xếp chồng

9

Làm thế nào tôi có thể tìm thấy một giá trị ô tối đa từ một raster xếp chồng.

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

hoạt động tốt nhưng

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

cho NA.

Tất nhiên không phải trong một raster xếp chồng lên nhau.

Đây là mã của tôi:

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)

raster r

— Không
nguồn

Bạn đang tìm kiếm tối đa của tất cả các lớp hoặc tối đa của mỗi lớp? Dù sao, bạn không sử dụng maxValueđúng cách. Theo trang trợ giúp, tốt hơn bạn nên sử dụng đối số bổ sung... Additional argument: layer number (for RasterStack or RasterBrick objects)

Tôi đang tìm kiếm tối đa của tất cả các lớp để có cùng một tỷ lệ như my.at <- seq (0, giá trị ô tối đa của tất cả các lớp, mức tăng). Cảm ơn, Nahm

— Nahm

Tôi đã nhận được nó với cellStats # geostat-cference.org/system/files/lewis_tutorAM.pdf Rad1998.max <- cellStats (RAD1998.all, 'max') Rad1998.all.max <- max (Rad1998.max) Rad1998.all .max

— Nahm

9

Ví dụ sau đây cho thấy hai cách để có được giá trị raster tối đa trong ngăn xếp. Việc đầu tiên sử dụng max()cũng cung cấp cho bạn một loạt các thông tin hữu ích khác. Phương thức thứ hai sử dụng maxValue(), chỉ cung cấp giá trị tối đa của cả hai trình quét trong ngăn xếp

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)

— Aaron
nguồn

0

max(MaxValue(raster)) sẽ cung cấp giá trị tối đa của tất cả các raster trong ngăn xếp.

— Umaomamaomao
nguồn