2015-06-19 51 views
0

我知道有一个类似的问题,但它不在rasterbrick中。我想找到这个光栅砖的最大值,我试过which.max,但我不明白输入。我想知道什么是最大值,并且还得到每个图层的所有最大值的列表。如何在rasterbrick中找到最大值?

> bshum2 
class  : RasterBrick 
dimensions : 5, 4, 20, 2920 (nrow, ncol, ncell, nlayers) 
resolution : 1, 1 (x, y) 
extent  : 14, 18, 0, 5 (xmin, xmax, ymin, ymax) 
coord. ref. : NA 
data source : in memory 
names  : X1979.01.01.00.00.00, X1979.01.01.03.00.00, 
min values :   0.011909009,   0.011799878,   
max values :   0.01536649,   0.01594792,   time  : 1979-01-01 00:00:00, 1979-12-31 21:00:00 (min, max) 

> which.max(bshum2) 
class  : RasterLayer 
dimensions : 5, 4, 20 (nrow, ncol, ncell) 
resolution : 1, 1 (x, y) 
extent  : 14, 18, 0, 5 (xmin, xmax, ymin, ymax) 
coord. ref. : NA 
data source : in memory 
names  : layer 
values  : 172, 1132 (min, max) 
+1

你阅读'which.min'应用于栅格*对象('which.min'?)帮助页面? –

+0

你想要的最大值,或只是位置? –

+0

@CarlWitthoft我想要的价值,而不是位置。是的,我读了帮助页面,但我不明白。 –

回答

2

这是你可以得到这些:

library(raster) 
b <- brick(system.file("external/rlogo.grd", package="raster")) 
minValue(b) 
# [1] 0 0 0 
maxValue(b) 
# [1] 255 255 255 
+0

这只是部分适用于我的情况。在多层rasterbrick中的maxValue返回每个图层的最大值,所以我做的是 a = maxValue(bshum2) 摘要(a) –

+1

您要求'所有图层的最大值'。但下一步确实微不足道:'min(minValue(b))' – RobertH

相关问题