大型数据集的栅格R提取功能

Question

我正在尝试为3公顷的面积为1平方米的网格（多边形）计算摘要统计信息（在这种情况下为平均值）。我的输入栅格数据集是从无人机以0.05 m的空间分辨率收集的。总共有20个覆盖带，覆盖3公顷，我想提取1 m网格的平均值。首先，我堆叠了所有栅格（速度很快），然后使用raster :: extract函数以及1 m2网格大小的shapefile来使用以下代码总结信息

VI_extract_mean <- raster::extract(final_stacked, data, fun=mean, df=TRUE, na.rm=TRUE)

这里； final_stacked：堆叠的栅格文件 数据：1平方米的网格形状文件

我有一台64 GB的RAM计算机。分析运行了8个多小时，但没有结果。问题是我有大量类似的数据集。加快过程的任何技巧。

Answer 1

您应该能够使用clusterR（）并利用计算机或环境中的多个内核。 MartinŠiklar的博客文章很好地满足了您的需求：https://www.gis-blog.com/increasing-the-speed-of-raster-processing-with-r-part-33-cluster/