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

---

TrecFeature.h

Go to the documentation of this file.

#ifndef Impala_Core_Trec_TrecFeature_h
#define Impala_Core_Trec_TrecFeature_h

#include <vector>
#include "Core/Array/Element/E1Sqrt.h"
#include "Core/Array/Pattern/PatM9N4PixOp.h"
#include "Core/Array/ColorSpace.h"
#include "Core/Array/Scale.h"
#include "Core/Array/ProjectRange.h"
#include "Core/Array/RecGauss.h"
#include "Core/Array/GaussDerivative.h"

namespace Impala
{
namespace Core
{
namespace Trec
{


using namespace Impala::Core::Array;

bool
IsWfeat(int f)
{
        int i = f % 4;
        return i >= 2;
}

template<class DstArrayT, class SrcArrayT>
class M9N4poConFeat
{
public:
    typedef Pattern::TagTransInVar TransVarianceCategory;

    typedef Pattern::TagCallValue CallCategory;

        typedef typename DstArrayT::StorType DstStorT;
        typedef typename DstArrayT::ArithType DstArithT;
        typedef typename SrcArrayT::ArithType SrcArithT;

        M9N4poConFeat(SrcArithT scale)
        {
                _scale = scale;
        }

    void
        DoIt(DstStorT* Cl, DstStorT* Cll, DstStorT* Clw, DstStorT* Cllw,
                 const SrcArithT& E, const SrcArithT& Ex, const SrcArithT& Ey,
                 const SrcArithT& El, const SrcArithT& Elx, const SrcArithT& Ely,
                 const SrcArithT& Ell, const SrcArithT& Ellx, const SrcArithT& Elly)
        {
                *Cl = El / E;
                *Cll = Ell / E;
                SrcArithT tmpx = Elx*E - El*Ex;
                SrcArithT tmpy = Ely*E - El*Ey;
                *Clw = ( Element::E1Sqrt(tmpx*tmpx + tmpy*tmpy) / (E*E) ) * _scale;
                tmpx = Ellx*E - Ell*Ex;
                tmpy = Elly*E - Ell*Ey;
                *Cllw = ( Element::E1Sqrt(tmpx*tmpx + tmpy*tmpy) / (E*E) ) * _scale;
        }

        SrcArithT _scale;
};

template<class ArrayT>
void
MyRestrict(ArrayT* im, int border)
{
        im->mBW += border;
        im->mBH += border;
        im->mCW -= 2*border;
        im->mCH -= 2*border;
}

std::vector<Array2dScalarReal64*>
BuildInvCwFeatureList(Array2dVec3Real64* im, double precision, int border)
{
        Array2dVec3Real64* imO = 0;
        ColorSpace(imO, im, Element::RGB, Element::OOO);

        const int nrScales = 3;
        //double scales[nrScales] = {1.0, 2.0, 3.5};
        double scales[nrScales] = {2.0, 4.0, 8.0};

        std::vector<Array2dScalarReal64*> imList;
        Array2dScalarReal64* ch = 0;
        Array2dScalarReal64* res = 0;
        Array2dScalarReal64* E = 0;
        Array2dScalarReal64* Ex = 0;
        Array2dScalarReal64* Ey = 0;
        Array2dScalarReal64* El = 0;
        Array2dScalarReal64* Elx = 0;
        Array2dScalarReal64* Ely = 0;
        Array2dScalarReal64* Ell = 0;
        Array2dScalarReal64* Ellx = 0;
        Array2dScalarReal64* Elly = 0;
        //for all scales
        for (int i=0 ; i<nrScales ; i++)
    {
                ProjectRange(ch, imO, 1);
                RecGauss(E, ch, scales[i], scales[i], 0, 0, precision);         
                RecGauss(Ex, ch, scales[i], scales[i], 1, 0, precision);                
                RecGauss(Ey, ch, scales[i], scales[i], 0, 1, precision);                

                ProjectRange(ch, imO, 2);
                RecGauss(El, ch, scales[i], scales[i], 0, 0, precision);                
                RecGauss(Elx, ch, scales[i], scales[i], 1, 0, precision);               
                RecGauss(Ely, ch, scales[i], scales[i], 0, 1, precision);               

                ProjectRange(ch, imO, 3);
                RecGauss(Ell, ch, scales[i], scales[i], 0, 0, precision);               
                RecGauss(Ellx, ch, scales[i], scales[i], 1, 0, precision);              
                RecGauss(Elly, ch, scales[i], scales[i], 0, 1, precision);              

                M9N4poConFeat<Array2dScalarReal64, Array2dScalarReal64>
                        mnpoConFeat(scales[i]);
                std::vector<Array2dScalarReal64*> fList;
                fList = Pattern::PatM9N4PixOp<Array2dScalarReal64>
            (E, Ex, Ey, El, Elx, Ely, Ell, Ellx, Elly, mnpoConFeat);
                if (border != 0)
        {
                        MyRestrict(fList[0], border);
                        MyRestrict(fList[1], border);
                        MyRestrict(fList[2], border);
                        MyRestrict(fList[3], border);
                }

                imList.push_back(fList[0]); // Cl
                imList.push_back(fList[1]); // Cll
                imList.push_back(fList[2]); // Clw
                imList.push_back(fList[3]); // Cllw
        }
        delete imO;
        delete ch;
        delete res;
        delete E;
        delete Ex;
        delete Ey;
        delete El;
        delete Elx;
        delete Ely;
        delete Ell;
        delete Ellx;
        delete Elly;
        return imList;
}

std::vector<Array2dScalarReal64*>
BuildInvECwFeatureList(Array2dVec3Real64* im, double precision, int border)
{
        Array2dVec3Real64* imO = 0;
        ColorSpace(imO, im, Element::RGB, Element::OOO);

        const int nrScales = 3;
        double scales[nrScales] = {1.0, 2.0, 3.5};
        //double scales[nrScales] = {2.0, 4.0, 8.0};

        std::vector<Array2dScalarReal64*> imList;
        Array2dScalarReal64* ch = 0;
        Array2dScalarReal64* res = 0;
        Array2dScalarReal64* E = 0;
        Array2dScalarReal64* Ex = 0;
        Array2dScalarReal64* Ey = 0;
        Array2dScalarReal64* El = 0;
        Array2dScalarReal64* Elx = 0;
        Array2dScalarReal64* Ely = 0;
        Array2dScalarReal64* Ell = 0;
        Array2dScalarReal64* Ellx = 0;
        Array2dScalarReal64* Elly = 0;
        //for all scales
        for (int i=0 ; i<nrScales ; i++)
    {
                E = 0;
                El = 0;
                Ell = 0;

                ProjectRange(ch, imO, 1);
                RecGauss(E, ch, scales[i], scales[i], 0, 0, precision);         
                RecGauss(Ex, ch, scales[i], scales[i], 1, 0, precision);                
                RecGauss(Ey, ch, scales[i], scales[i], 0, 1, precision);                

                ProjectRange(ch, imO, 2);
                RecGauss(El, ch, scales[i], scales[i], 0, 0, precision);                
                RecGauss(Elx, ch, scales[i], scales[i], 1, 0, precision);               
                RecGauss(Ely, ch, scales[i], scales[i], 0, 1, precision);               

                ProjectRange(ch, imO, 3);
                RecGauss(Ell, ch, scales[i], scales[i], 0, 0, precision);               
                RecGauss(Ellx, ch, scales[i], scales[i], 1, 0, precision);              
                RecGauss(Elly, ch, scales[i], scales[i], 0, 1, precision);              

                M9N4poConFeat<Array2dScalarReal64, Array2dScalarReal64>
                        mnpoConFeat(scales[i]);
                std::vector<Array2dScalarReal64*> fList;
                fList = Pattern::PatM9N4PixOp<Array2dScalarReal64>
            (E, Ex, Ey, El, Elx, Ely, Ell, Ellx, Elly, mnpoConFeat);
                if (border != 0)
        {
                        MyRestrict(E, border);
                        MyRestrict(El, border);
                        MyRestrict(Ell, border);
                        MyRestrict(fList[0], border);
                        MyRestrict(fList[1], border);
                        MyRestrict(fList[2], border);
                        MyRestrict(fList[3], border);
                }

                imList.push_back(E);
                imList.push_back(El);
                imList.push_back(Ell);
                imList.push_back(fList[0]); // Cl
                imList.push_back(fList[1]); // Cll
                imList.push_back(fList[2]); // Clw
                imList.push_back(fList[3]); // Cllw
        }
        delete imO;
        delete ch;
        delete res;
        //delete E;
        delete Ex;
        delete Ey;
        //delete El;
        delete Elx;
        delete Ely;
        //delete Ell;
        delete Ellx;
        delete Elly;
        return imList;
}

} // namespace Trec
} // namespace Core
} // namespace Impala

#endif

---