R:文本 xyz 到多个光栅
R: Text xyz to multiple raster
我有一个逗号分隔的文本文件,其中 header 个名称
lat,long,z1,z2,z3,z4,z5
这些列包含纬度、经度和多个高度值(z1、z2、z3、z4)。纬度和经度数据以 5 米的相等间隔间隔开。
给定投影信息为 UTM,区域 18,将数据转换为单独的栅格或每个具有 z 值的 R 栅格堆栈的最有效方法是什么。
可以访问样本数据here。
你可以这样做:
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
我有一个逗号分隔的文本文件,其中 header 个名称
lat,long,z1,z2,z3,z4,z5
这些列包含纬度、经度和多个高度值(z1、z2、z3、z4)。纬度和经度数据以 5 米的相等间隔间隔开。
给定投影信息为 UTM,区域 18,将数据转换为单独的栅格或每个具有 z 值的 R 栅格堆栈的最有效方法是什么。
可以访问样本数据here。
你可以这样做:
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