From 4b6b2ef3d5498dd06ca66bb007178fcb9371c8bb Mon Sep 17 00:00:00 2001
From: Philipp Arras <parras@mpa-garching.mpg.de>
Date: Tue, 26 Oct 2021 15:46:42 +0200
Subject: [PATCH] Add compulsory fits header entries
---
resolve/fits.py | 16 +++++++++++++++-
1 file changed, 15 insertions(+), 1 deletion(-)
diff --git a/resolve/fits.py b/resolve/fits.py
index 4cd012d1..21554802 100644
--- a/resolve/fits.py
+++ b/resolve/fits.py
@@ -36,10 +36,24 @@ def field2fits(field, file_name, overwrite, direction=None):
h["CDELT2"] = dom.distances[1] * 180 / np.pi
h["CRPIX2"] = dom.shape[1] / 2
h["CUNIT2"] = "deg"
+ # Compulsory header entries
+ h["SIMPLE"] = "T"
+ h["BITPIX"] = 64
+ h["NAXIS"] = 2
+ h["NAXIS1"] = dom.shape[0]
+ h["NAXIS2"] = dom.shape[1]
+ # /Compulsory header entries
+
h["DATE-MAP"] = Time(time.time(), format="unix").iso.split()[0]
if direction is not None:
h["EQUINOX"] = direction.equinox
- hdu = pyfits.PrimaryHDU(field.val[:, :].T, header=h)
+ # physical value = bzero + bscale * array value
+ h["BSCALE"] = 1.
+ h["BZERO"] = 0.
+ # physical unit
+ h["BUNIT"] = "Jy/sr"
+
+ hdu = pyfits.PrimaryHDU(field.val.T, header=h)
hdulist = pyfits.HDUList([hdu])
base, ext = splitext(file_name)
hdulist.writeto(base + ext, overwrite=overwrite)
--
GitLab