Add unit test for HpxNDMap.from_wcs_tiles()
@@ -369,3 +369,15 @@ def test_hpx_nd_map_to_wcs_tiles():
assert_allclose(tiles[0].data, 1)
assert_allclose(tiles[32].data, 1)
+
+def test_from_wcs_tiles():
+ geom = HpxGeom.create(nside=8)
+
+ wcs_geoms = geom.to_wcs_tiles(nside_tiles=4)
+
+ wcs_tiles = [Map.from_geom(geom, data=1) for geom in wcs_geoms]
+
+ m = HpxNDMap.from_wcs_tiles(wcs_tiles=wcs_tiles)
+
+ assert_allclose(m.data, 1)
+
