def using_quadratic_loop(data):
    return np.array([[i, j] for i in range(1, np.shape(data)[0] - 1) for j in
        range(1, np.shape(data)[1] - 1) if np.all(data[i - 1:i + 2, j - 1:j +
        2] == data[i, j])]).T
Получить положение элементов с установкой некоторых критериев в numpy