Add position and amplitude as GaussianPSF constructor parameters
Implement photutils.datasets.make_gaussian_image
Plot each example image in the docs