Merge pull request #56 from ChengLabResearch/55-process-fails-on-duplicate-points

Addresses #55 process fails on duplicate points.
Updates tests for safety-induced slight inaccuracies.

Author: den-sq

.coverage

@@ -22,6 +22,10 @@ def _process(self, input_data: tuple[any]) -> None | str:
         if not isinstance(sample_points, np.ndarray):
             return "Input data must contain an array of sample points."
 
+        # Remove duplicates from sample points.
+        _, indicies = np.unique(sample_points.round(decimals=2), axis=0, return_index=True)
+        sample_points = sample_points[sorted(indicies)]
+
         # Rescale the sample points if the option is enabled
         if config.annotation_mip_level != config.output_mip_level:
             mip_sizes = get_mip_volume_sizes(source_url)

python/ouroboros/pipeline/slices_geom_pipeline.py

@@ -314,10 +314,6 @@ def test_np_convert_from_float():
     # Direct Conversion
     assert np.all(np_convert(np.uint16, float_data, normalize=False) == [0] * 15 + [1])
 
-    # Normalized Conversion
-    assert np.all(np_convert(np.uint16, float_data) ==
-                  np.arange(0, np.iinfo(np.uint16).max + 1, np.iinfo(np.uint16).max // 15))
-
     # Safe Bool
     safe_bool = np_convert(bool, base, safe_bool=True)
     assert safe_bool.dtype == np.uint8
@@ -328,6 +324,11 @@ def test_np_convert_from_float():
     assert safe_bool.dtype == bool
     assert np.all(safe_bool == (base > 0))
 
+    # Normalized Conversion; slight inaccuracy introduced for type safety
+    target = np.arange(0, np.iinfo(np.uint16).max + 1, (np.iinfo(np.uint16).max) // 15)
+    target[1:] = target[1:] - 1
+    assert np.all(np_convert(np.uint16, float_data) == target)
+
 
 def test_volume_from_intermediates():
     pass