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

---

MatMahalanobis.h

Go to the documentation of this file.

#ifndef Impala_Core_Matrix_MatMahalanobis_h
#define Impala_Core_Matrix_MatMahalanobis_h

#include "Core/Matrix/MatSet.h"

namespace Impala
{
namespace Core
{
namespace Matrix
{


template<class ArrayT>
inline ArrayT*
MatMahalanobis(ArrayT* a, ArrayT* b, ArrayT* c)
{
    typedef typename ArrayT::ArithType ArithT;
    typedef typename ArrayT::StorType StorT;

    int nr = MatNrRow(a);
    int nc = MatNrCol(b);
    ArrayT* m = MatCreate<ArrayT>(nr, nc);
    MatSet(m, 0);

    for (int i=0 ; i<nr ; i++)   //for each pixel
    {
        StorT* miPtr = MatE(m, i, 0);
        StorT* aiPtr = MatE(a, i, 0);
        for (int j=0 ; j<nc ; j++)   //for each concept
        {
            double sum = 0;
            for (int k=0 ; k<MatNrCol(a) ; k++) //for each filter
            {
                double ckj = *MatE(c, k, j);
                double bkj = *MatE(b, k, j);
                if (ckj > 0)  //is this correct?
                {
                    sum += ( (aiPtr[k] - bkj) * (aiPtr[k] - bkj) ) / (ckj*ckj);
                }
            }   
            miPtr[j] = sqrt( sum ); 
        }
    }
    return m;
}

} // namespace Matrix
} // namespace Core
} // namespace Impala

#endif

---