使用矩形网格数据的球体积积分

Question

我将矩形网格中的周期性数据存储为Numpy数组。我想将数据集成在半径r和中心a的范围内。

• 我可以使用任何scipy功能吗？
• 如果没有，最有效的方法是什么？