From c1c5b3d5e899f8e385c90aae4ff881ecad3c7422 Mon Sep 17 00:00:00 2001 From: Aaron LI Date: Mon, 10 Oct 2016 00:29:31 +0800 Subject: 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 --- bin/healpix2hpx | 5 ----- bin/hpx2healpix | 8 -------- 2 files changed, 13 deletions(-) (limited to 'bin') 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) -- cgit v1.2.2