mitkNavigationDataTransformFilter.cpp

Go to the documentation of this file.

/*===================================================================

The Medical Imaging Interaction Toolkit (MITK)

Copyright (c) German Cancer Research Center,
Division of Medical and Biological Informatics.
All rights reserved.

This software is distributed WITHOUT ANY WARRANTY; without
even the implied warranty of MERCHANTABILITY or FITNESS FOR
A PARTICULAR PURPOSE.

See LICENSE.txt or http://www.mitk.org for details.

===================================================================*/

#include "mitkNavigationDataTransformFilter.h"


mitk::NavigationDataTransformFilter::NavigationDataTransformFilter()
: mitk::NavigationDataToNavigationDataFilter()
{
  m_Rigid3DTransform = NULL;
  m_Precompose = false;
}


mitk::NavigationDataTransformFilter::~NavigationDataTransformFilter()
{
  m_Rigid3DTransform = NULL;
}

void mitk::NavigationDataTransformFilter::GenerateData()
{

  // only update data if m_Transform was set
  if(m_Rigid3DTransform.IsNull())
  {
    itkExceptionMacro("Invalid parameter: Transform was not set!  Use SetRigid3DTransform() before updating the filter.");
    return;
  }
  else
  {
    this->CreateOutputsForAllInputs(); // make sure that we have the same number of outputs as inputs

    /* update outputs with tracking data from tools */
    for (unsigned int i = 0; i < this->GetNumberOfIndexedOutputs() ; ++i)
    {
      mitk::NavigationData* output = this->GetOutput(i);
      assert(output);
      const mitk::NavigationData* input = this->GetInput(i);
      assert(input);

      if (input->IsDataValid() == false)
      {
        output->SetDataValid(false);
        continue;
      }

      // Get input position/orientation as float numbers
      const NavigationData::PositionType    pInF = input->GetPosition();
      const NavigationData::OrientationType oInF = input->GetOrientation();

      // Cast the input NavigationData to double precision
      TransformType::OutputVectorType pInD;
      FillVector3D(pInD, pInF[0], pInF[1], pInF[2]);
      TransformType::VersorType oInD;
      oInD.Set(oInF.x(), oInF.y(), oInF.z(), oInF.r());

      TransformType::Pointer composedTransform = TransformType::New();
      // SetRotation+SetOffset defines the Tip-to-World coordinate frame
      // transformation ("World" is used in the generic sense)
      composedTransform->SetRotation(oInD);
      composedTransform->SetOffset(pInD);
      // If !m_Precompose: The resulting transform is Tip-to-UserWorld
      // If m_Precompose:  The resulting transform is UserTip-to-World
      composedTransform->Compose(m_Rigid3DTransform, m_Precompose);

      // Transformed position/orientation as double numbers
      const TransformType::OutputVectorType  pOutD = composedTransform->GetOffset();
      const TransformType::VersorType        oOutD = composedTransform->GetVersor();

      // Cast to transformed NavigationData back to float precision
      NavigationData::OrientationType oOutF(oOutD.GetX(), oOutD.GetY(), oOutD.GetZ(), oOutD.GetW());
      NavigationData::PositionType pOutF;
      FillVector3D(pOutF, pOutD[0], pOutD[1], pOutD[2]);

      output->SetOrientation(oOutF);
      output->SetPosition(pOutF);
      output->SetDataValid(true); // operation was successful, therefore data of output is valid.
    }
  }
}