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

---

InvQuasi.h

Go to the documentation of this file.

#ifndef Impala_Core_Array_InvQuasi_h
#define Impala_Core_Array_InvQuasi_h

#include <vector>
#include "Core/Array/InvCompDer.h"
#include "Core/Array/InvOpponentDer.h"
#include "Core/Array/InvSphericalDer.h"
#include "Core/Array/InvHsiDer.h"
#include "Core/Array/Arrays.h"
#include "Core/Array/Add.h"
#include "Core/Array/AddSquared.h"
#include "Core/Array/Sqrt.h"

namespace Impala
{
namespace Core
{
namespace Array
{


inline std::vector<Array2dScalarReal64*>
InvQuasi(Array2dVec3Real64* src, double sigma, double precision)
{
    std::vector<Array2dScalarReal64*> quasiInvariants;

    // get needed Derivatives
    std::vector<Array2dScalarReal64*> rgbDerList = InvCompDer(src, sigma, precision);
    std::vector<Array2dScalarReal64*> opponentDerList = InvOpponentDer(rgbDerList);
    std::vector<Array2dScalarReal64*> sphericalDerList = InvSphericalDer(rgbDerList);
    std::vector<Array2dScalarReal64*> hsiDerList = InvHsiDer(rgbDerList);

    //specular quasi invariant & variant
    //sp_inv=sqrt(O1_x.^2+O1_y.^2+O2_x.^2+O2_y.^2+eps);
    Array2dScalarReal64* temp = 0;
    Array2dScalarReal64* image = 0;
    AddSquared(image, opponentDerList[0], opponentDerList[1]);
    AddSquared(temp, opponentDerList[2], opponentDerList[3]);
    Add(image, image, temp);
    Sqrt(image, image);
    quasiInvariants.push_back(image);
 
    //sp_var=sqrt(O3_x.^2+O3_y.^2+eps);
    temp = 0;
    image = 0;
    AddSquared(image, opponentDerList[4], opponentDerList[5]);
    Sqrt(image, image);
    quasiInvariants.push_back(image);

    //shadow-shading quasi invariant & variant
    //ss_inv=sqrt(theta_x.^2+theta_y.^2+phi_x.^2+phi_y.^2+eps);
    temp = 0;
    image = 0;
    AddSquared(image, sphericalDerList[0], sphericalDerList[1]);
    AddSquared(temp, sphericalDerList[2], sphericalDerList[3]);
    Add(image, image, temp);
    Sqrt(image, image);
    quasiInvariants.push_back(image);

    //ss_var=sqrt(r_x.^2+r_y.^2+eps);
    temp = 0;
    image = 0;
    AddSquared(image, sphericalDerList[4], sphericalDerList[5]);
    Sqrt(image, image);
    quasiInvariants.push_back(image);
    
    //specular shadow-shading quasi invariant & variant
    //spss_inv=sqrt(h_x.^2+h_y.^2+eps);
    temp = 0;
    image = 0;
    AddSquared(image, hsiDerList[0], hsiDerList[1]);
    Sqrt(image, image);
    quasiInvariants.push_back(image);

    //spss_var=sqrt(i_x.^2+i_y.^2+s_x.^2+s_y.^2+eps);
    temp = 0;
    image = 0;
    AddSquared(image, hsiDerList[2], hsiDerList[3]);
    AddSquared(temp, hsiDerList[4], hsiDerList[5]);
    Add(image, image, temp);
    Sqrt(image, image);
    quasiInvariants.push_back(image);

    int i;
    for(i=0 ; i<rgbDerList.size() ; ++i)
        delete rgbDerList[i];
    for(i=0 ; i<opponentDerList.size() ; ++i)
        delete opponentDerList[i];
    for(i=0 ; i<sphericalDerList.size() ; ++i)
        delete sphericalDerList[i];
    for(i=0 ; i<hsiDerList.size() ; ++i)
        delete hsiDerList[i];

    return quasiInvariants;
}

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

#endif

---