itkMeasurementVectorTraits.h

Go to the documentation of this file.

/*=========================================================================
 *
 *  Copyright Insight Software Consortium
 *
 *  Licensed under the Apache License, Version 2.0 (the "License");
 *  you may not use this file except in compliance with the License.
 *  You may obtain a copy of the License at
 *
 *         http://www.apache.org/licenses/LICENSE-2.0.txt
 *
 *  Unless required by applicable law or agreed to in writing, software
 *  distributed under the License is distributed on an "AS IS" BASIS,
 *  WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
 *  See the License for the specific language governing permissions and
 *  limitations under the License.
 *
 *=========================================================================*/
#ifndef itkMeasurementVectorTraits_h
#define itkMeasurementVectorTraits_h

#include "itkVariableLengthVector.h"
#include "vnl/vnl_vector_fixed.h"
#include "itkRGBPixel.h"
#include "itkMatrix.h"
#include "itkVariableSizeMatrix.h"
#include "itkNumericTraits.h"
#include "itkNumericTraitsStdVector.h"
#include "itkSize.h"
#include <vector>

namespace itk
{
namespace Statistics
{
class MeasurementVectorTraits
{
public:

  typedef IdentifierType  InstanceIdentifier;

  typedef InstanceIdentifier                                     AbsoluteFrequencyType;
  typedef NumericTraits< AbsoluteFrequencyType >::RealType       RelativeFrequencyType;
  typedef NumericTraits< AbsoluteFrequencyType >::AccumulateType TotalAbsoluteFrequencyType;
  typedef NumericTraits< RelativeFrequencyType >::AccumulateType TotalRelativeFrequencyType;

  typedef unsigned int MeasurementVectorLength;

  template< typename TVectorType >
  static bool IsResizable(const TVectorType &)
  {
    // Test whether the vector type is resizable or not
    //
    // If the default constructor creates a vector of
    // length zero, we assume that it is resizable,
    // otherwise that is a pretty useless measurement vector.
    TVectorType             m;
    MeasurementVectorLength len = NumericTraits<TVectorType>::GetLength(m);

    return ( len == 0 );
  }

  template< typename TValue1, unsigned int VLength, typename TValue2, unsigned int VLength2 >
  static MeasurementVectorLength Assert(const FixedArray< TValue1, VLength > &,
                                        const FixedArray< TValue2, VLength2 > &,
                                        const char *errMsg = "Length Mismatch")
  {
    if ( VLength != VLength2 )
      {
      itkGenericExceptionMacro(<< errMsg);
      }
    return 0;
  }

  template< typename TValue1, unsigned int VLength, typename TValue2, unsigned int VLength2 >
  static MeasurementVectorLength Assert(const FixedArray< TValue1, VLength > *,
                                        const FixedArray< TValue2, VLength2 > *,
                                        const char *errMsg = "Length Mismatch")
  {
    if ( VLength != VLength2 )
      {
      itkGenericExceptionMacro(<< errMsg);
      }
    return 0;
  }

  template< typename TValue1, typename TValue2 >
  static MeasurementVectorLength Assert(const Array< TValue1 > & a,
                                        const Array< TValue2 > & b, const char *errMsg = "Length Mismatch")
  {
    if ( b.Size() != a.Size() )
      {
      itkGenericExceptionMacro(<< errMsg);
      }
    return 0;
  }

  template< typename TValue1, typename TValue2 >
  static MeasurementVectorLength Assert(const Array< TValue1 > *a,
                                        const Array< TValue2 > *b, const char *errMsg = "Length Mismatch")
  {
    if ( b->Size() != a->Size() )
      {
      itkGenericExceptionMacro(<< errMsg);
      }
    return 0;
  }

  template< typename TValue1, typename TValue2 >
  static MeasurementVectorLength Assert(const VariableLengthVector< TValue1 > & a,
                                        const VariableLengthVector< TValue2 > & b,
                                        const char *errMsg = "Length Mismatch")
  {
    if ( b.Size() != a.Size() )
      {
      itkGenericExceptionMacro(<< errMsg);
      }
    return 0;
  }

  template< typename TValue1, typename TValue2 >
  static MeasurementVectorLength Assert(const VariableLengthVector< TValue1 > *a,
                                        const VariableLengthVector< TValue2 > *b,
                                        const char *errMsg = "Length Mismatch")
  {
    if ( b->Size() != a->Size() )
      {
      itkGenericExceptionMacro(<< errMsg);
      }
    return 0;
  }

  template< typename TValue1, typename TValue2 >
  static MeasurementVectorLength Assert(const std::vector< TValue1 > & a,
                                        const std::vector< TValue2 > & b, const char *errMsg = "Length Mismatch")
  {
    if ( b.size() != a.size() )
      {
      itkGenericExceptionMacro(<< errMsg);
      }
    return 0;
  }

  template< typename TValue1, typename TValue2 >
  static MeasurementVectorLength Assert(const std::vector< TValue1 > *a,
                                        const std::vector< TValue2 > *b, const char *errMsg = "Length Mismatch")
  {
    if ( b->size() != a->size() )
      {
      itkGenericExceptionMacro(<< errMsg);
      }
    return 0;
  }

  template< typename TValue1, unsigned int VLength, typename TValue2 >
  static MeasurementVectorLength Assert(const FixedArray< TValue1, VLength > &,
                                        const Array< TValue2 > & b, const char *errMsg = "Length Mismatch")
  {
    if ( b.Size() == 0 )
      {
      return VLength;
      }
    if ( b.Size() != 0 )
      {
      if ( b.Size() != VLength )
        {
        itkGenericExceptionMacro(<< errMsg);
        }
      }
    return 0;
  }

  template< typename TValue1, unsigned int VLength, typename TValue2 >
  static MeasurementVectorLength Assert(const FixedArray< TValue1, VLength > *,
                                        const Array< TValue2 > *b, const char *errMsg = "Length Mismatch")
  {
    if ( b->Size() == 0 )
      {
      return VLength;
      }
    else if ( b->Size() != VLength )
      {
      itkGenericExceptionMacro(<< errMsg);
      }
    return 0;
  }

  template< typename TValue1, unsigned int VLength, typename TValue2 >
  static MeasurementVectorLength Assert(const FixedArray< TValue1, VLength > &,
                                        const VariableLengthVector< TValue2 > & b,
                                        const char *errMsg = "Length Mismatch")
  {
    if ( b.Size() == 0 )
      {
      return VLength;
      }
    if ( b.Size() != 0 )
      {
      if ( b.Size() != VLength )
        {
        itkGenericExceptionMacro(<< errMsg);
        }
      }
    return 0;
  }

  template< typename TValue1, unsigned int VLength, typename TValue2 >
  static MeasurementVectorLength Assert(const FixedArray< TValue1, VLength > *,
                                        const VariableLengthVector< TValue2 > *b,
                                        const char *errMsg = "Length Mismatch")
  {
    if ( b->Size() == 0 )
      {
      return VLength;
      }
    else if ( b->Size() != VLength )
      {
      itkGenericExceptionMacro(<< errMsg);
      }
    return 0;
  }

  template< typename TValue1, unsigned int VLength, typename TValue2 >
  static MeasurementVectorLength Assert(const FixedArray< TValue1, VLength > &,
                                        const std::vector< TValue2 > & b, const char *errMsg = "Length Mismatch")
  {
    if ( b.size() == 0 )
      {
      return VLength;
      }
    if ( b.size() != 0 )
      {
      if ( b.size() != VLength )
        {
        itkGenericExceptionMacro(<< errMsg);
        }
      }
    return 0;
  }

  template< typename TValue1, unsigned int VLength, typename TValue2 >
  static MeasurementVectorLength Assert(const FixedArray< TValue1, VLength > *,
                                        const std::vector< TValue2 > *b, const char *errMsg = "Length Mismatch")
  {
    if ( b->size() == 0 )
      {
      return VLength;
      }
    else if ( b->size() != VLength )
      {
      itkGenericExceptionMacro(<< errMsg);
      }
    return 0;
  }

  template< typename TValue1, unsigned int VLength >
  static MeasurementVectorLength Assert(const FixedArray< TValue1, VLength > &,
                                        const MeasurementVectorLength l, const char *errMsg = "Length Mismatch")
  {
    if ( l == 0 )
      {
      return VLength;
      }
    else if ( l != VLength )
      {
      itkGenericExceptionMacro(<< errMsg);
      }
    return 0;
  }

  template< typename TValue1, unsigned int VLength >
  static MeasurementVectorLength Assert(const FixedArray< TValue1, VLength > *,
                                        const MeasurementVectorLength l, const char *errMsg = "Length Mismatch")
  {
    if ( l == 0 )
      {
      return VLength;
      }
    else if ( l != VLength )
      {
      itkGenericExceptionMacro(<< errMsg);
      }
    return 0;
  }

  template< typename TValue >
  static MeasurementVectorLength Assert(const Array< TValue > & a,
                                        const MeasurementVectorLength l, const char *errMsg = "Length Mismatch")
  {
    if ( ( ( l != 0 ) && ( a.Size() != l ) ) || ( a.Size() == 0 ) )
      {
      itkGenericExceptionMacro(<< errMsg);
      }
    else if ( l == 0 )
      {
      return a.Size();
      }
    return 0;
  }

  template< typename TValue >
  static MeasurementVectorLength Assert(const Array< TValue > *a,
                                        const MeasurementVectorLength l, const char *errMsg = "Length Mismatch")
  {
    if ( ( ( l != 0 ) && ( a->Size() != l ) ) || ( a->Size() == 0 ) )
      {
      itkGenericExceptionMacro(<< errMsg);
      }
    else if ( l == 0 )
      {
      return a->Size();
      }
    return 0;
  }

  template< typename TValue >
  static MeasurementVectorLength Assert(const VariableLengthVector< TValue > & a,
                                        const MeasurementVectorLength l, const char *errMsg = "Length Mismatch")
  {
    if ( ( ( l != 0 ) && ( a.Size() != l ) ) || ( a.Size() == 0 ) )
      {
      itkGenericExceptionMacro(<< errMsg);
      }
    else if ( l == 0 )
      {
      return a.Size();
      }
    return 0;
  }

  template< typename TValue >
  static MeasurementVectorLength Assert(const VariableLengthVector< TValue > *a,
                                        const MeasurementVectorLength l, const char *errMsg = "Length Mismatch")
  {
    if ( ( ( l != 0 ) && ( a->Size() != l ) ) || ( a->Size() == 0 ) )
      {
      itkGenericExceptionMacro(<< errMsg);
      }
    else if ( l == 0 )
      {
      return a->Size();
      }
    return 0;
  }

  template< typename TValue >
  static MeasurementVectorLength Assert(const std::vector< TValue > & a,
                                        const MeasurementVectorLength l, const char *errMsg = "Length Mismatch")
  {
    if ( ( ( l != 0 ) && ( a.size() != l ) ) || ( a.size() == 0 ) )
      {
      itkGenericExceptionMacro(<< errMsg);
      }
    else if ( l == 0 )
      {
      return static_cast<MeasurementVectorLength>( a.size() );
      }
    return 0;
  }

  template< typename TValue >
  static MeasurementVectorLength Assert(const std::vector< TValue > *a,
                                        const MeasurementVectorLength l, const char *errMsg = "Length Mismatch")
  {
    if ( ( ( l != 0 ) && ( a->size() != l ) ) || ( a->size() == 0 ) )
      {
      itkGenericExceptionMacro(<< errMsg);
      }
    else if ( l == 0 )
      {
      return a->size();
      }
    return 0;
  }

  template< typename TArrayType >
  static void  Assign(TArrayType & m, const TArrayType & v)
  {
    m = v;
  }

  template< typename TValue, unsigned int VLength >
  static void  Assign(FixedArray< TValue, VLength > & m, const TValue & v)
  {
    m[0] = v;
  }
};

template< typename TMeasurementVector >
class MeasurementVectorTraitsTypes
{
public:
  typedef typename TMeasurementVector::ValueType ValueType;
};

template< typename T >
class MeasurementVectorTraitsTypes< std::vector< T > >
{
public:
  typedef T ValueType;
};

template< typename TPixelType >
class MeasurementVectorPixelTraits
{
public:
  /* type of the vector that matches this pixel type */
  typedef TPixelType MeasurementVectorType;
};


template< >
class MeasurementVectorPixelTraits< char >
{
public:
  typedef FixedArray< char, 1 > MeasurementVectorType;
};

template< >
class MeasurementVectorPixelTraits< unsigned char >
{
public:
  typedef FixedArray< unsigned char, 1 > MeasurementVectorType;
};

template< >
class MeasurementVectorPixelTraits< signed char >
{
public:
  typedef FixedArray< signed char, 1 > MeasurementVectorType;
};

template< >
class MeasurementVectorPixelTraits< unsigned short >
{
public:
  typedef FixedArray< unsigned short, 1 > MeasurementVectorType;
};

template< >
class MeasurementVectorPixelTraits< signed short >
{
public:
  typedef FixedArray< signed short, 1 > MeasurementVectorType;
};

template< >
class MeasurementVectorPixelTraits< unsigned int >
{
public:
  typedef FixedArray< unsigned int, 1 > MeasurementVectorType;
};

template< >
class MeasurementVectorPixelTraits< signed int >
{
public:
  typedef FixedArray< signed int, 1 > MeasurementVectorType;
};

template< >
class MeasurementVectorPixelTraits< unsigned long >
{
public:
  typedef FixedArray< unsigned long, 1 > MeasurementVectorType;
};

template< >
class MeasurementVectorPixelTraits< signed long >
{
public:
  typedef FixedArray< signed long, 1 > MeasurementVectorType;
};

template< >
class MeasurementVectorPixelTraits< unsigned long long >
{
public:
  typedef FixedArray< unsigned long long, 1 > MeasurementVectorType;
};

template< >
class MeasurementVectorPixelTraits< signed long long >
{
public:
  typedef FixedArray< signed long long, 1 > MeasurementVectorType;
};

template< >
class MeasurementVectorPixelTraits< float >
{
public:
  typedef FixedArray< float, 1 > MeasurementVectorType;
};

template< >
class MeasurementVectorPixelTraits< double >
{
public:
  typedef FixedArray< double, 1 > MeasurementVectorType;
};


} // namespace Statistics
} // namespace itk

#endif  // itkMeasurementVectorTraits_h