cmtkImagePairNonrigidRegistration.h

Go to the documentation of this file.

/*
//
//  Copyright 2004-2010 SRI International
//
//  Copyright 1997-2009 Torsten Rohlfing
//
//  This file is part of the Computational Morphometry Toolkit.
//
//  http://www.nitrc.org/projects/cmtk/
//
//  The Computational Morphometry Toolkit is free software: you can
//  redistribute it and/or modify it under the terms of the GNU General Public
//  License as published by the Free Software Foundation, either version 3 of
//  the License, or (at your option) any later version.
//
//  The Computational Morphometry Toolkit is distributed in the hope that it
//  will be useful, but WITHOUT ANY WARRANTY; without even the implied
//  warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
//  GNU General Public License for more details.
//
//  You should have received a copy of the GNU General Public License along
//  with the Computational Morphometry Toolkit.  If not, see
//  <http://www.gnu.org/licenses/>.
//
//  $Revision: 2752 $
//
//  $LastChangedDate: 2011-01-17 11:33:31 -0800 (Mon, 17 Jan 2011) $
//
//  $LastChangedBy: torstenrohlfing $
//
*/

#ifndef __cmtkImagePairNonrigidRegistration_h_included_
#define __cmtkImagePairNonrigidRegistration_h_included_

#include <cmtkconfig.h>

#include <Base/cmtkSplineWarpXform.h>
#include <Registration/cmtkImagePairRegistration.h>

#include <string.h>

namespace
cmtk
{


class ImagePairNonrigidRegistration : 
  public ImagePairRegistration 
{
public:
  typedef ImagePairNonrigidRegistration Self;

  typedef ImagePairRegistration Superclass;

protected:
  SplineWarpXform::SmartPtr InitialWarpXform;

  SplineWarpXform::SmartPtr InverseWarpXform;

  cmtkGetSetMacro(bool,MatchFltToRefHistogram);

  cmtkGetSetMacroDefault(bool,RepeatMatchFltToRefHistogram,true);

  cmtkGetSetMacro(int,RefineGrid);

  cmtkGetSetMacro(bool,DelayRefineGrid);

  cmtkGetSetMacro(Types::Coordinate,GridSpacing);

  cmtkGetSetMacro(bool,ExactGridSpacing);

  unsigned int IgnoreEdge;

  const char* RestrictToAxes;

  cmtkGetSetMacro(bool,FastMode);

  cmtkGetSetMacro(bool,AdaptiveFixParameters);

  cmtkGetSetMacro(float,AdaptiveFixThreshFactor);

  cmtkGetSetMacro(float,JacobianConstraintWeight);

  cmtkGetSetMacro(float,GridEnergyWeight);

  cmtkGetSetMacro(float,RelaxWeight);

  cmtkGetSetMacro(float,InverseConsistencyWeight);

  cmtkGetSetMacro(float,LandmarkErrorWeight);

  bool m_RelaxToUnfold;

  ImagePairNonrigidRegistration();

  virtual ~ImagePairNonrigidRegistration() {};

  virtual CallbackResult InitRegistration ();

  virtual void EnterResolution( CoordinateVector::SmartPtr&, Functional::SmartPtr&, const int, const int );

  virtual int DoneResolution( CoordinateVector::SmartPtr&, Functional::SmartPtr&, const int, const int );

  SplineWarpXform::SmartPtr GetTransformation() const
  {
    return SplineWarpXform::SmartPtr::DynamicCastFrom( this->m_Xform );
  }

  const UniformVolume::SmartPtr GetReformattedFloatingImage( Interpolators::InterpolationEnum interpolator = Interpolators::LINEAR ) const;

private:
  MatchedLandmarkList::SmartPtr m_MatchedLandmarks;

  int RefinedGridAtLevel;

  int RefineGridCount;

  bool RelaxationStep;

  bool RefineDelayed;

  SplineWarpXform::SmartPtr MakeWarpXform( const UniformVolume::CoordinateVectorType& size, const AffineXform* initialAffine ) const;

  class LevelParameters
    : public Superclass::LevelParameters
  {
  public:
    LevelParameters( const Types::Coordinate resolution ) : m_Resolution( resolution ) {}

    Types::Coordinate m_Resolution;
  };

  virtual Functional* MakeFunctional( const int level, const Superclass::LevelParameters* levelParameters );
};


} // namespace cmtk

#endif // __cmtkImagePairNonrigidRegistration_h_included_