zxhMetricPhaseFFD.h

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

  Program:   ZXH Registration Software
  Author:    Xiahai Zhuang
  Module:    $RCSfle: zxhMetricPhaseFFD.h    $
  Language:  C++
  Date:      $Date: From  2009-02 $
  Version:   $Revision: 2.0 $
  revision:
=========================================================================*/
#ifndef zxhMetricPhaseFFD_h
#define zxhMetricPhaseFFD_h

#include "zxh.h"
#include "zxhMetricBase.h"
#include "zxhMetricPhase.h"
#include "zxhMetricFFDBase.h"
#include "zxhImageProcessPhase.h"
#include <math.h> 


class zxhMetricPhaseFFD :   public zxhMetricFFDBase, public zxhMetricPhase
{
public:
    zxhMetricPhaseFFD(void);
    virtual ~zxhMetricPhaseFFD(void);

    virtual std::string GetMetricType() {return zxhMetricPhase::GetMetricType()+"_"+zxhMetricFFDBase::GetMetricType(); } ;

    virtual zxhMetricBase* Clone( zxhMetricBase * & pRet );

    virtual int ComputeConstanceForFFDGradient() ;

    virtual int GetValueByFFDGrid(zxhlfloat&fValue, int aiCtrPnt[ImageDimensionMax]);

    //GetDiffByFFDGrid
    //GetDiffByFFDGridOnFFDPhysicalVector

    virtual int  GetValue( zxhlfloat & fValue )
    {
        //zxhTransformBase * pTransformOrig = m_pTransform ;
        //m_pTransform = 0 ; This is commented because the line-search optimization
        int iret = zxhMetricPhase::GetValue( fValue ) ;
        if( m_fWeightFFDBendingEnergy > 0 )
        {
            zxhlfloat bending = this->GetBendingEnergy( (zxhTransformFFDBase*)m_pTransform )  ;
            if( glbVerboseOutput>0) std::cout<<"success: \t\t bending energy "<<bending<<"\n" ;
                fValue = fValue - m_fWeightFFDBendingEnergy*bending; //minimize bending energy
        }
        return iret ;
    };

protected:
    virtual void    SetFFDMetricBase()  {m_pMetricBase=this;};
};

#endif