如何按变量名拆分手动堆叠的terra包SpatRaster？

你通过c()手动堆叠带变量名的栅格生成SpatRaster后，尝试用lyrs或subds参数提取图层报错，是因为这两个参数仅用于从文件读取栅格时指定子数据集或图层，不适用于已加载到内存的SpatRaster对象。

要按变量名（varnames）拆分SpatRaster，可通过以下两种方法实现：

方法一：直接用逻辑索引提取 #

利用varnames()获取所有图层的变量名，生成逻辑向量筛选目标图层：

# 提取所有变量名为"mean"的图层
mean_layers <- r_stack[[varnames(r_stack) == "mean"]]

# 同理提取max或min图层
max_layers <- r_stack[[varnames(r_stack) == "max"]]
min_layers <- r_stack[[varnames(r_stack) == "min"]]

方法二：先获取索引再提取 #

先通过which()找到匹配变量名的图层索引，再用索引提取：

# 获取变量名为"mean"的图层索引
mean_idx <- which(varnames(r_stack) == "mean")
# 提取对应图层
mean_layers <- r_stack[[mean_idx]]

验证示例 #

针对你提供的测试代码，运行上述方法后，mean_layers会仅包含变量名为mean的r1栅格，查看其信息：

mean_layers

输出如下：

class       : SpatRaster 
dimensions  : 10, 10, 1  (nrow, ncol, nlyr)
resolution  : 1, 1  (x, y)
extent      : 0, 10, 0, 10  (xmin, xmax, ymin, ymax)
coord. ref. :  
sources     : memory  
varnames    : mean 
names       : lyr.1 
min values  :     1 
max values  :   100 

内容的提问来源于stack exchange，提问作者jd_fc