def test_Basic(self):
"""Test the case of a region with only one intensity value and some basic return values."""
# Create label image with only one region
label_image = scipy.zeros(2 * 2 * 2, dtype=scipy.int8).reshape(2, 2, 2)
# Create original image with only one intensity value
original_image = scipy.zeros(2 * 2 * 2,
dtype=scipy.int8).reshape(2, 2, 2)
# Initialize object
statistics = LabelImageStatistics(label_image, original_image)
# Check created intensity distribution
intensity_distributions = statistics.get_intensity_distributions()
self.assertEqual(len(intensity_distributions), 1)
self.assertEqual(intensity_distributions[0], 0)
intensity_distribution_histogram = statistics.get_intensity_distribution_histogram(
)
self.assertEqual(
len(intensity_distribution_histogram[0]),
statistics.get_intensity_distribution_histogram_width())
self.assertEqual(
len(intensity_distribution_histogram[1]),
statistics.get_intensity_distribution_histogram_width() + 1)
self.assertEqual(
intensity_distribution_histogram[0]
[statistics.get_intensity_distribution_histogram_width() / 2], 1)
self.assertEqual(intensity_distribution_histogram[0].max(), 1)
self.assertEqual(intensity_distribution_histogram[0].min(), 0)
# Check created size distribution
sizes = statistics.get_sizes()
self.assertEqual(len(sizes), 1)
self.assertEqual(sizes[0], 2 * 2 * 2)
sizes_histogram = statistics.get_size_histogram()
self.assertEqual(len(sizes_histogram[0]),
statistics.get_size_histogram_width())
self.assertEqual(len(sizes_histogram[1]),
statistics.get_size_histogram_width() + 1)
self.assertEqual(
sizes_histogram[0][statistics.get_size_histogram_width() / 2], 1)
self.assertEqual(sizes_histogram[0].max(), 1)
self.assertEqual(sizes_histogram[0].min(), 0)
# Check other statistics values
self.assertTrue(statistics.labels_are_consecutive())
self.assertEqual(statistics.get_min_intensity(), 0)
self.assertEqual(statistics.get_max_intensity(), 0)
self.assertEqual(statistics.get_min_label(), 0)
self.assertEqual(statistics.get_max_label(), 0)
self.assertEqual(statistics.get_label_count(), 1)
def test_Basic(self):
"""Test the case of a region with only one intensity value and some basic return values."""
# Create label image with only one region
label_image = scipy.zeros(2*2*2, dtype=scipy.int8).reshape(2,2,2)
# Create original image with only one intensity value
original_image = scipy.zeros(2*2*2, dtype=scipy.int8).reshape(2,2,2)
# Initialize object
statistics = LabelImageStatistics(label_image, original_image)
# Check created intensity distribution
intensity_distributions = statistics.get_intensity_distributions()
self.assertEqual(len(intensity_distributions), 1)
self.assertEqual(intensity_distributions[0], 0)
intensity_distribution_histogram = statistics.get_intensity_distribution_histogram()
self.assertEqual(len(intensity_distribution_histogram[0]), statistics.get_intensity_distribution_histogram_width())
self.assertEqual(len(intensity_distribution_histogram[1]), statistics.get_intensity_distribution_histogram_width() + 1)
self.assertEqual(intensity_distribution_histogram[0][statistics.get_intensity_distribution_histogram_width()/2], 1)
self.assertEqual(intensity_distribution_histogram[0].max(), 1)
self.assertEqual(intensity_distribution_histogram[0].min(), 0)
# Check created size distribution
sizes = statistics.get_sizes()
self.assertEqual(len(sizes), 1)
self.assertEqual(sizes[0], 2*2*2)
sizes_histogram = statistics.get_size_histogram()
self.assertEqual(len(sizes_histogram[0]), statistics.get_size_histogram_width())
self.assertEqual(len(sizes_histogram[1]), statistics.get_size_histogram_width() + 1)
self.assertEqual(sizes_histogram[0][statistics.get_size_histogram_width()/2], 1)
self.assertEqual(sizes_histogram[0].max(), 1)
self.assertEqual(sizes_histogram[0].min(), 0)
# Check other statistics values
self.assertTrue(statistics.labels_are_consecutive())
self.assertEqual(statistics.get_min_intensity(), 0)
self.assertEqual(statistics.get_max_intensity(), 0)
self.assertEqual(statistics.get_min_label(), 0)
self.assertEqual(statistics.get_max_label(), 0)
self.assertEqual(statistics.get_label_count(), 1)
