Tibeuleu/FOC_Reduction
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

fix WCS computation, cdelt should not be sorted

Browse Source
This commit is contained in:
Thibault Barnouin
2025-04-01 17:02:31 +02:00
parent 84bbdb9634
commit 075dc4fd15
2 changed files with 4 additions and 4 deletions
Download Patch File Download Diff File Expand all files Collapse all files

6
package/lib/fits.py
View File

@@ -16,7 +16,7 @@ from astropy.io import fits
from astropy.wcs import WCS from astropy.wcs import WCS
from .convex_hull import clean_ROI from .convex_hull import clean_ROI
from .utils import wcs_PA from .utils import wcs_PA, princ_angle
def get_obs_data(infiles, data_folder="", compute_flux=False): def get_obs_data(infiles, data_folder="", compute_flux=False):
@@ -72,13 +72,13 @@ def get_obs_data(infiles, data_folder="", compute_flux=False):
for key in keys: for key in keys:
header.remove(key, ignore_missing=True) header.remove(key, ignore_missing=True)
new_cdelt = np.linalg.eigvals(wcs.wcs.cd) new_cdelt = np.linalg.eigvals(wcs.wcs.cd)
new_cdelt.sort() # new_cdelt.sort()
new_wcs.wcs.pc = wcs.wcs.cd.dot(np.diag(1.0 / new_cdelt)) new_wcs.wcs.pc = wcs.wcs.cd.dot(np.diag(1.0 / new_cdelt))
new_wcs.wcs.cdelt = new_cdelt new_wcs.wcs.cdelt = new_cdelt
for key, val in new_wcs.to_header().items(): for key, val in new_wcs.to_header().items():
header[key] = val header[key] = val
try: try:
_ = header["ORIENTAT"] header["ORIENTAT"] = princ_angle(float(header["ORIENTAT"]))
except KeyError: except KeyError:
header["ORIENTAT"] = wcs_PA(new_wcs.wcs.pc[1, 0], np.diag(new_wcs.wcs.pc).mean()) header["ORIENTAT"] = wcs_PA(new_wcs.wcs.pc[1, 0], np.diag(new_wcs.wcs.pc).mean())

2
package/lib/utils.py
View File

@@ -177,4 +177,4 @@ def wcs_PA(PC21, PC22):
orient = np.arccos(PC22) * 180.0 / np.pi orient = np.arccos(PC22) * 180.0 / np.pi
elif (abs(PC21) < abs(PC22)) and (PC22 < 0): elif (abs(PC21) < abs(PC22)) and (PC22 < 0):
orient = -np.arccos(PC22) * 180.0 / np.pi orient = -np.arccos(PC22) * 180.0 / np.pi
return orient return princ_angle(orient)