r - 根据栅格图层名称上的特定模式从栅格堆栈中选择栅格图层

Question

我在 r 中有一个栅格堆栈，其中包含 499 个图层。每一层都有一个特定的名称模式。

[1] "Sacharovi_PA1_RUN1_GAM"               "Sacharovi_PA1_RUN1_GLM"              
[3] "Sacharovi_PA1_RUN1_RF"                "Sacharovi_PA1_RUN1_CTA"              
[5] "Sacharovi_PA1_RUN1_MAXENT.Phillips"   "Sacharovi_PA1_RUN2_GAM"              
[7] "Sacharovi_PA1_RUN2_GLM"               "Sacharovi_PA1_RUN2_RF"               
[9] "Sacharovi_PA1_RUN2_CTA"                 

我想要做的是将每个模型(GAM、RF、GLM、CTA、MAXENT)保存到不同的堆栈。如何选择我想要的模型(所有包含“GAM”/“GLM”/“RF”等词的模型)？

到目前为止，我一直在尝试这样处理光栅堆栈的一个子集:

  result <- subset(my_stack, grep("GAM"))

但不会起作用。你能帮我解决这个问题吗？

Answer 1

您似乎没有正确完成 grep - 您必须告诉它搜索名称，见下文。还要确保您通过指定栅格包访问正确的函数。 Subset 是一个非常通用的函数名称，它可能已被您在 raster 之后添加的另一个库所取代 - 或者如果您根本没有加载该库，R 将默认为 base::subset。

我也更喜欢使用 value = TRUE 进行调试，尽管这两种方式都可以。

all_GAM <- raster::subset(my_stack, grep('_GAM', names(my_stack), value = T))