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)
masks = annulus_aperture.to_mask(method='exact')
fig, ax = plt.subplots()
ax.imshow(masks[0].multiply(data), origin='lower')