flat_idx = np.ravel_multi_index((Px, Py), dims=a.shape)
unique_idx, idx_idx = np.unique(flat_idx, return_inverse=True)
deltas = np.bincount(idx_idx, weights=x)
Одновременное назначение одинаковых индексов массива в Python/Numpy