Substitute astropy.io.fits for pyfits.
In create_fitsheader, default to ICRS.
Implement create_fitsheader_for function, and deprecate the fromdata keyword in create_fitsheader.
Simplify projection_scale_gnomonic.