Method TestCase6

Modules/ImageStatistics/Testing/mitkImageStatisticsCalculatorTest.cpp:308–336

}

void mitkImageStatisticsCalculatorTestSuite::TestCase6()
{
  /*****************************
   * quarter pixel (black) + whole pixel (white) + half pixel (gray) in x-direction
   * -> mean of 191.5 expected
   ******************************/
  MITK_INFO << std::endl << "Test case 6:-----------------------------------------------------------------------------------";

  std::string filename = this->GetTestDataFilePath("ImageStatisticsTestData/testimage.nrrd");
  m_TestImage = mitk::IOUtil::Load<mitk::Image>(filename);
  m_Geometry = m_TestImage->GetSlicedGeometry()->GetPlaneGeometry(0);

  mitk::Point2D pnt1; pnt1[0] = 11.0; pnt1[1] = 3.5;
  mitk::Point2D pnt2; pnt2[0] = 9.25; pnt2[1] = 3.5;
  mitk::Point2D pnt3; pnt3[0] = 9.25; pnt3[1] = 4.5;
  mitk::Point2D pnt4; pnt4[0] = 11.0; pnt4[1] = 4.5;
  std::vector<mitk::Point2D> points{ pnt1,pnt2,pnt3,pnt4 };
        auto figure = GeneratePlanarPolygon(m_Geometry, points);

  mitk::PlanarFigureMaskGenerator::Pointer planFigMaskGen = mitk::PlanarFigureMaskGenerator::New();
  planFigMaskGen->SetInputImage(m_TestImage);
        planFigMaskGen->SetPlanarFigure(figure.GetPointer());

  mitk::ImageStatisticsContainer::Pointer statisticsContainer;
  CPPUNIT_ASSERT_NO_THROW(statisticsContainer = ComputeStatistics(m_TestImage, planFigMaskGen.GetPointer()));
  auto statisticsObjectTimestep0 = statisticsContainer->GetStatistics(1, 0);

  this->VerifyStatistics(statisticsObjectTimestep0, 191.5, 89.802561210691536, 128.63499999046327);
}

void mitkImageStatisticsCalculatorTestSuite::TestCase7()
{