diff options
| -rw-r--r-- | acispy/image.py | 22 | ||||
| -rwxr-xr-x | bin/make_expmap.py | 20 |
2 files changed, 25 insertions, 17 deletions
diff --git a/acispy/image.py b/acispy/image.py new file mode 100644 index 0000000..0bcaf67 --- /dev/null +++ b/acispy/image.py @@ -0,0 +1,22 @@ +# Copyright (c) 2017 Weitian LI <liweitianux@live.com> +# MIT license + +""" +FITS image utilities +""" + +import subprocess + + +def get_xygrid(image): + """ + Get the ``xygrid`` of the input image. + """ + subprocess.check_call(["punlearn", "get_sky_limits"]) + subprocess.check_call([ + "get_sky_limits", "image=%s" % image, "verbose=0" + ]) + xygrid = subprocess.check_output([ + "pget", "get_sky_limits", "xygrid" + ]).decode("utf-8").strip() + return xygrid diff --git a/bin/make_expmap.py b/bin/make_expmap.py index ff1fbb0..1df4766 100755 --- a/bin/make_expmap.py +++ b/bin/make_expmap.py @@ -51,28 +51,13 @@ from acispy.manifest import get_manifest from acispy.pfiles import setup_pfiles from acispy.acis import ACIS from acispy.header import write_keyword +from acispy.image import get_xygrid logging.basicConfig(level=logging.INFO) logger = logging.getLogger(__name__) -def get_xygrid(image): - """ - Get the ``xygrid`` of the input image for later ``mkexpmap``'s use. - """ - logger.info("Get the 'xygrid' of image: %s" % image) - subprocess.check_call(["punlearn", "get_sky_limits"]) - subprocess.check_call([ - "get_sky_limits", "image=%s" % image, "verbose=0" - ]) - xygrid = subprocess.check_output([ - "pget", "get_sky_limits", "xygrid" - ]).decode("utf-8").strip() - logger.info("xygrid: %s" % xygrid) - return xygrid - - def make_aspect_histogram(outfile, asol, evtfile, chip, clobber=False): """ Create the aspect histogram for each chip, detailing the aspect @@ -211,6 +196,7 @@ def main(): logger.info("msk: %s" % msk) xygrid = get_xygrid(args.infile) + logger.info("%s:xygrid: %s" % (args.infile, xygrid)) expmaps = [] for c in chips: @@ -237,7 +223,7 @@ def main(): logger.info("Update keyword 'DETNAM' to %s" % detnam) write_keyword(args.outfile, keyword="DETNAM", value=detnam) - # Add created exposure map and exposure-corrected image to manifest + logger.info("Add created exposure map to manifest ...") key = "expmap" manifest.setpath(key, args.outfile) logger.info("Added '%s' to manifest: %s" % (key, manifest.get(key))) |