微信公众号搜"智元新知"关注
微信扫一扫可直接关注哦!

结合或替换索引与 dplyr 过滤绘制光栅

如何解决结合或替换索引与 dplyr 过滤绘制光栅

我有一个 data.frame/tibble,其中一个变量是一个光栅砖列表,例如 a

library(dplyr)
library(raster)

x <- matrix(
  data = rnorm(100,1),nrow = 10,ncol = 10
) %>%
  raster

y <- matrix(
  data = rnorm(100,ncol = 10
) %>%
  raster

z <- matrix(
  data = rnorm(100,ncol = 10
) %>%
  raster

a <- tribble(
  ~thing,~map,"b",brick(x,y),"c",brick(z,"d",z),"e",brick(y,)

a
#> # A tibble: 4 x 2
#>   thing map       
#>   <chr> <list>    
#> 1 b     <RstrBrck>
#> 2 c     <RstrBrck>
#> 3 d     <RstrBrck>
#> 4 e     <RstrBrck>

我想从砖块中取出一些栅格来绘制,例如:

plot(a$map[[3]][[2]])

但我想使用 dplyr::filter 选择要绘制的适当行,而不是索引行,我可以使用 pull 到达适当的列,但我无法理解围绕如何获得砖对象内的适当光栅以对其进行绘制,而不是将其分配为新对象然后索引新对象。我想要做的是如下所示(失败):


a %>%
  filter(thing == "d") %>%
  pull(map)[[2]] %>%
  plot
#> Error in pull(map): object 'map' not found

解决方案?

reprex package (v0.3.0) 于 2021 年 2 月 3 日创建

解决方法

试试这个

a %>%
  filter(thing == "d") %>%
  pull(map) %>%
  .[[1]] %>% 
  .[[2]] %>% 
  plot

a %>%
  filter(thing == "d") %>%
  .$map %>% 
  .[[1]] %>% 
  .[[2]] %>% 
  plot

请注意,pull(map) 等价于 .$map

版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。