Question

我的问题非常类似于一个遗憾的未回答的问题recently posted on Stackoverflow。我正在处理一个由{12}层组成的RasterStack对象（一年中每个月一个），我想复制这些层十次，最后得到一个由{120}层组成的RasterStack，每个第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, ...

当然，它可以像最后一行代码一样手动完成，但这对我来说似乎很不方便。如何以更好的方式实现我的目标的任何建议都将受到高度赞赏！

干杯，
弗洛里安

Answer 1

你可能在这里mget，因为它需要character对象名称向量并返回对象，所以你可以这样做：

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

nlayers( big.stack )
#[1] 144

如果您不喜欢使用replicate()，或者使用list将mget()堆叠为list，请使用list将它们放入rasterLayers stack()是ll <- replicate( 12 , s ) big.stack2 <- stack( ll ) identical( big.stack , big.stack2 ) #[1] TRUE ...

的有效输入

{{1}}

多次堆叠现有的RasterStack

1 个答案: