使用箭头赋值 function 作为 R purrr map2

Question

我有一个栅格列表：

filelist <- as.list(c("rasterA", "rasterB") # name them
rasterlist <- lapply(filelist, function(x) raster::raster(x)) # read them in

以及与这些栅格对应的 CRSes 列表：

crslist <- as.list(c("crsA.Rds", "crsB.Rds")) # name them
crslist %<>% lapply(function(x) readRDS(x)) # read them in

要在您使用的栅格中填充@crs插槽（例如）：

raster::crs(rasterA) <- crsA

但是我不知道如何在purrr::map2调用中使用这个箭头赋值。 R 备忘单中的示例是map2(x, y, sum) ，但是 sum 适用于向量，是无方向的，并且在其括号内包含它的所有术语。

阅读 map2 帮助，我尝试了公式格式：

purrr::map2(rasterlist, crslist, ~ raster::crs(.x) <- .y)

但没有骰子：

as_mapper(.f, ...) 中的错误：object 未找到“.y”

有谁知道我会怎么做？ 目前我使用的是一个简单的循环，但我正试图强迫自己学习映射并不断碰到这些类型的墙。

for (i in length(rasterlist)) {
 crs(rasterlist[[i]]) <- crslist[[i]]
}

谢谢！

起始条件：

for循环后：

在 raster_assigner function 之后改为：

希望这是可以挽救的，因为我还有其他类似的用例。 干杯!

Answer 1

我们可能需要

purrr::map2(rasterlist, crslist, ~ {crs(.x) <- .y;.x})

它与raster无关——因为下面的简单示例显示了相同的错误

> map2(list(1, 2, 3), list(3, 4, 5), ~ names(.x) <- .y)
Error in as_mapper(.f, ...) : object '.y' not found
> map2(list(1, 2, 3), list(3, 4, 5), ~ {names(.x) <- .y; .x})
[[1]]
3 
1 

[[2]]
4 
2 

[[3]]
5 
3 

然而，这不会在Map中返回错误（但这是不正确的，因为返回的 output 没有名称，除非我们用{}包装它并返回x

> Map(function(x, y) names(x) <- y, list(1, 2, 3), list(3, 4, 5))
[[1]]
[1] 3

[[2]]
[1] 4

[[3]]
[1] 5

Answer 2

可以做一个专用的 function：

raster_assigner <- function(raster, crs_input){
  raster::crs(raster) <- crs_input
  raster
}
purrr::map2(rasterlist, crslist, raster_assigner)