00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017 #ifndef __itkImageFunction_h
00018 #define __itkImageFunction_h
00019
00020 #include "itkFunctionBase.h"
00021 #include "itkPoint.h"
00022 #include "itkIndex.h"
00023 #include "itkContinuousIndex.h"
00024 #include "itkImageBase.h"
00025
00026 namespace itk
00027 {
00028
00029
00055 template <
00056 class TInputImage,
00057 class TOutput,
00058 class TCoordRep = float
00059 >
00060 class ITK_EXPORT ImageFunction :
00061 public FunctionBase< Point<TCoordRep,
00062 ::itk::GetImageDimension<TInputImage>::ImageDimension>,
00063 TOutput >
00064 {
00065 public:
00067 itkStaticConstMacro(ImageDimension, unsigned int,
00068 TInputImage::ImageDimension);
00069
00071 typedef ImageFunction Self;
00072 typedef FunctionBase< Point<TCoordRep,
00073 itkGetStaticConstMacro(ImageDimension)>,
00074 TOutput > Superclass;
00075 typedef SmartPointer<Self> Pointer;
00076 typedef SmartPointer<const Self> ConstPointer;
00077
00079 itkTypeMacro(ImageFunction, FunctionBase);
00080
00082 typedef TInputImage InputImageType;
00083
00085 typedef typename InputImageType::PixelType InputPixelType;
00086
00088 typedef typename InputImageType::ConstPointer InputImageConstPointer;
00089
00091 typedef TOutput OutputType;
00092
00094 typedef TCoordRep CoordRepType;
00095
00097 typedef typename InputImageType::IndexType IndexType;
00098
00100 typedef ContinuousIndex<TCoordRep,itkGetStaticConstMacro(ImageDimension)>
00101 ContinuousIndexType;
00102
00104 typedef Point<TCoordRep,itkGetStaticConstMacro(ImageDimension)> PointType;
00105
00110 virtual void SetInputImage( const InputImageType * ptr );
00111
00113 const InputImageType * GetInputImage() const
00114 { return m_Image.GetPointer(); }
00115
00118 virtual TOutput Evaluate( const PointType& point ) const = 0;
00119
00122 virtual TOutput EvaluateAtIndex( const IndexType & index ) const = 0;
00123
00126 virtual TOutput EvaluateAtContinuousIndex(
00127 const ContinuousIndexType & index ) const = 0;
00128
00132 virtual bool IsInsideBuffer( const IndexType & index ) const
00133 {
00134 for ( unsigned int j = 0; j < ImageDimension; j++ )
00135 {
00136 if ( index[j] < m_StartIndex[j] ) { return false; };
00137 if ( index[j] > m_EndIndex[j] ) { return false; };
00138 }
00139 return true;
00140 }
00142
00146 virtual bool IsInsideBuffer( const ContinuousIndexType & index ) const
00147 {
00148 for ( unsigned int j = 0; j < ImageDimension; j++ )
00149 {
00150 if ( index[j] < m_StartContinuousIndex[j] ) { return false; };
00151 if ( index[j] > m_EndContinuousIndex[j] ) { return false; };
00152 }
00153 return true;
00154 }
00156
00160 virtual bool IsInsideBuffer( const PointType & point ) const
00161 {
00162 ContinuousIndexType index;
00163 m_Image->TransformPhysicalPointToContinuousIndex( point, index );
00164 return this->IsInsideBuffer( index );
00165 }
00167
00169 void ConvertPointToNearestIndex( const PointType & point,
00170 IndexType & index ) const
00171 {
00172 ContinuousIndexType cindex;
00173 m_Image->TransformPhysicalPointToContinuousIndex( point, cindex );
00174 this->ConvertContinuousIndexToNearestIndex( cindex, index );
00175 }
00177
00179 void ConvertPointToContinuousIndex( const PointType & point,
00180 ContinuousIndexType & cindex ) const
00181 {
00182 m_Image->TransformPhysicalPointToContinuousIndex( point, cindex );
00183 }
00184
00186 inline void ConvertContinuousIndexToNearestIndex( const ContinuousIndexType & cindex,
00187 IndexType & index ) const
00188 {
00189 index.CopyWithRound( cindex );
00190 }
00191
00192 itkGetConstReferenceMacro(StartIndex, IndexType);
00193 itkGetConstReferenceMacro(EndIndex, IndexType);
00194
00195 itkGetConstReferenceMacro(StartContinuousIndex, ContinuousIndexType);
00196 itkGetConstReferenceMacro(EndContinuousIndex, ContinuousIndexType);
00197
00198 protected:
00199 ImageFunction();
00200 ~ImageFunction() {}
00201 void PrintSelf(std::ostream& os, Indent indent) const;
00202
00204 InputImageConstPointer m_Image;
00205
00207 IndexType m_StartIndex;
00208 IndexType m_EndIndex;
00209 ContinuousIndexType m_StartContinuousIndex;
00210 ContinuousIndexType m_EndContinuousIndex;
00211
00212 private:
00213 ImageFunction(const Self&);
00214 void operator=(const Self&);
00215
00216 };
00217
00218 }
00219
00220
00221 #define ITK_TEMPLATE_ImageFunction(_, EXPORT, x, y) namespace itk { \
00222 _(3(class EXPORT ImageFunction< ITK_TEMPLATE_3 x >)) \
00223 namespace Templates { typedef ImageFunction< ITK_TEMPLATE_3 x > ImageFunction##y; } \
00224 }
00225
00226 #if ITK_TEMPLATE_EXPLICIT
00227 # include "Templates/itkImageFunction+-.h"
00228 #endif
00229
00230 #if ITK_TEMPLATE_TXX
00231 # include "itkImageFunction.txx"
00232 #endif
00233
00234 #endif
00235