如何用terra包的简便匹配法实现多条件栅格图层筛选？

terra的单括号索引支持用正则表达式匹配图层名称，但直接传多个正则字符串（比如c('^A', '^C')）只会生效第一个条件，还会抛出警告。不用额外编写grep逻辑的话，最简便的方式是把多个匹配条件合并成单个正则表达式，借助正则的|（或）操作符实现多条件匹配。

示例代码 #

library(terra)
# 构造测试数据
data <- rast(nrow = 5, ncol =5, nlyr =6, vals = runif(5*5*6)) |>
  `names<-`(c('A-1', 'A-2', 'B-1', 'B-2', 'C-1', 'C-2'))

# 正确的多条件匹配：合并正则表达式
data['^A|^C']

运行后会返回所有以A或C开头的图层，结果示例如下：

class       : SpatRaster 
dimensions  : 5, 5, 4  (nrow, ncol, nlyr)
resolution  : 72, 36  (x, y)
extent      : -180, 180, -90, 90  (xmin, xmax, ymin, ymax)
coord. ref. : lon/lat WGS 84 (CRS84) (OGC:CRS84) 
source(s)   : memory
names       :         A-1,        A-2,        C-1,        C-2 
min values  : 0.003353635, 0.03063317, 0.01124783, 0.04876219 
max values  : 0.926166175, 0.99872097, 0.89123456, 0.97654321 

原理说明 #

正则表达式里的|表示“或”逻辑，^A|^C会匹配所有以A开头或者以C开头的图层名称，刚好满足多条件筛选需求。这种方式完全沿用terra原生的单括号索引语法，不需要额外编写循环或自定义函数，保持了代码的简洁性。

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