00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017 #ifndef __itkDistanceToCentroidMembershipFunction_h
00018 #define __itkDistanceToCentroidMembershipFunction_h
00019
00020 #include <vnl/vnl_vector.h>
00021 #include <vnl/vnl_transpose.h>
00022 #include <vnl/vnl_matrix.h>
00023 #include <vnl/algo/vnl_matrix_inverse.h>
00024 #include <vnl/algo/vnl_determinant.h>
00025
00026
00027 #include "itkMembershipFunctionBase.h"
00028
00029
00030 namespace itk{
00031 namespace Statistics{
00032
00042 template< class TVector >
00043 class ITK_EXPORT DistanceToCentroidMembershipFunction :
00044 public MembershipFunctionBase< TVector >
00045 {
00046 public:
00048 typedef DistanceToCentroidMembershipFunction Self;
00049 typedef MembershipFunctionBase< TVector > Superclass ;
00050 typedef SmartPointer<Self> Pointer;
00051 typedef SmartPointer<const Self> ConstPointer;
00052
00054 itkTypeMacro(DistanceToCentroidMembershipFunction, MembershipFunctionBase);
00055 itkNewMacro(Self);
00056
00058 typedef TVector MeasurementVectorType ;
00059
00061 void SetCentroid(const vnl_vector<double> ¢roid) ;
00062
00064 const vnl_vector<double> & GetCentroid() const;
00065
00067 itkSetMacro( NumberOfSamples, int );
00068
00070 itkGetMacro( NumberOfSamples, int );
00071
00075 double Evaluate(const MeasurementVectorType &measurement) const;
00076
00077 protected:
00078 DistanceToCentroidMembershipFunction(void) ;
00079 virtual ~DistanceToCentroidMembershipFunction(void) {}
00080 void PrintSelf(std::ostream& os, Indent indent) const;
00081
00082 private:
00083 vnl_vector< double > m_Centroid;
00084
00085
00086 int m_NumberOfSamples;
00087
00088 itkStaticConstMacro(VectorDimension, unsigned int, TVector::Dimension);
00089 typedef vnl_matrix_fixed<double,1,itkGetStaticConstMacro(VectorDimension)> ColumnVectorType;
00090
00091 mutable ColumnVectorType m_TempVec;
00092 mutable ColumnVectorType m_TempMat;
00093 };
00094
00095 }
00096 }
00097
00098 #ifndef ITK_MANUAL_INSTANTIATION
00099 #include "itkDistanceToCentroidMembershipFunction.txx"
00100 #endif
00101
00102 #endif