R中的并行马哈拉诺比斯光栅分类

Question

我正在寻找使用 Mahalanobis距离上课对栅格进行分类的最有效方法。我的课程已经通过均值和协方差进行了总结。

到目前为止，我正在将我的6波段输入堆栈图像作为数据帧进行调整，然后我使用stats::mahalanobis 我尝试过以下解决方案：

• 使用spatial.tools::rasterEngine
• raster::clusterR

是否有专用的包装或工具用于图像分类，并行实现马哈拉诺比斯距离？ 我有大约5000个8000×8000像素的光栅，因此性能至关重要。 欢迎任何想法！