img = np.zeros((256, 256))
for i in range(data.shape[0]):
    img[data[i, 0], data[i, 1]] += 1
img = ndi.gaussian_filter(img, (10, 10))
Эффективный метод расчета плотности неравномерно расположенных точек