zxhMetricPhase.h

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

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

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

#define PHASEBINN 68









class zxhMetricPhase :  public zxhMetricBase
{
public:
    zxhMetricPhase(void);
    virtual ~zxhMetricPhase(void);

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

    virtual zxhMetricBase* Clone( zxhMetricBase * & pRet ); 

    /* * * * * * * * * * * *
     *  for phase registration V. Grau et al.
     * * * * * * * * * * * */
    virtual int ComputePhaseMetric4D(zxhlfloat &fValue, 
                        float xfTestFrom,float xfTestTo,float yfTestFrom,float yfTestTo,
                        float zfTestFrom=0,float zfTestTo=0,float tfTestFrom=0,float tfTestTo=0, bool nonrigid=false ) ;
    virtual int  GetValue( zxhlfloat & fValue ) ;

    virtual void UpdatePhaseProcess( int processimage =0) ;

    virtual void SetPhaseScale( float s )                                       {m_procPhaseTest.SetScale(s);m_procPhaseRef.SetScale(s);};
    virtual zxhImageProcessPhaseT<short, float> * GetProcPhaseTest()            { return &m_procPhaseTest ; } ;
    virtual zxhImageProcessPhaseT<short, float> * GetProcPhaseRef()             { return &m_procPhaseRef ; } ;
    virtual void SetSimilarityWeighting(float wi,float wf, float wo)        
    { m_fPhaseWeightings[0]=wi; m_fPhaseWeightings[1]=wf; m_fPhaseWeightings[2]=wo;}
    virtual void SetIsEchoDataUsingMask(bool b)                                 {m_bIsEchoDataUsingMask=b;};

protected:
    virtual void EvaluateSimilarityOnPixel( float *coTest, float *coRef,
                                            float *OrientTest2Ref, float *OrientRef,
                                            zxhlfloat & SumGrayAB, zxhlfloat & SumGrayA2, zxhlfloat & SumGrayB2, zxhlfloat & SumGrayA, zxhlfloat & SumGrayB,
                                            zxhlfloat & SumPhaseAB,zxhlfloat & SumPhaseA2,zxhlfloat & SumPhaseB2,zxhlfloat & SumPhaseA,zxhlfloat & SumPhaseB,
                                            zxhlfloat & SumOrientCos, zxhlfloat & SumPhaseCos, zxhlfloat& Wrt  ) ;
    virtual void EvaluateSimilarityOnPixel( float *coTest, float *coRef,
                                            float *OrientTest2Ref, float *OrientRef,
                                            zxhlfloat & SumGrayAB, zxhlfloat & SumGrayA2, zxhlfloat & SumGrayB2, zxhlfloat & SumGrayA, zxhlfloat & SumGrayB,
                                            zxhlfloat & SumPhaseAB,zxhlfloat & SumPhaseA2,zxhlfloat & SumPhaseB2,zxhlfloat & SumPhaseA,zxhlfloat & SumPhaseB,
                                            zxhlfloat & SumOrientCos, zxhlfloat & SumPhaseCos,
                                            zxhlfloat *pfJointHistogramGray, zxhlfloat *pfJointHistogramPhase, zxhlfloat& Wrt  ) ;

    virtual int GetValueFromIntensityPhase( zxhlfloat & fValue, zxhlfloat SumOrientCos, zxhlfloat SumPhaseCos, zxhlfloat Wrt,
                                                zxhlfloat SumGrayAB, zxhlfloat SumGrayA2, zxhlfloat SumGrayB2, zxhlfloat SumGrayA, zxhlfloat SumGrayB,
                                                zxhlfloat SumPhaseAB,zxhlfloat SumPhaseA2,zxhlfloat SumPhaseB2,zxhlfloat SumPhaseA,zxhlfloat SumPhaseB );
    virtual int GetValueFromIntensityPhase( zxhlfloat & fValue, zxhlfloat SumOrientCos, zxhlfloat SumPhaseCos, zxhlfloat Wrt,  
                                            zxhlfloat SumGrayAB, zxhlfloat SumGrayA2, zxhlfloat SumGrayB2, zxhlfloat SumGrayA, zxhlfloat SumGrayB,
                                            zxhlfloat SumPhaseAB,zxhlfloat SumPhaseA2,zxhlfloat SumPhaseB2,zxhlfloat SumPhaseA,zxhlfloat SumPhaseB,
                                            zxhlfloat *pfJointHistogramGray, zxhlfloat *pfJointHistogramPhase);
    zxhImageProcessPhaseT<short, float> m_procPhaseRef ;
    zxhImageProcessPhaseT<short, float> m_procPhaseTest ;
    zxhImageDataT<float>                m_imgFeatureIntensityTest;
    zxhImageDataT<float>                m_imgFeatureIntensityRef;

    int m_iNumberOfHistogramBins;
    float m_afImageTestGraySpan[2];
    float m_afImageRefGraySpan[2];
    float m_fBinSizeImageIntensity[2] ;
    float m_afImageTestPhaseSpan[2];
    float m_afImageRefPhaseSpan[2];
    float m_fBinSizeImagePhase[2] ;
    zxhlfloat m_afCJointHistogramIntensity[PHASEBINN][PHASEBINN] ;
    zxhlfloat m_afCJointHistogramPhase[PHASEBINN][PHASEBINN] ;

    float                               m_fPhaseWeightings[3];
    bool                                m_bIsEchoDataUsingMask;

};

#endif