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

itkDisplacementFieldJacobianDeterminantFilter.h

Go to the documentation of this file.
00001 /*=========================================================================
00002 
00003   Program:   Insight Segmentation & Registration Toolkit
00004   Module:    $RCSfile: itkDisplacementFieldJacobianDeterminantFilter.h,v $
00005   Language:  C++
00006   Date:      $Date: 2008-10-14 19:20:33 $
00007   Version:   $Revision: 1.2 $
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 __itkDisplacementFieldJacobianDeterminantFilter_h
00018 #define __itkDisplacementFieldJacobianDeterminantFilter_h
00019 
00020 #include "itkConstNeighborhoodIterator.h"
00021 #include "itkNeighborhoodIterator.h"
00022 #include "itkImageToImageFilter.h"
00023 #include "itkImage.h"
00024 #include "itkVector.h"
00025 #include "vnl/vnl_matrix.h"
00026 #include "vnl/vnl_det.h"
00027 
00028 namespace itk
00029 {
00111 template < typename TInputImage,
00112            typename TRealType = float,
00113            typename TOutputImage = Image< TRealType,
00114                                           ::itk::GetImageDimension<TInputImage>::ImageDimension >
00115 >
00116 class ITK_EXPORT DisplacementFieldJacobianDeterminantFilter :
00117     public ImageToImageFilter< TInputImage, TOutputImage >
00118 {
00119 public:
00121   typedef DisplacementFieldJacobianDeterminantFilter      Self;
00122   typedef ImageToImageFilter< TInputImage, TOutputImage > Superclass;
00123   typedef SmartPointer<Self>                              Pointer;
00124   typedef SmartPointer<const Self>                        ConstPointer;
00125 
00127   itkNewMacro(Self);
00128 
00130   itkTypeMacro(DisplacementFieldJacobianDeterminantFilter, ImageToImageFilter);
00131 
00134   typedef typename TOutputImage::PixelType OutputPixelType;
00135   typedef typename TInputImage::PixelType  InputPixelType;
00136 
00138   typedef TInputImage                       InputImageType;
00139   typedef TOutputImage                      OutputImageType;
00140   typedef typename InputImageType::Pointer  InputImagePointer;
00141   typedef typename OutputImageType::Pointer OutputImagePointer;
00142 
00144   itkStaticConstMacro(ImageDimension, unsigned int,
00145                       TOutputImage::ImageDimension);
00146 
00148   itkStaticConstMacro(VectorDimension, unsigned int,
00149                       InputPixelType::Dimension);
00150 
00152   typedef TRealType RealType;
00153   typedef Vector<TRealType, ::itk::GetVectorDimension<InputPixelType>::VectorDimension>
00154                     RealVectorType;
00155   typedef Image<RealVectorType, ::itk::GetImageDimension<TInputImage>::ImageDimension>
00156                     RealVectorImageType;
00157 
00158 
00161   typedef ConstNeighborhoodIterator<RealVectorImageType> ConstNeighborhoodIteratorType;
00162   typedef typename ConstNeighborhoodIteratorType::RadiusType RadiusType;
00163 
00165   typedef typename Superclass::OutputImageRegionType OutputImageRegionType;
00166 
00175   virtual void GenerateInputRequestedRegion() throw(InvalidRequestedRegionError);
00176 
00180   void SetUseImageSpacingOn()
00181     { this->SetUseImageSpacing(true); }
00182 
00186   void SetUseImageSpacingOff()
00187     { this->SetUseImageSpacing(false); }
00188 
00191   void SetUseImageSpacing(bool);
00192   itkGetMacro(UseImageSpacing, bool);
00194 
00197   void SetDerivativeWeights(TRealType data[]);
00198   itkGetVectorMacro(DerivativeWeights, const TRealType, itk::GetImageDimension<TInputImage>::ImageDimension);
00200 
00201 protected:
00202   DisplacementFieldJacobianDeterminantFilter();
00203   virtual ~DisplacementFieldJacobianDeterminantFilter() {}
00204 
00208   void BeforeThreadedGenerateData ();
00209 
00222   void ThreadedGenerateData(const OutputImageRegionType& outputRegionForThread,
00223                             int threadId );
00224 
00225   void PrintSelf(std::ostream& os, Indent indent) const;
00226 
00227   typedef typename InputImageType::Superclass ImageBaseType;
00228 
00230   itkGetConstObjectMacro( RealValuedInputImage, ImageBaseType );
00231 
00233   itkGetConstReferenceMacro( NeighborhoodRadius, RadiusType );
00234   itkSetMacro( NeighborhoodRadius, RadiusType );
00236 
00237 
00238   virtual TRealType EvaluateAtNeighborhood(const ConstNeighborhoodIteratorType &it) const;
00239 
00241   TRealType m_DerivativeWeights[itk::GetImageDimension<TInputImage>::ImageDimension];
00242 
00244   TRealType m_HalfDerivativeWeights[itk::GetImageDimension<TInputImage>::ImageDimension];
00245 
00246 private:
00247   bool m_UseImageSpacing;
00248   int  m_RequestedNumberOfThreads;
00249 
00250   typename ImageBaseType::ConstPointer m_RealValuedInputImage;
00251 
00252   DisplacementFieldJacobianDeterminantFilter(const Self&); //purposely not implemented
00253   void operator=(const Self&); //purposely not implemented
00254 
00255   RadiusType    m_NeighborhoodRadius;
00256 };
00257 
00258 } // end namespace itk
00259 
00260 #ifndef ITK_MANUAL_INSTANTIATION
00261 #include "itkDisplacementFieldJacobianDeterminantFilter.txx"
00262 #endif
00263 
00264 #endif
00265 

Generated at Wed Nov 5 21:11:25 2008 for ITK by doxygen 1.5.1 written by Dimitri van Heesch, © 1997-2000