Thank you Melissa for offering to make a demo that would be quite helpful!
What I need:
I need to find the PSFs for each source in a FITS image so I can pass it as a parameter to pysersic’s FitSingle class. It needs to be graphable; and similar to the below screenshot (from the PSF in the examples folder on Pysersic Readthedocs Github). I have to do this locally, so I can’t use any lsst resources that are notebook-specific (butler).
What I have tried/the information I have:
The getPSF function in the 12a notebook tutorial uses a CoaddPSF object as shown below, so perhaps the CoaddPsf section in the FITS header should be used for reconstructing the PSF (but the PsfexPsf also looks promising).
I have tried printing out the CoaddPSFs (there are two included in the header; it seems as though one is metadata and the other one contains the actual PSF..) and their columns, and I can easily share the screenshots of those if helpful, just not sure how comfortable Rubin is with that information being public.
It appears as though someone tried to do something similar in a previous Rubin post but the code didn’t work for the psfex PSF header in my image and left me with a 2D array in which the same 3-column 1D array was repeated 255 times.
Thank you for your time :]