mitkToFTestingCommon.h

Go to the documentation of this file.

/*============================================================================
 
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.
 
============================================================================*/
#ifndef mitkToFTestingCommon_h
#define mitkToFTestingCommon_h
 
#include <MitkToFProcessingExports.h>
#include "mitkNumericTypes.h"
#include <mitkToFProcessingCommon.h>
#include <mitkSurface.h>
#include <mitkPointSet.h>
#include <itksys/SystemTools.hxx>
 
#include <vtkSmartPointer.h>
#include <vtkPolyData.h>
 
namespace mitk
{
class MITKTOFPROCESSING_EXPORT ToFTestingCommon
{
public:
 
    static bool PointSetsEqual(mitk::PointSet::Pointer pointSet1, mitk::PointSet::Pointer pointSet2)
    {
        bool pointSetsEqual = true;
        if (pointSet1->GetSize()==pointSet2->GetSize())
        {
            for (int i=0; i<pointSet1->GetSize(); i++)
            {
                mitk::Point3D expectedPoint = pointSet1->GetPoint(i);
                mitk::Point3D resultPoint = pointSet2->GetPoint(i);
                if (!mitk::Equal(expectedPoint,resultPoint))
                {
                    std::cout << std::endl;
                    std::cout << std::setprecision(12) << "expected: " << expectedPoint;
                    std::cout << std::endl;
                    std::cout << std::setprecision(12) << "resultPoint: " << resultPoint;
                    std::cout << std::endl;
                    pointSetsEqual = false;
                }
            }
        }
        else
        {
            pointSetsEqual = false;
            MITK_INFO<<"Point sets have different size: "<<pointSet1->GetSize()<<" vs. "<<pointSet2->GetSize();
        }
        return pointSetsEqual;
    }
 
    static bool VtkPolyDatasEqual( vtkSmartPointer<vtkPolyData> poly1, vtkSmartPointer<vtkPolyData> poly2  )
    {
        return PointSetsEqual(VtkPolyDataToMitkPointSet(poly1), VtkPolyDataToMitkPointSet(poly2));
    }
 
    static mitk::PointSet::Pointer VtkPolyDataToMitkPointSet( vtkSmartPointer<vtkPolyData> poly )
    {
        mitk::PointSet::Pointer result = mitk::PointSet::New();
        int numberOfPoints = poly->GetNumberOfPoints();
        for (int i=0; i<numberOfPoints; i++)
        {
          double* currentPoint = poly->GetPoint(i);
          mitk::Point3D point;
          point[0] = currentPoint[0];
          point[1] = currentPoint[1];
          point[2] = currentPoint[2];
          result->InsertPoint(i,point);
        }
        return result;
    }
};
 
}
#endif