Adapt `psf.py` to new array utils api in `astropy.ndata.utils`

This PR  modifies `psf.py` to se the new array utility functions in `astropy.nddata.utils`, which are included in astropy/astropy#3201. It also includes the changes of #216. 
@larrybradley @cdeil