Go to the documentation of this file.00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017 #ifndef __itkMutualInformationImageToImageMetric_h
00018 #define __itkMutualInformationImageToImageMetric_h
00019
00020 #include "itkImageToImageMetric.h"
00021 #include "itkCovariantVector.h"
00022 #include "itkPoint.h"
00023
00024 #include "itkIndex.h"
00025 #include "itkKernelFunction.h"
00026 #include "itkCentralDifferenceImageFunction.h"
00027
00028 namespace itk
00029 {
00030
00089 template <class TFixedImage,class TMovingImage >
00090 class ITK_EXPORT MutualInformationImageToImageMetric :
00091 public ImageToImageMetric< TFixedImage, TMovingImage >
00092 {
00093 public:
00094
00096 typedef MutualInformationImageToImageMetric Self;
00097 typedef ImageToImageMetric< TFixedImage, TMovingImage > Superclass;
00098 typedef SmartPointer<Self> Pointer;
00099 typedef SmartPointer<const Self> ConstPointer;
00100
00102 itkNewMacro(Self);
00103
00105 itkTypeMacro(MutualInformationImageToImageMetric, ImageToImageMetric);
00106
00108 typedef typename Superclass::TransformType TransformType;
00109 typedef typename Superclass::TransformPointer TransformPointer;
00110 typedef typename Superclass::TransformJacobianType TransformJacobianType;
00111 typedef typename Superclass::InterpolatorType InterpolatorType;
00112 typedef typename Superclass::MeasureType MeasureType;
00113 typedef typename Superclass::DerivativeType DerivativeType;
00114 typedef typename Superclass::ParametersType ParametersType;
00115 typedef typename Superclass::FixedImageType FixedImageType;
00116 typedef typename Superclass::MovingImageType MovingImageType;
00117 typedef typename Superclass::FixedImageConstPointer FixedImageConstPointer;
00118 typedef typename Superclass::MovingImageConstPointer MovingImageCosntPointer;
00119
00121 typedef typename FixedImageType::IndexType FixedImageIndexType;
00122 typedef typename FixedImageIndexType::IndexValueType FixedImageIndexValueType;
00123 typedef typename MovingImageType::IndexType MovingImageIndexType;
00124 typedef typename TransformType::InputPointType FixedImagePointType;
00125 typedef typename TransformType::OutputPointType MovingImagePointType;
00126
00128 itkStaticConstMacro(MovingImageDimension, unsigned int,
00129 MovingImageType::ImageDimension);
00130
00132 void GetDerivative(
00133 const ParametersType& parameters,
00134 DerivativeType & Derivative ) const;
00135
00137 MeasureType GetValue( const ParametersType& parameters ) const;
00138
00140 void GetValueAndDerivative( const ParametersType& parameters,
00141 MeasureType& Value, DerivativeType& Derivative ) const;
00142
00147 void SetNumberOfSpatialSamples( unsigned int num );
00148
00150 itkGetConstReferenceMacro( NumberOfSpatialSamples, unsigned int );
00151
00157 itkSetClampMacro( MovingImageStandardDeviation, double,
00158 NumericTraits<double>::NonpositiveMin(), NumericTraits<double>::max() );
00159 itkGetConstReferenceMacro( MovingImageStandardDeviation, double );
00161
00167 itkSetClampMacro( FixedImageStandardDeviation, double,
00168 NumericTraits<double>::NonpositiveMin(), NumericTraits<double>::max() );
00169 itkGetConstMacro( FixedImageStandardDeviation, double );
00171
00174 itkSetObjectMacro( KernelFunction, KernelFunction );
00175 itkGetObjectMacro( KernelFunction, KernelFunction );
00177
00187 void ReinitializeSeed();
00188 void ReinitializeSeed(int);
00190
00191 protected:
00192 MutualInformationImageToImageMetric();
00193 virtual ~MutualInformationImageToImageMetric() {};
00194 void PrintSelf(std::ostream& os, Indent indent) const;
00195
00196 private:
00197 MutualInformationImageToImageMetric(const Self&);
00198 void operator=(const Self&);
00199
00203 class SpatialSample
00204 {
00205 public:
00206 SpatialSample():FixedImageValue(0.0),MovingImageValue(0.0)
00207 { FixedImagePointValue.Fill( 0.0 ); }
00208 ~SpatialSample(){};
00209
00210 FixedImagePointType FixedImagePointValue;
00211 double FixedImageValue;
00212 double MovingImageValue;
00213 };
00214
00216 typedef std::vector<SpatialSample> SpatialSampleContainer;
00217
00220 mutable SpatialSampleContainer m_SampleA;
00221
00224 mutable SpatialSampleContainer m_SampleB;
00225
00226 unsigned int m_NumberOfSpatialSamples;
00227 double m_MovingImageStandardDeviation;
00228 double m_FixedImageStandardDeviation;
00229 typename KernelFunction::Pointer m_KernelFunction;
00230 double m_MinProbability;
00231
00237 virtual void SampleFixedImageDomain( SpatialSampleContainer& samples ) const;
00238
00242 void CalculateDerivatives( const FixedImagePointType& , DerivativeType& ) const;
00243
00244 typedef typename Superclass::CoordinateRepresentationType
00245 CoordinateRepresentationType;
00246 typedef CentralDifferenceImageFunction< MovingImageType,
00247 CoordinateRepresentationType > DerivativeFunctionType;
00248
00249 typename DerivativeFunctionType::Pointer m_DerivativeCalculator;
00250
00251 bool m_ReseedIterator;
00252 int m_RandomSeed;
00253
00254
00255 };
00256
00257 }
00258
00259 #ifndef ITK_MANUAL_INSTANTIATION
00260 #include "itkMutualInformationImageToImageMetric.txx"
00261 #endif
00262
00263 #endif
00264