def write_moc(moc_object):
"""
Write MOC to a buffer and return it.
"""
with closing(StringIO()) as f:
write_moc_fits(moc_object, f)
return f.getvalue()
def jsa_tile_moc(in_, out, max_order):
"""Create image MOC representation of a JSA tile."""
if in_.lower().endswith('.fits') or in_.lower().endswith('.fit'):
(data, tile_number, nside) = read_jsa_tile_fits(in_)
good = np.logical_not(np.isnan(data))
else:
(data, tile_number, nside) = read_jsa_tile_sdf(in_)
try:
good = data != np.finfo(data.dtype).min
except ValueError:
good = data != np.iinfo(data.dtype).min
(tile_order, map_suborder) = get_tile_params(data, nside)
pixels = get_pixel_set(good, tile_number, map_suborder)
moc = MOC(tile_order + map_suborder, pixels, moctype='image')
moc.normalize(max_order=max_order)
write_moc_fits(moc, out, clobber=True)
def jsa_tile_moc(in_, out, max_order):
"""Create image MOC representation of a JSA tile."""
if in_.lower().endswith('.fits') or in_.lower().endswith('.fit'):
(data, tile_number, nside) = read_jsa_tile_fits(in_)
good = np.logical_not(np.isnan(data))
else:
(data, tile_number, nside) = read_jsa_tile_sdf(in_)
try:
good = data != np.finfo(data.dtype).min
except ValueError:
good = data != np.iinfo(data.dtype).min
(tile_order, map_suborder) = get_tile_params(data, nside)
pixels = get_pixel_set(good, tile_number, map_suborder)
moc = MOC(tile_order + map_suborder, pixels, moctype='image')
moc.normalize(max_order=max_order)
write_moc_fits(moc, out, clobber=True)
