00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017 #ifndef __itkDeformationFieldJacobianDeterminantFilter_h
00018 #define __itkDeformationFieldJacobianDeterminantFilter_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
00029
00030 #include "itkDisplacementFieldJacobianDeterminantFilter.h"
00031
00032 namespace itk
00033 {
00073 template < typename TInputImage,
00074 typename TRealType = float,
00075 typename TOutputImage = Image< TRealType,
00076 ::itk::GetImageDimension<TInputImage>::ImageDimension >
00077 >
00078 class ITK_EXPORT DeformationFieldJacobianDeterminantFilter :
00079 public DisplacementFieldJacobianDeterminantFilter< TInputImage,TRealType,TOutputImage>
00080 {
00081 public:
00083 typedef DeformationFieldJacobianDeterminantFilter Self;
00084 typedef DisplacementFieldJacobianDeterminantFilter<TInputImage,TRealType,TOutputImage> Superclass;
00085 typedef SmartPointer<Self> Pointer;
00086 typedef SmartPointer<const Self> ConstPointer;
00087
00089 itkNewMacro(Self);
00090
00092 itkTypeMacro(DeformationFieldJacobianDeterminantFilter, DisplacementFieldJacobianDeterminantFilter);
00093
00096 typedef typename TOutputImage::PixelType OutputPixelType;
00097 typedef typename TInputImage::PixelType InputPixelType;
00098
00100 typedef TInputImage InputImageType;
00101 typedef TOutputImage OutputImageType;
00102 typedef typename InputImageType::Pointer InputImagePointer;
00103 typedef typename OutputImageType::Pointer OutputImagePointer;
00104
00106 itkStaticConstMacro(ImageDimension, unsigned int,
00107 TOutputImage::ImageDimension);
00108
00110 itkStaticConstMacro(VectorDimension, unsigned int,
00111 InputPixelType::Dimension);
00112
00114 typedef TRealType RealType;
00115 typedef Vector<TRealType, ::itk::GetVectorDimension<InputPixelType>::VectorDimension> RealVectorType;
00116 typedef Image<RealVectorType, ::itk::GetImageDimension<TInputImage>::ImageDimension> RealVectorImageType;
00117
00118
00121 typedef ConstNeighborhoodIterator<RealVectorImageType> ConstNeighborhoodIteratorType;
00122 typedef typename ConstNeighborhoodIteratorType::RadiusType RadiusType;
00123
00125 typedef typename Superclass::OutputImageRegionType OutputImageRegionType;
00126
00127 void PrintSelf(std::ostream& os, Indent indent) const;
00128 virtual TRealType EvaluateAtNeighborhood(const ConstNeighborhoodIteratorType &it) const;
00129 protected:
00130 DeformationFieldJacobianDeterminantFilter();
00131 virtual ~DeformationFieldJacobianDeterminantFilter() {}
00132
00133 private:
00134 };
00135
00136 }
00137
00138 #ifndef ITK_MANUAL_INSTANTIATION
00139 #include "itkDeformationFieldJacobianDeterminantFilter.txx"
00140 #endif
00141
00142 #endif
00143