Experimental rotation.setter for PolygonPixelRegion

Author: dhomeier

regions/shapes/polygon.py

@@ -70,6 +70,7 @@ def __init__(self, vertices, meta=None, visual=None,
         self.visual = visual or RegionVisual()
         self.origin = origin
         self.vertices = vertices + origin
+        self._rotation = 0.0 * u.degree
 
     @property
     def area(self):
@@ -301,6 +302,24 @@ def rotate(self, center, angle):
         vertices = self.vertices.rotate(center, angle)
         return self.copy(vertices=vertices)
 
+    @property
+    def rotation(self):
+        """
+        Rotation angle to apply in-place rotations (operating on this instance).
+        Since `.setter` will apply the rotation directly on the vertices, this
+        value will always be reset to 0.
+        """
+        return self._rotation
+
+    @rotation.setter
+    def rotation(self, angle):
+        self.vertices = self.vertices.rotate(self.centroid, angle - self._rotation)
+        self._rotation = 0.0 * u.degree
+        if hasattr(self, '_mpl_selector'):
+            self._mpl_selector.verts = list(zip(self.vertices.x, self.vertices.y))
+            if getattr(self, '_mpl_selector_callback', None) is not None:
+                self._mpl_selector_callback(self)
+
 
 class RegularPolygonPixelRegion(PolygonPixelRegion):
     """

regions/shapes/tests/test_polygon.py

@@ -115,6 +115,22 @@ def test_origin(self):
         reg2 = PolygonPixelRegion(relverts, origin=origin)
         assert_equal(reg1.vertices, reg2.vertices)
 
+    def test_rotation(self):
+        """Test 'in-place' rotation of polygon instance, including full rotation"""
+        self.reg.rotation = 90 * u.deg
+        assert_allclose(self.reg.vertices.x, [8/3., 8/3., -1/3.], rtol=1e-9)
+        assert_allclose(self.reg.vertices.y, [4/3., 10/3., 4/3.], rtol=1e-9)
+        assert_allclose(self.reg.rotation, 0 * u.deg, rtol=1e-9)
+
+        self.reg.rotation = 90 * u.deg
+        assert_allclose(self.reg.vertices.x, [7/3., 1/3., 7/3.], rtol=1e-9)
+        assert_allclose(self.reg.vertices.y, [3, 3, 0], rtol=1e-9)
+
+        self.reg.rotation = 180 * u.deg
+        assert_allclose(self.reg.vertices.x, [1, 3, 1], rtol=1e-9)
+        assert_allclose(self.reg.vertices.y, [1, 1, 4], rtol=1e-9)
+        assert_allclose(self.reg.rotation, 0 * u.deg, rtol=1e-9)
+
 
 class TestPolygonSkyRegion(BaseTestSkyRegion):