正如@Brian L的评论中所述，要获得每个单元格中的累积平方英尺值，您需要将reduce_C_function更改为sum。这意味着Matplotlib将使用Python sum函数（默认值是均值函数）来确定该bin的颜色，基于C（在本例中为SQFT）。下面是一个使用随机生成数据的示例。

import numpy as np
import matplotlib.pyplot as plt

plt.close("all")

N = 158                             # OP said 25k datapoints; sqrt(25000) ~ 158
latitude = np.random.standard_normal(N**2) 
longitude = 12 * np.random.standard_normal(N**2) 

var = 100
SQFT = 1_000*np.exp(-(latitude**2 + longitude**2)/2/var)/2/np.pi/var + 1

fig, ax = plt.subplots()
p = ax.hexbin(longitude, latitude, C=SQFT, gridsize=50, reduce_C_function=sum)
fig.colorbar(p, ax=ax, label="SQFT")
ax.set_xlabel("Longitude")
ax.set_ylabel("Latitude")
fig.show()