Commit 1429a6dd authored by Leo Pound Singer's avatar Leo Pound Singer

All tools now output multi-resolution HEALPix files

parent 4fa2c8e8
......@@ -10,6 +10,9 @@ Changelog
- Un-pin Numpy version now that Numpy 1.15.4 is out.
- The ``bayestar-localize-lvalert`` and ``ligo-skymap-from-samples`` tools will
now natively output multi-resolution HEALPix files.
0.0.17 (2018-10-24)
===================
......
......@@ -49,7 +49,7 @@ def parser():
'--no-tag', action='store_true',
help='Do not set lvem tag for GraceDB entry')
parser.add_argument(
'-o', '--output', metavar='FILE.fits[.gz]', default='bayestar.fits.gz',
'-o', '--output', metavar='FILE.fits', default='bayestar.fits',
help='Name for uploaded file')
parser.add_argument(
'graceid', metavar='G123456', nargs='*',
......@@ -67,7 +67,7 @@ def main(args=None):
import sys
import tempfile
import urllib.parse
from ..bayestar import localize, rasterize
from ..bayestar import localize
from ..io import fits
from ..io import events
from ..util.file import rename
......@@ -133,12 +133,12 @@ def main(args=None):
try:
# perform sky localization
log.info("starting sky localization")
sky_map = rasterize(localize(
sky_map = localize(
event, opts.waveform, opts.f_low, opts.min_distance,
opts.max_distance, opts.prior_distance_power, opts.cosmology,
mcmc=opts.mcmc, chain_dump=chain_dump,
enable_snr_series=opts.enable_snr_series,
f_high_truncate=opts.f_high_truncate))
f_high_truncate=opts.f_high_truncate)
sky_map.meta['objid'] = str(graceid)
sky_map.meta['url'] = '{}/{}'.format(base_url, graceid)
log.info("sky localization complete")
......
......@@ -66,7 +66,6 @@ def main(args=None):
# Late imports
from .. import io
from ..bayestar import rasterize
from .. import version
from astropy.table import Table
from astropy.time import Time
......@@ -124,7 +123,7 @@ def main(args=None):
skypost.multiprocess = args.j
log.info('making skymap')
hpmap = rasterize(skypost.as_healpix())
hpmap = skypost.as_healpix()
hpmap.meta.update(io.fits.metadata_for_version_module(version))
hpmap.meta['creator'] = _parser.prog
hpmap.meta['origin'] = 'LIGO/Virgo'
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment