在特定 window 中整合热图值

Integrating heatmap values within specific window

我正在绘制二维函数,使用:

import numpy as np
import matplotlib.pyplot as plt

# Build x,y mesh
nx = 1000; ny = 1000 
x = np.linspace(-0.5, 0.5, nx)
y = np.linspace(-0.5, 0.5, ny)
X,Y = np.meshgrid(x, y)

# Evaluate the Gaussian on grid
func = np.exp(-np.pi*((X)**2 + (Y)**2))
plt.imshow(func); plt.colorbar()

此时,我在网格的 xy 向量的整个范围内集成函数,包括正确的 dxdy 递增。值如标题所示:

# Integrate the function
intval = np.trapz(np.trapz(func, y, dx=x[1]-x[0]), x, dx=x[1]-x[0])
plt.title('Integrated value = %s' % intval)

如何在 window 的 部分 上集成相同的功能?比如,在 x[400:600] 和 y[400:600]?

之间

只传递数组的那部分怎么样?

cut = slice(400, 600)
intval = np.trapz(np.trapz(func[cut, cut], y[cut], dx=x[1]-x[0]),
                  x[cut], dx=x[1]-x[0])