def apply_hpxmoc(self, lon=None, lat=None):
try:
moc = MOC()
moc.read(self.filename, filetype="fits")
except:
raise Exception(f"Unable to find/open Healpix MOC file: {self.filename}")
# get the moc nside at the max resolution of the MOC
nside = hp.order2nside(moc.order)
#get the healpix pixels indices of the targets at the max resolution
idx = hp.ang2pix(nside, lon, lat, lonlat=True, nest=True )
m_mask = moc.contains(idx)
return m_mask
def test_contains(self):
m = MOC()
m.add(0, (10, 11))
m.add(1, (36, 37))
m.add(2, (128, 129))
m.add(3, (448, 499))
self.assertEqual(m.contains(0, 10), True)
self.assertEqual(m.contains(0, 11), True)
self.assertEqual(m.contains(0, 0, True), False)
self.assertEqual(m.contains(0, 0, False), False)
self.assertEqual(m.contains(1, 40), True)
self.assertEqual(m.contains(2, 160), True)
self.assertEqual(m.contains(0, 7, True), True)
self.assertEqual(m.contains(0, 7, False), False)
def test_contains(self):
m = MOC()
m.add(0, (10, 11))
m.add(1, (36, 37))
m.add(2, (128, 129))
m.add(3, (448, 499))
self.assertEqual(m.contains(0, 10), True)
self.assertEqual(m.contains(0, 11), True)
self.assertEqual(m.contains(0, 0, True), False)
self.assertEqual(m.contains(0, 0, False), False)
self.assertEqual(m.contains(1, 40), True)
self.assertEqual(m.contains(2, 160), True)
self.assertEqual(m.contains(0, 7, True), True)
self.assertEqual(m.contains(0, 7, False), False)
