def neuron_neighbours(neurons, center, radius):
    center_x, center_y = center
    nx, ny = neurons.shape

    y, x = np.ogrid[-center_x:nx - center_x, -center_y:ny - center_y]
    mask = (x * x + y * y) <= radius ** 2
    neurons[mask] = 1

    return neurons
Function find all neighbours neurons by radius and coords.
dataset/ETHPy150Open itdxer/neupy/neupy/algorithms/competitive/sofm.py/neuron_neighbours