ENH: Add a Python test

Author: dzenanz

Modules/Filtering/ImageStatistics/test/Input/DzZ_Seeds.seg.nrrd.sha512

@@ -0,0 +1 @@
+4106f7a97659761a7fd42594a4f19aa9dad312445a1c9ff63f397f0c32cb952ee5d6ba6c6582951d883e6f40587091437a872c5fce966670d57f4984125c7a5d

Modules/Filtering/ImageStatistics/test/Input/DzZ_T1.seg.nrrd.sha512

@@ -0,0 +1 @@
+eb93426d1ee5f00d722979fd9dbfd8701f27bbb4256522e3d25c4db0f88f33b6a2db1ce572613e38c1a91787ca911e2399d206a137f7e987421b357bd5925b5d

Modules/Filtering/ImageStatistics/wrapping/test/CMakeLists.txt

@@ -0,0 +1,12 @@
+set(test_input_dir ${itk-module_SOURCE_DIR}/test/Input)
+
+# let's make sure 3D uchar images are wrapped
+list(FIND ITK_WRAP_IMAGE_DIMS 3 wrap_3_index)
+if(ITK_WRAP_PYTHON AND ITK_WRAP_unsigned_char AND wrap_3_index GREATER -1)
+  itk_python_add_test(NAME OrientImageFilterTest
+    TEST_DRIVER_ARGS
+    COMMAND ${CMAKE_CURRENT_SOURCE_DIR}/itkLabelOverlapMeasuresImageFilterTest.py
+      DATA{${test_input_dir}/DzZ_T1.seg.nrrd}
+      DATA{${test_input_dir}/DzZ_Seeds.seg.nrrd}
+    )
+endif()

Modules/Filtering/ImageStatistics/wrapping/test/itkLabelOverlapMeasuresImageFilterTest.py

@@ -0,0 +1,34 @@
+# ==========================================================================
+#
+#   Copyright NumFOCUS
+#
+#   Licensed under the Apache License, Version 2.0 (the "License");
+#   you may not use this file except in compliance with the License.
+#   You may obtain a copy of the License at
+#
+#          https://www.apache.org/licenses/LICENSE-2.0.txt
+#
+#   Unless required by applicable law or agreed to in writing, software
+#   distributed under the License is distributed on an "AS IS" BASIS,
+#   WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
+#   See the License for the specific language governing permissions and
+#   limitations under the License.
+#
+# ==========================================================================*/
+
+
+import itk
+from sys import argv
+
+itk.auto_progress(2)
+
+ref = itk.imread(argv[1], itk.UC)
+seg = itk.imread(argv[2], itk.UC)
+
+lom_filter = itk.LabelOverlapMeasuresImageFilter[itk.Image[itk.UC, 3]].New()
+lom_filter.SetTargetImage(seg)
+lom_filter.SetSourceImage(ref)
+lom_filter.UpdateLargestPossibleRegion()
+lsm = lom_filter.GetLabelSetMeasures()
+for label, measure in lsm.items():
+    print(f"Label: {label}, i: {measure.m_Intersection}, u: {measure.m_Union}")