import matplotlib.pyplot as plt
from astropy.visualization import simple_norm
from photutils.datasets import make_100gaussians_image

data = make_100gaussians_image()
norm = simple_norm(data, 'sqrt', percent=99.5)
fig, ax = plt.subplots()
ax.imshow(data, norm=norm, origin='lower')
ax.set_title('Data')