lsst-sphgeom 30.0.0rc2__tar.gz → 30.0.1__tar.gz

{lsst_sphgeom-30.0.0rc2/python/lsst_sphgeom.egg-info → lsst_sphgeom-30.0.1}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: lsst-sphgeom
-Version: 30.0.0rc2
+Version: 30.0.1
 Summary: A spherical geometry library.
 Author-email: Rubin Observatory Data Management <dm-admin@lists.lsst.org>
 License-Expression: BSD-3-Clause OR GPL-3.0-or-later

{lsst_sphgeom-30.0.0rc2 → lsst_sphgeom-30.0.1}/python/lsst/sphgeom/_healpixPixelization.py

@@ -150,7 +150,17 @@ class HealpixPixelization(PixelizationABC):
                 raise ValueError("Invalid circle radius.")
         elif isinstance(region, ConvexPolygon):
             vertices = np.array([[v.x(), v.y(), v.z()] for v in region.getVertices()])
-            pixels = hpg.query_polygon_vec(nside, vertices)
+            retry = False
+            try:
+                pixels = hpg.query_polygon_vec(nside, vertices)
+            except RuntimeError:
+                # This happens if the polygon is determined to have a
+                # degenerate corner (which can happen with rounding).
+                # In this case, redo with a bounding circle.
+                retry = True
+
+            if retry:
+                pixels = self._interior_pixels_from_region(nside, region.getBoundingCircle())
         elif isinstance(region, Box):
             pixels = hpg.query_box(
                 nside,

{lsst_sphgeom-30.0.0rc2 → lsst_sphgeom-30.0.1}/python/lsst/sphgeom/version.py

@@ -1,2 +1,2 @@
 __all__ = ["__version__"]
-__version__ = "30.0.0rc2"
+__version__ = "30.0.1"

{lsst_sphgeom-30.0.0rc2 → lsst_sphgeom-30.0.1/python/lsst_sphgeom.egg-info}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: lsst-sphgeom
-Version: 30.0.0rc2
+Version: 30.0.1
 Summary: A spherical geometry library.
 Author-email: Rubin Observatory Data Management <dm-admin@lists.lsst.org>
 License-Expression: BSD-3-Clause OR GPL-3.0-or-later

{lsst_sphgeom-30.0.0rc2 → lsst_sphgeom-30.0.1}/tests/test_HealpixPixelization.py

@@ -215,6 +215,51 @@ class HealpixPixelizationTestCase(unittest.TestCase):
         )
         self._check_interior(h, ellipse, [pix])
 
+    def test_interior_nearly_degenerate_polygon(self):
+        """Check interior pixels with nearly degenerate polygon.
+
+        See DM-53933.
+        """
+        poly = ConvexPolygon(
+            [
+                UnitVector3d(-0.016144196513206758, -0.9235039286924489, -0.3832490816799893),
+                UnitVector3d(-0.01988138785117867, -0.9228560783513854, -0.38463149775728545),
+                UnitVector3d(-0.022348863683228713, -0.9209177038528237, -0.3891158066983548),
+                UnitVector3d(-0.024841110338799998, -0.918545770328335, -0.3945333788782152),
+                UnitVector3d(-0.026637426968850634, -0.9164294757774608, -0.39930873194901123),
+                UnitVector3d(-0.02509225617228684, -0.9149335621264763, -0.4028237276709774),
+                UnitVector3d(-0.0221015078570923, -0.911998223315482, -0.4095982959191202),
+                UnitVector3d(-0.020553220684181986, -0.9104571830599119, -0.4130923418972051),
+                UnitVector3d(-0.017371565494301272, -0.9072522623419655, -0.4202279871541906),
+                UnitVector3d(-0.014347349160349979, -0.9041551644174409, -0.4269632211670476),
+                UnitVector3d(-0.014303105070771129, -0.9041094310633312, -0.42706153871271785),
+                UnitVector3d(-0.012802279210872009, -0.9025514012895127, -0.4303917630221838),
+                UnitVector3d(-0.007862275539177377, -0.9015265103248066, -0.43265244227315125),
+                UnitVector3d(-0.0019484762556141095, -0.900471339489277, -0.4349109911219402),
+                UnitVector3d(0.0032282853511799436, -0.8996943229589627, -0.436508537613075),
+                UnitVector3d(0.006969847869770022, -0.9003377930402403, -0.4351359323752772),
+                UnitVector3d(0.022099717156107892, -0.9027944227880381, -0.42950417074160424),
+                UnitVector3d(0.033589177231523826, -0.9045049275054283, -0.4251383342999184),
+                UnitVector3d(0.03732521827078148, -0.9050320420971739, -0.4237025263772424),
+                UnitVector3d(0.039795546158565454, -0.9070052460628258, -0.41923477684998206),
+                UnitVector3d(0.04228779372661228, -0.9093772751891867, -0.41381724694752126),
+                UnitVector3d(0.04409395578698435, -0.9114581951867685, -0.4090228373696683),
+                UnitVector3d(0.04255005806475871, -0.9130743444969162, -0.4055671756691021),
+                UnitVector3d(0.04109145924730725, -0.91458405813035, -0.4023027374885082),
+                UnitVector3d(0.03791490987091337, -0.9178181221381234, -0.3951864803916365),
+                UnitVector3d(0.03335712714052637, -0.9223387473238281, -0.38493965404208763),
+                UnitVector3d(0.03180762744041508, -0.9238417869862389, -0.3814506881035668),
+                UnitVector3d(0.030248610691484885, -0.9253347742243898, -0.3779425580195124),
+                UnitVector3d(0.02530895801451361, -0.9263964395908428, -0.3756981412751856),
+                UnitVector3d(0.01939515964342646, -0.9274517060282444, -0.37343963470379626),
+                UnitVector3d(0.014215536919345694, -0.9281937989444676, -0.37182548340738003),
+                UnitVector3d(0.010475126504663994, -0.9276710642902235, -0.3732514811803906),
+            ],
+        )
+
+        h = HealpixPixelization(6)
+        self._check_interior(h, poly, [28890, 28891, 28912, 28913, 28914])
+
     def _check_interior(self, pixelization, region, check_pixels):
         """Check the interior from a region.