diff --git a/Modules/ImageStatistics/Testing/mitkImageStatisticsContainerTest.cpp b/Modules/ImageStatistics/Testing/mitkImageStatisticsContainerTest.cpp
new file mode 100644
index 0000000000..6578f14cb6
--- /dev/null
+++ b/Modules/ImageStatistics/Testing/mitkImageStatisticsContainerTest.cpp
@@ -0,0 +1,298 @@
+/*============================================================================
+
+The Medical Imaging Interaction Toolkit (MITK)
+
+Copyright (c) German Cancer Research Center (DKFZ)
+All rights reserved.
+
+Use of this source code is governed by a 3-clause BSD license that can be
+found in the LICENSE file.
+
+============================================================================*/
+// Testing
+#include "mitkTestingMacros.h"
+#include "mitkTestFixture.h"
+
+//MITK includes
+#include "mitkIOUtil.h"
+#include "mitkStandaloneDataStorage.h"
+#include "mitkStatisticsToImageRelationRule.h"
+#include "mitkArbitraryTimeGeometry.h"
+#include "mitkImageStatisticsContainerNodeHelper.h"
+#include "mitkImageStatisticsContainerManager.h"
+#include "mitkImageStatisticsContainer.h"
+
+class mitkImageStatisticsContainerTestSuite : public mitk::TestFixture
+{
+ CPPUNIT_TEST_SUITE(mitkImageStatisticsContainerTestSuite);
+ MITK_TEST(StatisticNamesIO);
+ MITK_TEST(InvalidKey);
+ MITK_TEST(TimeSteps);
+ MITK_TEST(PrintSelf);
+ MITK_TEST(InternalClone);
+ MITK_TEST(StatisticNames);
+ MITK_TEST(Reset);
+ CPPUNIT_TEST_SUITE_END();
+
+private:
+ mitk::ImageStatisticsContainer::Pointer m_statisticsContainer;
+ mitk::ImageStatisticsContainer::Pointer m_statisticsContainer2;
+ mitk::ImageStatisticsContainer::Pointer m_statisticsContainer3;
+
+ mitk::ImageStatisticsContainer::ImageStatisticsObject m_statisticsObject;
+ mitk::ImageStatisticsContainer::ImageStatisticsObject m_statisticsObject2;
+ mitk::ImageStatisticsContainer::ImageStatisticsObject m_statisticsObject3;
+
+ mitk::ArbitraryTimeGeometry::Pointer m_timeGeometry;
+
+ mitk::Geometry3D::Pointer m_Geometry1;
+ mitk::Geometry3D::Pointer m_Geometry2;
+ mitk::Geometry3D::Pointer m_Geometry3;
+ mitk::Geometry3D::Pointer m_Geometry3_5;
+ mitk::Geometry3D::Pointer m_Geometry4;
+ mitk::Geometry3D::Pointer m_Geometry5;
+
+ mitk::TimePointType m_Geometry1MinTP;
+ mitk::TimePointType m_Geometry2MinTP;
+ mitk::TimePointType m_Geometry3MinTP;
+ mitk::TimePointType m_Geometry3_5MinTP;
+ mitk::TimePointType m_Geometry4MinTP;
+ mitk::TimePointType m_Geometry5MinTP;
+
+ mitk::TimePointType m_Geometry1MaxTP;
+ mitk::TimePointType m_Geometry2MaxTP;
+ mitk::TimePointType m_Geometry3MaxTP;
+ mitk::TimePointType m_Geometry3_5MaxTP;
+ mitk::TimePointType m_Geometry4MaxTP;
+ mitk::TimePointType m_Geometry5MaxTP;
+
+ std::vector<std::string> estimatedDefaultStatisticNames
+ {
+ "Mean",
+ "Median",
+ "StandardDeviation",
+ "RMS",
+ "Max",
+ "MaxPosition",
+ "Min",
+ "MinPosition",
+ "#Voxel",
+ "Volume [mm^3]",
+ "Skewness",
+ "Kurtosis",
+ "Uniformity",
+ "Entropy",
+ "MPP",
+ "UPP"
+ };
+
+public:
+ void setUp() override
+ {
+ m_statisticsContainer = mitk::ImageStatisticsContainer::New();
+ m_statisticsContainer2 = mitk::ImageStatisticsContainer::New();
+ m_statisticsContainer3 = mitk::ImageStatisticsContainer::New();
+
+ m_Geometry1 = mitk::Geometry3D::New();
+ m_Geometry2 = mitk::Geometry3D::New();
+ m_Geometry3 = mitk::Geometry3D::New();
+ m_Geometry3_5 = mitk::Geometry3D::New();
+ m_Geometry4 = mitk::Geometry3D::New();
+ m_Geometry5 = mitk::Geometry3D::New();
+
+ m_Geometry1MinTP = 1;
+ m_Geometry2MinTP = 2;
+ m_Geometry3MinTP = 3;
+ m_Geometry3_5MinTP = 3.5;
+ m_Geometry4MinTP = 4;
+ m_Geometry5MinTP = 5;
+
+ m_Geometry1MaxTP = 1.9;
+ m_Geometry2MaxTP = 2.9;
+ m_Geometry3MaxTP = 3.9;
+ m_Geometry3_5MaxTP = 3.9;
+ m_Geometry4MaxTP = 4.9;
+ m_Geometry5MaxTP = 5.9;
+
+ m_timeGeometry = mitk::ArbitraryTimeGeometry::New();
+ m_timeGeometry->ClearAllGeometries();
+ m_timeGeometry->AppendNewTimeStep(m_Geometry1, m_Geometry1MinTP, m_Geometry1MaxTP);
+ m_timeGeometry->AppendNewTimeStep(m_Geometry2, m_Geometry2MinTP, m_Geometry2MaxTP);
+ m_timeGeometry->AppendNewTimeStep(m_Geometry3, m_Geometry3MinTP, m_Geometry3MaxTP);
+ m_timeGeometry->AppendNewTimeStep(m_Geometry4, m_Geometry4MinTP, m_Geometry4MaxTP);
+ m_timeGeometry->AppendNewTimeStep(m_Geometry5, m_Geometry5MinTP, m_Geometry5MaxTP);
+ }
+
+ void tearDown() override
+ {
+ m_statisticsContainer->Reset();
+ m_statisticsContainer2->Reset();
+ m_statisticsContainer3->Reset();
+
+ m_statisticsObject.Reset();
+ m_statisticsObject2.Reset();
+ m_statisticsObject3.Reset();
+
+ m_timeGeometry->ClearAllGeometries();
+ }
+
+ void CreateNodeRelationImage(mitk::BaseData::Pointer statistics, mitk::BaseData::ConstPointer image)
+ {
+ auto rule = mitk::StatisticsToImageRelationRule::New();
+ rule->Connect(statistics, image);
+ }
+
+ void StatisticNamesIO()
+ {
+ auto defaultStatisticNames = mitk::ImageStatisticsContainer::ImageStatisticsObject::GetDefaultStatisticNames();
+ for (size_t i = 0; i < estimatedDefaultStatisticNames.size(); i++)
+ {
+ CPPUNIT_ASSERT_EQUAL_MESSAGE("Estimated default statistics names are not correct.", defaultStatisticNames.at(i), estimatedDefaultStatisticNames.at(i));
+ }
+
+ std::string testName = "testName";
+ double realTypeTest = 3.14;
+ m_statisticsObject.AddStatistic(testName, realTypeTest);
+ auto customStatisticNames = m_statisticsObject.GetCustomStatisticNames();
+ size_t one = 1;
+ CPPUNIT_ASSERT_EQUAL_MESSAGE("Amount of estimated custom statistics names is not correct.", customStatisticNames.size(), one);
+ CPPUNIT_ASSERT_EQUAL_MESSAGE("Estimated custom statistics names are not correct.", customStatisticNames.front(), testName);
+
+ auto allStatisticNames = m_statisticsObject.GetAllStatisticNames();
+ estimatedDefaultStatisticNames.push_back(testName);
+ for (size_t i = 0; i < estimatedDefaultStatisticNames.size(); i++)
+ {
+ CPPUNIT_ASSERT_EQUAL_MESSAGE("Estimated statistics names are not correct.", allStatisticNames.at(i), estimatedDefaultStatisticNames.at(i));
+ }
+
+ auto existingStatisticNames = m_statisticsObject.GetExistingStatisticNames();
+ CPPUNIT_ASSERT_EQUAL_MESSAGE("Amount of estimated existing statistics names is not correct.", existingStatisticNames.size(), one);
+ CPPUNIT_ASSERT_EQUAL_MESSAGE("Estimated existing statistics names are not correct.", existingStatisticNames.front(), testName);
+
+ m_statisticsObject.Reset();
+ CPPUNIT_ASSERT_MESSAGE("Custom statistics names were not removed correctly.", m_statisticsObject.GetCustomStatisticNames().empty());
+ }
+
+ void InvalidKey()
+ {
+ CPPUNIT_ASSERT_THROW_MESSAGE("Exception should have been thrown because key does not exists.", m_statisticsObject.GetValueNonConverted("Test"), mitk::Exception);
+ }
+
+ void TimeSteps()
+ {
+ m_statisticsContainer->SetTimeGeometry(m_timeGeometry);
+
+ m_statisticsContainer->SetStatisticsForTimeStep(0, m_statisticsObject);
+ CPPUNIT_ASSERT_EQUAL_MESSAGE("Previous added time step was not saved correctly.", m_statisticsContainer->TimeStepExists(0), true);
+ m_statisticsContainer->SetStatisticsForTimeStep(1, m_statisticsObject);
+ CPPUNIT_ASSERT_EQUAL_MESSAGE("Previous added time step was not saved correctly.", m_statisticsContainer->TimeStepExists(1), true);
+ m_statisticsContainer->SetStatisticsForTimeStep(2, m_statisticsObject);
+ CPPUNIT_ASSERT_EQUAL_MESSAGE("Previous added time step was not saved correctly.", m_statisticsContainer->TimeStepExists(2), true);
+ m_statisticsContainer->SetStatisticsForTimeStep(3, m_statisticsObject);
+ CPPUNIT_ASSERT_EQUAL_MESSAGE("Previous added time step was not saved correctly.", m_statisticsContainer->TimeStepExists(3), true);
+ m_statisticsContainer->SetStatisticsForTimeStep(4, m_statisticsObject);
+ CPPUNIT_ASSERT_EQUAL_MESSAGE("Previous added time step was not saved correctly.", m_statisticsContainer->TimeStepExists(4), true);
+
+ CPPUNIT_ASSERT_THROW_MESSAGE("Out of timeStep geometry bounds timeStep was added but no exception was thrown.", m_statisticsContainer->SetStatisticsForTimeStep(42, m_statisticsObject), mitk::Exception);
+
+ CPPUNIT_ASSERT_EQUAL_MESSAGE("Statistics container does not contain the right amount of timeSteps.", static_cast<int> (m_statisticsContainer->GetNumberOfTimeSteps()), 5);
+ }
+
+ void PrintSelf()
+ {
+ m_statisticsContainer->SetTimeGeometry(m_timeGeometry);
+
+ m_statisticsObject.AddStatistic("Test", 4.2);
+
+ m_statisticsContainer->SetStatisticsForTimeStep(0, m_statisticsObject);
+ m_statisticsContainer->SetStatisticsForTimeStep(1, m_statisticsObject);
+ m_statisticsContainer->SetStatisticsForTimeStep(2, m_statisticsObject);
+ m_statisticsContainer->SetStatisticsForTimeStep(3, m_statisticsObject);
+ m_statisticsContainer->SetStatisticsForTimeStep(4, m_statisticsObject);
+
+ std::stringstream print;
+ m_statisticsContainer->Print(print);
+
+ CPPUNIT_ASSERT_EQUAL_MESSAGE("Data was changed.", static_cast<int> (print.str().length()), 1237);
+ }
+
+ void InternalClone()
+ {
+ auto clone = m_statisticsContainer->Clone();
+
+ CPPUNIT_ASSERT_EQUAL_MESSAGE("Internal clone was not cloned correctly.", m_statisticsContainer->GetGeometry(), clone->GetGeometry());
+ CPPUNIT_ASSERT_EQUAL_MESSAGE("Internal clone was not cloned correctly.", m_statisticsContainer->GetNumberOfTimeSteps(), clone->GetNumberOfTimeSteps());
+ }
+
+ void StatisticNames()
+ {
+ m_statisticsContainer->SetTimeGeometry(m_timeGeometry);
+
+ m_statisticsObject.AddStatistic("Test", 4.2);
+ m_statisticsObject2.AddStatistic("Test2", 4.2);
+ m_statisticsObject2.AddStatistic("Test3", 4.2);
+
+ m_statisticsContainer->SetStatisticsForTimeStep(0, m_statisticsObject);
+ m_statisticsContainer->SetStatisticsForTimeStep(1, m_statisticsObject);
+ m_statisticsContainer->SetStatisticsForTimeStep(2, m_statisticsObject);
+ m_statisticsContainer->SetStatisticsForTimeStep(3, m_statisticsObject);
+ m_statisticsContainer->SetStatisticsForTimeStep(4, m_statisticsObject);
+
+ m_statisticsContainer->SetStatisticsForTimeStep(0, m_statisticsObject2);
+ m_statisticsContainer->SetStatisticsForTimeStep(1, m_statisticsObject2);
+ m_statisticsContainer->SetStatisticsForTimeStep(2, m_statisticsObject2);
+ m_statisticsContainer->SetStatisticsForTimeStep(3, m_statisticsObject2);
+ m_statisticsContainer->SetStatisticsForTimeStep(4, m_statisticsObject2);
+
+ m_statisticsContainer->SetStatisticsForTimeStep(0, m_statisticsObject3);
+ m_statisticsContainer->SetStatisticsForTimeStep(1, m_statisticsObject3);
+ m_statisticsContainer->SetStatisticsForTimeStep(2, m_statisticsObject3);
+ m_statisticsContainer->SetStatisticsForTimeStep(3, m_statisticsObject3);
+ m_statisticsContainer->SetStatisticsForTimeStep(4, m_statisticsObject3);
+
+ CPPUNIT_ASSERT_EQUAL_MESSAGE("Estimated statistics names are not correct.", mitk::GetAllStatisticNames(m_statisticsContainer).size(), estimatedDefaultStatisticNames.size() + 1);
+
+ std::vector<mitk::ImageStatisticsContainer::ConstPointer> containers;
+
+ containers.push_back(m_statisticsContainer.GetPointer());
+ containers.push_back(m_statisticsContainer2.GetPointer());
+ containers.push_back(m_statisticsContainer3.GetPointer());
+
+ CPPUNIT_ASSERT_EQUAL_MESSAGE("Estimated statistics names are not correct.", mitk::GetAllStatisticNames(containers).size(), estimatedDefaultStatisticNames.size() + 1);
+ }
+
+ void Reset()
+ {
+ m_statisticsContainer->SetTimeGeometry(m_timeGeometry);
+
+ m_statisticsObject.AddStatistic("Test", 4.2);
+
+ m_statisticsContainer->SetStatisticsForTimeStep(0, m_statisticsObject);
+ m_statisticsContainer->SetStatisticsForTimeStep(1, m_statisticsObject);
+ m_statisticsContainer->SetStatisticsForTimeStep(2, m_statisticsObject);
+ m_statisticsContainer->SetStatisticsForTimeStep(3, m_statisticsObject);
+ m_statisticsContainer->SetStatisticsForTimeStep(4, m_statisticsObject);
+
+ m_statisticsContainer->Reset();
+
+ m_statisticsObject.Reset();
+
+ m_timeGeometry->ClearAllGeometries();
+
+ CPPUNIT_ASSERT_EQUAL_MESSAGE("Statistic object for time step was not deleted correctly.", m_statisticsContainer->GetStatisticsForTimeStep(0).GetExistingStatisticNames().size(), size_t());
+ CPPUNIT_ASSERT_EQUAL_MESSAGE("Statistic object for time step was not deleted correctly.", m_statisticsContainer->GetStatisticsForTimeStep(0).GetCustomStatisticNames().size(), size_t());
+ CPPUNIT_ASSERT_EQUAL_MESSAGE("Statistic object for time step was not deleted correctly.", m_statisticsContainer->GetStatisticsForTimeStep(1).GetExistingStatisticNames().size(), size_t());
+ CPPUNIT_ASSERT_EQUAL_MESSAGE("Statistic object for time step was not deleted correctly.", m_statisticsContainer->GetStatisticsForTimeStep(1).GetCustomStatisticNames().size(), size_t());
+ CPPUNIT_ASSERT_EQUAL_MESSAGE("Statistic object for time step was not deleted correctly.", m_statisticsContainer->GetStatisticsForTimeStep(2).GetExistingStatisticNames().size(), size_t());
+ CPPUNIT_ASSERT_EQUAL_MESSAGE("Statistic object for time step was not deleted correctly.", m_statisticsContainer->GetStatisticsForTimeStep(2).GetCustomStatisticNames().size(), size_t());
+ CPPUNIT_ASSERT_EQUAL_MESSAGE("Statistic object for time step was not deleted correctly.", m_statisticsContainer->GetStatisticsForTimeStep(3).GetExistingStatisticNames().size(), size_t());
+ CPPUNIT_ASSERT_EQUAL_MESSAGE("Statistic object for time step was not deleted correctly.", m_statisticsContainer->GetStatisticsForTimeStep(3).GetCustomStatisticNames().size(), size_t());
+ CPPUNIT_ASSERT_EQUAL_MESSAGE("Statistic object for time step was not deleted correctly.", m_statisticsContainer->GetStatisticsForTimeStep(4).GetExistingStatisticNames().size(), size_t());
+ CPPUNIT_ASSERT_EQUAL_MESSAGE("Statistic object for time step was not deleted correctly.", m_statisticsContainer->GetStatisticsForTimeStep(4).GetCustomStatisticNames().size(), size_t());
+
+ CPPUNIT_ASSERT_EQUAL_MESSAGE("Statistic was not deleted correctly.", m_statisticsObject.HasStatistic("Test"), false);
+
+ CPPUNIT_ASSERT_EQUAL_MESSAGE("Time geometry was not cleared correctly.", m_timeGeometry->IsValid(), false);
+ }
+};
+MITK_TEST_SUITE_REGISTRATION(mitkImageStatisticsContainer)