|| ||           || Doxygen's  

---

HxExportExtraIdentMaskMoments.h

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

#ifndef HxExportExtraIdentMaskMoments_h
#define HxExportExtraIdentMaskMoments_h

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


template<class ResultT, class ImValT, class ExtraValT>
class HxExportExtraIdentMaskMoments
{
public:
    typedef HxTagTransVar       TransVarianceCategory;

    typedef HxTag1Phase         PhaseCategory;


                        HxExportExtraIdentMaskMoments(HxTagList& tags);

                        ~HxExportExtraIdentMaskMoments();

    void                doIt(const ImValT& imV, const ExtraValT& extraV,
                             int x, int y, int z);

    static HxString     className();

private:
    HxTagList&          _tags;
    HxValueList*        _valList;
    ExtraValT           _maskVal;
    int                 _order;
    int                 _number;
    ResultT*            _sums;
};

template<class ResultT, class ImValT, class ExtraValT>

HxExportExtraIdentMaskMoments<ResultT,ImValT,ExtraValT>::HxExportExtraIdentMaskMoments(HxTagList& tags) : _tags(tags)
{
    _maskVal = HxGetTag<int>(tags, "maskVal");
    _order = HxGetTag<int>(tags, "order");
    _valList = HxGetTag<HxValueList*>(tags, "valList");
    _number = ((_order+1)*(_order+2))/2;
    _sums = new ResultT[_number];
    ResultT *ptr = _sums;
    for (int i=0 ; i<_number; i++)
        *ptr++ = HxScalarInt(0);
}

template<class ResultT, class ImValT, class ExtraValT>

HxExportExtraIdentMaskMoments<ResultT,ImValT,ExtraValT>    ::~HxExportExtraIdentMaskMoments()
{
    HxValueListBackInserter res = std::back_inserter(*_valList);
    for (int i=0 ; i<_number ; i++)
        *res++ = HxValue(_sums[i]);
    delete _sums;
}

template<class ResultT, class ImValT, class ExtraValT>
inline void

HxExportExtraIdentMaskMoments<ResultT,ImValT,ExtraValT>    ::doIt(const ImValT& imV, const ExtraValT& extraV, int x, int y, int)
{
    if (extraV == _maskVal) {
        int i=0;
        for (int q=0 ; q<=_order ; q++) {
            for (int p=0 ; p<=_order ; p++) {
                if (p+q <= _order) {
                    HxScalarDouble f = ::pow(x, p) * ::pow(y, q);
                    _sums[i] += ResultT(f) * ResultT(imV);
                    i++;
                }
            }
        }
    }
}

template<class ResultT, class ImValT, class ExtraValT>
inline HxString
HxExportExtraIdentMaskMoments<ResultT,ImValT,ExtraValT>::className()
{
    return HxString("identMaskMoments");
}


#endif

---