00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017 #ifndef __itkBloxBoundaryPointImageToBloxBoundaryProfileImageFilter_h
00018 #define __itkBloxBoundaryPointImageToBloxBoundaryProfileImageFilter_h
00019
00020 #include "itkBloxBoundaryPointImage.h"
00021 #include "itkBloxBoundaryProfileImage.h"
00022 #include "itkImageToImageFilter.h"
00023 #include "itkSize.h"
00024
00025 namespace itk
00026 {
00027
00054 template< typename TSourceImage >
00055 class ITK_EXPORT BloxBoundaryPointImageToBloxBoundaryProfileImageFilter :
00056 public ImageToImageFilter<TSourceImage,BloxBoundaryProfileImage< ::itk::GetImageDimension<TSourceImage>::ImageDimension> >
00057 {
00058 public:
00060 itkStaticConstMacro(NDimensions, unsigned int, TSourceImage::ImageDimension);
00061
00063 typedef BloxBoundaryPointImageToBloxBoundaryProfileImageFilter Self;
00064 typedef ImageToImageFilter<TSourceImage,BloxBoundaryProfileImage<itkGetStaticConstMacro(NDimensions)> > Superclass;
00065 typedef SmartPointer<Self> Pointer;
00066 typedef SmartPointer<const Self> ConstPointer;
00067
00069 itkNewMacro(Self);
00070
00072 itkTypeMacro(itkBloxBoundaryPointImageToBloxBoundaryProfileImageFilter, ImageToImageFilter);
00073
00075 typedef BloxBoundaryPointImage<itkGetStaticConstMacro(NDimensions)> BoundaryPointImageType;
00076 typedef typename BoundaryPointImageType::Pointer BoundaryPointImagePointer;
00077 typedef typename BoundaryPointImageType::RegionType BoundaryPointImageRegionType;
00078 typedef typename BoundaryPointImageType::PixelType BoundaryPointImagePixelType;
00079 typedef typename BoundaryPointImageType::ConstPointer BoundaryPointImageConstPointer;
00080
00082 typedef TSourceImage SourceImageType;
00083 typedef typename SourceImageType::Pointer SourceImagePointer;
00084 typedef typename SourceImageType::RegionType SourceImageRegionType;
00085 typedef typename SourceImageType::PixelType SourceImagePixelType;
00086 typedef typename SourceImageType::ConstPointer SourceImageConstPointer;
00087
00089 typedef BloxBoundaryProfileImage<itkGetStaticConstMacro(NDimensions)> OutputImageType;
00090 typedef typename OutputImageType::Pointer OutputImagePointer;
00091 typedef typename OutputImageType::RegionType OutputImageRegionType;
00092 typedef typename OutputImageType::PixelType OutputImagePixelType;
00093
00095 typedef typename BloxBoundaryProfileImage<itkGetStaticConstMacro(NDimensions)>::IndexType IndexType;
00096
00098 typedef typename BloxBoundaryProfileImage<itkGetStaticConstMacro(NDimensions)>::PixelType PixelType;
00099
00101 typedef Point<double, itkGetStaticConstMacro(NDimensions)> PositionType;
00102
00104 typedef typename PositionType::VectorType VectorType;
00105
00107 void SetInput1( const SourceImageType * image1 );
00108
00110 void SetInput2( const BoundaryPointImageType * image2 );
00111
00113 double FindAccumulatorMaximum();
00114
00116 double FindAccumulatorMinimum();
00117
00119 void FindBoundaryProfiles();
00120
00122 bool AddSplatToAccumulatorAndNormalizer(int binNumber, double weight, double sourcePixelValue);
00123
00125 void NormalizeSplatAccumulator();
00126
00128 int FitProfile();
00129
00131 void Initialize(double setUniqueAxis, double setSymmetricAxes, unsigned int numberOfBins,
00132 unsigned int splatMethod, unsigned int spaceDimension);
00133
00134 protected:
00135 BloxBoundaryPointImageToBloxBoundaryProfileImageFilter();
00136 ~BloxBoundaryPointImageToBloxBoundaryProfileImageFilter();
00137
00138 void PrintSelf(std::ostream& os, Indent indent) const;
00139
00141 void GenerateData();
00142
00143 private:
00144 BloxBoundaryPointImageToBloxBoundaryProfileImageFilter(const Self&);
00145 void operator=(const Self&);
00146
00148 double m_UniqueAxis;
00149
00151 double m_SymmetricAxes;
00152
00154 unsigned int m_NumberOfBins;
00155
00157 unsigned int m_SplatMethod;
00158
00160 unsigned long int m_NumBoundaryProfiles;
00161
00163 double * m_Accumulator;
00164
00166 double * m_Normalizer;
00167
00169 double * m_NormalizedAccumulator;
00170
00172 double * m_FinalParameters;
00173
00175 unsigned int m_SpaceDimension;
00176 };
00177
00178 }
00179
00180 #ifndef ITK_MANUAL_INSTANTIATION
00181 #include "itkBloxBoundaryPointImageToBloxBoundaryProfileImageFilter.txx"
00182 #endif
00183
00184 #endif