gpt4 book ai didi

r - 多次堆叠现有的 RasterStack

转载 作者:行者123 更新时间:2023-12-03 01:30:41 25 4
gpt4 key购买 nike

我的问题与一个不幸的未解答问题非常相似 recently posted on Stackoverflow 。我正在处理一个由十二层组成的 RasterStack 对象(一年中的每个月一层),我想将这些层复制十次,最终得到一个 RasterStack > 由 120 层组成,每 12 层相似(即,第 1 层与第 13 层相同,与第 25 层相同,依此类推)。

出于复制目的,我们以 raster 包中的示例为例:

library(raster)

file <- system.file("external/test.grd", package = "raster")
s <- stack(file, file, file, file, file, file, file, file, file, file, file, file)

stack(s, s, s, s, s, s, s, s, s, s)

class : RasterStack
dimensions : 115, 80, 9200, 120 (nrow, ncol, ncell, nlayers)
resolution : 40, 40 (x, y)
extent : 178400, 181600, 329400, 334000 (xmin, xmax, ymin, ymax)
coord. ref. : +init=epsg:28992 +towgs84=565.237,50.0087,465.658,-0.406857,0.350733,-1.87035,4.0812 +proj=sterea +lat_0=52.15616055555555 +lon_0=5.38763888888889 +k=0.9999079 +x_0=155000 +y_0=463000 +ellps=bessel +units=m +no_defs
names : test.1.1, test.2.1, test.3.1, test.4.1, test.5.1, test.6.1, test.7.1, test.8.1, test.9.1, test.10.1, test.11.1, test.12.1, test.1.2, test.2.2, test.3.2, ...
min values : 128.434, 128.434, 128.434, 128.434, 128.434, 128.434, 128.434, 128.434, 128.434, 128.434, 128.434, 128.434, 128.434, 128.434, 128.434, ...
max values : 1805.78, 1805.78, 1805.78, 1805.78, 1805.78, 1805.78, 1805.78, 1805.78, 1805.78, 1805.78, 1805.78, 1805.78, 1805.78, 1805.78, 1805.78, ...

当然,也可以像最后一行代码那样手动完成,但这对我来说似乎很不方便。任何有关如何以更好的方式实现我的目标的建议将不胜感激!

最佳答案

您可能喜欢此处的 mget,因为它采用对象名称的字符 向量并返回对象,因此您可以执行以下操作:

big.stack <- stack( mget( rep( "s" , 12 ) ) )

nlayers( big.stack )
#[1] 144

如果您不喜欢使用 mget(),或者使用 replicate() 将它们放入 list 中,然后堆叠list 作为 rasterLayersliststack() 的有效输入...

ll <- replicate( 12 , s )

big.stack2 <- stack( ll )

identical( big.stack , big.stack2 )
#[1] TRUE

关于r - 多次堆叠现有的 RasterStack,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/21041499/

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