cmtkElasticRegistration.h

Go to the documentation of this file.

/*
//
//  Copyright 1997-2009 Torsten Rohlfing
//
//  Copyright 2004-2010 SRI International
//
//  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 __cmtkElasticRegistration_h_included_
#define __cmtkElasticRegistration_h_included_

#include <cmtkconfig.h>

#include <Registration/cmtkVoxelRegistration.h>
#include <Base/cmtkSplineWarpXform.h>
#include <Base/cmtkMatchedLandmarkList.h>

#include <string.h>

namespace
cmtk
{


class ElasticRegistration : 
  public VoxelRegistration 
{
protected:
  SplineWarpXform::SmartPtr InitialWarpXform;

  SplineWarpXform::SmartPtr InverseWarpXform;

  bool ForceSwitchVolumes;
  
  cmtkGetSetMacro(bool,MatchFltToRefHistogram);

  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,RigidityConstraintWeight);

  cmtkGetSetMacro(UniformVolume::SmartPtr,RigidityConstraintMap);

  cmtkGetSetMacro(float,GridEnergyWeight);

  cmtkGetSetMacro(float,RelaxWeight);

  cmtkGetSetMacro(float,InverseConsistencyWeight);

  cmtkGetSetMacro(float,LandmarkErrorWeight);

  cmtkGetSetMacro(MatchedLandmarkList::SmartPtr,LandmarkList);

  bool m_RelaxToUnfold;

  virtual void SetForceOutside( const bool flag = true, const Types::DataItem value = 0 )
  {
    this->m_ForceOutsideFlag = flag;
    this->m_ForceOutsideValue = value;
  }

  ElasticRegistration();

  virtual ~ElasticRegistration() {};

  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 );

private:
  int RefinedGridAtLevel;

  int RefineGridCount;

  typedef VoxelRegistration Superclass;

  bool RelaxationStep;

  bool RefineDelayed;

  bool m_ForceOutsideFlag;

  Types::DataItem m_ForceOutsideValue;

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

  Functional* MakeFunctional( UniformVolume::SmartPtr& refVolume, UniformVolume::SmartPtr& fltVolume, UniformVolume::SmartPtr& rigidityMap ) const;
};


} // namespace cmtk

#endif // __cmtkElasticRegistration_h_included_