Method TestCase9

Modules/ImageStatistics/Testing/mitkImageStatisticsCalculatorTest.cpp:398–426

}

void mitkImageStatisticsCalculatorTestSuite::TestCase9()
{
  /*****************************
   * whole pixel (gray) + half pixel (white) in y-direction
   * -> mean of 191.5 expected
   ******************************/
  MITK_INFO << std::endl << "Test case 9:-----------------------------------------------------------------------------------";

  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.5; pnt1[1] = 10.5;
  mitk::Point2D pnt2; pnt2[0] = 11.5; pnt2[1] = 12.0;
  mitk::Point2D pnt3; pnt3[0] = 12.5; pnt3[1] = 12.0;
  mitk::Point2D pnt4; pnt4[0] = 12.5; pnt4[1] = 10.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::TestCase10()
{