---

HxSftMoments.h

/*
 *  Copyright (c) 1998, University of Amsterdam, The Netherlands.
 *  All rights reserved.
 *
 *
 *  Author(s):
 *  Dennis Koelma (koelma@wins.uva.nl)
 *  Edo Poll (poll@wins.uva.nl)
 */

#ifndef HxSftMoments_h
#define HxSftMoments_h

#include "HxSampleFunctorTem.h"

/*  Generates all moments up to order N.
 *
 *  Mpq = Sum[ x^p * y^q * f(x,y) ]
 *        x,y
 *
 *  The order in the resulting list is determined by:
 *  for (int q=0 ; q<=N ; q++)
 *      for (int p=0 ; p<=N ; p++)
 *          if (p+q <= N)
 *              [Mpq];
 *  For N=1 the order is: m00, m10, m01
 *  For N=2 the order is: m00, m10, m20, m01, m11, m02
 *  For N=3 the order is: m00, m10, m20, m30, m01, m11, m21, m02, m12, m03.
 */
template <class ArgType, class ResType>
class HxSftMoments : public HxSampleFunctorTem<ArgType,ResType>
{
public:
                        HxSftMoments(int order, HxString name);
    virtual             ~HxSftMoments();

    virtual void        init();
    virtual void        next(ArgType pixV, ArgType maskV, HxPoint p);
    virtual void        result(HxValueListBackInserter res);

private:
    int                 _order;
    int                 _number;
    ResType*            _sums;
};

#ifdef INC_TEMPLATE_SRC
#include "HxSftMoments.c"
#endif

#endif

---