gpt4 book ai didi

r - Terra 等效于 raster::stack()?

转载 作者:行者123 更新时间:2023-12-05 01:05:00 27 4
gpt4 key购买 nike

基本上是标题。我知道您可以使用 rast() 读取栅格文件夹,但我只想堆叠两个单独读取的栅格。谢谢

最佳答案

请注意,使用 raster::stack,您可以在多个参数 (stack(x1,x2,x3)) 或列表 ( stack(list(x1,x2, x3))).

terra 的 c 不再适用。你需要区分:

  • 在提供单独的参数时使用 c(x1, x2, x3)
  • 在以列表形式提供参数时使用 rast(list(x1,x2,x3))
library(terra)
#> terra 1.5.21

x <- rast(xmin=-110, xmax=-80, ymin=40, ymax=70, ncols=30, nrows=30)
values(x) <- 1:ncell(x)

many_rasters <- list(x,x)

## this works
rast( many_rasters)
#> class : SpatRaster
#> dimensions : 30, 30, 2 (nrow, ncol, nlyr)
#> resolution : 1, 1 (x, y)
#> extent : -110, -80, 40, 70 (xmin, xmax, ymin, ymax)
#> coord. ref. : lon/lat WGS 84
#> sources : memory
#> memory
#> names : lyr.1, lyr.1
#> min values : 1, 1
#> max values : 900, 900

## just using c creates a list
c(many_rasters)
#> [[1]]
#> class : SpatRaster
#> dimensions : 30, 30, 1 (nrow, ncol, nlyr)
#> resolution : 1, 1 (x, y)
#> extent : -110, -80, 40, 70 (xmin, xmax, ymin, ymax)
#> coord. ref. : lon/lat WGS 84
#> source : memory
#> name : lyr.1
#> min value : 1
#> max value : 900
#>
#> [[2]]
#> class : SpatRaster
#> dimensions : 30, 30, 1 (nrow, ncol, nlyr)
#> resolution : 1, 1 (x, y)
#> extent : -110, -80, 40, 70 (xmin, xmax, ymin, ymax)
#> coord. ref. : lon/lat WGS 84
#> source : memory
#> name : lyr.1
#> min value : 1
#> max value : 900

关于r - Terra 等效于 raster::stack()?,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/71213802/

27 4 0
Copyright 2021 - 2024 cfsdn All Rights Reserved 蜀ICP备2022000587号
广告合作:1813099741@qq.com 6ren.com