00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017 #ifndef __itkWarpVectorImageFilter_h
00018 #define __itkWarpVectorImageFilter_h
00019
00020 #include "itkImageToImageFilter.h"
00021 #include "itkVectorInterpolateImageFunction.h"
00022 #include "itkVectorLinearInterpolateImageFunction.h"
00023 #include "itkPoint.h"
00024 #include "itkFixedArray.h"
00025
00026 namespace itk
00027 {
00028
00082 template <
00083 class TInputImage,
00084 class TOutputImage,
00085 class TDeformationField
00086 >
00087 class ITK_EXPORT WarpVectorImageFilter :
00088 public ImageToImageFilter<TInputImage, TOutputImage>
00089 {
00090 public:
00092 typedef WarpVectorImageFilter Self;
00093 typedef ImageToImageFilter<TInputImage,TOutputImage> Superclass;
00094 typedef SmartPointer<Self> Pointer;
00095 typedef SmartPointer<const Self> ConstPointer;
00096
00098 itkNewMacro(Self);
00099
00101 itkTypeMacro( WarpVectorImageFilter, ImageToImageFilter );
00102
00104 typedef typename TOutputImage::RegionType OutputImageRegionType;
00105
00107 typedef typename Superclass::InputImageType InputImageType;
00108 typedef typename Superclass::InputImagePointer InputImagePointer;
00109 typedef typename Superclass::OutputImageType OutputImageType;
00110 typedef typename Superclass::OutputImagePointer OutputImagePointer;
00111 typedef typename Superclass::InputImageConstPointer InputImageConstPointer;
00112 typedef typename OutputImageType::IndexType IndexType;
00113 typedef typename OutputImageType::SizeType SizeType;
00114 typedef typename OutputImageType::PixelType PixelType;
00115 typedef typename OutputImageType::SpacingType SpacingType;
00116 typedef typename OutputImageType::PixelType::ValueType ValueType;
00117
00119 itkStaticConstMacro(ImageDimension, unsigned int,
00120 TOutputImage::ImageDimension );
00121
00123 itkStaticConstMacro(PixelDimension, unsigned int,
00124 PixelType::Dimension );
00125
00127 typedef TDeformationField DeformationFieldType;
00128 typedef typename DeformationFieldType::Pointer DeformationFieldPointer;
00129 typedef typename DeformationFieldType::PixelType DisplacementType;
00130
00132 typedef double CoordRepType;
00133 typedef VectorInterpolateImageFunction<InputImageType,CoordRepType> InterpolatorType;
00134 typedef typename InterpolatorType::Pointer InterpolatorPointer;
00135 typedef VectorLinearInterpolateImageFunction<InputImageType,CoordRepType>
00136 DefaultInterpolatorType;
00137
00139 typedef Point<CoordRepType,itkGetStaticConstMacro(ImageDimension)> PointType;
00140
00142 void SetDeformationField( DeformationFieldType * field );
00143
00145 DeformationFieldType * GetDeformationField(void);
00146
00148 itkSetObjectMacro( Interpolator, InterpolatorType );
00149
00151 itkGetObjectMacro( Interpolator, InterpolatorType );
00152
00154 itkSetMacro(OutputSpacing, SpacingType);
00155 virtual void SetOutputSpacing( const double* values);
00156
00158 itkGetConstReferenceMacro(OutputSpacing, SpacingType);
00159
00161 itkSetMacro(OutputOrigin, PointType);
00162 virtual void SetOutputOrigin( const double* values);
00163
00165 itkGetConstReferenceMacro(OutputOrigin, PointType);
00166
00168 itkSetMacro( EdgePaddingValue, PixelType );
00169
00171 itkGetMacro( EdgePaddingValue, PixelType );
00172
00178 virtual void GenerateOutputInformation();
00179
00186 virtual void GenerateInputRequestedRegion();
00187
00190 virtual void BeforeThreadedGenerateData();
00191
00192 protected:
00193 WarpVectorImageFilter();
00194 ~WarpVectorImageFilter() {};
00195 void PrintSelf(std::ostream& os, Indent indent) const;
00196
00200 void ThreadedGenerateData(const OutputImageRegionType& outputRegionForThread,
00201 int threadId );
00202
00203 private:
00204 WarpVectorImageFilter(const Self&);
00205 void operator=(const Self&);
00206
00207 PixelType m_EdgePaddingValue;
00208 SpacingType m_OutputSpacing;
00209 PointType m_OutputOrigin;
00210
00211 InterpolatorPointer m_Interpolator;
00212
00213 };
00214
00215 }
00216
00217 #ifndef ITK_MANUAL_INSTANTIATION
00218 #include "itkWarpVectorImageFilter.txx"
00219 #endif
00220
00221 #endif