栅格计算中calc与stack组合相比直接调用目标函数有什么优势？

使用stack与calc组合的方式相比直接调用mean或算术运算计算栅格平均值，优势主要有以下几点：

• 代码更简洁易维护，可扩展性更强
  如果后续需要处理的图层数量变多（比如几十上百个），不需要手动逐个罗列L1到Ln的变量名，可以配合list.files()等函数批量读取指定路径下的所有栅格文件直接生成栅格栈，大幅降低手动输入变量的出错概率，后续增减图层也不需要修改统计部分的代码。
• 缺失值处理更灵活便捷
  用(L1+L2+……+L12)/12这种算术运算的方式时，只要有一个图层的对应像元值为NA，最终结果的该位置就会返回NA，无法灵活忽略缺失值。而calc配合mean(na.rm = T)可以自动仅统计该像元位置有有效值的图层计算平均值，不需要额外编写复杂的缺失值判断逻辑。
• 内存利用效率更高
  全球高分辨率栅格数据的单图层体量通常很大，直接算术运算需要一次性将所有图层全量加载到内存中，非常容易出现内存不足报错。而stack支持栅格数据延迟读取，calc运算时会自动分块处理数据，不需要一次性把全部图层加载到内存，更适合处理大体量的栅格数据集。
• 自带数据一致性校验
  调用stack堆叠多个栅格时，工具会自动校验所有图层的投影、分辨率、空间范围是否完全一致，不符合条件会直接报错提醒，避免误用参数不匹配的栅格运算得到错误结果还难以排查问题。
• 统计功能扩展性更强
  除了平均值之外，如果后续需要计算中位数、标准差、自定义统计量，只需要修改calc传入的函数参数即可，不需要重写整个运算逻辑，代码复用性更高。

内容的提问来源于stack exchange，提问作者Polly Wu