00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017 #ifndef __itkFastMarchingExtensionImageFilter_h
00018 #define __itkFastMarchingExtensionImageFilter_h
00019
00020 #include "itkFastMarchingImageFilter.h"
00021 #include "itkImage.h"
00022 #include "itkLevelSet.h"
00023
00024 namespace itk
00025 {
00026
00052 template <
00053 class TLevelSet,
00054 class TAuxValue,
00055 unsigned int VAuxDimension = 1,
00056 class TSpeedImage = Image<float,::itk::GetImageDimension<TLevelSet>::ImageDimension>
00057 >
00058 class ITK_EXPORT FastMarchingExtensionImageFilter :
00059 public FastMarchingImageFilter<TLevelSet,TSpeedImage>
00060 {
00061 public:
00063 typedef FastMarchingExtensionImageFilter Self;
00064 typedef FastMarchingImageFilter<TLevelSet,TSpeedImage> Superclass;
00065 typedef SmartPointer<Self> Pointer;
00066 typedef SmartPointer<const Self> ConstPointer;
00067
00069 itkNewMacro(Self);
00070
00072 itkTypeMacro(FastMarchingExtensionImageFilter, FastMarchingImageFilter);
00073
00075 typedef typename Superclass::LevelSetType LevelSetType;
00076 typedef typename Superclass::SpeedImageType SpeedImageType;
00077 typedef typename Superclass::LevelSetImageType LevelSetImageType;
00078
00080 itkStaticConstMacro(SetDimension, unsigned int,Superclass::SetDimension);
00081
00083 itkStaticConstMacro(AuxDimension, unsigned int,VAuxDimension);
00084
00086 typedef AuxVarTypeDefault<TAuxValue,
00087 itkGetStaticConstMacro(AuxDimension),
00088 itkGetStaticConstMacro(SetDimension)>
00089 AuxVarType;
00090 typedef typename AuxVarType::AuxValueType AuxValueType;
00091 typedef typename AuxVarType::AuxValueVectorType AuxValueVectorType;
00092 typedef typename AuxVarType::AuxValueContainer AuxValueContainer;
00093 typedef typename AuxVarType::AuxImageType AuxImageType;
00094 typedef typename AuxVarType::AuxImagePointer AuxImagePointer;
00096
00098 typedef Index<itkGetStaticConstMacro(SetDimension)> IndexType;
00099
00101 AuxImageType * GetAuxiliaryImage( unsigned int idx );
00102
00104 void SetAuxiliaryAliveValues( AuxValueContainer * values )
00105 {
00106 m_AuxAliveValues = values;
00107 }
00108
00110 AuxValueContainer * GetAuxiliaryAliveValues(void)
00111 {
00112 return m_AuxAliveValues.GetPointer();
00113 }
00114
00116 void SetAuxiliaryTrialValues( AuxValueContainer * values )
00117 {
00118 m_AuxTrialValues = values;
00119 }
00120
00122 typename AuxValueContainer::Pointer GetAuxiliaryTrialValues()
00123 {
00124 return m_AuxTrialValues;
00125 }
00126
00127 #ifdef ITK_USE_CONCEPT_CHECKING
00128
00129 itkConceptMacro(AuxValueHasNumericTraitsCheck,
00130 (Concept::HasNumericTraits<TAuxValue>));
00131
00133 #endif
00134
00135 protected:
00136 FastMarchingExtensionImageFilter();
00137 ~FastMarchingExtensionImageFilter(){};
00138 void PrintSelf( std::ostream& os, Indent indent ) const;
00139
00140 virtual void Initialize( LevelSetImageType * );
00141 virtual double UpdateValue( const IndexType & index,
00142 const SpeedImageType * speed, LevelSetImageType * output);
00143
00145 virtual void GenerateOutputInformation();
00146 virtual void EnlargeOutputRequestedRegion( DataObject *output );
00148
00149 private:
00150 FastMarchingExtensionImageFilter(const Self&);
00151 void operator=(const Self&);
00152
00153 typename AuxValueContainer::Pointer m_AuxAliveValues;
00154 typename AuxValueContainer::Pointer m_AuxTrialValues;
00155
00156 };
00157
00158 }
00159
00160
00161 #ifndef ITK_MANUAL_INSTANTIATION
00162 #include "itkFastMarchingExtensionImageFilter.txx"
00163 #endif
00164
00165 #endif
00166