Merge branch 'develop' into copilot/fix-raster-arithmetic-result

Author: b-j-karl

src/rastr/raster.py

@@ -271,8 +271,12 @@ def __truediv__(self, other: float | Self | ArrayLike) -> Self:
         else:
             return NotImplemented
 
-    def __rtruediv__(self, other: float | ArrayLike) -> Self:
-        return self / other
+    def __rtruediv__(self, other: float) -> Self:
+        if not isinstance(other, float | int):
+            return NotImplemented
+        cls = self.__class__
+        new_arr = other / self.arr
+        return cls(arr=new_arr, raster_meta=self.raster_meta)
 
     def __sub__(self, other: float | Self | ArrayLike) -> Self:
         cls = self.__class__

tests/rastr/test_raster.py

@@ -911,15 +911,17 @@ def test_right_div_float(self):
                 transform=Affine(1.0, 0.0, 0.0, 0.0, 1.0, 0.0),
             )
             raster1 = Raster(
-                arr=np.array([[1, 2], [3, 4]]),
+                arr=np.array([[1, 2], [3, 4]], dtype=float),
                 raster_meta=raster_meta,
             )
 
             # Act
             result = 2.0 / raster1
 
             # Assert
-            np.testing.assert_array_equal(result.arr, np.array([[0.5, 1], [1.5, 2]]))
+            np.testing.assert_array_equal(
+                result.arr, np.array([[2.0 / 1, 2.0 / 2], [2.0 / 3, 2.0 / 4]])
+            )
 
         def test_shape_mismatch(self):
             # Arrange