Main Page   Groups   Namespace List   Class Hierarchy   Alphabetical List   Compound List   File List   Namespace Members   Compound Members   File Members   Concepts

itkKernelTransform.h

Go to the documentation of this file.
00001 /*=========================================================================
00002 
00003   Program:   Insight Segmentation & Registration Toolkit
00004   Module:    $RCSfile: itkKernelTransform.h,v $
00005   Language:  C++
00006   Date:      $Date: 2005/01/22 02:47:25 $
00007   Version:   $Revision: 1.37 $
00008 
00009   Copyright (c) Insight Software Consortium. All rights reserved.
00010   See ITKCopyright.txt or http://www.itk.org/HTML/Copyright.htm for details.
00011 
00012      This software is distributed WITHOUT ANY WARRANTY; without even 
00013      the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR 
00014      PURPOSE.  See the above copyright notices for more information.
00015 
00016 =========================================================================*/
00017 #ifndef __itkKernelTransform_h
00018 #define __itkKernelTransform_h
00019 
00020 #include "itkTransform.h"
00021 #include "itkPoint.h"
00022 #include "itkVector.h"
00023 #include "itkMatrix.h"
00024 #include "itkPointSet.h"
00025 #include <deque>
00026 #include <math.h>
00027 #include "vnl/vnl_matrix_fixed.h"
00028 #include "vnl/vnl_matrix.h"
00029 #include "vnl/vnl_vector.h"
00030 #include "vnl/vnl_vector_fixed.h"
00031 #include "vnl/algo/vnl_svd.h"
00032 #include "vnl/vnl_sample.h"
00033 
00034 namespace itk
00035 {
00036 
00062 template <class TScalarType, // probably only float and double make sense here
00063           unsigned int NDimensions>   // Number of dimensions
00064 class ITK_EXPORT KernelTransform : 
00065                     public Transform<TScalarType, NDimensions,NDimensions>
00066 {
00067 public:
00069   typedef KernelTransform Self;
00070   typedef Transform<TScalarType, NDimensions, NDimensions >   Superclass;
00071   typedef SmartPointer<Self>        Pointer;
00072   typedef SmartPointer<const Self>  ConstPointer;
00073   
00075   itkTypeMacro( KernelTransform, Transform );
00076 
00078   itkNewMacro( Self );
00079 
00081   itkStaticConstMacro(SpaceDimension, unsigned int, NDimensions);
00082 
00084   typedef typename Superclass::ScalarType  ScalarType;
00085 
00087   typedef typename Superclass::ParametersType  ParametersType;
00088 
00090   typedef typename Superclass::JacobianType  JacobianType;
00091 
00093   typedef typename Superclass::InputPointType   InputPointType;
00094   typedef typename Superclass::OutputPointType  OutputPointType;
00095   
00097   typedef typename Superclass::InputVectorType   InputVectorType;
00098   typedef typename Superclass::OutputVectorType  OutputVectorType;
00099   
00102   typedef DefaultStaticMeshTraits<TScalarType,
00103                                   NDimensions,
00104                                   NDimensions,
00105                                   TScalarType,
00106                                   TScalarType> PointSetTraitsType;
00107   typedef PointSet<InputPointType, NDimensions, PointSetTraitsType> PointSetType;
00108   typedef typename PointSetType::Pointer                        PointSetPointer;
00109   typedef typename PointSetType::PointsContainer                PointsContainer;
00110   typedef typename PointSetType::PointsContainerIterator        PointsIterator;
00111   typedef typename PointSetType::PointsContainerConstIterator   PointsConstIterator;
00112     
00114   typedef itk::VectorContainer<unsigned long,InputVectorType> VectorSetType;
00115   typedef typename VectorSetType::Pointer        VectorSetPointer;
00116   
00118   itkGetObjectMacro( SourceLandmarks, PointSetType);
00119   
00121   virtual void SetSourceLandmarks(PointSetType *);
00122   
00124   itkGetObjectMacro( TargetLandmarks, PointSetType);
00125   
00127   virtual void SetTargetLandmarks(PointSetType *);
00128   
00131   itkGetObjectMacro( Displacements, VectorSetType );
00132   
00134   void ComputeWMatrix(void);
00135   
00137   virtual OutputPointType TransformPoint(const InputPointType& thisPoint) const;
00138   
00140   typedef vnl_matrix_fixed<TScalarType, NDimensions, NDimensions> IMatrixType;
00141 
00142 
00144   virtual const JacobianType & GetJacobian(const InputPointType  &point ) const;
00145 
00150   virtual void SetParameters(const ParametersType &);
00151 
00153   virtual void UpdateParameters(void) const;
00154 
00156   virtual const ParametersType& GetParameters(void) const;
00157 
00168   itkSetClampMacro(Stiffness, double, 0.0, NumericTraits<double>::max());
00169   itkGetMacro(Stiffness, double);
00170 
00171   /* Supply a Reinitialize seed API */
00172   static void ReinitializeSeed(int seed);
00173   static void ReinitializeSeed();
00174   
00175 
00176 protected:
00177   KernelTransform();
00178   virtual ~KernelTransform();
00179   void PrintSelf(std::ostream& os, Indent indent) const;
00180 
00181 public:  
00183   typedef vnl_matrix_fixed<TScalarType, NDimensions, NDimensions> GMatrixType;
00184   
00186   typedef vnl_matrix<TScalarType> LMatrixType;
00187   
00189   typedef vnl_matrix<TScalarType> KMatrixType;
00190   
00192   typedef vnl_matrix<TScalarType> PMatrixType;
00193   
00195   typedef vnl_matrix<TScalarType> YMatrixType;
00196 
00198   typedef vnl_matrix<TScalarType> WMatrixType;
00199   
00201   typedef vnl_matrix<TScalarType> DMatrixType;
00202   
00204   typedef vnl_matrix_fixed<TScalarType,NDimensions,NDimensions> AMatrixType;
00205   
00207   typedef vnl_vector_fixed<TScalarType,NDimensions> BMatrixType;
00208   
00210   typedef vnl_matrix_fixed<TScalarType, 1, NDimensions> RowMatrixType;
00211   
00213   typedef vnl_matrix_fixed<TScalarType, NDimensions, 1> ColumnMatrixType;
00214 
00216   PointSetPointer m_SourceLandmarks;
00217   
00219   PointSetPointer m_TargetLandmarks;
00220   
00221 protected:
00228   virtual const GMatrixType & ComputeG(const InputVectorType & landmarkVector) const;
00229 
00236   virtual const GMatrixType & ComputeReflexiveG(PointsIterator) const;
00237 
00238   
00241   virtual void ComputeDeformationContribution( const InputPointType & inputPoint,
00242                                                      OutputPointType & result ) const;
00243 
00245   void ComputeK();
00246   
00248   void ComputeL();
00249   
00251   void ComputeP();
00252   
00254   void ComputeY();
00255   
00257   void ComputeD();
00258 
00263   void ReorganizeW(void);
00264 
00266   double m_Stiffness;
00267 
00270   VectorSetPointer m_Displacements;
00271 
00273   LMatrixType m_LMatrix;
00274 
00276   KMatrixType m_KMatrix;
00277 
00279   PMatrixType m_PMatrix;
00280 
00282   YMatrixType m_YMatrix;
00283   
00285   WMatrixType m_WMatrix;
00286 
00292   DMatrixType m_DMatrix;
00293 
00295   AMatrixType m_AMatrix;
00296 
00298   BMatrixType m_BVector;
00299 
00303   mutable GMatrixType m_GMatrix;
00304 
00306   bool m_WMatrixComputed;
00307 
00309   IMatrixType m_I;
00310 
00311  private:
00312   KernelTransform(const Self&); //purposely not implemented
00313   void operator=(const Self&); //purposely not implemented
00314 
00315 };
00316 
00317 } // end namespace itk
00318 
00319 #ifndef ITK_MANUAL_INSTANTIATION
00320 #include "itkKernelTransform.txx"
00321 #endif
00322 
00323 #endif // __itkKernelTransform_h

Generated at Wed Mar 30 00:00:35 2005 for ITK by doxygen 1.3.9.1 written by Dimitri van Heesch, © 1997-2000