Minor change

add_fits_site_info
Cees Bassa 2022-03-29 14:57:36 +02:00
parent bc24bbd88c
commit 91b6304c3a
1 changed files with 2 additions and 0 deletions

View File

@ -188,6 +188,7 @@ def add_wcs(fname, w, rmsx, rmsy):
# Read fits
hdu = fits.open(fname)
# Keywords to add
whdr = {"CRPIX1": w.wcs.crpix[0], "CRPIX2": w.wcs.crpix[1],
"CRVAL1": w.wcs.crval[0], "CRVAL2": w.wcs.crval[1],
"CD1_1": w.wcs.cd[0, 0], "CD1_2": w.wcs.cd[0, 1],
@ -199,6 +200,7 @@ def add_wcs(fname, w, rmsx, rmsy):
hdr = hdu[0].header
for k, v in whdr.items():
hdr[k] = v
hdu = fits.PrimaryHDU(header=hdr, data=hdu[0].data)
hdu.writeto(fname, overwrite=True, output_verify="ignore")