Why Features On Wcs Projected Subplot Are In The Wrong Place In Matplotlib?

I have a fits file about an astronomical object. I can plot it like this: from astropy.io import fits from astropy.wcs import WCS hdul = fits.open(fitsfilename)[0] wcs = WCS(hdul

Solution 1:

To plot in world coordinates, you need to specify the transform, e.g.:

ax = fig.gca()
ax.scatter([34], [3.2], transform=ax.get_transform('world'))

Generally, you can ignore those FITS header warnings, since none of the FITS header cards there are related to the WCS (afaik).