itk::PSMMeanCurvatureAttribute< TNumericType, VDimension > Class Template Reference

Inheritance diagram for itk::PSMMeanCurvatureAttribute< TNumericType, VDimension >:

Collaboration diagram for itk::PSMMeanCurvatureAttribute< TNumericType, VDimension >:

Public Types
typedef TNumericType	NumericType
typedef PSMMeanCurvatureAttribute	Self
typedef PSMContainerArrayAttribute< TNumericType, VDimension >	Superclass
typedef SmartPointer< Self >	Pointer
typedef SmartPointer< const Self >	ConstPointer
typedef WeakPointer< const Self >	ConstWeakPointer
typedef PSMParticleSystem< VDimension >	ParticleSystemType
typedef ParticleSystemType::PointType	PointType
typedef vnl_vector_fixed< TNumericType, VDimension >	VnlVectorType
Public Types inherited from itk::PSMContainerArrayAttribute< TNumericType, VDimension >
typedef TNumericType	DataType
typedef PSMContainerArrayAttribute	Self
typedef PSMAttribute< VDimension >	Superclass
typedef SmartPointer< Self >	Pointer
typedef SmartPointer< const Self >	ConstPointer
typedef WeakPointer< const Self >	ConstWeakPointer
Public Types inherited from itk::PSMAttribute< VDimension >
typedef PSMAttribute	Self
typedef DataObject	Superclass
typedef SmartPointer< Self >	Pointer
typedef SmartPointer< const Self >	ConstPointer
typedef WeakPointer< const Self >	ConstWeakPointer

Public Member Functions
	itkNewMacro (Self)
	itkTypeMacro (PSMMeanCurvatureAttribute, PSMContainerArrayAttribute)
virtual void	PositionAddEventCallback (Object *o, const EventObject &e)
virtual void	PositionSetEventCallback (Object *o, const EventObject &e)
virtual void	DomainAddEventCallback (Object *o, const EventObject &e)
void	ComputeMeanCurvature (const ParticleSystemType *system, unsigned int idx, unsigned int dom)
virtual void	ComputeCurvatureStatistics (const ParticleSystemType *, unsigned int d)
double	GetMeanCurvature (int d)
double	GetCurvatureStandardDeviation (int d)
Public Member Functions inherited from itk::PSMContainerArrayAttribute< TNumericType, VDimension >
	itkNewMacro (Self)
	itkTypeMacro (PSMContainerArrayAttribute, PSMAttribute)
virtual void	PositionRemoveEventCallback (Object *, const EventObject &)
void	ZeroAllValues ()
Public Member Functions inherited from itk::PSMAttribute< VDimension >
	itkNewMacro (Self)
	itkTypeMacro (PSMAttribute, DataObject)
virtual void	EventCallback (Object *, const EventObject &)
virtual void	EventWithIndexCallback (Object *, const EventObject &)
virtual void	TransformSetEventCallback (Object *, const EventObject &)
virtual void	PrefixTransformSetEventCallback (Object *, const EventObject &)
virtual void	NeighborhoodSetEventCallback (Object *, const EventObject &)

Protected Member Functions
void	PrintSelf (std::ostream &os, Indent indent) const
Protected Member Functions inherited from itk::PSMContainerArrayAttribute< TNumericType, VDimension >
void	PrintSelf (std::ostream &os, Indent indent) const
Protected Member Functions inherited from itk::PSMAttribute< VDimension >
void	PrintSelf (std::ostream &os, Indent indent) const

Additional Inherited Members
Public Attributes inherited from itk::PSMAttribute< VDimension >
DefinedCallbacksStruct	m_DefinedCallbacks

Detailed Description

template<class TNumericType, unsigned int VDimension>
class itk::PSMMeanCurvatureAttribute< TNumericType, VDimension >

Definition at line 34 of file itkPSMMeanCurvatureAttribute.h.

Member Typedef Documentation

template<class TNumericType, unsigned int VDimension>

typedef TNumericType itk::PSMMeanCurvatureAttribute< TNumericType, VDimension >::NumericType

Standard class typedefs

Definition at line 39 of file itkPSMMeanCurvatureAttribute.h.

template<class TNumericType, unsigned int VDimension>

typedef PSMParticleSystem<VDimension> itk::PSMMeanCurvatureAttribute< TNumericType, VDimension >::ParticleSystemType

Numeric types.

Definition at line 47 of file itkPSMMeanCurvatureAttribute.h.

Member Function Documentation

template<class TNumericType , unsigned int VDimension>

void itk::PSMMeanCurvatureAttribute< TNumericType, VDimension >::ComputeCurvatureStatistics ( const ParticleSystemType *  system, unsigned int  d  )

virtual

Compute the mean and std deviation of the curvature on the image surface.

Definition at line 31 of file itkPSMMeanCurvatureAttribute.hxx.

{
  typedef PSMImageDomainWithCurvature<TNumericType, VDimension> DomainType;
  typedef typename DomainType::ImageType ImageType;

  // Loop through a zero crossing image, project all the zero crossing points
  // to the surface, and use those points to comput curvature stats.  
  typedef itk::ZeroCrossingImageFilter<ImageType, ImageType> ZeroCrossingImageFilterType;
  typename ZeroCrossingImageFilterType::Pointer zc = ZeroCrossingImageFilterType::New() ;
  zc->SetInput( dynamic_cast<const DomainType *>(system->GetDomain(d))->GetImage() );
  zc->Update();
  
  itk::ImageRegionConstIteratorWithIndex<ImageType> it(zc->GetOutput(),
                                                       zc->GetOutput()->GetRequestedRegion());
  std::vector<double> datalist;
  m_MeanCurvatureList[d] = 0.0;
  m_CurvatureStandardDeviationList[d] = 0.0;
  const DomainType *domain = static_cast<const DomainType *>(system->GetDomain(d));
  
  for (; ! it.IsAtEnd(); ++it)
    {
      if (it.Get() == 1.0)
        {
          // Find closest pixel location to surface.
          PointType pos;
          domain->GetImage()->TransformIndexToPhysicalPoint(it.GetIndex(), pos);
          
          // Project point to surface.
          // Make sure constraints are enabled
          //      bool c = domain->GetConstraintsEnabled();
          //      domain->EnableConstraints();
          domain->ApplyConstraints(pos);
          //      domain->SetConstraintsEnabled(c);
          
          // Compute curvature at point.
          double mc = domain->GetCurvature(pos);      
          m_MeanCurvatureList[d] += mc;
          datalist.push_back(mc);
        }
    }
  double n = static_cast<double>(datalist.size());
  m_MeanCurvatureList[d] /= n;
  
  // Compute std deviation using point list
  for (unsigned int i = 0; i < datalist.size(); i++)
    {
      m_CurvatureStandardDeviationList[d] += (datalist[i] - m_MeanCurvatureList[d]) * (datalist[i] - m_MeanCurvatureList[d]);
    }
  m_CurvatureStandardDeviationList[d] = sqrt(m_CurvatureStandardDeviationList[d] / (n-1));
}

template<class TNumericType, unsigned int VDimension>

virtual void itk::PSMMeanCurvatureAttribute< TNumericType, VDimension >::DomainAddEventCallback ( Object *  , const EventObject &    )

inlinevirtual

Callbacks that may be defined by a subclass. If a subclass defines one of these callback methods, the corresponding flag in m_DefinedCallbacks should be set to true so that the ParticleSystem will know to register the appropriate event with this method.

Reimplemented from itk::PSMContainerArrayAttribute< TNumericType, VDimension >.

Definition at line 72 of file itkPSMMeanCurvatureAttribute.h.

  {
    Superclass::DomainAddEventCallback(o, e);
    m_MeanCurvatureList.push_back(0.0);
    m_CurvatureStandardDeviationList.push_back(0.0);
    const ParticleDomainAddEvent &event = dynamic_cast<const ParticleDomainAddEvent &>(e);
    const ParticleSystemType *ps= dynamic_cast<const ParticleSystemType *>(o);
    this->ComputeCurvatureStatistics(ps, event.GetDomainIndex());
  }

template<class TNumericType, unsigned int VDimension>

itk::PSMMeanCurvatureAttribute< TNumericType, VDimension >::itkNewMacro

(

Self

)

Method for creation through the object factory.

template<class TNumericType, unsigned int VDimension>

itk::PSMMeanCurvatureAttribute< TNumericType, VDimension >::itkTypeMacro	(	PSMMeanCurvatureAttribute< TNumericType, VDimension >	,
		PSMContainerArrayAttribute
	)

Run-time type information (and related methods).

---

The documentation for this class was generated from the following files:

• C:/Users/Brig/Documents/ShapeWorksStudio/src/ParticleShapeworks/include/itkPSMMeanCurvatureAttribute.h
• C:/Users/Brig/Documents/ShapeWorksStudio/src/ParticleShapeworks/include/itkPSMMeanCurvatureAttribute.hxx