Filter out negative (missing) neighbors when computing healpix envelope.

lsst · Oct 29, 2024 · e18d690 · e18d690
1 parent 27dc372
commit e18d690
Showing 1 changed file with 2 additions and 2 deletions.
diff --git a/python/lsst/sphgeom/_healpixPixelization.py b/python/lsst/sphgeom/_healpixPixelization.py
@@ -89,9 +89,9 @@ def envelope(self, region: Region, maxRanges: int = 0):
 
         # Dilate the high resolution pixels by one to ensure that the full
         # region is completely covered at high resolution.
-        neighbors = hpg.neighbors(self._nside_highres, pixels_highres)
+        neighbors = hpg.neighbors(self._nside_highres, pixels_highres).ravel()
         # Shift back to the original resolution and uniquify
-        pixels = np.unique(np.right_shift(neighbors.ravel(), self._bit_shift))
+        pixels = np.unique(np.right_shift(neighbors[neighbors >= 0], self._bit_shift))
 
         return RangeSet(pixels)