Skip to content
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

Jump to bottom

Data Only #331

Open
wants to merge 1 commit into
base: main
Choose a base branch
from
Open

Data Only #331

Changes from all commits
Commits
Show all changes
1 commit
Select commit
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
25 changes: 24 additions & 1 deletion EclipsingBinaries/IRAF_Reduction.py
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -1,7 +1,7 @@
"""
Author: Kyle Koeller
Created: 11/08/2022
Last Edited: 06/17/2023
Last Edited: 09/09/2023

This program is meant to automatically do the data reduction of the raw images from the
Ball State University Observatory (BSUO) and SARA data. The new calibrated images are placed into a new folder as to
Expand Down Expand Up @@ -337,9 +337,12 @@ def bias(files, calibrated_data, path, pipeline):

list_of_words = file_name.split(".")
new_fname = "{}.fits".format(list_of_words[0])
# original_header = new_ccd.header # save the original header

# Save the result
new_ccd.write(calibrated_data / new_fname, overwrite=overwrite)
# fits.writeto(calibrated_data / new_fname, new_ccd.data.astype('float32'), header=original_header,
# overwrite=overwrite)

add_header(calibrated_data, new_fname, "BIAS", "None", None, None, None, trim_region, overscan_region)

Expand All @@ -361,7 +364,10 @@ def bias(files, calibrated_data, path, pipeline):

fname = "zero.fits"
combined_bias.meta['combined'] = True
# original_header = combined_bias.header # save the original header
combined_bias.write(calibrated_data / fname, overwrite=overwrite)
# this only outputs the data image and not the mask or uncertainty and in single precision
# fits.writeto(calibrated_data / fname, combined_bias.data.astype('float32'), header=original_header, overwrite=overwrite)

add_header(calibrated_data, fname, "BIAS", "None", None, None, None, trim_region, overscan_region)

Expand Down Expand Up @@ -409,9 +415,12 @@ def dark(files, zero, calibrated_path, overscan_region, trim_region):
# new file name that uses the number from the original image
list_of_words = file_name.split(".")
new_fname = "{}.fits".format(list_of_words[0])
# original_header = sub_ccd.header # save the original header

# Save the result
sub_ccd.write(calibrated_path / new_fname, overwrite=overwrite)
# fits.writeto(calibrated_path / new_fname, sub_ccd.data.astype('float32'), header=original_header,
# overwrite=overwrite)

add_header(calibrated_path, new_fname, "DARK", "None", None, None, None, trim_region, overscan_region)

Expand All @@ -433,7 +442,11 @@ def dark(files, zero, calibrated_path, overscan_region, trim_region):

fname = "master_dark.fits"
combined_dark.meta['combined'] = True
# original_header = combined_dark.header # save the original header

combined_dark.write(calibrated_path / fname, overwrite=overwrite)
# this only outputs the data image and not the mask or uncertainty and in single precision
# fits.writeto(calibrated_path / new_fname, combined_dark.data.astype('float32'), header=original_header, overwrite=overwrite)

add_header(calibrated_path, fname, "DARK", "None", None, None, None, trim_region, overscan_region)

Expand Down Expand Up @@ -462,9 +475,12 @@ def flat(files, zero, combined_dark, calibrated_path, overscan_region, trim_regi
# new file name with the filter and number from the original file
list_of_words = file_name.split(".")
new_fname = "{}.fits".format(list_of_words[0])
original_header = final_ccd.header # save the original header

# Save the result
final_ccd.write(calibrated_path / new_fname, overwrite=overwrite)
# fits.writeto(calibrated_path / new_fname, final_ccd.data.astype('float32'), header=original_header,
# overwrite=overwrite)

add_header(calibrated_path, new_fname, "FLAT", "None", None, None, None, trim_region, overscan_region)

Expand All @@ -488,8 +504,12 @@ def flat(files, zero, combined_dark, calibrated_path, overscan_region, trim_regi

combined_flats.meta['combined'] = True
flat_file_name = 'master_flat_{}.fits'.format(filt.replace("Empty/", ""))
# original_header = combined_flats.header # save the original header

combined_flats.write(calibrated_path / flat_file_name, overwrite=overwrite)
# this only outputs the data image and not the mask or uncertainty and in single precision
# fits.writeto(calibrated_path / new_fname, combined_flats.data.astype('float32'),
# header=original_header, overwrite=overwrite)

add_header(calibrated_path, flat_file_name, "FLAT", "None", None, None, None, trim_region, overscan_region)

Expand All @@ -516,7 +536,10 @@ def science_images(files, calibrated_data, zero, combined_dark, trim_region, ove
new_fname = "{}.fits".format(list_of_words[0])

all_reds.append(reduced)
original_header = reduced.header # save the original header
reduced.write(calibrated_data / new_fname, overwrite=overwrite)
# this only outputs the data image and not the mask or uncertainty and in single precision
# fits.writeto(calibrated_data / new_fname, reduced.data.astype('float32'), header=original_header, overwrite=overwrite)

hjd = light.header["JD-HELIO"]
ra = light.header["RA"]
Expand Down
Loading