Function testItkImageIsCenterBased

Modules/Core/test/mitkGeometry3DTest.cpp:239–291

#include <itkImage.h>

int testItkImageIsCenterBased()
{
  MITK_TEST_OUTPUT(<< "Testing whether itk::Image coordinates are center-based.");
  typedef itk::Image<int, 3> ItkIntImage3D;
  ItkIntImage3D::Pointer itkintimage = ItkIntImage3D::New();
  ItkIntImage3D::SizeType size;
  size.Fill(10);
  mitk::Point3D origin;
  mitk::FillVector3D(origin, 2, 3, 7);
  itkintimage->Initialize();
  itkintimage->SetRegions(size);
  itkintimage->SetOrigin(origin);
  std::cout << "[PASSED]" << std::endl;

  MITK_TEST_OUTPUT(<< " Testing itk::Image::TransformPhysicalPointToContinuousIndex(origin)==(0,0,0)");
  mitk::Point3D globalOrigin;
  mitk::FillVector3D(globalOrigin, 0, 0, 0);

  auto originContinuousIndex = itkintimage->TransformPhysicalPointToContinuousIndex<mitk::ScalarType, mitk::ScalarType>(origin);
  MITK_TEST_CONDITION_REQUIRED(originContinuousIndex == globalOrigin, "");

  MITK_TEST_OUTPUT(<< " Testing itk::Image::TransformPhysicalPointToIndex(origin)==(0,0,0)");
  itk::Index<3> itkindex = itkintimage->TransformPhysicalPointToIndex(origin);
  itk::Index<3> globalOriginIndex;
  for (unsigned int i = 0; i < 3; ++i)
    globalOriginIndex[i] = static_cast<itk::IndexValueType>(globalOrigin[i]);
  MITK_TEST_CONDITION_REQUIRED(itkindex == globalOriginIndex, "");

  MITK_TEST_OUTPUT(<< " Testing itk::Image::TransformPhysicalPointToIndex(origin-0.5*spacing)==(0,0,0)");
  mitk::Vector3D halfSpacingStep = itkintimage->GetSpacing() * 0.5;
  mitk::Matrix3D rotation;
  mitk::Point3D originOffCenter = origin - halfSpacingStep;
  itkindex = itkintimage->TransformPhysicalPointToIndex(originOffCenter);
  MITK_TEST_CONDITION_REQUIRED(itkindex == globalOriginIndex, "");

  MITK_TEST_OUTPUT(
    << " Testing itk::Image::TransformPhysicalPointToIndex(origin+0.5*spacing-eps, itk::Index)==(0,0,0)");
  originOffCenter = origin + halfSpacingStep;
  originOffCenter -= mitk::Vector(0.0001);
  itkindex = itkintimage->TransformPhysicalPointToIndex(originOffCenter);
  MITK_TEST_CONDITION_REQUIRED(itkindex == globalOriginIndex, "");

  MITK_TEST_OUTPUT(<< " Testing itk::Image::TransformPhysicalPointToIndex(origin+0.5*spacing, itk::Index)==(1,1,1)");
  originOffCenter = origin + halfSpacingStep;
  itk::Index<3> global111;
  mitk::FillVector3D(global111, 1, 1, 1);
  itkindex = itkintimage->TransformPhysicalPointToIndex(originOffCenter);
  MITK_TEST_CONDITION_REQUIRED(itkindex == global111, "");

  MITK_TEST_OUTPUT(<< "=> Yes, itk::Image coordinates are center-based.");

  return EXIT_SUCCESS;
}

int testGeometry3D(bool imageGeometry)
{