using AstroPy.wcs directly