star psf
>>> import matplotlib.pyplot as plt
>>> from astropy.visualization import simple_norm
>>> nrows = 5
>>> ncols = 5
>>> fig, ax = plt.subplots(nrows=nrows, ncols=ncols, figsize=(20, 20),
... squeeze=True)
>>> ax = ax.ravel()
>>> for i in range(nrows*ncols):
... norm = simple_norm(stars[i], 'log', percent=99.)
... ax[i].imshow(stars[i], norm=norm, origin='lower', cmap='viridis')