I have a question regarding creating WCS using the wcs.utils.fit_wcs_from_points function. The data we have has some obvious distortion. Therefore we would like to write the WCS of the data with a SIP projection, however even when I set the sip_degree to some number in the fit_wcs_from_points function, the resulting header doesn’t have the SIP projection keywords or any SIP coefficients in the header.
How might I go about creating a wcs object with a distortion?
Thanks in advance.