|| ||           || Doxygen's  

---

HxNgbOpticalFlow.h

/*
 *  Copyright (c) 2002, University of Amsterdam, The Netherlands.
 *  All rights reserved.
 *
 *  Author(s):
 *  Leon Todoran            (todoran@science.uva.nl)
 */

#ifndef HxNgbOpticalFlow_h
#define HxNgbOpticalFlow_h

#include "HxClassName.h"
#include "HxSizes.h"
#include "HxCategories.h"
#include "HxCnum.h"
#include "HxTagList.h"


template<class DstT, class SrcT>
class HxNgbOpticalFlow
{
public:

    typedef HxTagLoop           IteratorCategory;

    typedef HxTag1Phase         PhaseCategory;


                            HxNgbOpticalFlow(HxTagList& tags)
                            {
                                _size = HxGetTag(tags, "ngbSize", HxSizes(5, 5, 1));
                                _factor = HxScalarInt(_size.x() * _size.y());
                                _result = HxVec2Double(0,0);
                                
                                //initialize the weighting factors
//w
/*                              _w2[0]  = 0.0625; _w2[1]  = 0.0625; _w2[2]  = 0.0625; _w2[3]  = 0.0625; _w2[4]  = 0.0625;
                                _w2[5]  = 0.0625; _w2[6]  = 0.2500; _w2[7]  = 0.2500; _w2[8]  = 0.2500; _w2[9]  = 0.0625;
                                _w2[10] = 0.0625; _w2[11] = 0.2500; _w2[12] = 0.3750; _w2[13] = 0.2500; _w2[14] = 0.0625;
                                _w2[15] = 0.0625; _w2[16] = 0.2500; _w2[17] = 0.2500; _w2[18] = 0.2500; _w2[19] = 0.0625;
                                _w2[20] = 0.0625; _w2[21] = 0.0625; _w2[22] = 0.0625; _w2[23] = 0.0625; _w2[24] = 0.0625;
*/
//wij*wij
/*
                                _w2[0]  = 0.00390625; _w2[1]  = 0.00390625; _w2[2]  = 0.00390625; _w2[3]  = 0.00390625; _w2[4]  = 0.00390625;
                                _w2[5]  = 0.00390625; _w2[6]  = 0.0625; _w2[7]  = 0.0625; _w2[8]  = 0.0625; _w2[9]  = 0.00390625;
                                _w2[10] = 0.00390625; _w2[11] = 0.0625; _w2[12] = 0.140625; _w2[13] = 0.0625; _w2[14] = 0.00390625;
                                _w2[15] = 0.00390625; _w2[16] = 0.0625; _w2[17] = 0.0625; _w2[18] = 0.0625; _w2[19] = 0.00390625;
                                _w2[20] = 0.00390625; _w2[21] = 0.00390625; _w2[22] = 0.00390625; _w2[23] = 0.00390625; _w2[24] = 0.00390625;
*/
/*
                            0.0195    0.0547    0.0625    0.0547    0.0195
                            0.0547    0.1953    0.2266    0.1953    0.0547
                            0.0625    0.2266    0.2734    0.2266    0.0625
                            0.0547    0.1953    0.2266    0.1953    0.0547
                            0.0195    0.0547    0.0625    0.0547    0.0195
*/
//w*w
                                _w2[0]  = 0.0195; _w2[1]  = 0.0547; _w2[2]  = 0.0625; _w2[3]  = 0.0547; _w2[4]  = 0.0195;
                                _w2[5]  = 0.0547; _w2[6]  = 0.1953; _w2[7]  = 0.2266; _w2[8]  = 0.1953; _w2[9]  = 0.0547;
                                _w2[10] = 0.0625; _w2[11] = 0.2266; _w2[12] = 0.2734; _w2[13] = 0.2266; _w2[14] = 0.0625;
                                _w2[15] = 0.0547; _w2[16] = 0.1953; _w2[17] = 0.2266; _w2[18] = 0.1953; _w2[19] = 0.0547;
                                _w2[20] = 0.0195; _w2[21] = 0.0547; _w2[22] = 0.0625; _w2[23] = 0.0547; _w2[24] = 0.0195;
                            }

    HxSizes                 size()
                            {
                                return _size;
                            }

    void                    init(int x, int y, const SrcT& v1, const SrcT& v2,
                                 const SrcT& v3)
                            {
                                _result = HxVec2Double(0,0);
                                _IxIxw2 = HxScalarDouble(0);
                                _IxIyw2 = HxScalarDouble(0);
                                _IxItw2 = HxScalarDouble(0);
                                _IyItw2 = HxScalarDouble(0);
                                _IyIyw2 = HxScalarDouble(0);
                                
                                _idx = 0;
                            }

    void                    next(int x, int y, const SrcT& v1, const SrcT& v2,
                                 const SrcT& v3)
                            {
                                _IxIxw2 += v1*v1*_w2[_idx];
                                _IxIyw2 += v1*v2*_w2[_idx];
                                _IxItw2 += v1*v3*_w2[_idx];
                                _IyItw2 += v2*v3*_w2[_idx];
                                _IyIyw2 += v2*v2*_w2[_idx];

                                _idx++;
                            }

    DstT                    result() const
                            {
                                SrcT det = (_IxIxw2*_IyIyw2-_IxIyw2*_IxIyw2);
//                              if(det.x()==0)
                                if(fabs(det.x()) <1E-11)
                                    return _result;
                                DstT tmp( ((_IyItw2*_IxIyw2-_IxItw2*_IyIyw2)/det).x() ,
                                          ((_IxItw2*_IxIyw2-_IyItw2*_IxIxw2)/det).x() 
                                          );

                                return tmp;
                            }

    static HxString         className()
                            {
                                return HxString("opticalFlow");
                            }

private:
    HxSizes                 _size;
    DstT                    _result;
    DstT                    _factor;
    SrcT                    _IxIyw2;
    SrcT                    _IxItw2;
    SrcT                    _IyItw2;
    SrcT                    _IxIxw2;
    SrcT                    _IyIyw2;

    double                  _w2[5*5];
    int                     _idx;
};

#endif

---