img_scalar.h

Go to the documentation of this file.

#ifndef IMG_SCALAR_H_
#define IMG_SCALAR_H_

#include "img/img_base.h"

namespace img {

template<typename ScalarType>
inline ScalarType clampValue(ScalarType value, ScalarType minval=ScalarType(0.0), ScalarType maxval=ScalarType(255.0))
{
  STATIC_FLOAT_OR_DOUBLE_TYPECHECK(ScalarType);
  if(value < minval) return minval;
  if(value > maxval) return maxval;
  return value;
}

template<typename ScalarType>
inline int valueAsInt(ScalarType value)
{
  STATIC_FLOAT_OR_DOUBLE_TYPECHECK(ScalarType);
  return static_cast<int>(floor(static_cast<double>(value)+0.5f));
}

template<typename ScalarType>
inline ScalarType max3(ScalarType a, ScalarType b, ScalarType c)
{
  STATIC_FLOAT_OR_DOUBLE_TYPECHECK(ScalarType);
  if ( b<=a && c<=a ) return a;
  if ( a<=b && c<=b ) return b;
  return c;
}

template<typename ScalarType>
inline ScalarType min3(ScalarType a, ScalarType b, ScalarType c)
{
  STATIC_FLOAT_OR_DOUBLE_TYPECHECK(ScalarType);
  if ( a<=b && a<=c ) return a;
  if ( b<=a && b<=c ) return b;
  return c;
}

template<typename ScalarType>
inline ScalarType median3(ScalarType a, ScalarType b, ScalarType c)
{
  STATIC_FLOAT_OR_DOUBLE_TYPECHECK(ScalarType);
  if ( (b<=a && a<=c) || (c<=a && a<=b) ) return a;
  if ( (a<=b && b<=c) || (c<=b && b<=a) ) return b;
  return c;
}

template<typename ScalarType>
inline bool almostEqual(ScalarType a, ScalarType b, ScalarType EPSILON=ScalarType(10e-5))
{
  STATIC_FLOAT_OR_DOUBLE_TYPECHECK(ScalarType);
  ScalarType diff=a-b;
  if (diff<0)
    diff *= ScalarType(-1);
  return diff < EPSILON;
}

} //end namespace img

#endif /*IMG_SCALAR_H_*/