neighbors = np.concatenate([x[:, :-1].flatten(), x[:, +1:].flatten(), x[+1:,
    :].flatten(), x[:-1, :].flatten()])
centers = np.concatenate([x[:, +1:].flatten(), x[:, :-1].flatten(), x[:-1,
    :].flatten(), x[+1:, :].flatten()])
Подсчет ячеек соседних областей numpy