ITK
4.1.0
Insight Segmentation and Registration Toolkit
|
00001 /*========================================================================= 00002 * 00003 * Copyright Insight Software Consortium 00004 * 00005 * Licensed under the Apache License, Version 2.0 (the "License"); 00006 * you may not use this file except in compliance with the License. 00007 * You may obtain a copy of the License at 00008 * 00009 * http://www.apache.org/licenses/LICENSE-2.0.txt 00010 * 00011 * Unless required by applicable law or agreed to in writing, software 00012 * distributed under the License is distributed on an "AS IS" BASIS, 00013 * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. 00014 * See the License for the specific language governing permissions and 00015 * limitations under the License. 00016 * 00017 *=========================================================================*/ 00018 #ifndef __itkIterativeInverseDisplacementFieldImageFilter_h 00019 #define __itkIterativeInverseDisplacementFieldImageFilter_h 00020 00021 00022 #include "itkWarpVectorImageFilter.h" 00023 #include "itkImageRegionIterator.h" 00024 #include "itkTimeProbe.h" 00025 00026 namespace itk 00027 { 00051 template< class TInputImage, class TOutputImage > 00052 class ITK_EXPORT IterativeInverseDisplacementFieldImageFilter: 00053 public ImageToImageFilter< TInputImage, TOutputImage > 00054 { 00055 public: 00057 typedef IterativeInverseDisplacementFieldImageFilter Self; 00058 typedef ImageToImageFilter< TInputImage, TOutputImage > Superclass; 00059 typedef SmartPointer< Self > Pointer; 00060 typedef SmartPointer< const Self > ConstPointer; 00061 00063 itkNewMacro(Self); 00064 00066 itkTypeMacro(IterativeInverseDisplacementFieldImageFilter, ImageToImageFilter); 00067 00069 typedef TInputImage InputImageType; 00070 typedef typename InputImageType::ConstPointer InputImageConstPointer; 00071 typedef typename InputImageType::Pointer InputImagePointer; 00072 typedef typename InputImageType::PointType InputImagePointType; 00073 typedef typename InputImageType::RegionType InputImageRegionType; 00074 typedef typename InputImageType::SpacingType InputImageSpacingType; 00075 typedef TOutputImage OutputImageType; 00076 typedef typename OutputImageType::Pointer OutputImagePointer; 00077 typedef typename OutputImageType::PixelType OutputImagePixelType; 00078 typedef typename OutputImageType::PointType OutputImagePointType; 00079 typedef typename OutputImageType::IndexType OutputImageIndexType; 00080 typedef typename OutputImagePixelType::ValueType OutputImageValueType; 00081 00082 typedef TimeProbe TimeType; 00083 00084 typedef ImageRegionConstIterator< InputImageType > InputConstIterator; 00085 typedef ImageRegionIterator< InputImageType > InputIterator; 00086 typedef ImageRegionIterator< OutputImageType > OutputIterator; 00087 00088 typedef WarpVectorImageFilter< TOutputImage, TInputImage, TOutputImage > VectorWarperType; 00089 00090 typedef VectorLinearInterpolateImageFunction< TInputImage, double > FieldInterpolatorType; 00091 typedef typename FieldInterpolatorType::Pointer FieldInterpolatorPointer; 00092 typedef typename FieldInterpolatorType::OutputType FieldInterpolatorOutputType; 00093 00094 itkSetMacro(NumberOfIterations, unsigned int); 00095 itkGetConstMacro(NumberOfIterations, unsigned int); 00096 00097 // If the error (in mm) between forward and backward mapping is smaller than 00098 // the StopValue, 00099 // the algorithm stops. 00100 // This value can be used to speed up the calculation. 00101 itkSetMacro(StopValue, double); 00102 itkGetConstMacro(StopValue, double); 00103 00104 #ifdef ITK_USE_CONCEPT_CHECKING 00105 00106 itkConceptMacro( OutputHasNumericTraitsCheck, 00107 ( Concept::HasNumericTraits< OutputImageValueType > ) ); 00108 00110 #endif 00111 protected: 00112 IterativeInverseDisplacementFieldImageFilter(); 00113 ~IterativeInverseDisplacementFieldImageFilter() {} 00114 00115 void PrintSelf(std::ostream & os, Indent indent) const; 00116 00117 void GenerateData(); 00118 00119 unsigned int m_NumberOfIterations; 00120 00121 double m_StopValue; 00122 double m_Time; 00123 private: 00124 IterativeInverseDisplacementFieldImageFilter(const Self &); //purposely not implemented 00125 void operator=(const Self &); //purposely not implemented 00126 }; 00127 } // end namespace itk 00128 00129 #ifndef ITK_MANUAL_INSTANTIATION 00130 #include "itkIterativeInverseDisplacementFieldImageFilter.hxx" 00131 #endif 00132 00133 #endif 00134