def test_hpxgeom_read_write(tmp_path, nside, nested, frame, region, axes):
geom0 = HpxGeom(nside, nested, frame, region=region, axes=axes)
hdu_bands = geom0.to_bands_hdu(hdu_bands="TEST_BANDS")
hdu_prim = fits.PrimaryHDU()
hdu_prim.header.update(geom0.to_header())
hdulist = fits.HDUList([hdu_prim, hdu_bands])
hdulist.writeto(tmp_path / "tmp.fits")
with fits.open(tmp_path / "tmp.fits", memmap=False) as hdulist:
geom1 = HpxGeom.from_header(hdulist[0].header, hdulist["TEST_BANDS"])
assert_allclose(geom0.nside, geom1.nside)
assert_allclose(geom0.npix, geom1.npix)
assert_allclose(geom0.nest, geom1.nest)
assert geom0.frame == geom1.frame
def test_hpxgeom_from_header():
pars = {
"HPX_REG": "DISK(110.,75.,2.)",
"EXTNAME": "SKYMAP",
"NSIDE": 2**6,
"ORDER": 6,
"PIXTYPE": "HEALPIX",
"ORDERING": "RING",
"COORDSYS": "CEL",
"TTYPE1": "PIX",
"TFORM1": "K",
"TTYPE2": "CHANNEL1",
"TFORM2": "D",
"INDXSCHM": "EXPLICIT",
}
header = fits.Header()
header.update(pars)
hpx = HpxGeom.from_header(header)
assert hpx.frame == "icrs"
assert not hpx.nest
assert_allclose(hpx.nside, np.array([64]))
