00001 #ifndef Impala_Core_Array_Trait_M9poHUEw_h 00002 #define Impala_Core_Array_Trait_M9poHUEw_h 00003 00004 #include "Core/Array/Pattern/Categories.h" 00005 #include "Core/Array/Element/E1Sqrt.h" 00006 00007 namespace Impala 00008 { 00009 namespace Core 00010 { 00011 namespace Array 00012 { 00013 namespace Trait 00014 { 00015 00016 00019 template<class DstArrayT, class SrcArrayT> 00020 class M9poHUEw 00021 { 00022 public: 00024 typedef Pattern::TagTransInVar TransVarianceCategory; 00025 00027 typedef Pattern::TagCallValue CallCategory; 00028 00029 typedef typename DstArrayT::ArithType DstArithT; 00030 typedef typename SrcArrayT::ArithType SrcArithT; 00031 00033 M9poHUEw() 00034 { 00035 } 00036 00038 DstArithT 00039 DoIt(const SrcArithT& R, const SrcArithT& G, const SrcArithT& B, 00040 const SrcArithT& Rx, const SrcArithT& Gx, const SrcArithT& Bx, 00041 const SrcArithT& Ry, const SrcArithT& Gy, const SrcArithT& By) 00042 { 00043 SrcArithT div = (B*B+G*G+R*R-G*R-B*R-B*G); 00044 SrcArithT f_Obx = (G-B)*Rx+(B-R)*Gx+(R-G)*Bx; 00045 SrcArithT f_Oby = (G-B)*Ry+(B-R)*Gy+(R-G)*By; 00046 00047 SrcArithT f_Hue = .5*sqrt(3.)*Element::E1Sqrt(f_Obx*f_Obx + f_Oby*f_Oby) / div; 00048 return f_Hue; 00049 } 00050 00051 }; 00052 00053 } // namespace Trait 00054 } // namespace Array 00055 } // namespace Core 00056 } // namespace Impala 00057 00058 #endif