import matplotlib.pyplot as plt
from photutils.aperture import CircularAnnulus, CircularAperture
from photutils.datasets import make_100gaussians_image

data = make_100gaussians_image()
positions = [(145.1, 168.3), (84.5, 224.1), (48.3, 200.3)]
aperture = CircularAperture(positions, r=5)
annulus_aperture = CircularAnnulus(positions, r_in=10, r_out=15)
masks2 = annulus_aperture.to_mask(method='center')
fig, ax = plt.subplots()
ax.imshow(masks2[0], origin='lower')