Add  Skycoord.to_pixels_wcs

This is a fairly trivial add on to #2885 that adds a convinence method to go from SkyCoord to pixel coordinates.

Not that there is *not* an inverse method, because it would make most sense as a classmethod, and `astropy.wcs.utils.pixel_to_skycoord` only uses the `SkyCoord` object explicitly.  It wouldn't be too hard to change that, though... What do you think, @astrofrog?
local imports do significantly slow down things. If, e.g., you are repeatedly converting a set of reference coordinates to pixels in order to do astrometry on a frame, this could be costly. That said, I would have feared a circular import just going through `wcs` at the top. Great that that doesn't happen!