fix get_stats when an empty mask is passed

Author: vschaffn

geoutils/raster/raster.py

@@ -1918,7 +1918,7 @@ def _statistics(self, band: int = 1, counts: tuple[int, int] | None = None) -> d
             "90th percentile": np.nanpercentile(mdata, 90),
             "LE90": linear_error(mdata, interval=90),
             "NMAD": nmad(data),
-            "RMSE": np.sqrt(np.ma.mean(np.square(data))),
+            "RMSE": np.sum(np.sqrt(np.ma.mean(np.square(data)))),
             "Standard deviation": np.ma.std(data),
             "Valid count": valid_count,
             "Total count": data.size,
@@ -1932,7 +1932,7 @@ def _statistics(self, band: int = 1, counts: tuple[int, int] | None = None) -> d
                     "Valid inlier count": valid_inlier_count,
                     "Total inlier count": counts[1],
                     "Percentage inlier points": (valid_inlier_count / counts[0]) * 100,
-                    "Percentage valid inlier points": (valid_inlier_count / counts[1]) * 100,
+                    "Percentage valid inlier points": (valid_inlier_count / counts[1]) * 100 if counts[1] != 0 else 0,
                 }
             )
 

tests/test_raster/test_raster.py

@@ -2002,11 +2002,13 @@ def test_stats(self, example: str, caplog) -> None:
             "Percentage valid points",
         ]
 
+        stat_types = (int, float, np.integer, np.floating)
+
         # Full stats
         stats = raster.get_stats()
         for name in expected_stats:
             assert name in stats
-            assert stats.get(name) is not None
+            assert isinstance(stats.get(name), stat_types)
 
         # With mask
         inlier_mask = raster.get_mask()
@@ -2018,10 +2020,21 @@ def test_stats(self, example: str, caplog) -> None:
             "Percentage valid inlier points",
         ]:
             assert name in stats_masked
-            assert stats_masked.get(name) is not None
+            assert isinstance(stats_masked.get(name), stat_types)
             stats_masked.pop(name)
         assert stats_masked == stats
 
+        # Empty mask
+        empty_mask = np.ones_like(inlier_mask)
+        stats_masked = raster.get_stats(inlier_mask=empty_mask)
+        for name in [
+            "Valid inlier count",
+            "Total inlier count",
+            "Percentage inlier points",
+            "Percentage valid inlier points",
+        ]:
+            assert stats_masked.get(name) == 0
+
         # Single stat
         for name in expected_stats:
             stat = raster.get_stats(stats_name=name)