我有一个球体上的二维质量分布，作为经纬度网格单元中的质量矩阵给出。我需要将这些质量投影到同一个球体上的另一个网格，具有不同的极点位置和不同的单元格大小，再次作为网格单元格中的质量。

投影必须是：

1. 快速地

2. 单调，即不会产生实质性的人为浓度最大值/最小值（质量/细胞面积）

3. 严格质量守恒（数值精度）

我找不到任何现成的工具，尽管我相信应该有一些可用的东西（C、fortran 或其他）。有什么线索吗？

谢谢！