nbmorph 0.3.0__tar.gz → 0.3.1__tar.gz

{nbmorph-0.3.0/src/nbmorph.egg-info → nbmorph-0.3.1}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: nbmorph
-Version: 0.3.0
+Version: 0.3.1
 Summary: A small package with Numba-accelerated morphological operations.
 Author-email: Marius Causemann <mariusca@simula.no>
 License: MIT License

{nbmorph-0.3.0 → nbmorph-0.3.1}/pyproject.toml

@@ -2,7 +2,7 @@
 
 [project]
 name = "nbmorph"
-version = "0.3.0"
+version = "0.3.1"
 authors = [
   { name="Marius Causemann", email="mariusca@simula.no" },
 ]

{nbmorph-0.3.0 → nbmorph-0.3.1}/src/nbmorph/box_kernel.py

@@ -19,7 +19,7 @@ def min27(v0, v1, v2, v3, v4, v5, v6, v7, v8, v9,
     return  min((
         min((v0, v1, v2, v3, v4, v5, v6, v7, v8)),
         min((v9, v10, v11, v12, v13, v14, v15, v16, v17)),
-        max((v18, v19, v20, v21, v22, v23, v24, v25, v26))
+        min((v18, v19, v20, v21, v22, v23, v24, v25, v26))
         ))
 
 @numba.njit(inline="always")
@@ -64,7 +64,8 @@ def choose_op27(opname,
                 v10, v11, v12, v13, v14, v15, v16, v17,v18,
                 v19, v20, v21, v22, v23, v24, v25, v26,
                 onlyzero=False):
-    if onlyzero and v0 >0: return v0
+    # v13 is the centre voxel: the 27 values arrive in z-1/z/z+1 slice order.
+    if onlyzero and v13 >0: return v13
     match opname:
         case "min": 
             return min27(v0, v1, v2, v3, v4, v5, v6, v7, v8,

{nbmorph-0.3.0 → nbmorph-0.3.1/src/nbmorph.egg-info}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: nbmorph
-Version: 0.3.0
+Version: 0.3.1
 Summary: A small package with Numba-accelerated morphological operations.
 Author-email: Marius Causemann <mariusca@simula.no>
 License: MIT License

{nbmorph-0.3.0 → nbmorph-0.3.1}/tests/test_minmax.py

@@ -24,3 +24,21 @@ def test_minimum_diamond():
     labels[1, 0, 1] = 1
     result = nbmorph.minimum_diamond(labels)
     assert result[1, 1, 1] == 1
+
+
+def test_minimum_box_footprint_is_full_3x3x3():
+    """A single zero voxel must spread to all 27 voxels of the box footprint."""
+    labels = np.ones((5, 5, 5), dtype=np.uint8)
+    labels[2, 2, 2] = 0
+    result = nbmorph.minimum_box(labels)
+    assert (result == 0).sum() == 27
+    assert (result[1:4, 1:4, 1:4] == 0).all()
+
+
+def test_maximum_box_footprint_is_full_3x3x3():
+    """A single high voxel must spread to all 27 voxels of the box footprint."""
+    labels = np.ones((5, 5, 5), dtype=np.uint8)
+    labels[2, 2, 2] = 5
+    result = nbmorph.maximum_box(labels)
+    assert (result == 5).sum() == 27
+    assert (result[1:4, 1:4, 1:4] == 5).all()