|| || || || || || || ||

---

PixStat.h

Go to the documentation of this file.

#ifndef Impala_Core_Array_PixStat_h
#define Impala_Core_Array_PixStat_h

#include "Core/Array/Pattern/PatInOutOp.h"
#include "Core/Array/Trait/ExportStat.h"
#include "Core/Array/Element/E1Cast.h"

namespace Impala
{
namespace Core
{
namespace Array
{


template<class ArrayT, class ArithT>
inline void
PixStat(ArrayT* im, ArithT* minVal, ArithT* maxVal, ArithT* mean,
        ArithT* variance)
{
    Trait::ExportStat<ArithT, ArithT> pixOp;
    Pattern::PatInOutOp(im, pixOp);
    if (minVal)
        *minVal = pixOp.mMinVal;
    if (maxVal)
        *maxVal = pixOp.mMaxVal;
    ArithT num = Element::E1Cast(pixOp.mNum, ArithT());
    ArithT num1 = num - Element::E1Cast(1, ArithT());
    if (mean)
        *mean = pixOp.mSum / num;
    if (variance)
        *variance = (pixOp.mNum == 1) ? Element::E1Cast(0, ArithT())
            : (pixOp.mSumSqr - ((pixOp.mSum * pixOp.mSum) / num)) / (num1);
}

template<class ArrayT, class ArithT>
inline void
PixStatBasic(ArrayT* im, ArithT* nr, ArithT* sum, ArithT* sumSqr, ArithT* sumAbs)
{
    Trait::ExportStat<ArithT, ArithT> pixOp;
    Pattern::PatInOutOp(im, pixOp);
    if (nr)
        *nr = Element::E1Cast(pixOp.mNum, ArithT());
    if (sum)
        *sum = pixOp.mSum;
    if (sumSqr)
        *sumSqr = pixOp.mSumSqr;
    if (sumAbs)
        *sumAbs = pixOp.mSumAbs;
}

} // namespace Array
} // namespace Core
} // namespace Impala

#endif

---