aboutsummaryrefslogtreecommitdiffstats
path: root/bin
diff options
context:
space:
mode:
authorAaron LI <aaronly.me@outlook.com>2016-10-10 00:29:31 +0800
committerAaron LI <aaronly.me@outlook.com>2016-10-10 00:29:31 +0800
commitc1c5b3d5e899f8e385c90aae4ff881ecad3c7422 (patch)
treef8e1a627bb437ba04ef78558507f532f31181763 /bin
parentcbbd8eec608ac6d1b037e33ecbe0b76c4bfa6ccc (diff)
downloadfg21sim-c1c5b3d5e899f8e385c90aae4ff881ecad3c7422.tar.bz2
utils: Preseve the dtype when read/write FITS files
* utils/fits.py: hack the dtype to ignore the byteorder (FITS data use big endianness, e.g., dtype(">f4")) * utils/healpix.py: explicit convert the dtype and log the dtype * bin/healpix2hpx, bin/hpx2healpix: remove the --float argument * other minor fixes/updates
Diffstat (limited to 'bin')
-rwxr-xr-xbin/healpix2hpx5
-rwxr-xr-xbin/hpx2healpix8
2 files changed, 0 insertions, 13 deletions
diff --git a/bin/healpix2hpx b/bin/healpix2hpx
index c684e09..670e27a 100755
--- a/bin/healpix2hpx
+++ b/bin/healpix2hpx
@@ -28,8 +28,6 @@ def main():
parser.add_argument("outfile", help="output FITS image in HPX projection")
parser.add_argument("-C", "--clobber", action="store_true",
help="overwrite the existing output file")
- parser.add_argument("-F", "--float", action="store_true",
- help="use float (single precision) instead of double")
parser.add_argument("-l", "--log", dest="loglevel", default=None,
choices=["DEBUG", "INFO", "WARNING",
"ERROR", "CRITICAL"],
@@ -68,9 +66,6 @@ def main():
hpx_data, hpx_header = healpix2hpx(args.infile,
append_history=history,
append_comment=comments)
- if args.float:
- logger.info("HPX FITS images: use single-precision float numbers")
- hpx_data = hpx_data.astype(np.float32)
hdu = fits.PrimaryHDU(data=hpx_data, header=hpx_header)
hdu.writeto(args.outfile, clobber=args.clobber, checksum=True)
logger.info("HPX FITS images write to: %s" % args.outfile)
diff --git a/bin/hpx2healpix b/bin/hpx2healpix
index ac65fbd..ca485ac 100755
--- a/bin/hpx2healpix
+++ b/bin/hpx2healpix
@@ -13,9 +13,6 @@ import sys
import argparse
import logging
-import numpy as np
-from astropy.io import fits
-
import fg21sim
from fg21sim.configs import configs
from fg21sim.utils import hpx2healpix, write_fits_healpix, setup_logging
@@ -28,8 +25,6 @@ def main():
parser.add_argument("outfile", help="output HEALPix data file")
parser.add_argument("-C", "--clobber", action="store_true",
help="overwrite the existing output file")
- parser.add_argument("-F", "--float", action="store_true",
- help="use float (single precision) instead of double")
parser.add_argument("-l", "--log", dest="loglevel", default=None,
choices=["DEBUG", "INFO", "WARNING",
"ERROR", "CRITICAL"],
@@ -68,9 +63,6 @@ def main():
hp_data, hp_header = hpx2healpix(args.infile,
append_history=history,
append_comment=comments)
- if args.float:
- logger.info("HEALPix data: use single-precision float numbers")
- hp_data = hp_data.astype(np.float32)
write_fits_healpix(args.outfile, hpmap=hp_data, header=hp_header,
clobber=args.clobber)
logger.info("HEALPix data write to FITS file: %s" % args.outfile)