R中光栅砖中的单独层
Separate layers in raster brick in R
R 新手:
我有一个包含 5 个栅格层的栅格砖。我如何打破这块砖,以便每一层都可以分配给它自己的对象,我可以单独操作每一层?
这是砖的信息,以备不时之需。
>>eco.brick
class : RasterBrick
dimensions : 1920, 1440, 2764800, 5 (nrow, ncol, ncell, nlayers)
resolution : 0.008333333, 0.008333333 (x, y)
extent : -71.00014, -59.00014, -29.99181, -13.99181 (xmin, xmax, ymin, ymax)
crs : +proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +towgs84=0,0,0
source : /tmp/RtmpdeDrSf/raster/r_tmp_2021-05-05_043302_2341_32265.grd
names : Central.Andes, Amazonian.Orinocan.Lowland, Eastern.Highlands, Gran.Chaco, Monte.Patagonian
min values : 0, 0, 0, 0, 0
max values : 1, 1, 1, 1, 1
示例数据
library(raster)
b <- brick(system.file("external/rlogo.grd", package="raster"))
你可以像这样列出图层
x <- as.list(b)
您还可以将每一层分配给一个变量
r1 <- b[[1]]
r2 <- b[[2]]
R 新手:
我有一个包含 5 个栅格层的栅格砖。我如何打破这块砖,以便每一层都可以分配给它自己的对象,我可以单独操作每一层?
这是砖的信息,以备不时之需。
>>eco.brick
class : RasterBrick
dimensions : 1920, 1440, 2764800, 5 (nrow, ncol, ncell, nlayers)
resolution : 0.008333333, 0.008333333 (x, y)
extent : -71.00014, -59.00014, -29.99181, -13.99181 (xmin, xmax, ymin, ymax)
crs : +proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +towgs84=0,0,0
source : /tmp/RtmpdeDrSf/raster/r_tmp_2021-05-05_043302_2341_32265.grd
names : Central.Andes, Amazonian.Orinocan.Lowland, Eastern.Highlands, Gran.Chaco, Monte.Patagonian
min values : 0, 0, 0, 0, 0
max values : 1, 1, 1, 1, 1
示例数据
library(raster)
b <- brick(system.file("external/rlogo.grd", package="raster"))
你可以像这样列出图层
x <- as.list(b)
您还可以将每一层分配给一个变量
r1 <- b[[1]]
r2 <- b[[2]]