python-3.x - 如何将 3d 直方图的边缘值转换为直方图的中点值？

Question

如果为您提供与三维函数的 bin边缘相对应的hist, (x, y, z) = numpy.histogramdd()位置 ( x , y , z ) 的直方图值，如何计算（插值）中点的直方图值，即 ( x+d/2 , y +d/2 , z+d/2 ) 其中d是所有三个方向上的固定箱大小？

Answer 1

实际上，np.histogramdd为您提供 bin 边界x, y, z，但计数是“在中点”，而不是边界（严格来说，它们都不是，它们d x d x d位于以中点为中心的立方体上）。

如果——由于某种原因——你在边界处有值并想要插值，并且你对线性插值没问题：

np.lib.stride_tricks.as_strided(hist, (2, 2, 2, *map((-1).__add__, hist.shape)), 2 * hist.strides).mean(axis=(0, 1, 2))

这只是取 8 个最近邻居的平均值。