Function TEST

Testing/Unit/sitkLabelStatisticsTest.cxx:25–67

#include <memory>

TEST(LabelStatistics, Simple)
{
  // By using the same image, the label min/max values should equal the label itself.
  itk::simple::Image intensityImage = itk::simple::ReadImage(dataFinder.GetFile("Input/2th_cthead1.png"));
  itk::simple::Image labelImage = itk::simple::ReadImage(dataFinder.GetFile("Input/2th_cthead1.png"));

  itk::simple::LabelStatisticsImageFilter lsFilter;

  EXPECT_TRUE(lsFilter.GetUseHistograms());
  lsFilter.UseHistogramsOff();
  EXPECT_FALSE(lsFilter.GetUseHistograms());
  lsFilter.UseHistogramsOn();
  EXPECT_TRUE(lsFilter.GetUseHistograms());
  lsFilter.SetUseHistograms(false);
  EXPECT_FALSE(lsFilter.GetUseHistograms());
  lsFilter.SetUseHistograms(true);
  EXPECT_TRUE(lsFilter.GetUseHistograms());

  try
  {
    lsFilter.Execute(intensityImage, labelImage);
  }
  catch (itk::ExceptionObject & e)
  {
    std::cout << "LabelStatistics failed: " << e.what() << std::endl;
  }

  std::vector<int64_t> labels = lsFilter.GetLabels();
  for (std::vector<int64_t>::const_iterator i = labels.begin(); i != labels.end(); ++i)
  {
    // By using the same image, the label min/max/mean values should equal the label itself.
    ASSERT_EQ(lsFilter.GetMinimum(*i), *i);
    ASSERT_EQ(lsFilter.GetMaximum(*i), *i);
    ASSERT_EQ(lsFilter.GetMean(*i), *i);
    ASSERT_EQ(lsFilter.GetMedian(*i), *i);
    // By using the same image, the label variance values should equal to Zero.
    ASSERT_EQ(lsFilter.GetSigma(*i), 0.0);
    ASSERT_EQ(lsFilter.GetVariance(*i), 0.0);
  }

  ASSERT_EQ(lsFilter.GetSum(0), 0);
  ASSERT_EQ(lsFilter.GetCount(0), 33390u);
}


TEST(LabelStatistics, Commands)