import matplotlib.pyplot as plt
from astropy.visualization import simple_norm
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)

norm = simple_norm(data, 'sqrt', percent=99)
fig, ax = plt.subplots()
ax.imshow(data, norm=norm, origin='lower')
ax.set_xlim(0, 170)
ax.set_ylim(130, 250)

ap_patches = aperture.plot(color='white', lw=2,
                           label='Photometry aperture')
ann_patches = annulus_aperture.plot(color='red', lw=2,
                                    label='Background annulus')
handles = (ap_patches[0], ann_patches[0])
ax.legend(loc=(0.17, 0.05), facecolor='#458989', labelcolor='white',
          handles=handles, prop={'weight': 'bold', 'size': 11})