如何在terra::mosaic()中使用自定义函数？

Question

我有数千个相邻和（部分）重叠的geotiff，我想将它们镶嵌到一个geotiff中。我需要使用“模式”来计算新的 geotiff 的值，以防像素重叠。

我在R的terra包中找到了mosaic()函数，但是“mode”不在“fun”参数列表中。有没有办法提供自定义功能？

顺便说一下，我担心是否可以使用此功能来处理如此大的数据（即 6000 张 512x512 像素的图像），马赛克后的数据约为。 10000 x 10000 像素。因此，如果有人提出替代解决方案（例如 gdal_calc_py？），我们将不胜感激。

Answer 1

您（尚）无法使用自定义函数，但使用“terra”版本 1.7.57，您现在可以使用

fun="modal"

。

这是目前的开发版本。您可以从 R-Universe 安装此版本的 terra

install.packages('terra', repos='https://rspatial.r-universe.dev')

你的问题似乎不是特别大，“terra”被设计用于处理大数据。