r - 从 R 中的光栅文件中获取数据

Question

我是光栅文件的新手，但它们似乎是打开包含所有天气数据的大型政府文件的最佳方式，所以我想弄清楚如何使用它们。作为引用，我正在下载位于 here 的文件(只是一些磨坊天气的东西)。当我使用R的光栅包导入文件时是这样的

> r <- raster("/path/to/file.grb")

一切正常。当我输入

时，我什至可以获得一些元数据

> r
class       : RasterLayer 
band        : 1  (of  37  bands)
dimensions  : 224, 464, 103936  (nrow, ncol, ncell)
resolution  : 0.125, 0.125  (x, y)
extent      : -125.0005, -67.0005, 25.0005, 53.0005  (xmin, xmax, ymin, ymax)
coord. ref. : +proj=longlat +a=6371200 +b=6371200 +no_defs 
data source : /path/to/file.grb 
names       : NLDAS_MOS0125_H.A20140629.0100.002 

此时我所做的就是以一种非常明显的方式为栅格编制索引。

> r[100,100]
267.1

所以，我想我可以“索引”它，但我不知道数字 267.1 是什么意思。这肯定不是牢房里的全部。应该有一堆变量，包括但不限于土壤水分、地表径流和蒸发。

如何使用 R 以相同的方式访问此信息？

Answer 1

# create two rasters
r1 <- raster(matrix(ncol = 10, nrow = 10, runif(100))) 
r2 <- raster(matrix(ncol = 10, nrow = 10, runif(100)))

# creates a raster stack -- the stack (or brick function) allows you to
# to use multilayer band rasters
# http://www.inside-r.org/packages/cran/raster/docs/stack
st_r <- stack(r1, r2)

# extract values -- will create a matrix with 100 rows and two columns
vl <- getValues(st_r)