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

---

InvWCxy.h

Go to the documentation of this file.

#ifndef Impala_Core_Array_InvWCxy_h
#define Impala_Core_Array_InvWCxy_h

#include <vector>
#include "Core/Array/Arrays.h"
#include "Core/Array/InvCompDer.h"
#include "Core/Array/Rgb2Ooo.h"
#include "Core/Array/Add.h"
#include "Core/Array/Sub.h"
#include "Core/Array/Div.h"
#include "Core/Array/Norm2Sqr.h"
#include "Core/Array/Pattern/PatM9N4PixOp.h"
#include "Core/Array/Trait/M9N4poCxy.h"

namespace Impala
{
namespace Core
{
namespace Array
{


/*
input is an rgb image output is a list of:
- gradient (Ex*Ex + Ey*Ey)
- Wx
- W 45
- Wy
- W 135
- Clx
- Cl 45
- Cly
- Cl 135
- Cllx
- Cll 45
- Clly
- Cll 135

perhaps it would be better to return gradient through another parameter??
*/

inline std::vector<Array2dScalarReal64*>
//InvWiccest(Array2dVec3UInt8* rgb, double derivativeSigma)
InvWCxy(Array2dVec3UInt8* rgb, double derivativeSigma)
{
    std::vector<Array2dScalarReal64*> result;
    result.resize(13);
    for(int i=0 ; i<13 ; ++i)
        result[i] = new Array2dScalarReal64(rgb->CW(), rgb->CH(), rgb->BW(),
                                            rgb->BH());

    // convert to ooo
    Array2dVec3Real64* ooo = new Array2dVec3Real64(rgb->mCW, rgb->mCH, 0, 0);
    Rgb2Ooo(ooo, rgb);
    std::vector<Array2dScalarReal64*> listDer = InvCompDer(ooo, derivativeSigma,
                                                               3);
    delete ooo;

    // compute Wx and Wy
    Div(result[1], listDer[1], listDer[0]);
    Div(result[3], listDer[2], listDer[0]);

    // compute {Cl,Cll}{x,y}
    Trait::M9N4poCxy<Array2dScalarReal64, Array2dScalarReal64> mnpoConFeat;
    Pattern::PatM9N4PixOp(result[5], result[7], result[9], result[11],
                          listDer[0], listDer[1], listDer[2],
                          listDer[3], listDer[4], listDer[5],
                          listDer[6], listDer[7], listDer[8], mnpoConFeat);

    // compute the rest
    Norm2Sqr(result[0], listDer[1], listDer[2]); //gradient
    Add(result[ 2], result[1], result[3]); //W 45
    Sub(result[ 4], result[1], result[3]); //W 135
    Add(result[ 6], result[5], result[7]); //Cl 45
    Sub(result[ 8], result[5], result[7]); //Cl 135
    Add(result[10], result[9], result[11]); //Cll 45
    Sub(result[12], result[9], result[11]); //Cll 135

    for(i=0 ; i<listDer.size() ; i++)
        delete listDer[i];

    return result;
}

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

#endif

---