R：将xyz文本复制到多个栅格

Question

我有一个逗号分隔的文本文件，标题名称为

lat,long,z1,z2,z3,z4,z5

这些列包含lat，long和multiple height值（z1，z2，z3，z4）。纬度和长度数据以5米的相等间隔隔开。

在投影信息为UTM，区域18的情况下，将数据转换为具有z值的单个栅格或R栅格堆栈的最有效方法是什么？

可以访问样本数据here。

Answer 1

你可以这样做：

x <- data.frame(lon=1:11, lat=5:15, z1=runif(11), z2=11:1)
library(raster)
b <- rasterFromXYZ(x)

#> b
#class       : RasterBrick 
#dimensions  : 11, 11, 121, 2  (nrow, ncol, ncell, nlayers)
#resolution  : 1, 1  (x, y)
#extent      : 0.5, 11.5, 4.5, 15.5  (xmin, xmax, ymin, ymax)
#coord. ref. : NA 
#data source : in memory
#names       :         z1,         z2 
#min values  : 0.08495835, 1.00000000