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

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)