在 python 中重新创建 NCL bin_avg 函数

Recreate NCL bin_avg function in python

我正在尝试在 python 中重新创建 NCL 的 bin_avg 函数。该函数接受这些点的纬度、经度和观测值的列表,并计算某个网格中每个框内包含的观测值的平均值。

我能想到的最接近的方法是使用 scipy.spatial.cKDTree.query_ball_point,但这会找到半径内的所有点,而不是网格框内的所有点...

看来我要找的东西在 scipy.stats.binned_statistic