Skip to content

fixing headers for LSST compatibility #50

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Open
martin-millon wants to merge 2 commits into
base: main
Choose a base branch
from
Open

fixing headers for LSST compatibility #50

Changes from all commits
Commits
Show all changes
2 commits
Select commit Hold shift + click to select a range
b8d4701
fixing headers
martin-millon Jul 7, 2025
a2f28e4
Merge branch 'main' into lsst_fix
martin-millon Jul 7, 2025
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
18 changes: 9 additions & 9 deletions lightcurver/processes/frame_importation.py
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -86,20 +86,20 @@ def process_new_frame(fits_file, user_config):
logger.info(f' file {fits_file}: background estimated.')

# before we write, let's keep as much as we can from our previous header
new_header = wcs.to_header()
new_header = header.copy()
# then copy non-WCS entries from the original header
for key in header:
if key not in new_header and not key.startswith('WCSAXES') and not key.startswith(
'CRPIX') and not key.startswith('CRVAL') and not key.startswith('CDELT') and not key.startswith(
'CTYPE') and not key.startswith('CUNIT') and not key.startswith('CD') and key not in ['COMMENT',
'HISTORY']:
if key.strip(): # some headers are weird
new_header[key] = header[key]
# for key in header:
# if key not in new_header and not key.startswith('WCSAXES') and not key.startswith(
# 'CRPIX') and not key.startswith('CRVAL') and not key.startswith('CDELT') and not key.startswith(
# 'CTYPE') and not key.startswith('CUNIT') and not key.startswith('CD') and key not in ['COMMENT',
# 'HISTORY']:
# if key.strip(): # some headers are weird
# new_header[key] = header[key]
# now we can write the file
write_file = user_config['workdir'] / copied_image_relpath
logger.info(f' Writing file: {write_file}')
fits.writeto(write_file, cutout_data_sub.astype(np.float32),
header=new_header, overwrite=True)
header=new_header, overwrite=True, output_verify="silentfix")
# and find sources
# (do plots if toggle set)
do_plot = user_config.get('source_extraction_do_plots', False)
Expand Down