mitkArray.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 mitkArray_h
#define mitkArray_h
 
#include <itkFixedArray.h>
 
#include "mitkEqual.h"
#include "mitkNumericConstants.h"
 
namespace mitk
{
  template <typename ArrayType, typename TCoordRep, unsigned int NVectorDimension>
  void FillArray(itk::FixedArray<TCoordRep, NVectorDimension> &toArray, const ArrayType &array)
  {
    for (unsigned short int var = 0; var < NVectorDimension; ++var)
    {
      toArray[var] = array[var];
    }
  }
 
  template <typename ArrayType, typename TCoordRep, unsigned int NVectorDimension>
  itk::FixedArray<TCoordRep, NVectorDimension> FillArray(const ArrayType &array)
  {
    itk::FixedArray<TCoordRep, NVectorDimension> vectorOrPoint;
 
    mitk::FillArray(vectorOrPoint, array);
 
    return vectorOrPoint;
  }
 
  template <typename ArrayType, typename TCoordRep, unsigned int NVectorDimension>
  void ToArray(ArrayType &array, const itk::FixedArray<TCoordRep, NVectorDimension> &vectorOrPoint)
  {
    for (unsigned short int var = 0; var < NVectorDimension; ++var)
    {
      array[var] = vectorOrPoint[var];
    }
  }
 
  template <typename ArrayType, typename TCoordRep, unsigned int NVectorDimension>
  ArrayType ToArray(const itk::FixedArray<TCoordRep, NVectorDimension> &vectorOrPoint)
  {
    ArrayType result;
 
    mitk::ToArray(result, vectorOrPoint);
 
    return result;
  }
 
  // The FillVector3D and FillVector4D methods are implemented for all common array types here
 
  template <class Tout>
  inline void FillVector3D(Tout &out, mitk::ScalarType x, mitk::ScalarType y, mitk::ScalarType z)
  {
    out[0] = x;
    out[1] = y;
    out[2] = z;
  }
 
  template <class Tout>
  inline void FillVector4D(Tout &out, mitk::ScalarType x, mitk::ScalarType y, mitk::ScalarType z, mitk::ScalarType t)
  {
    out[0] = x;
    out[1] = y;
    out[2] = z;
    out[3] = t;
  }
 
  template <typename TArrayType1, typename TArrayType2>
  inline bool EqualArray(
    TArrayType1 &arrayType1, TArrayType2 &arrayType2, int size, ScalarType eps = mitk::eps, bool verbose = false)
  {
    bool isEqual = true;
    for (int var = 0; var < size; ++var)
    {
      isEqual = isEqual && Equal(arrayType1[var], arrayType2[var], eps);
    }
 
    ConditionalOutputOfDifference(arrayType1, arrayType2, eps, verbose, isEqual);
 
    return isEqual;
  }
}
 
#endif