从R中的栅格文件中获取数据

Alex Chumbley 发表于 Dev

亚历克斯·丘比

我是栅格文件的新手，但是它们似乎是打开包含所有天气数据的大型政府文件的最佳方法，因此，我试图弄清楚如何使用它们。作为参考，我正在下载位于此处的文件（只是一些有关工厂天气的信息）。当我使用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以相同的方式访问此信息？

用户1738753

# 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)

本文收集自互联网，转载请注明来源。

如有侵权，请联系[email protected] 删除。