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)