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

---

M9N4poFHS.h

Go to the documentation of this file.

#ifndef Impala_Core_Array_Trait_M9N4poFHS_h
#define Impala_Core_Array_Trait_M9N4poFHS_h

#include "Core/Array/Pattern/Categories.h"
#include "Core/Array/Element/E1Sqrt.h"

namespace Impala
{
namespace Core
{
namespace Array
{
namespace Trait
{


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

    typedef Pattern::TagCallValue CallCategory;

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

    M9N4poFHS()
    {
        mC = ::sqrt(3.0);
    }

    void
    DoIt(DstStorT* FHx, DstStorT* FHy, DstStorT* FSx, DstStorT* FSy,
         const SrcArithT& R, const SrcArithT& G, const SrcArithT& B,
         const SrcArithT& Rx, const SrcArithT& Gx, const SrcArithT& Bx,
         const SrcArithT& Ry, const SrcArithT& Gy, const SrcArithT& By)
    {
        SrcArithT saturation = Element::E1Sqrt(2*(R*R + G*G + B*B - R*G - R*B - G*B));
        
        *FHx = (saturation == 0) ? 0 : (R*(Bx-Gx) + G*(Rx-Bx) + B*(Gx-Rx)) / saturation;
        *FHy = (saturation == 0) ? 0 : (R*(By-Gy) + G*(Ry-By) + B*(Gy-Ry)) / saturation;
        saturation *= mC; //because S divides by sqrt(3)*saturation
        *FSx = (saturation == 0) ? 0 : (R*(2.*Rx-Gx-Bx)+G*(2.*Gx-Rx-Bx)+B*(2.*Bx-Rx-Gx)) / saturation;
        *FSy = (saturation == 0) ? 0 : (R*(2.*Ry-Gy-By)+G*(2.*Gy-Ry-By)+B*(2.*By-Ry-Gy)) / saturation;
    }

    SrcArithT mC;
};

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

#endif

---