zxhMetricJE.h

/*=========================================================================

  Program:   ZXH Registration Software
  Author:    Xiahai Zhuang
  Module:    $RCSfle: zxhMetricJE.h    $
  Language:  C++
  Date:      $Date: From  2007-07 $
  Version:   $Revision: 2.0 $
  revision:
=========================================================================*/
#ifndef zxhMetricJE_h
#define zxhMetricJE_h

#include "zxhMetricBase.h"
class zxhMetricJE :
    public zxhMetricBase
{
public:
    zxhMetricJE(void);
    ~zxhMetricJE(void);

    virtual std::string GetMetricType() {return "_JE_"; } ;

    const static int m_stc_iPaddedBins =4;

    virtual zxhMetricBase* Clone(zxhMetricBase*&pRet);

    virtual bool SetImage(zxhImageData*pImageTest,zxhImageData*pImageRef);

    virtual int     GetValue(zxhlfloat& fValue); //

    void    SetNumberOfHistogramBins(int i=64)          {m_iNumberOfHistogramBins=i;};

    int GetNumberOfHistogramBins()          {return m_iNumberOfHistogramBins;};

    virtual std::string GetPrintString();

    virtual bool        SetMetricFromStream(std::ifstream & ifs);

    virtual bool PreComputeImageGraylevel();
protected:
    int m_iNumberOfHistogramBins;
    PixelTypeDefault m_aiImageTestGraySpan[2];

    PixelTypeDefault m_aiImageRefGraySpan[2];

    float   m_fEntropyT;

    float   m_fEntropyR;

    float   m_fEntropyJ;

    float   m_fTotalFreq; 

    float   *m_afHistogram;

    inline float    BSplinei(int iOrd,float u)
    { // -1,0,1,2
        float v;
        switch(iOrd)
        {
        case -1:
            v=1.0f-u;
            return v*v*v/6.0f;
        case 0:
            v=u*u;
            return 0.5f*u*v-v+2.0f/3.0f;
        case 1:
            v=u*u;
            return (-3.0f*u*v+3.0f*v+3.0f*u+1)/6;
        case 2:return u*u*u/6;
        }
        return -1;
    }
 
    int ComputeHistogramFromPDF4D(float * &afPDF,float&fTotalFreq,
                        float xfTestFrom,float xfTestTo,float yfTestFrom,float yfTestTo,
                        float zfTestFrom,float zfTestTo,float tfTestFrom,float tfTestTo);
};

#endif