Visualizing images in sky coordinates using WCS in a notebook?

Is there a recommended way to visualize images processed with the Stack (e.g., the masked image of an Exposure) in sky coordinates using the WCS in a notebook? Is this possible with afw.display or another tool?

If you use one of the notebook magic commands that gives you interactive plots, like %matplotlib ipympl, you can hover over the afw.display-ed images and see the sky coordinates. If you want to plot objects at certain coordinates on top of the image, you can use the pixelToSky and/or skyToPixel methods available with the image’s WCS.

Thanks @mrawls. So it seems that one can get the transformation for specific coordinates in the image (either with the hover technique you mentioned, or by overplotting points with known sky coordinates into the pixel coordinates of the image), but maybe it is not possible to directly plot the image in another coordinate system with the current tooling?

@bechtol Do you have an example of what you have in mind?

I was imagining something like this:
Same image, WCS-based axes.

But are you wanting to warp the image to another projection?

When using Firefly for image visualization, driven by afw.display from a Rubin/LSST stack image, or “natively” against a data service, the WCS is known to the display, and sky coordinates (in equatorial or Galactic systems) for the location of the mouse are displayed continuously, or can be captured for a selected point. In addition, a coordinate grid can be plotted over the image in any of several systems.