Add new functions to make images with Gaussian sources

This PR adds two new functions for making images with Gaussian sources.

`make_gaussian_sources` extends the existing `make_gaussian_image` by allowing elliptical Gaussians, including optional Gaussian or Poisson noise, and allowing an optional `seed` to reproduce the random noise images.

`make_random_gaussians` generates random 2D elliptical Gaussians using `make_gaussian_image`.

Here are the docs for the new functions:
* http://www.stsci.edu/~lbradley/photutils/make_gaussians/api/photutils.datasets.make_gaussian_sources.html#photutils.datasets.make_gaussian_sources
* http://www.stsci.edu/~lbradley/photutils/make_gaussians/api/photutils.datasets.make_random_gaussians.html#photutils.datasets.make_random_gaussians


@cdeil Do you mind if I remove `make_gaussian_image` since that functionality (circular Gaussians) is covered by the `make_gaussian_sources`?
I should add that most optical/NIR images I work with are "background limited" where the Poisson source noise is small compared to the background.  So this is a good approximation in that case.  I'm sure gamma-ray images are much different!