zxhMetricMIFFD2.h

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

  Program:   ZXH Registration Software
  Author:    Xiahai Zhuang
  Module:    $RCSfle: zxhMetricMIFFD2.h    $
  Language:  C++
  Date:      $Date: From  2008-10 $
  Version:   $Revision: 2.0 $
  revision:
=========================================================================*/
#ifndef zxhMetricMIFFD2_h
#define zxhMetricMIFFD2_h

#include "zxh.h"
#include <time.h>
#include "zxhMetricFFD2Base.h"
#include "zxhMetricMIConstraint.h"
#include "zxhImageGipl.h"
#include <math.h>
class zxhMetricMIFFD2 : public zxhMetricMIConstraint,
                        public zxhMetricFFD2Base // interface for FFD2
{
public:
    zxhMetricMIFFD2(void);
    ~zxhMetricMIFFD2(void);

    virtual std::string GetMetricType() {return zxhMetricMIConstraint::GetMetricType()+"_"+zxhMetricFFD2Base::GetMetricType(); } ;

    virtual zxhMetricBase* Clone( zxhMetricBase * & pRet ); 
    virtual int GetValue( zxhlfloat & fvalue ) ;

    virtual int ComputeHistogramFromPDF4D(zxhlfloat * &afPDF,zxhlfloat&fTotalFreq,
                                                float xfTestFrom,float xfTestTo,float yfTestFrom,float yfTestTo,
                                                float zfTestFrom,float zfTestTo,float tfTestFrom,float tfTestTo);
    virtual int ComputeConstanceForFFDGradient() ;
    
    virtual void SetAlwaysUpdateGetValue( bool b )              {m_bAlwaysUpdateGetValue=b;};
    virtual bool GetAlwaysUpdateGetValue( void )                {return m_bAlwaysUpdateGetValue;};

protected:
    virtual void    SetFFDMetricBase()  {m_pMetricBase=this;};
    virtual void    EvaluateOnPixelSource2TargetForce(  zxhlfloat *afForce, int iLengthOfHistogramDim,
                                                        const float * gridCoTest, const float * gridCoRef ) const ;

    /***   methods for ComputeConstanceForFFDGradient    ***/



    virtual bool    EvaluateAllPixelOnTargetForceForMetric();


    bool    m_bAlwaysUpdateGetValue ;
    //virtual bool  ComputeAllControlPointDerivativeBySumLocalEstimation();
//  virtual int GetDiffByFFDGrid(float * afPD, // output partial differentiation
//                                  zxhTransformFFDBase * pFFDGrid,
//                                  int * aiCtrPnt,
//                                  zxhTransformFFDBase * pDifferential )
//  {return this->GetDerivativeFFDGrid(afPD,pFFDGrid,aiCtrPnt,pDifferential); } ;
//  ///
//  virtual int GetDerivativeFFDGrid(float * afPD, // output partial differentiation
//                                   zxhTransformFFDBase * pFFDGrid,
//                                   int * aiCtrPnt,
//                                   zxhTransformFFDBase * pDifferential);
};

#endif