Hi!
I am writing again with a new difficulty that I found when using gnuastro on my files. I am trying to run the subtracting the PSF tutorial on some HST images. When I run the astscript-psf-select-stars script on any of my individual HST frames, the following error appears. I attached a minimal reproducible jupyter-notebook that downloads and executes the code to this email for reference:
astquery: gal_wcs_read_fitsptr: (warning) wcsfix status for cdfix: -1
astquery: gal_wcs_read_fitsptr: (warning) wcsfix status for obsfix: -1
astquery: gal_wcs_read_fitsptr: (warning) wcsfix status for unitfix: -1
astquery: gal_wcs_read_fitsptr: (warning) wcsfix status for spcfix: -1
astquery: gal_wcs_read_fitsptr: (warning) wcsfix status for celfix: 5
astquery: gal_wcs_read_fitsptr: (warning) wcsfix status for cylfix: -1
astquery: gal_wcs_read_fitsptr: WCSLIB warning: wcsset error 5: Invalid parameter value
astquery: ./mastDownload/HST/j9en0pnaq/j9en0pnaq_flc.fits (hdu 1): contains no WCS to derive the sky coverage
I am using Gnuastro 0.19 on a MacOS 12.6.2. I am guessing that WCSlib can't deal correctly with the distorted frames of HST, and that the PSF subtraction should be run in the final mosaics, or at least the undistorted frames. That would be less than ideal, since the PSF is better defined in the original frames before drizzling. But before going to the workaround, do you think that there is a solution to this error?
Thank you!
--
==================================================================
Ex Luna, Scientia
Alejandro S. Borlaff, PhD - NASA / BAERI Researcher
NASA Space Science and Astrobiology Division (Code ST)
NASA Ames Research Center, N245-6, Moffett Field, California (94035) USA
Webpage :
https://borlaff.github.io/GitHub:
https://github.com/BorlaffTwitter: @asborlaff
==================================================================