用圆而不是矩形的binning 2D数据 - 来自熊猫DF

Question

我有一个X，Y数据的数据框，需要将其汇成圆圈。即一定尺寸的圆圈和间距以某个点为中心。因此，例如，在此采样/binning之后将遗漏一些数据。怎么可能？

我尝试了np.histogram2d并创建蒙版/广播。面具太慢了，我似乎无法广播圆圈。仅通过以下答案判断该点是否在圆圈的上述网格之内：将2D数据纳入x，y 中的重叠圆。

如果有一种方法可以将边缘或其他方法输入直方图2D并使边缘循环，请告诉我。干杯

Answer 1

这可以做到的唯一方法是将您的观点和网格循环循环：这样：

def inside_circle(x, y, x0, y0, r):
    return (x - x0)*(x - x0) + (y - y0)*(y - y0) < r*r


x_bins = np.linspace(-9, 9, 30)
y_bins = np.linspace(-9, 9, 30)

h = df['upmm']
w = df['anode_entrance']

histo = np.zeros((32,32))

for i in range(0, len(h)):
    for j in range(0, len(x_bins)):
        for k in range(0, len(y_bins)):
            if inside_circle(h[i], w[i], x_bins[j], y_bins[k], 0.01):
                histo[j][k] = histo[j][k] + 1
            
plt.imshow(histo, cmap='hot', interpolation='nearest')
plt.show()