aboutsummaryrefslogtreecommitdiffstats
path: root/astro
diff options
context:
space:
mode:
authorAaron LI <aly@aaronly.me>2017-09-12 19:47:52 +0800
committerAaron LI <aly@aaronly.me>2017-09-12 19:47:52 +0800
commitf1560262d88a255fe4c63b8b4e7f171c4af12111 (patch)
tree6fe632c7ac66c0592004057154abd75854f1ce13 /astro
parent97bc6f61444c312c5326b5d6f5ed8771c8dfbeba (diff)
downloadatoolbox-f1560262d88a255fe4c63b8b4e7f171c4af12111.tar.bz2
taper_sky.py: Convert input image to be float32
Diffstat (limited to 'astro')
-rwxr-xr-xastro/oskar/taper_sky.py6
1 files changed, 3 insertions, 3 deletions
diff --git a/astro/oskar/taper_sky.py b/astro/oskar/taper_sky.py
index 282fe66..f70b47c 100755
--- a/astro/oskar/taper_sky.py
+++ b/astro/oskar/taper_sky.py
@@ -147,8 +147,8 @@ def main():
args.outfile_taper)
with fits.open(args.infile) as f:
- image = f[0].data
- header = f[0].header
+ image = f[0].data.astype(np.float32)
+ header = f[0].header.copy(strip=True)
logger.info("Read sky image from file: %s" % args.infile)
L = args.Router
@@ -189,7 +189,7 @@ def main():
header["OBJECT"] = ("Tukey Window", "Taper window")
header["Rinner"] = (args.Rinner, "[pixel] inner radius")
header["Router"] = (args.Router, "[pixel] outer radius")
- hdu = fits.PrimaryHDU(data=taper, header=header)
+ hdu = fits.PrimaryHDU(data=taper.astype(np.float32), header=header)
hdu.writeto(args.outfile_taper)
logger.info("Wrote tapered to file: %s" % args.outfile_taper)