Method IsSubGeometry

Modules/Core/src/DataManagement/mitkBaseGeometry.cpp:1065–1181

}

bool mitk::IsSubGeometry(const mitk::BaseGeometry& testGeo,
  const mitk::BaseGeometry& referenceGeo,
  ScalarType coordinateEps,
  ScalarType directionEps,
  bool verbose)
{
  bool result = true;

  // Compare spacings (must be equal)
  const auto testedSpacing = testGeo.GetSpacing();
  if (!mitk::Equal(testedSpacing, referenceGeo.GetSpacing(), coordinateEps))
  {
    if (verbose)
    {
      MITK_INFO << "[( Geometry3D )] Spacing differs.";
      MITK_INFO << "testedGeometry is " << setprecision(12) << testedSpacing << " : referenceGeometry is "
        << referenceGeo.GetSpacing() << " and tolerance is " << coordinateEps;
    }
    result = false;
  }

  // Compare ImageGeometry Flag (must be equal)
  if (referenceGeo.GetImageGeometry() != testGeo.GetImageGeometry())
  {
    if (verbose)
    {
      MITK_INFO << "[( Geometry3D )] GetImageGeometry is different.";
      MITK_INFO << "referenceGeo is " << referenceGeo.GetImageGeometry() << " : testGeo is "
        << testGeo.GetImageGeometry();
    }
    result = false;
  }

  // Compare IndexToWorldTransform Matrix (must be equal -> same axis directions)
  if (!Equal(*(testGeo.GetIndexToWorldTransform()), *(referenceGeo.GetIndexToWorldTransform()), directionEps, verbose))
  {
    result = false;
  }

  //check if the geometry is within or equal to the bounds of reference geomentry.
  for (int i = 0; i<8; ++i)
  {
    auto testCorner = testGeo.GetCornerPoint(i);
    bool isInside = false;
    mitk::Point3D testCornerIndex;
    referenceGeo.WorldToIndex(testCorner, testCornerIndex);

    std::bitset<sizeof(int)> bs(i);
    //To regard the coordinateEps, we subtract or add it to the index elements
    //depending on whether it was constructed by a lower or an upper bound value
    //(see implementation of BaseGeometry::GetCorner()).
    if (bs.test(0))
    {
      testCornerIndex[2] -= coordinateEps;
    }
    else
    {
      testCornerIndex[2] += coordinateEps;