def generate(ysize, xsize, ycenter, xcenter, radius):
x = np.arange(xsize)[None,:]
y = np.arange(ysize)[:,None]
return ((xcenter - x) ** 2 + (ycenter - y) ** 2 <= radius ** 2) * 1
generate(10,8,4,3,2)
array([[0, 0, 0, 0, 0, 0, 0, 0],
[0, 0, 0, 0, 0, 0, 0, 0],
[0, 0, 0, 1, 0, 0, 0, 0],
[0, 0, 1, 1, 1, 0, 0, 0],
[0, 1, 1, 1, 1, 1, 0, 0],
[0, 0, 1, 1, 1, 0, 0, 0],
[0, 0, 0, 1, 0, 0, 0, 0],
[0, 0, 0, 0, 0, 0, 0, 0],
[0, 0, 0, 0, 0, 0, 0, 0],
[0, 0, 0, 0, 0, 0, 0, 0]])