Deformably register two images using the PDEs and NCC. More...

#include <otbNCCRegistrationFilter.h>

Inheritance diagram for otb::NCCRegistrationFilter< TFixedImage, TMovingImage, TDeformationField >:

Collaboration diagram for otb::NCCRegistrationFilter< TFixedImage, TMovingImage, TDeformationField >:

Public Types
typedef NCCRegistrationFilter	Self
typedef itk::PDEDeformableRegistrationFilter < TFixedImage, TMovingImage, TDeformationField >	Superclass
typedef itk::SmartPointer< Self >	Pointer
typedef itk::SmartPointer < const Self >	ConstPointer
typedef Superclass::TimeStepType	TimeStepType
typedef Superclass::FixedImageType	FixedImageType
typedef Superclass::FixedImagePointer	FixedImagePointer
typedef Superclass::MovingImageType	MovingImageType
typedef Superclass::MovingImagePointer	MovingImagePointer
typedef Superclass::DeformationFieldType	DeformationFieldType
typedef Superclass::DeformationFieldPointer	DeformationFieldPointer
typedef Superclass::FiniteDifferenceFunctionType	FiniteDifferenceFunctionType
typedef NCCRegistrationFunction < FixedImageType, MovingImageType, DeformationFieldType >	NCCRegistrationFunctionType
typedef NCCRegistrationFunctionType::RadiusType	RadiusType
typedef FixedImageType::ConstPointer	FixedImageConstPointer
typedef MovingImageType::ConstPointer	MovingImageConstPointer
typedef Superclass::OutputImageType	OutputImageType
typedef PDEDeformableRegistrationFunction < FixedImageType, MovingImageType, DeformationFieldType >	PDEDeformableRegistrationFunctionType
typedef Superclass::InputImageType	InputImageType
typedef Superclass::PixelType	PixelType
typedef OutputImageType	UpdateBufferType
enum	FilterStateType
typedef TDeformationField::PixelType	OutputPixelType
typedef TDeformationField::PixelType	InputPixelType
typedef NumericTraits < OutputPixelType >::ValueType	OutputPixelValueType
typedef NumericTraits < InputPixelType >::ValueType	InputPixelValueType
typedef FiniteDifferenceFunctionType::NeighborhoodScalesType	NeighborhoodScalesType
typedef Superclass::OutputImagePointer	OutputImagePointer
typedef Superclass::OutputImageRegionType	OutputImageRegionType
typedef Superclass::OutputImagePixelType	OutputImagePixelType
typedef InputImageType::Pointer	InputImagePointer
typedef InputImageType::ConstPointer	InputImageConstPointer
typedef InputImageType::RegionType	InputImageRegionType
typedef InputImageType::PixelType	InputImagePixelType
typedef DataObject::Pointer	DataObjectPointer
typedef std::vector < DataObjectPointer >	DataObjectPointerArray
typedef DataObjectPointerArray::size_type	DataObjectPointerArraySizeType
Public Member Functions
virtual const char *	GetNameOfClass () const
virtual double	GetMetric () const
virtual void	SetNCCRadius (RadiusType radius)
virtual RadiusType	GetNCCRadius () const
void	SetFixedImage (const FixedImageType *ptr)
const FixedImageType *	GetFixedImage (void) const
void	SetMovingImage (const MovingImageType *ptr)
const MovingImageType *	GetMovingImage (void) const
void	SetInitialDeformationField (const DeformationFieldType *ptr)
DeformationFieldType *	GetDeformationField ()
virtual std::vector < SmartPointer< DataObject > >::size_type	GetNumberOfValidRequiredInputs () const
virtual void	SetSmoothDeformationField (bool _arg)
virtual bool	GetSmoothDeformationField () const
virtual void	SmoothDeformationFieldOn ()
virtual void	SmoothDeformationFieldOff ()
virtual void	SetStandardDeviations (double data[])
virtual void	SetStandardDeviations (double value)
const double *	GetStandardDeviations (void)
virtual void	SetSmoothUpdateField (bool _arg)
virtual bool	GetSmoothUpdateField () const
virtual void	SmoothUpdateFieldOn ()
virtual void	SmoothUpdateFieldOff ()
virtual void	SetUpdateFieldStandardDeviations (double data[])
virtual void	SetUpdateFieldStandardDeviations (double value)
const double *	GetUpdateFieldStandardDeviations (void)
virtual void	StopRegistration ()
virtual void	SetMaximumError (double _arg)
virtual double	GetMaximumError () const
virtual void	SetMaximumKernelWidth (unsigned int _arg)
virtual unsigned int	GetMaximumKernelWidth () const
virtual const unsigned int &	GetElapsedIterations ()
virtual const FiniteDifferenceFunctionType::Pointer &	GetDifferenceFunction () const
virtual void	SetDifferenceFunction (FiniteDifferenceFunctionType *_arg)
virtual void	SetNumberOfIterations (unsigned int _arg)
virtual const unsigned int &	GetNumberOfIterations ()
virtual void	SetUseImageSpacing (bool _arg)
virtual void	UseImageSpacingOn ()
virtual void	UseImageSpacingOff ()
virtual const bool &	GetUseImageSpacing ()
virtual void	SetMaximumRMSError (double _arg)
virtual const double &	GetMaximumRMSError ()
virtual void	SetRMSChange (double _arg)
virtual const double &	GetRMSChange ()
void	SetStateToInitialized ()
void	SetStateToUninitialized ()
virtual void	SetState (FilterStateType _arg)
virtual const FilterStateType &	GetState ()
virtual void	SetManualReinitialization (bool _arg)
virtual const bool &	GetManualReinitialization ()
virtual void	ManualReinitializationOn ()
virtual void	ManualReinitializationOff ()
virtual void	SetInPlace (bool _arg)
virtual bool	GetInPlace () const
virtual void	InPlaceOn ()
virtual void	InPlaceOff ()
virtual bool	CanRunInPlace () const
virtual void	SetInput (const InputImageType *image)
virtual void	SetInput (unsigned int, const TDeformationField *image)
const InputImageType *	GetInput (void)
const InputImageType *	GetInput (unsigned int idx)
virtual void	PushBackInput (const InputImageType *image)
virtual void	PopBackInput ()
virtual void	PushFrontInput (const InputImageType *image)
virtual void	PopFrontInput ()
OutputImageType *	GetOutput (void)
OutputImageType *	GetOutput (unsigned int idx)
virtual void	GraftOutput (DataObject *output)
virtual void	GraftNthOutput (unsigned int idx, DataObject *output)
virtual DataObjectPointer	MakeOutput (unsigned int idx)
DataObjectPointerArray &	GetInputs ()
DataObjectPointerArraySizeType	GetNumberOfInputs () const
DataObjectPointerArray &	GetOutputs ()
DataObjectPointerArraySizeType	GetNumberOfOutputs () const
virtual void	SetAbortGenerateData (bool _arg)
virtual const bool &	GetAbortGenerateData ()
virtual void	AbortGenerateDataOn ()
virtual void	AbortGenerateDataOff ()
virtual void	SetProgress (float _arg)
virtual const float &	GetProgress ()
void	UpdateProgress (float amount)
virtual void	Update ()
virtual void	UpdateLargestPossibleRegion ()
virtual void	UpdateOutputInformation ()
virtual void	PropagateRequestedRegion (DataObject *output)
virtual void	UpdateOutputData (DataObject *output)
virtual void	EnlargeOutputRequestedRegion (DataObject *)
virtual void	ResetPipeline ()
virtual void	SetReleaseDataFlag (bool flag)
virtual bool	GetReleaseDataFlag () const
void	ReleaseDataFlagOn ()
void	ReleaseDataFlagOff ()
virtual void	SetReleaseDataBeforeUpdateFlag (bool _arg)
virtual const bool &	GetReleaseDataBeforeUpdateFlag ()
virtual void	ReleaseDataBeforeUpdateFlagOn ()
virtual void	ReleaseDataBeforeUpdateFlagOff ()
virtual void	SetNumberOfThreads (int _arg)
virtual const int &	GetNumberOfThreads ()
MultiThreader *	GetMultiThreader ()
virtual void	PrepareOutputs ()
virtual LightObject::Pointer	CreateAnother () const
virtual void	DebugOn () const
virtual void	DebugOff () const
bool	GetDebug () const
void	SetDebug (bool debugFlag) const
virtual unsigned long	GetMTime () const
virtual void	Modified () const
virtual void	Register () const
virtual void	UnRegister () const
virtual void	SetReferenceCount (int)
unsigned long	AddObserver (const EventObject &event, Command *)
unsigned long	AddObserver (const EventObject &event, Command *) const
Command *	GetCommand (unsigned long tag)
void	InvokeEvent (const EventObject &)
void	InvokeEvent (const EventObject &) const
void	RemoveObserver (unsigned long tag)
void	RemoveAllObservers ()
bool	HasObserver (const EventObject &event) const
MetaDataDictionary &	GetMetaDataDictionary (void)
const MetaDataDictionary &	GetMetaDataDictionary (void) const
void	SetMetaDataDictionary (const MetaDataDictionary &rhs)
virtual void	Delete ()
void	Print (std::ostream &os, Indent indent=0) const
virtual int	GetReferenceCount () const
Static Public Member Functions
static Pointer	New ()
static void	SetGlobalWarningDisplay (bool flag)
static bool	GetGlobalWarningDisplay ()
static void	GlobalWarningDisplayOn ()
static void	GlobalWarningDisplayOff ()
static void	BreakOnError ()
Static Public Attributes
static const unsigned int	ImageDimension = Superclass::ImageDimension
static const unsigned int	InputImageDimension
static const unsigned int	OutputImageDimension
Protected Types
typedef UpdateBufferType::RegionType	ThreadRegionType
typedef ImageToImageFilterDetail::ImageRegionCopier < itkGetStaticConstMacro(OutputImageDimension), itkGetStaticConstMacro(InputImageDimension)>	InputToOutputRegionCopierType
typedef ImageToImageFilterDetail::ImageRegionCopier < itkGetStaticConstMacro(InputImageDimension), itkGetStaticConstMacro(OutputImageDimension)>	OutputToInputRegionCopierType
typedef int	InternalReferenceCountType
Protected Member Functions
void	PrintSelf (std::ostream &os, itk::Indent indent) const
virtual void	InitializeIteration ()
virtual void	ApplyUpdate (TimeStepType dt)
virtual void	GenerateInputRequestedRegion ()
virtual bool	Halt ()
virtual void	CopyInputToOutput ()
virtual void	SmoothDeformationField ()
virtual void	SmoothUpdateField ()
virtual void	PostProcessOutput ()
virtual void	Initialize ()
virtual void	GenerateOutputInformation ()
virtual void	ApplyUpdate (TimeStepType dt)
virtual void	ApplyUpdate (TimeStepType dt)=0
virtual UpdateBufferType *	GetUpdateBuffer ()
virtual TimeStepType	CalculateChange ()
virtual void	AllocateUpdateBuffer ()
virtual void	ThreadedApplyUpdate (TimeStepType dt, const ThreadRegionType &regionToProcess, int threadId)
virtual TimeStepType	ThreadedCalculateChange (const ThreadRegionType &regionToProcess, int threadId)
virtual void	GenerateData ()
virtual bool	ThreadedHalt (void *)
virtual TimeStepType	ResolveTimeStep (const TimeStepType timeStepList, const bool valid, int size)
virtual void	SetElapsedIterations (unsigned int _arg)
virtual void	AllocateOutputs ()
virtual void	ReleaseInputs ()
const DataObject *	GetInput (unsigned int idx) const
void	PushBackInput (const DataObject *input)
void	PushFrontInput (const DataObject *input)
virtual void	CallCopyOutputRegionToInputRegion (InputImageRegionType &destRegion, const OutputImageRegionType &srcRegion)
virtual void	CallCopyInputRegionToOutputRegion (OutputImageRegionType &destRegion, const InputImageRegionType &srcRegion)
const DataObject *	GetOutput (unsigned int idx) const
virtual void	ThreadedGenerateData (const OutputImageRegionType &outputRegionForThread, int threadId)
virtual void	BeforeThreadedGenerateData ()
virtual void	AfterThreadedGenerateData ()
virtual int	SplitRequestedRegion (int i, int num, OutputImageRegionType &splitRegion)
virtual void	SetNthInput (unsigned int num, DataObject *input)
virtual void	AddInput (DataObject *input)
virtual void	RemoveInput (DataObject *input)
virtual void	SetNumberOfRequiredInputs (unsigned int _arg)
virtual const unsigned int &	GetNumberOfRequiredInputs ()
void	SetNumberOfInputs (unsigned int num)
virtual void	SetNthOutput (unsigned int num, DataObject *output)
virtual void	AddOutput (DataObject *output)
virtual void	RemoveOutput (DataObject *output)
virtual void	SetNumberOfRequiredOutputs (unsigned int _arg)
virtual const unsigned int &	GetNumberOfRequiredOutputs ()
void	SetNumberOfOutputs (unsigned int num)
virtual void	GenerateOutputRequestedRegion (DataObject *output)
virtual void	PropagateResetPipeline ()
virtual void	CacheInputReleaseDataFlags ()
virtual void	RestoreInputReleaseDataFlags ()
bool	PrintObservers (std::ostream &os, Indent indent) const
virtual void	PrintHeader (std::ostream &os, Indent indent) const
virtual void	PrintTrailer (std::ostream &os, Indent indent) const
Static Protected Member Functions
static ITK_THREAD_RETURN_TYPE	ThreaderCallback (void *arg)
Protected Attributes
unsigned int	m_NumberOfIterations
unsigned int	m_ElapsedIterations
bool	m_ManualReinitialization
double	m_RMSChange
double	m_MaximumRMSError
bool	m_Updating
TimeStamp	m_OutputInformationMTime
InternalReferenceCountType	m_ReferenceCount
SimpleFastMutexLock	m_ReferenceCountLock
Private Member Functions
	NCCRegistrationFilter (const Self &)
void	operator= (const Self &)

Detailed Description

template<class TFixedImage, class TMovingImage, class TDeformationField>
class otb::NCCRegistrationFilter< TFixedImage, TMovingImage, TDeformationField >

Deformably register two images using the PDEs and NCC.

NCCRegistrationFilter implements the ####### register two images by computing the deformation field which will map a moving image onto a fixed image.

A deformation field is represented as a image whose pixel type is some vector type with at least N elements, where N is the dimension of the fixed image. The vector type must support element access via operator []. It is assumed that the vector elements behave like floating point scalars.

This class is templated over the fixed image type, moving image type and the deformation field type.

The input fixed and moving images are set via methods SetFixedImage and SetMovingImage respectively. An initial deformation field maybe set via SetInitialDeformationField or SetInput. If no initial field is set, a zero field is used as the initial condition.

The algorithm has one parameter: the number of iteration to be performed.

The output deformation field can be obtained via methods GetOutput or GetDeformationField.

This class make use of the finite difference solver hierarchy. Update for each iteration is computed in NCCRegistrationFunction.

Warning:: This filter assumes that the fixed image type, moving image type and deformation field type all have the same number of dimensions.

See also:: NCCRegistrationFunction

Examples:: DisparityMap/NCCRegistrationFilterExample.cxx.

Definition at line 65 of file otbNCCRegistrationFilter.h.

Member Typedef Documentation

typedef DataObject::Pointer itk::ImageSource< TDeformationField >::DataObjectPointer [inherited]

Smart Pointer type to a DataObject.

Reimplemented from itk::ProcessObject.

Definition at line 62 of file itkImageSource.h.

typedef std::vector<DataObjectPointer> itk::ProcessObject::DataObjectPointerArray [inherited]

STL Array of SmartPointers to DataObjects

Definition at line 103 of file itkProcessObject.h.

typedef DataObjectPointerArray::size_type itk::ProcessObject::DataObjectPointerArraySizeType [inherited]

Size type of an std::vector

Definition at line 112 of file itkProcessObject.h.

template<class TFixedImage, class TMovingImage, class TDeformationField>

typedef Superclass::DeformationFieldType otb::NCCRegistrationFilter< TFixedImage, TMovingImage, TDeformationField >::DeformationFieldType

Deformation field type.

Reimplemented from itk::PDEDeformableRegistrationFilter< TFixedImage, TMovingImage, TDeformationField >.

Definition at line 97 of file otbNCCRegistrationFilter.h.

template<class TFixedImage, class TMovingImage, class TDeformationField>

typedef Superclass::FiniteDifferenceFunctionType otb::NCCRegistrationFilter< TFixedImage, TMovingImage, TDeformationField >::FiniteDifferenceFunctionType

FiniteDifferenceFunction type.

Reimplemented from itk::PDEDeformableRegistrationFilter< TFixedImage, TMovingImage, TDeformationField >.

Definition at line 103 of file otbNCCRegistrationFilter.h.

template<class TFixedImage, class TMovingImage, class TDeformationField>

typedef Superclass::FixedImageType otb::NCCRegistrationFilter< TFixedImage, TMovingImage, TDeformationField >::FixedImageType

FixedImage image type.

Reimplemented from itk::PDEDeformableRegistrationFilter< TFixedImage, TMovingImage, TDeformationField >.

Definition at line 88 of file otbNCCRegistrationFilter.h.

typedef Superclass::InputImageType itk::DenseFiniteDifferenceImageFilter< TDeformationField , TDeformationField >::InputImageType [inherited]

Convenient typedefs

Reimplemented from itk::FiniteDifferenceImageFilter< TDeformationField, TDeformationField >.

Definition at line 76 of file itkDenseFiniteDifferenceImageFilter.h.

typedef ImageToImageFilterDetail::ImageRegionCopier<itkGetStaticConstMacro(OutputImageDimension), itkGetStaticConstMacro(InputImageDimension)> itk::ImageToImageFilter< TDeformationField , TDeformationField >::InputToOutputRegionCopierType [protected, inherited]

Typedef for the region copier function object that converts an input region to an output region.

Definition at line 161 of file itkImageToImageFilter.h.

typedef int itk::LightObject::InternalReferenceCountType [protected, inherited]

Define the type of the reference count according to the target. This allows the use of atomic operations

Definition at line 137 of file itkLightObject.h.

template<class TFixedImage, class TMovingImage, class TDeformationField>

typedef Superclass::MovingImageType otb::NCCRegistrationFilter< TFixedImage, TMovingImage, TDeformationField >::MovingImageType

MovingImage image type.

Reimplemented from itk::PDEDeformableRegistrationFilter< TFixedImage, TMovingImage, TDeformationField >.

Definition at line 92 of file otbNCCRegistrationFilter.h.

template<class TFixedImage, class TMovingImage, class TDeformationField>

typedef NCCRegistrationFunction<FixedImageType, MovingImageType, DeformationFieldType> otb::NCCRegistrationFilter< TFixedImage, TMovingImage, TDeformationField >::NCCRegistrationFunctionType

NCCRegistrationFilterFunction type.

Definition at line 107 of file otbNCCRegistrationFilter.h.

typedef Superclass::OutputImageRegionType itk::InPlaceImageFilter< TDeformationField , TDeformationField >::OutputImageRegionType [inherited]

Superclass typedefs.

Reimplemented from itk::ImageToImageFilter< TDeformationField, TDeformationField >.

Definition at line 78 of file itkInPlaceImageFilter.h.

template<class TFixedImage , class TMovingImage , class TDeformationField >

typedef Superclass::OutputImageType itk::PDEDeformableRegistrationFilter< TFixedImage, TMovingImage, TDeformationField >::OutputImageType [inherited]

Types inherithed from the superclass

Reimplemented from itk::DenseFiniteDifferenceImageFilter< TDeformationField, TDeformationField >.

Definition at line 104 of file itkPDEDeformableRegistrationFilter.h.

typedef TDeformationField ::PixelType itk::FiniteDifferenceImageFilter< TDeformationField , TDeformationField >::OutputPixelType [inherited]

The pixel type of the output image will be used in computations.

Definition at line 143 of file itkFiniteDifferenceImageFilter.h.

typedef NumericTraits< OutputPixelType >::ValueType itk::FiniteDifferenceImageFilter< TDeformationField , TDeformationField >::OutputPixelValueType [inherited]

Extract value type in case the pixel is of vector type

Definition at line 148 of file itkFiniteDifferenceImageFilter.h.

typedef ImageToImageFilterDetail::ImageRegionCopier<itkGetStaticConstMacro(InputImageDimension), itkGetStaticConstMacro(OutputImageDimension)> itk::ImageToImageFilter< TDeformationField , TDeformationField >::OutputToInputRegionCopierType [protected, inherited]

Typedef for the region copier function object that converts an output region to an input region.

Definition at line 166 of file itkImageToImageFilter.h.

template<class TFixedImage , class TMovingImage , class TDeformationField >

typedef PDEDeformableRegistrationFunction<FixedImageType,MovingImageType, DeformationFieldType> itk::PDEDeformableRegistrationFilter< TFixedImage, TMovingImage, TDeformationField >::PDEDeformableRegistrationFunctionType [inherited]

PDEDeformableRegistrationFilterFunction type.

Definition at line 112 of file itkPDEDeformableRegistrationFilter.h.

typedef Superclass::PixelType itk::DenseFiniteDifferenceImageFilter< TDeformationField , TDeformationField >::PixelType [inherited]

The pixel type of the output image will be used in computations. Inherited from the superclass.

Reimplemented from itk::FiniteDifferenceImageFilter< TDeformationField, TDeformationField >.

Definition at line 90 of file itkDenseFiniteDifferenceImageFilter.h.

template<class TFixedImage, class TMovingImage, class TDeformationField>

typedef NCCRegistrationFilter otb::NCCRegistrationFilter< TFixedImage, TMovingImage, TDeformationField >::Self

Standard class typedefs.

Reimplemented from itk::PDEDeformableRegistrationFilter< TFixedImage, TMovingImage, TDeformationField >.

Definition at line 71 of file otbNCCRegistrationFilter.h.

typedef UpdateBufferType::RegionType itk::DenseFiniteDifferenceImageFilter< TDeformationField , TDeformationField >::ThreadRegionType [protected, inherited]

The type of region used for multithreading

Definition at line 140 of file itkDenseFiniteDifferenceImageFilter.h.

template<class TFixedImage, class TMovingImage, class TDeformationField>

typedef Superclass::TimeStepType otb::NCCRegistrationFilter< TFixedImage, TMovingImage, TDeformationField >::TimeStepType

Inherit types from superclass.

Reimplemented from itk::DenseFiniteDifferenceImageFilter< TDeformationField, TDeformationField >.

Definition at line 82 of file otbNCCRegistrationFilter.h.

typedef OutputImageType itk::DenseFiniteDifferenceImageFilter< TDeformationField , TDeformationField >::UpdateBufferType [inherited]

The container type for the update buffer.

Definition at line 96 of file itkDenseFiniteDifferenceImageFilter.h.

Member Function Documentation

virtual void itk::ProcessObject::AbortGenerateDataOn ( ) [virtual, inherited]

Turn on and off the AbortGenerateData flag.

void itk::ProcessObject::AddInput ( DataObject * input ) [protected, virtual, inherited]

Adds an input to the first null position in the input list. Expands the list memory if necessary

Reimplemented in otb::PersistentDescriptorsListSampleGenerator< TInputImage, TVectorData, TFunctionType, TListSample, TLabelListSample >, otb::DescriptorsListSampleGenerator< TInputImage, TVectorData, TListSample, TLabelListSample, TOutputPrecision, TCoordRep >, otb::PersistentObjectDetectionClassifier< TInputImage, TOutputVectorData, TLabel, TFunctionType >, and otb::ObjectDetectionClassifier< TInputImage, TOutputVectorData, TLabel, TFunctionPrecision, TCoordRep >.

Definition at line 144 of file itkProcessObject.cxx.

{
  DataObjectPointerArraySizeType idx;
  
  this->Modified();
  
  for (idx = 0; idx < m_Inputs.size(); ++idx)
    {
    if (!m_Inputs[idx])
      {
      m_Inputs[idx] = input;
      return;
      }
    }
  
  this->SetNumberOfInputs( static_cast<int>( m_Inputs.size() + 1 ) );
  m_Inputs[ static_cast<int>( m_Inputs.size() ) - 1] = input;
}

unsigned long itk::Object::AddObserver	(	const EventObject &	event,
		Command *	cmd
	)			`[inherited]`

Allow people to add/remove/invoke observers (callbacks) to any ITK object. This is an implementation of the subject/observer design pattern. An observer is added by specifying an event to respond to and an itk::Command to execute. It returns an unsigned long tag which can be used later to remove the event or retrieve the command. The memory for the Command becomes the responsibility of this object, so don't pass the same instance of a command to two different objects

Definition at line 389 of file itkObject.cxx.

Referenced by otb::StreamingImageVirtualWriter< TInputImage >::GenerateData(), itk::ProgressAccumulator::RegisterInternalFilter(), otb::StreamingImageFileWriter< TInputImage >::UpdateOutputData(), and otb::WriterWatcherBase::WriterWatcherBase().

{
  if (!this->m_SubjectImplementation)
    {
    this->m_SubjectImplementation = new SubjectImplementation;
    }
  return this->m_SubjectImplementation->AddObserver(event,cmd);
}

void itk::ProcessObject::AddOutput ( DataObject * output ) [protected, virtual, inherited]

Adds an output to the first null position in the output list. Expands the list memory if necessary

Definition at line 390 of file itkProcessObject.cxx.

References itk::DataObject::ConnectSource().

{
  unsigned int idx;
  
  for (idx = 0; idx < m_Outputs.size(); ++idx)
    {
    if ( m_Outputs[idx].IsNull() )
      {
      m_Outputs[idx] = output;

      if (output)
        {
        output->ConnectSource(this, idx);
        }
      this->Modified();
  
      return;
      }
    }
  
  this->SetNumberOfOutputs( static_cast<int>( m_Outputs.size() ) + 1);
  m_Outputs[ static_cast<int>( m_Outputs.size() ) - 1] = output;
  if (output)
    {
    output->ConnectSource(this, static_cast<int>( m_Outputs.size() ) - 1 );
    }
  this->Modified();
}

virtual void itk::ImageSource< TDeformationField >::AfterThreadedGenerateData ( void ) [inline, protected, virtual, inherited]

If an imaging filter needs to perform processing after all processing threads have completed, the filter can can provide an implementation for AfterThreadedGenerateData(). The execution flow in the default GenerateData() method will be: 1) Allocate the output buffer 2) Call BeforeThreadedGenerateData() 3) Spawn threads, calling ThreadedGenerateData() in each thread. 4) Call AfterThreadedGenerateData() Note that this flow of control is only available if a filter provides a ThreadedGenerateData() method and NOT a GenerateData() method.

Definition at line 264 of file itkImageSource.h.

{};

virtual void itk::InPlaceImageFilter< TDeformationField , TDeformationField >::AllocateOutputs ( ) [protected, virtual, inherited]

The GenerateData method normally allocates the buffers for all of the outputs of a filter. Since InPlaceImageFilter's can use an overwritten version of the input for its output, the output buffer should not be allocated. When possible, we graft the input to the filter to the output. If an InPlaceFilter has multiple outputs, then it would need to override this method to graft one of its outputs and allocate the remaining. If a filter is threaded (i.e. it provides an implementation of ThreadedGenerateData()), this method is called automatically. If an InPlaceFilter is not threaded (i.e. it provides an implementation of GenerateData()), then this method (or equivalent) must be called in GenerateData().

Reimplemented from itk::ImageSource< TDeformationField >.

virtual void itk::DenseFiniteDifferenceImageFilter< TDeformationField , TDeformationField >::AllocateUpdateBuffer ( ) [protected, virtual, inherited]

This method allocates storage in m_UpdateBuffer. It is called from Superclass::GenerateData().

Implements itk::FiniteDifferenceImageFilter< TDeformationField, TDeformationField >.

Reimplemented in itk::DiffeomorphicDemonsRegistrationFilter< TFixedImage, TMovingImage, TDeformationField >, and itk::FastSymmetricForcesDemonsRegistrationFilter< TFixedImage, TMovingImage, TDeformationField >.

virtual void itk::FiniteDifferenceImageFilter< TDeformationField , TDeformationField >::ApplyUpdate ( TimeStepType dt ) [protected, pure virtual, inherited]

This method is defined by a subclass to apply changes to the output from an update buffer and a time step value "dt".

Parameters:

dt

Time step value.

virtual void itk::DenseFiniteDifferenceImageFilter< TDeformationField , TDeformationField >::ApplyUpdate ( TimeStepType dt ) [protected, virtual, inherited]

This method applies changes from the m_UpdateBuffer to the output using the ThreadedApplyUpdate() method and a multithreading mechanism. "dt" is the time step to use for the update of each pixel.

template<class TFixedImage , class TMovingImage , class TDeformationField >

void otb::NCCRegistrationFilter< TFixedImage, TMovingImage, TDeformationField >::ApplyUpdate ( TimeStepType dt ) [protected, virtual]

Apply update.

Definition at line 231 of file otbNCCRegistrationFilter.txx.

{
  // If we smooth the update buffer before applying it, then the are
  // approximating a viscuous problem as opposed to an elastic problem
  /*  if ( this->GetSmoothUpdateField() )
      {
      this->SmoothUpdateField();
      }
    */
  this->Superclass::ApplyUpdate(dt);

  NCCRegistrationFunctionType *drfp =
    dynamic_cast<NCCRegistrationFunctionType *>
      (this->GetDifferenceFunction().GetPointer());

  if (!drfp)
    {
    itkExceptionMacro(<<
                      "Could not cast difference function to NCCRegistrationFunction");
    }

//  this->SetRMSChange( drfp->GetRMSChange() );

}

virtual void itk::ImageSource< TDeformationField >::BeforeThreadedGenerateData ( void ) [inline, protected, virtual, inherited]

If an imaging filter needs to perform processing after the buffer has been allocated but before threads are spawned, the filter can can provide an implementation for BeforeThreadedGenerateData(). The execution flow in the default GenerateData() method will be: 1) Allocate the output buffer 2) Call BeforeThreadedGenerateData() 3) Spawn threads, calling ThreadedGenerateData() in each thread. 4) Call AfterThreadedGenerateData() Note that this flow of control is only available if a filter provides a ThreadedGenerateData() method and NOT a GenerateData() method.

Definition at line 252 of file itkImageSource.h.

{};

void itk::LightObject::BreakOnError ( ) [static, inherited]

This method is called when itkExceptionMacro executes. It allows the debugger to break on error.

Definition at line 149 of file itkLightObject.cxx.

{
  ;  
}

void itk::ProcessObject::CacheInputReleaseDataFlags ( ) [protected, virtual, inherited]

Cache the state of any ReleaseDataFlag's on the inputs. While the filter is executing, we need to set the ReleaseDataFlag's on the inputs to false in case the current filter is implemented using a mini-pipeline (which will try to release the inputs). After the filter finishes, we restore the state of the ReleaseDataFlag's before the call to ReleaseInputs().

Definition at line 1050 of file itkProcessObject.cxx.

{
  unsigned int idx;
  
  m_CachedInputReleaseDataFlags.resize( m_Inputs.size() );
  for (idx = 0; idx < m_Inputs.size(); ++idx)
    {
    if (m_Inputs[idx])
      {
      m_CachedInputReleaseDataFlags[idx]=m_Inputs[idx]->GetReleaseDataFlag();
      m_Inputs[idx]->ReleaseDataFlagOff();
      }
    else
      {
      m_CachedInputReleaseDataFlags[idx] = false;
      }
    }
}

virtual TimeStepType itk::DenseFiniteDifferenceImageFilter< TDeformationField , TDeformationField >::CalculateChange ( ) [protected, virtual, inherited]

This method populates an update buffer with changes for each pixel in the output using the ThreadedCalculateChange() method and a multithreading mechanism. Returns value is a time step to be used for the update.

Implements itk::FiniteDifferenceImageFilter< TDeformationField, TDeformationField >.

virtual void itk::ImageToImageFilter< TDeformationField , TDeformationField >::CallCopyInputRegionToOutputRegion	(	OutputImageRegionType &	destRegion,
		const InputImageRegionType &	srcRegion
	)			`[protected, virtual, inherited]`

This function calls the actual region copier to do the mapping from input image space to output image space. It uses a Function object used for dispatching to various routines to copy an input region (start index and size) to an output region. For most filters, this is a trivial copy because most filters require the input dimension to match the output dimension. However, some filters like itk::UnaryFunctorImageFilter can support output images of a higher dimension that the input.

This function object is used by the default implementation of GenerateOutputInformation(). It can also be used in routines like ThreadedGenerateData() where a filter may need to map an input region to an output region.

The default copier uses a "dispatch pattern" to call one of three overloaded functions depending on whether the input and output images are the same dimension, the input is a higher dimension that the output, or the input is of a lower dimension than the output. The use of an overloaded function is required for proper compilation of the various cases.

For the latter two cases, trivial implementations are used. If the input image is a higher dimension than the output, the first portion of the input region is copied to the output region. If the input region is a lower dimension than the output, the input region information is copied into the first portion of the output region and the rest of the output region is set to zero.

If a filter needs a different default behavior, it can override this method.

virtual void itk::ImageToImageFilter< TDeformationField , TDeformationField >::CallCopyOutputRegionToInputRegion	(	InputImageRegionType &	destRegion,
		const OutputImageRegionType &	srcRegion
	)			`[protected, virtual, inherited]`

This function calls the actual region copier to do the mapping from output image space to input image space. It uses a Function object used for dispatching to various routines to copy an output region (start index and size) to an input region. For most filters, this is a trivial copy because most filters require the input dimension to match the output dimension. However, some filters like itk::ExtractImageFilter can support output images of a lower dimension that the input.

This function object can be used by GenerateOutputInformation() to copy the input LargestPossibleRegion to the output LargestPossibleRegion and can also be used in GenerateData or ThreadedGenerateData() where a filter may need to map an output region to an input region.

The default copier uses a "dispatch pattern" to call one of three overloaded functions depending on whether the input and output images are the same dimension, the input is a higher dimension that the output, or the input is of a lower dimension than the output. The use of an overloaded function is required for proper compilation of the various cases.

For the latter two cases, trivial implementations are used. If the input image is a higher dimension than the output, the output region information is copied into the first portion of the input region and the rest of the input region is set to zero. If the input region is a lower dimension than the output, the first portion of the output region is copied to the input region.

If a filter needs a different default behavior, it can override this method. The ExtractImageFilter overrides this function object so that if the input image is a higher dimension than the output image, the filter can control "where" in the input image the output subimage is extracted (as opposed to mapping to first few dimensions of the input).

virtual bool itk::InPlaceImageFilter< TDeformationField , TDeformationField >::CanRunInPlace ( ) const [inline, virtual, inherited]

Can the filter run in place? To do so, the filter's first input and output must have the same dimension and pixel type. This method can be used in conjunction with the InPlace ivar to determine whether a particular use of the filter is really running in place. Some filters may be able to optimize their operation if the InPlace is true and CanRunInPlace is true. CanRunInPlace may also be overridded by InPlaceImageFilter subclasses to fine tune its behavior.

Definition at line 111 of file itkInPlaceImageFilter.h.

    {
    return (typeid(TInputImage) == typeid(TOutputImage));
    }

template<class TFixedImage , class TMovingImage , class TDeformationField >

void itk::PDEDeformableRegistrationFilter< TFixedImage, TMovingImage, TDeformationField >::CopyInputToOutput ( ) [protected, virtual, inherited]

A simple method to copy the data from the input to the output. If the input does not exist, a zero field is written to the output.

Reimplemented from itk::DenseFiniteDifferenceImageFilter< TDeformationField, TDeformationField >.

Definition at line 280 of file itkPDEDeformableRegistrationFilter.txx.

{

  typename Superclass::InputImageType::ConstPointer  inputPtr  = this->GetInput();
  
  if( inputPtr )
    {
    this->Superclass::CopyInputToOutput();
    }
  else
    {
    typename Superclass::PixelType zeros;
    for( unsigned int j = 0; j < ImageDimension; j++ )
      {
      zeros[j] = 0;
      }

    typename OutputImageType::Pointer output = this->GetOutput();
  
    ImageRegionIterator<OutputImageType> out(output, output->GetRequestedRegion());

    while( ! out.IsAtEnd() )
      {
      out.Value() =  zeros;
      ++out;
      }
    }
}

LightObject::Pointer itk::Object::CreateAnother ( void ) const [virtual, inherited]

Create an object from an instance, potentially deferring to a factory. This method allows you to create an instance of an object that is exactly the same type as the referring object. This is useful in cases where an object has been cast back to a base class.

Reimplemented from itk::LightObject.

Definition at line 231 of file itkObject.cxx.

References itk::SmartPointer< TObjectType >::GetPointer(), and itk::Object::New().

Referenced by itk::StreamingImageIOBase::GetActualNumberOfSplitsForWriting().

{
  return Object::New().GetPointer();
}

void itk::Object::DebugOff ( ) const [virtual, inherited]

Turn debugging output off.

Definition at line 253 of file itkObject.cxx.

{
  m_Debug = false;
}

void itk::Object::DebugOn ( ) const [virtual, inherited]

Turn debugging output on.

Definition at line 242 of file itkObject.cxx.

{
  m_Debug = true;
}

void itk::LightObject::Delete ( ) [virtual, inherited]

Delete an itk object. This method should always be used to delete an object when the new operator was used to create it. Using the C delete method will not work with reference counting.

Delete a itk object. This method should always be used to delete an object when the new operator was used to create it. Using the C++ delete method will not work with reference counting.

Definition at line 88 of file itkLightObject.cxx.

{
  this->UnRegister();
}

virtual void itk::ProcessObject::EnlargeOutputRequestedRegion ( DataObject * ) [inline, virtual, inherited]

Give the process object a chance to indictate that it will produce more output than it was requested to produce. For example, many imaging filters must compute the entire output at once or can only produce output in complete slices. Such filters cannot handle smaller requested regions. These filters must provide an implementation of this method, setting the output requested region to the size they will produce. By default, a process object does not modify the size of the output requested region.

Definition at line 224 of file itkProcessObject.h.

{};

virtual void itk::FiniteDifferenceImageFilter< TDeformationField , TDeformationField >::GenerateData ( void ) [protected, virtual, inherited]

This is the default, high-level algorithm for calculating finite difference solutions. It calls virtual methods in its subclasses to implement the major steps of the algorithm.

Reimplemented from itk::ImageSource< TDeformationField >.

template<class TFixedImage , class TMovingImage , class TDeformationField >

void otb::NCCRegistrationFilter< TFixedImage, TMovingImage, TDeformationField >::GenerateInputRequestedRegion ( void ) [protected, virtual]

Update the Input requested region.

Reimplemented from itk::PDEDeformableRegistrationFilter< TFixedImage, TMovingImage, TDeformationField >.

Definition at line 159 of file otbNCCRegistrationFilter.txx.

References itk::DataObjectError::SetDataObject().

{
  // get pointers to the input and output
  typename Superclass::FixedImagePointer fixedPtr =
    const_cast<TFixedImage *>(this->GetFixedImage());
  typename Superclass::MovingImagePointer movingPtr =
    const_cast<TMovingImage *>(this->GetMovingImage());
  typename TDeformationField::Pointer outputPtr = this->GetOutput();

  if (!fixedPtr || !movingPtr || !outputPtr)
    {
    return;
    }

  // get a copy of the input requested region (should equal the output
  // requested region)
  typename TDeformationField::RegionType requestedRegion;
  requestedRegion = outputPtr->GetRequestedRegion();

  // pad the input requested region by the operator radius
  requestedRegion.PadByRadius(this->GetNCCRadius());

  // crop the input requested region at the input's largest possible region
  if (requestedRegion.Crop(fixedPtr->GetLargestPossibleRegion()))
    {
    if (requestedRegion.Crop(movingPtr->GetLargestPossibleRegion()))
      {
      fixedPtr->SetRequestedRegion(requestedRegion);
      movingPtr->SetRequestedRegion(requestedRegion);
      return;
      }
    else
      {
      // Couldn't crop the region (requested region is outside the largest
      // possible region).  Throw an exception.

      // store what we tried to request (prior to trying to crop)
      movingPtr->SetRequestedRegion(requestedRegion);

      // build an exception
      itk::InvalidRequestedRegionError e(__FILE__, __LINE__);
      e.SetLocation(ITK_LOCATION);
      e.SetDescription(
        "Requested region is (at least partially) outside the largest possible region of the moving image.");
      e.SetDataObject(movingPtr);
      throw e;

      }
    }
  else
    {
    // Couldn't crop the region (requested region is outside the largest
    // possible region).  Throw an exception.

    // store what we tried to request (prior to trying to crop)
    fixedPtr->SetRequestedRegion(requestedRegion);

    // build an exception
    itk::InvalidRequestedRegionError e(__FILE__, __LINE__);
    e.SetLocation(ITK_LOCATION);
    e.SetDescription("Requested region is (at least partially) outside the largest possible region of the fixed image.");
    e.SetDataObject(fixedPtr);
    throw e;
    }
}

template<class TFixedImage , class TMovingImage , class TDeformationField >

void itk::PDEDeformableRegistrationFilter< TFixedImage, TMovingImage, TDeformationField >::GenerateOutputInformation ( void ) [protected, virtual, inherited]

By default the output deformation field has the same Spacing, Origin and LargestPossibleRegion as the input/initial deformation field. If the initial deformation field is not set, the output information is copied from the fixed image.

Reimplemented from itk::ProcessObject.

Definition at line 313 of file itkPDEDeformableRegistrationFilter.txx.

{

  typename DataObject::Pointer output;

  if( this->GetInput(0) )
    {
    // Initial deformation field is set.
    // Copy information from initial field.
    this->Superclass::GenerateOutputInformation();

    }
  else if( this->GetFixedImage() )
    {
    // Initial deforamtion field is not set. 
    // Copy information from the fixed image.
    for (unsigned int idx = 0; idx < 
           this->GetNumberOfOutputs(); ++idx )
      {
      output = this->GetOutput(idx);
      if (output)
        {
        output->CopyInformation(this->GetFixedImage());
        }  
      }

    }

}

void itk::ProcessObject::GenerateOutputRequestedRegion ( DataObject * output ) [protected, virtual, inherited]

Given one output whose requested region has been set, how should the requested regions for the remaining outputs of the process object be set? By default, all the outputs are set to the same requested region. If a filter needs to produce different requested regions for each output, for instance an image processing filter producing several outputs at different resolutions, then that filter may override this method and set the requested regions appropriatedly.

Note that a filter producing multiple outputs of different types is required to override this method. The default implementation can only correctly handle multiple outputs of the same type.

By default we set all the output requested regions to be the same.

Reimplemented in itk::MultiResolutionPyramidImageFilter< TInputImage, TOutputImage >, itk::RecursiveMultiResolutionPyramidImageFilter< TInputImage, TOutputImage >, itk::watershed::BoundaryResolver< TPixelType, TDimension >, itk::watershed::EquivalenceRelabeler< TScalarType, TImageDimension >, itk::watershed::Relabeler< TScalarType, TImageDimension >, itk::watershed::Segmenter< TInputImage >, itk::watershed::SegmentTreeGenerator< TScalarType >, itk::watershed::Relabeler< ScalarType, itkGetStaticConstMacro(ImageDimension)>, itk::watershed::Segmenter< InputImageType >, and itk::watershed::SegmentTreeGenerator< ScalarType >.

Definition at line 846 of file itkProcessObject.cxx.

{
  DataObjectPointerArraySizeType idx;
  for (idx = 0; idx < m_Outputs.size(); ++idx)
    {
    if (m_Outputs[idx] && m_Outputs[idx] != output)
      {
      m_Outputs[idx]->SetRequestedRegion(output);
      }
    }  
}

virtual const bool& itk::ProcessObject::GetAbortGenerateData ( ) [virtual, inherited]

Get the AbortGenerateData flag for the process object. Process objects may handle premature termination of execution in different ways.

Command * itk::Object::GetCommand ( unsigned long tag ) [inherited]

Get the command associated with the given tag. NOTE: This returns a pointer to a Command, but it is safe to asign this to a Command::Pointer. Since Command inherits from LightObject, at this point in the code, only a pointer or a reference to the Command can be used.

Definition at line 414 of file itkObject.cxx.

References itk::Object::GetCommand().

Referenced by itk::Object::GetCommand().

{
  if (this->m_SubjectImplementation)
    {
    return this->m_SubjectImplementation->GetCommand(tag);
    }
  return NULL;
}

bool itk::Object::GetDebug ( ) const [inherited]

Get the value of the debug flag.

Definition at line 264 of file itkObject.cxx.

{
  return m_Debug;
}

template<class TFixedImage , class TMovingImage , class TDeformationField >

DeformationFieldType* itk::PDEDeformableRegistrationFilter< TFixedImage, TMovingImage, TDeformationField >::GetDeformationField ( void ) [inline, inherited]

Get output deformation field.

Definition at line 135 of file itkPDEDeformableRegistrationFilter.h.

    { return this->GetOutput(); }

virtual const FiniteDifferenceFunctionType ::Pointer& itk::FiniteDifferenceImageFilter< TDeformationField , TDeformationField >::GetDifferenceFunction ( ) const [virtual, inherited]

This method returns a pointer to a FiniteDifferenceFunction object that will be used by the filter to calculate updates at image pixels.

Returns:: A FiniteDifferenceObject pointer.

virtual const unsigned int& itk::FiniteDifferenceImageFilter< TDeformationField , TDeformationField >::GetElapsedIterations ( ) [virtual, inherited]

Get the number of elapsed iterations of the filter.

template<class TFixedImage , class TMovingImage , class TDeformationField >

const PDEDeformableRegistrationFilter< TFixedImage, TMovingImage, TDeformationField >::FixedImageType * itk::PDEDeformableRegistrationFilter< TFixedImage, TMovingImage, TDeformationField >::GetFixedImage ( void ) const [inherited]

Get the fixed image.

Definition at line 83 of file itkPDEDeformableRegistrationFilter.txx.

References itk::ProcessObject::GetInput().

{
  return dynamic_cast< const FixedImageType * >
    ( this->ProcessObject::GetInput( 1 ) );
}

bool itk::Object::GetGlobalWarningDisplay ( ) [static, inherited]

Get the value of the global debug output control flag.

Definition at line 381 of file itkObject.cxx.

{
  return m_GlobalWarningDisplay;
}

const InputImageType* itk::ImageToImageFilter< TDeformationField , TDeformationField >::GetInput ( unsigned int idx ) [inherited]

Method used internally for getting an input.

Reimplemented from itk::ProcessObject.

DataObjectPointerArray& itk::ProcessObject::GetInputs ( ) [inline, inherited]

Return an array with all the inputs of this process object. This is useful for tracing back in the pipeline to construct graphs etc.

Definition at line 108 of file itkProcessObject.h.

Referenced by otb::PipelineMemoryPrintCalculator::EvaluateProcessObjectPrintRecursive().

    {return m_Inputs;}

const MetaDataDictionary & itk::Object::GetMetaDataDictionary ( void ) const [inherited]

Returns:: A constant reference to this objects MetaDataDictionary.

Reimplemented in otb::ImageMetadataInterfaceBase.

Definition at line 543 of file itkObject.cxx.

{
  if(m_MetaDataDictionary==NULL)
    {
    m_MetaDataDictionary=new MetaDataDictionary;
    }
  return *m_MetaDataDictionary;
}

MetaDataDictionary & itk::Object::GetMetaDataDictionary ( void ) [inherited]

Returns:: A reference to this objects MetaDataDictionary.

Warning:: This reference may be changed.

Reimplemented in otb::ImageMetadataInterfaceBase.

Definition at line 532 of file itkObject.cxx.

Referenced by otb::RenderingImageFilter< TInputImage, TOutputImage >::BeforeThreadedGenerateData(), itk::Brains2MaskImageIO::CanReadFile(), otb::DataNode< TPrecision, VDimension, TValuePrecision >::CopyFieldList(), otb::VectorImage< TPixel, VImageDimension >::CopyInformation(), otb::Image< TPixel, VImageDimension >::CopyInformation(), otb::VectorImage< TPixel, VImageDimension >::GetImageKeywordlist(), otb::Image< TPixel, VImageDimension >::GetImageKeywordlist(), itk::GDCMImageIO::GetPatientName(), itk::GDCMImageIO::GetValueFromTag(), otb::RAMDrivenAdaptativeStreamingManager< TImage >::PrepareStreaming(), otb::VectorImage< TPixel, VImageDimension >::PrintSelf(), otb::Image< TPixel, VImageDimension >::PrintSelf(), otb::JPEG2000ImageIO::ReadImageInformation(), itk::IPLCommonImageIO::ReadImageInformation(), otb::ImageMetadataInterfaceBase::SetImage(), itk::GDCMImageIO::Write(), and itk::PolygonGroupSpatialObjectXMLFileWriter::WriteFile().

{
  if(m_MetaDataDictionary==NULL)
    {
    m_MetaDataDictionary=new MetaDataDictionary;
    }
  return *m_MetaDataDictionary;
}

template<class TFixedImage , class TMovingImage , class TDeformationField >

double otb::NCCRegistrationFilter< TFixedImage, TMovingImage, TDeformationField >::GetMetric ( ) const [virtual]

Get the metric value. The metric value is computed over the the overlapping region between the two images. This is value is only available for the previous iteration and NOT the current iteration.

Returns:: metric

Definition at line 92 of file otbNCCRegistrationFilter.txx.

{

  NCCRegistrationFunctionType *drfp =
    dynamic_cast<NCCRegistrationFunctionType *>
      (this->GetDifferenceFunction().GetPointer());

  if (!drfp)
    {
    itkExceptionMacro(<<
                      "Could not cast difference function to NCCRegistrationFunction");
    }

  return drfp->GetEnergy();

}

template<class TFixedImage , class TMovingImage , class TDeformationField >

const PDEDeformableRegistrationFilter< TFixedImage, TMovingImage, TDeformationField >::MovingImageType * itk::PDEDeformableRegistrationFilter< TFixedImage, TMovingImage, TDeformationField >::GetMovingImage ( void ) const [inherited]

Get the moving image.

Definition at line 110 of file itkPDEDeformableRegistrationFilter.txx.

References itk::ProcessObject::GetInput().

{
  return dynamic_cast< const MovingImageType * >
    ( this->ProcessObject::GetInput( 2 ) );
}

unsigned long itk::Object::GetMTime ( void ) const [virtual, inherited]

Return this objects modified time.

Return the modification for this object.

Definition at line 286 of file itkObject.cxx.

Referenced by itk::BoundingBox< TPointIdentifier, VPointDimension, TCoordRep, TPointsContainer >::GetMTime(), itk::ResampleImageFilter< TInputImage, TOutputImage, TInterpolatorPrecisionType >::GetMTime(), itk::DeformationFieldSource< TOutputImage >::GetMTime(), itk::SpatialObject< TDimension >::GetMTime(), itk::TransformToDeformationFieldSource< TOutputImage, TTransformPrecisionType >::GetMTime(), itk::VectorResampleImageFilter< TInputImage, TOutputImage, TInterpolatorPrecisionType >::GetMTime(), itk::InverseDeformationFieldImageFilter< TInputImage, TOutputImage >::GetMTime(), itk::SpatialObject< ::itk::GetMeshDimension< TMesh >::PointDimension >::GetObjectMTime(), otb::Function::StandardRenderingFunction< TPixel, TRGBPixel, TPixelRepresentationFunction, TTransferFunction >::Initialize(), itk::VTKImageExportBase::PipelineModifiedCallback(), and itk::ProcessObject::UpdateOutputInformation().

{
  return m_MTime.GetMTime();
}

MultiThreader* itk::ProcessObject::GetMultiThreader ( ) [inline, inherited]

Return the multithreader used by this class.

Definition at line 280 of file itkProcessObject.h.

    {return m_Threader;}

template<class TFixedImage, class TMovingImage, class TDeformationField>

virtual const char* otb::NCCRegistrationFilter< TFixedImage, TMovingImage, TDeformationField >::GetNameOfClass ( ) const [virtual]

Run-time type information (and related methods).

Reimplemented from itk::PDEDeformableRegistrationFilter< TFixedImage, TMovingImage, TDeformationField >.

template<class TFixedImage , class TMovingImage , class TDeformationField >

NCCRegistrationFilter< TFixedImage, TMovingImage, TDeformationField >::RadiusType otb::NCCRegistrationFilter< TFixedImage, TMovingImage, TDeformationField >::GetNCCRadius ( ) const [virtual]

Get the radius.

Returns:: radius.

Definition at line 115 of file otbNCCRegistrationFilter.txx.

References itk::FiniteDifferenceFunction< TImageType >::GetRadius().

{

  NCCRegistrationFunctionType *drfp =
    dynamic_cast<NCCRegistrationFunctionType *>
      (this->GetDifferenceFunction().GetPointer());

  if (!drfp)
    {
    itkExceptionMacro(<<
                      "Could not cast difference function to NCCRegistrationFunction");
    }

  return drfp->GetRadius();

}

DataObjectPointerArraySizeType itk::ProcessObject::GetNumberOfInputs ( ) const [inline, inherited]

Get the size of the input vector. This is merely the size of the input vector, not the number of inputs that have valid DataObject's assigned. Use GetNumberOfValidRequiredInputs() to determine how many inputs are non-null.

Definition at line 118 of file itkProcessObject.h.

Referenced by otb::PipelineMemoryPrintCalculator::EvaluateProcessObjectPrintRecursive(), otb::WaveletFilterBank< TInputImage, TOutputImage, TWaveletOperator, Wavelet::INVERSE >::GenerateInputRequestedRegion(), otb::WaveletFilterBank< TInputImage, TOutputImage, TWaveletOperator, Wavelet::INVERSE >::GenerateOutputInformation(), otb::ModulusAndDirectionImageToImageFilter< TInputImage, TInputImageDirection, TOutputImage >::GetInput(), and otb::ModulusAndDirectionImageToImageFilter< TInputImage, TInputImageDirection, TOutputImage >::GetInputDirection().

    {return m_Inputs.size();}

template<class TFixedImage , class TMovingImage , class TDeformationField >

std::vector< SmartPointer< DataObject > >::size_type itk::PDEDeformableRegistrationFilter< TFixedImage, TMovingImage, TDeformationField >::GetNumberOfValidRequiredInputs ( ) const [virtual, inherited]

Get the number of valid inputs. For PDEDeformableRegistration, this checks whether the fixed and moving images have been set. While PDEDeformableRegistration can take a third input as an initial deformation field, this input is not a required input.

Reimplemented from itk::ProcessObject.

Definition at line 123 of file itkPDEDeformableRegistrationFilter.txx.

{
  typename std::vector<SmartPointer<DataObject> >::size_type num = 0;

  if (this->GetFixedImage())
    {
    num++;
    }

  if (this->GetMovingImage())
    {
    num++;
    }
  
  return num;
}

OutputImageType* itk::ImageSource< TDeformationField >::GetOutput ( void ) [inherited]

Get the output data of this process object. The output of this function is not valid until an appropriate Update() method has been called, either explicitly or implicitly. Both the filter itself and the data object have Update() methods, and both methods update the data. Here are three ways to use GetOutput() and make sure the data is valid. In these examples, image is a pointer to some Image object, and the particular ProcessObjects involved are filters. The same examples apply to non-image (e.g. Mesh) data as well.

   anotherFilter->SetInput( someFilter->GetOutput() );
   anotherFilter->Update();

In this situation, someFilter and anotherFilter are said to constitute a pipeline.

   image = someFilter->GetOutput();
   image->Update();

   someFilter->Update();
   image = someFilter->GetOutput();

(In the above example, the two lines of code can be in either order.)

Note that Update() is not called automatically except within a pipeline as in the first example. When streaming (using a StreamingImageFilter) is activated, it may be more efficient to use a pipeline than to call Update() once for each filter in turn.

For an image, the data generated is for the requested Region, which can be set using ImageBase::SetRequestedRegion(). By default, the largest possible region is requested.

For Filters which have multiple outputs of different types, the GetOutput() method assumes the output is of OutputImageType. For the GetOutput(unsigned int) method, a dynamic_cast is performed incase the filter has outputs of different types or image types. Derived classes should have names get methods for these outputs.

OutputImageType* itk::ImageSource< TDeformationField >::GetOutput ( unsigned int idx ) [inherited]

Method used internally for getting an output.

Reimplemented from itk::ProcessObject.

DataObjectPointerArray& itk::ProcessObject::GetOutputs ( ) [inline, inherited]

Return an array with all the outputs of this process object. This is useful for tracing forward in the pipeline to contruct graphs etc.

Definition at line 133 of file itkProcessObject.h.

Referenced by otb::PipelineMemoryPrintCalculator::EvaluateProcessObjectPrintRecursive().

    { return m_Outputs; }

virtual const float& itk::ProcessObject::GetProgress ( ) [virtual, inherited]

Get the execution progress of a process object. The progress is a floating number in [0,1] with 0 meaning no progress and 1 meaning the filter has completed execution.

Referenced by itk::WatershedMiniPipelineProgressCommand::Execute(), itk::LaplacianRecursiveGaussianImageFilter< TInputImage, TOutputImage >::ReportProgress(), and itk::XMLFilterWatcher::ShowProgress().

virtual int itk::LightObject::GetReferenceCount ( ) const [inline, virtual, inherited]

Gets the reference count on this object.

Definition at line 105 of file itkLightObject.h.

    {return static_cast<int>(m_ReferenceCount);}

template<class TFixedImage , class TMovingImage , class TDeformationField >

const double* itk::PDEDeformableRegistrationFilter< TFixedImage, TMovingImage, TDeformationField >::GetStandardDeviations ( void ) [inline, inherited]

Get the Gaussian smoothing standard deviations use for smoothing the deformation field.

Definition at line 162 of file itkPDEDeformableRegistrationFilter.h.

    { return (double *) m_StandardDeviations; }

virtual UpdateBufferType* itk::DenseFiniteDifferenceImageFilter< TDeformationField , TDeformationField >::GetUpdateBuffer ( ) [inline, protected, virtual, inherited]

Method to allow subclasses to get direct access to the update buffer

Definition at line 127 of file itkDenseFiniteDifferenceImageFilter.h.

    { return m_UpdateBuffer; }

template<class TFixedImage , class TMovingImage , class TDeformationField >

const double* itk::PDEDeformableRegistrationFilter< TFixedImage, TMovingImage, TDeformationField >::GetUpdateFieldStandardDeviations ( void ) [inline, inherited]

Get the Gaussian smoothing standard deviations used for smoothing the update field.

Definition at line 181 of file itkPDEDeformableRegistrationFilter.h.

    { return (double *) m_UpdateFieldStandardDeviations; }

virtual void itk::ImageSource< TDeformationField >::GraftNthOutput	(	unsigned int	idx,
		DataObject *	output
	)			`[virtual, inherited]`

Graft the specified data object onto this ProcessObject's idx'th output. This is similar to the GraftOutput method except it allows you to specify which output is affected. The specified index must be a valid output number (less than ProcessObject::GetNumberOfOutputs()). See the GraftOutput for general usage information.

virtual void itk::ImageSource< TDeformationField >::GraftOutput ( DataObject * output ) [virtual, inherited]

Graft the specified DataObject onto this ProcessObject's output. This method grabs a handle to the specified DataObject's bulk data to used as its output's own bulk data. It also copies the region ivars (RequestedRegion, BufferedRegion, LargestPossibleRegion) and meta-data (Spacing, Origin) from the specified data object into this filter's output data object. Most importantly, however, it leaves the Source ivar untouched so the original pipeline routing is intact. This method is used when a process object is implemented using a mini-pipeline which is defined in its GenerateData() method. The usage is:

    // setup the mini-pipeline to process the input to this filter
    firstFilterInMiniPipeline->SetInput( this->GetInput() );
   
    // setup the mini-pipeline to calculate the correct regions
    // and write to the appropriate bulk data block
    lastFilterInMiniPipeline->GraftOutput( this->GetOutput() );

    // execute the mini-pipeline
    lastFilterInMiniPipeline->Update();

    // graft the mini-pipeline output back onto this filter's output.
    // this is needed to get the appropriate regions passed back.
    this->GraftOutput( lastFilterInMiniPipeline->GetOutput() );

For proper pipeline execution, a filter using a mini-pipeline must implement the GenerateInputRequestedRegion(), GenerateOutputRequestedRegion(), GenerateOutputInformation() and EnlargeOutputRequestedRegion() methods as necessary to reflect how the mini-pipeline will execute (in other words, the outer filter's pipeline mechanism must be consistent with what the mini-pipeline will do).

template<class TFixedImage , class TMovingImage , class TDeformationField >

virtual bool itk::PDEDeformableRegistrationFilter< TFixedImage, TMovingImage, TDeformationField >::Halt ( ) [inline, protected, virtual, inherited]

Supplies the halting criteria for this class of filters. The algorithm will stop after a user-specified number of iterations.

Reimplemented from itk::FiniteDifferenceImageFilter< TDeformationField, TDeformationField >.

Reimplemented in itk::LevelSetMotionRegistrationFilter< TFixedImage, TMovingImage, TDeformationField >.

Definition at line 208 of file itkPDEDeformableRegistrationFilter.h.

    {

    if ( m_StopRegistrationFlag )
      {
      return true;
      }

    return this->Superclass::Halt();
    }

bool itk::Object::HasObserver ( const EventObject & event ) const [inherited]

Return true if an observer is registered for this event.

Definition at line 467 of file itkObject.cxx.

{
  if (this->m_SubjectImplementation)
    {
    return this->m_SubjectImplementation->HasObserver(event);
    }
  return false;
}

template<class TFixedImage , class TMovingImage , class TDeformationField >

void itk::PDEDeformableRegistrationFilter< TFixedImage, TMovingImage, TDeformationField >::Initialize ( void ) [protected, virtual, inherited]

This method is called before iterating the solution.

Reimplemented from itk::FiniteDifferenceImageFilter< TDeformationField, TDeformationField >.

Definition at line 401 of file itkPDEDeformableRegistrationFilter.txx.

{
  this->Superclass::Initialize();
  m_StopRegistrationFlag = false;
}

template<class TFixedImage , class TMovingImage , class TDeformationField >

void otb::NCCRegistrationFilter< TFixedImage, TMovingImage, TDeformationField >::InitializeIteration ( ) [protected, virtual]

Initialize the state of filter and equation before each iteration.

Reimplemented from itk::PDEDeformableRegistrationFilter< TFixedImage, TMovingImage, TDeformationField >.

Definition at line 59 of file otbNCCRegistrationFilter.txx.

{

  // call the superclass  implementation
  Superclass::InitializeIteration();

  // set the gradient selection flag
  NCCRegistrationFunctionType *drfp =
    dynamic_cast<NCCRegistrationFunctionType *>
      (this->GetDifferenceFunction().GetPointer());

  if (!drfp)
    {
    itkExceptionMacro(<<
                      "Could not cast difference function to NCCRegistrationFunction");
    }

  /*
   * Smooth the deformation field

  if ( this->GetSmoothDeformationField() )
    {
    this->SmoothDeformationField();
    }
   */
}

void itk::Object::InvokeEvent ( const EventObject & event ) [inherited]

Call Execute on all the Commands observing this event id.

Definition at line 446 of file itkObject.cxx.

Referenced by itk::VTKImageExportBase::UpdateDataCallback().

{
  if (this->m_SubjectImplementation)
    {
    this->m_SubjectImplementation->InvokeEvent(event,this);
    }
}

void itk::Object::InvokeEvent ( const EventObject & event ) const [inherited]

Call Execute on all the Commands observing this event id. The actions triggered by this call doesn't modify this object.

Definition at line 457 of file itkObject.cxx.

{
  if (this->m_SubjectImplementation)
    {
    this->m_SubjectImplementation->InvokeEvent(event,this);
    }
}

virtual DataObjectPointer itk::ImageSource< TDeformationField >::MakeOutput ( unsigned int idx ) [virtual, inherited]

Make a DataObject of the correct type to used as the specified output. Every ProcessObject subclass must be able to create a DataObject that can be used as a specified output. This method is automatically called when DataObject::DisconnectPipeline() is called. DataObject::DisconnectPipeline, disconnects a data object from being an output of its current source. When the data object is disconnected, the ProcessObject needs to construct a replacement output data object so that the ProcessObject is in a valid state. So DataObject::DisconnectPipeline eventually calls ProcessObject::MakeOutput. Note that MakeOutput always returns a SmartPointer to a DataObject. If a subclass of ImageSource has multiple outputs of different types, then that class must provide an implementation of MakeOutput().

Reimplemented from itk::ProcessObject.

void itk::Object::Modified ( void ) const [virtual, inherited]

Update the modification time for this object. Many filters rely on the modification time to determine if they need to recompute their data.

Make sure this object's modified time is greater than all others.

Definition at line 297 of file itkObject.cxx.

{
  m_MTime.Modified();
  InvokeEvent( ModifiedEvent() );
}

template<class TFixedImage, class TMovingImage, class TDeformationField>

static Pointer otb::NCCRegistrationFilter< TFixedImage, TMovingImage, TDeformationField >::New ( ) [static]

Method for creation through the object factory.

Reimplemented from itk::PDEDeformableRegistrationFilter< TFixedImage, TMovingImage, TDeformationField >.

virtual void itk::ImageToImageFilter< TDeformationField , TDeformationField >::PopBackInput ( ) [virtual, inherited]

Model a stack on the input list by providing a pop back

Reimplemented from itk::ProcessObject.

template<class TFixedImage , class TMovingImage , class TDeformationField >

void itk::PDEDeformableRegistrationFilter< TFixedImage, TMovingImage, TDeformationField >::PostProcessOutput ( ) [protected, virtual, inherited]

This method is called after the solution has been generated. In this case, the filter release the memory of the internal buffers.

Reimplemented from itk::FiniteDifferenceImageFilter< TDeformationField, TDeformationField >.

Definition at line 388 of file itkPDEDeformableRegistrationFilter.txx.

{
  this->Superclass::PostProcessOutput();
  m_TempField->Initialize();
}

void itk::ProcessObject::PrepareOutputs ( ) [virtual, inherited]

An opportunity to deallocate a ProcessObject's bulk data storage. Some filters may wish to reuse existing bulk data storage to avoid unnecessary deallocation/allocation sequences. The default implementation calls Initialize() on each output. DataObject::Initialize() frees its bulk data by default.

Reimplemented in itk::WatershedImageFilter< TInputImage >.

Definition at line 864 of file itkProcessObject.cxx.

{  
  unsigned int idx;
  
  if (this->GetReleaseDataBeforeUpdateFlag())
    {
    for (idx = 0; idx < m_Outputs.size(); idx++)
      {
      if (m_Outputs[idx])
        {
        m_Outputs[idx]->PrepareForNewData(); 
        }
      }
    }
}

void itk::LightObject::Print	(	std::ostream &	os,
		Indent	indent = `0`
	)			const `[inherited]`

Cause the object to print itself out.

Avoid DLL boundary problems. This function will be common to all itk objects. It just calls the header/self/trailer virtual print methods, which can be overriden by subclasses (any itk object).

Definition at line 135 of file itkLightObject.cxx.

References itk::Indent::GetNextIndent().

Referenced by itk::operator<<(), and itk::WeakPointer< ProcessObject >::Print().

{
  this->PrintHeader(os, indent); 
  this->PrintSelf(os, indent.GetNextIndent());
  this->PrintTrailer(os, indent);
}

void itk::LightObject::PrintHeader	(	std::ostream &	os,
		Indent	indent
	)			const `[protected, virtual, inherited]`

Define a default print header for all objects.

Definition at line 320 of file itkLightObject.cxx.

{
  os << indent << this->GetNameOfClass() << " (" << this << ")\n";
}

template<class TFixedImage , class TMovingImage , class TDeformationField >

void otb::NCCRegistrationFilter< TFixedImage, TMovingImage, TDeformationField >::PrintSelf	(	std::ostream &	os,
		itk::Indent	indent
	)			const `[protected, virtual]`

Methods invoked by Print() to print information about the object including superclasses. Typically not called by the user (use Print() instead) but used in the hierarchical print process to combine the output of several classes.

Chaining method to print an object's instance variables, as well as its superclasses.

Reimplemented from itk::PDEDeformableRegistrationFilter< TFixedImage, TMovingImage, TDeformationField >.

Definition at line 46 of file otbNCCRegistrationFilter.txx.

{
  Superclass::PrintSelf(os, indent);
  os << indent << "NCC Radius: " <<
  this->GetNCCRadius() << std::endl;
}

void itk::LightObject::PrintTrailer	(	std::ostream &	os,
		Indent	indent
	)			const `[protected, virtual, inherited]`

Define a default print trailer for all objects.

Definition at line 331 of file itkLightObject.cxx.

{
}

void itk::ProcessObject::PropagateRequestedRegion ( DataObject * output ) [virtual, inherited]

Send the requested region information back up the pipeline (to the filters that preceed this one).

check flag to avoid executing forever if there is a loop

Give the subclass a chance to indicate that it will provide more data then required for the output. This can happen, for example, when a source can only produce the whole output. Although this is being called for a specific output, the source may need to enlarge all outputs.

Give the subclass a chance to define how to set the requested regions for each of its outputs, given this output's requested region. The default implementation is to make all the output requested regions the same. A subclass may need to override this method if each output is a different resolution.

Give the subclass a chance to request a larger requested region on the inputs. This is necessary when, for example, a filter requires more data at the "internal" boundaries to produce the boundary values - such as an image filter that derives a new pixel value by applying some operation to a neighborhood of surrounding original values.

Now that we know the input requested region, propagate this through all the inputs.

Reimplemented in itk::StreamingImageFilter< TInputImage, TOutputImage >, and itk::VTKImageImport< TOutputImage >.

Definition at line 766 of file itkProcessObject.cxx.

{
  if (m_Updating)
    {
    return;
    }

  this->EnlargeOutputRequestedRegion( output );

  this->GenerateOutputRequestedRegion( output );
  
  this->GenerateInputRequestedRegion();

  m_Updating = true;
  DataObjectPointerArraySizeType idx;
  for (idx = 0; idx < m_Inputs.size(); ++idx)
    {
    if (m_Inputs[idx])
      {
      m_Inputs[idx]->PropagateRequestedRegion();
      }
    }
  m_Updating = false;
}

void itk::ProcessObject::PropagateResetPipeline ( ) [protected, virtual, inherited]

Called to allocate the input array. Copies old inputs. Propagate a call to ResetPipeline() up the pipeline. Called only from DataObject.

Propagate the ResetPipeline call

Definition at line 631 of file itkProcessObject.cxx.

{
  // 
  // Reset this object.
  //
  // Clear the updating flag.
  m_Updating = 0;

  //
  // Loop through the inputs
  //
  unsigned int idx;
  DataObject::Pointer input;
  for (idx = 0; idx < m_Inputs.size(); ++idx)
    {
    if (m_Inputs[idx])
      {
      input = m_Inputs[idx];

      input->PropagateResetPipeline();
      }
    }
}

virtual void itk::ImageToImageFilter< TDeformationField , TDeformationField >::PushBackInput ( const InputImageType * image ) [virtual, inherited]

Push/Pop the input of this process object. These methods allow a filter to model its input vector as a queue or stack. These routines may not be appropriate for all filters, especially filters with different types of inputs. These routines follow the semantics of STL.

The routines are useful for applications that need to process "rolling" sets of images. For instance, if an application has 10 images and they need to run a filter on images 1, 2, 3, 4, then run the filter on images 2, 3, 4, 5, then run the filter on images 3, 4, 5, 6, the application can accomplish this by popping an input off the front of the input list and push a new image onto the back of input list. Again, this only makes sense for filters that single type of input.

Other uses are also possible. For a single input filter, pushing and popping inputs allow the application to temporarily replace an input to a filter.

void itk::ImageToImageFilter< TDeformationField , TDeformationField >::PushBackInput ( const DataObject * input ) [inline, protected, virtual, inherited]

PushBackInput(), PushFronInput() in the public section force the input to be the type expected by an ImageToImageFilter. However, these methods end of "hiding" the versions from the superclass (ProcessObject) whose arguments are DataObjects. Here, we re-expose the versions from ProcessObject to avoid warnings about hiding methods from the superclass.

Reimplemented from itk::ProcessObject.

Definition at line 248 of file itkImageToImageFilter.h.

    { Superclass::PushBackInput(input); }

void itk::Object::Register ( ) const [virtual, inherited]

Increase the reference count (mark as used by another object).

Reimplemented from itk::LightObject.

Definition at line 309 of file itkObject.cxx.

Referenced by itk::ObjectFactoryBase::RegisterFactory().

{
  itkDebugMacro(<< "Registered, "
                << "ReferenceCount = " << (m_ReferenceCount+1));

  // call the parent
  Superclass::Register();
}

virtual void itk::InPlaceImageFilter< TDeformationField , TDeformationField >::ReleaseInputs ( ) [protected, virtual, inherited]

InPlaceImageFilter may transfer ownership of the input bulk data to the output object. Once the output object owns the bulk data (done in AllocateOutputs()), the input object must release its hold on the bulk data. ProcessObject::ReleaseInputs() only releases the input bulk data when the user has set the ReleaseDataFlag. InPlaceImageFilter::ReleaseInputs() also releases the input that it has overwritten.

See also:: ProcessObject::ReleaseInputs()

Reimplemented from itk::ProcessObject.

void itk::Object::RemoveAllObservers ( ) [inherited]

Remove all observers .

Definition at line 435 of file itkObject.cxx.

{
  if (this->m_SubjectImplementation)
    {
    this->m_SubjectImplementation->RemoveAllObservers();
    }
}

void itk::ProcessObject::RemoveInput ( DataObject * input ) [protected, virtual, inherited]

Remove an input.

Removes the first occurence of the given OutputObject from the inputs to this ProcessObject. If it's the last object on the list, shortens the list.

Definition at line 173 of file itkProcessObject.cxx.

{
  if (!input)
    {
    return;
    }
  
  // find the input in the list of inputs
  DataObjectPointerArray::iterator pos = 
    std::find(m_Inputs.begin(), m_Inputs.end(), input);

  if(pos == m_Inputs.end())
    {
    itkDebugMacro("tried to remove an input that was not in the list");
    return;
    }

  // Set the position in the m_Inputs containing input to 0
  *pos = 0;

  // if that was the last input, then shrink the list
  if (pos == m_Inputs.end() - 1 )
    {
    this->SetNumberOfInputs( static_cast<int>( m_Inputs.size() ) - 1);
    }

  this->Modified();
}

void itk::Object::RemoveObserver ( unsigned long tag ) [inherited]

Remove the observer with this tag value.

Definition at line 425 of file itkObject.cxx.

Referenced by otb::StreamingImageFileWriter< TInputImage >::UpdateOutputData().

{
  if (this->m_SubjectImplementation)
    {
    this->m_SubjectImplementation->RemoveObserver(tag);
    }
}

void itk::ProcessObject::ResetPipeline ( ) [virtual, inherited]

Reset the pipeline. If an exception is thrown during an Update(), the pipeline may be in an inconsistent state. This method clears the internal state of the pipeline so Update() can be called.

Definition at line 621 of file itkProcessObject.cxx.

{
  if (this->GetOutput(0))
    {
    this->GetOutput(0)->ResetPipeline();
    }
}

virtual TimeStepType itk::FiniteDifferenceImageFilter< TDeformationField , TDeformationField >::ResolveTimeStep	(	const TimeStepType *	timeStepList,
		const bool *	valid,
		int	size
	)			`[protected, virtual, inherited]`

Virtual method for resolving a single time step from a set of time steps returned from processing threads.

Returns:: Time step (dt) for the iteration update based on a list of time steps generated from the threaded calculated change method (one for each region processed).

Parameters:

	timeStepList	The set of time changes compiled from all the threaded calls to ThreadedGenerateData.
	valid	The set of flags indicating which of "list" elements are valid
	size	The size of "list" and "valid"

The default is to return the minimum value in the list.

void itk::ProcessObject::RestoreInputReleaseDataFlags ( ) [protected, virtual, inherited]

Restore the cached input ReleaseDataFlags.

Definition at line 1075 of file itkProcessObject.cxx.

{
  unsigned int idx;
  
  for (idx = 0;
       idx < m_Inputs.size() && idx < m_CachedInputReleaseDataFlags.size();
       ++idx)
    {
    if (m_Inputs[idx])
      {
      m_Inputs[idx]->SetReleaseDataFlag(m_CachedInputReleaseDataFlags[idx]);
      }
    }
}

virtual void itk::ProcessObject::SetAbortGenerateData ( bool _arg ) [virtual, inherited]

Set the AbortGenerateData flag for the process object. Process objects may handle premature termination of execution in different ways.

void itk::Object::SetDebug ( bool debugFlag ) const [inherited]

Set the value of the debug flag. A non-zero value turns debugging on.

Definition at line 275 of file itkObject.cxx.

{
  m_Debug = debugFlag;
}

virtual void itk::FiniteDifferenceImageFilter< TDeformationField , TDeformationField >::SetDifferenceFunction ( FiniteDifferenceFunctionType * _arg ) [virtual, inherited]

This method sets the pointer to a FiniteDifferenceFunction object that will be used by the filter to calculate updates at image pixels.

Returns:: A FiniteDifferenceObject pointer.

virtual void itk::FiniteDifferenceImageFilter< TDeformationField , TDeformationField >::SetElapsedIterations ( unsigned int _arg ) [protected, virtual, inherited]

Set the number of elapsed iterations of the filter.

template<class TFixedImage , class TMovingImage , class TDeformationField >

void itk::PDEDeformableRegistrationFilter< TFixedImage, TMovingImage, TDeformationField >::SetFixedImage ( const FixedImageType * ptr ) [inherited]

Set the fixed image.

Definition at line 69 of file itkPDEDeformableRegistrationFilter.txx.

References itk::ProcessObject::SetNthInput().

{
  this->ProcessObject::SetNthInput( 1, const_cast< FixedImageType * >( ptr ) );
}

void itk::Object::SetGlobalWarningDisplay ( bool val ) [static, inherited]

This is a global flag that controls whether any debug, warning or error messages are displayed.

Set the value of the global debug output control flag.

Definition at line 370 of file itkObject.cxx.

{
  m_GlobalWarningDisplay = val;
}

template<class TFixedImage , class TMovingImage , class TDeformationField >

void itk::PDEDeformableRegistrationFilter< TFixedImage, TMovingImage, TDeformationField >::SetInitialDeformationField ( const DeformationFieldType * ptr ) [inline, inherited]

Set initial deformation field.

Definition at line 131 of file itkPDEDeformableRegistrationFilter.h.

  { this->SetInput( ptr ); }

virtual void itk::InPlaceImageFilter< TDeformationField , TDeformationField >::SetInPlace ( bool _arg ) [virtual, inherited]

In place operation can be turned on and off. Asking for in-place operation, i.e. calling SetInplace(true) or InplaceOn(), will be effective only if CanRunInPlace also returns true. By default CanRunInPlace checks whether the input and output image type match.

virtual void itk::ImageToImageFilter< TDeformationField , TDeformationField >::SetInput ( const InputImageType * image ) [virtual, inherited]

Set/Get the image input of this process object.

virtual void itk::FiniteDifferenceImageFilter< TDeformationField , TDeformationField >::SetManualReinitialization ( bool _arg ) [virtual, inherited]

Require the filter to be manually reinitialized (by calling SetStateToUninitialized()

template<class TFixedImage , class TMovingImage , class TDeformationField >

virtual void itk::PDEDeformableRegistrationFilter< TFixedImage, TMovingImage, TDeformationField >::SetMaximumError ( double _arg ) [virtual, inherited]

Set/Get the desired maximum error of the Guassian kernel approximate.

See also:: GaussianOperator.

template<class TFixedImage , class TMovingImage , class TDeformationField >

virtual void itk::PDEDeformableRegistrationFilter< TFixedImage, TMovingImage, TDeformationField >::SetMaximumKernelWidth ( unsigned int _arg ) [virtual, inherited]

Set/Get the desired limits of the Gaussian kernel width.

See also:: GaussianOperator.

virtual void itk::FiniteDifferenceImageFilter< TDeformationField , TDeformationField >::SetMaximumRMSError ( double _arg ) [virtual, inherited]

Set/Get the maximum error allowed in the solution. This may not be defined for all solvers and its meaning may change with the application.

void itk::Object::SetMetaDataDictionary ( const MetaDataDictionary & rhs ) [inherited]

Returns:: Set the MetaDataDictionary

Reimplemented in otb::ImageMetadataInterfaceBase.

Definition at line 554 of file itkObject.cxx.

Referenced by otb::VectorImage< TPixel, VImageDimension >::CopyInformation(), and otb::Image< TPixel, VImageDimension >::CopyInformation().

{
  if(m_MetaDataDictionary==NULL)
    {
    m_MetaDataDictionary=new MetaDataDictionary;
    }
  *m_MetaDataDictionary=rhs;
}

template<class TFixedImage , class TMovingImage , class TDeformationField >

void itk::PDEDeformableRegistrationFilter< TFixedImage, TMovingImage, TDeformationField >::SetMovingImage ( const MovingImageType * ptr ) [inherited]

Set the moving image.

Definition at line 96 of file itkPDEDeformableRegistrationFilter.txx.

References itk::ProcessObject::SetNthInput().

{
  this->ProcessObject::SetNthInput( 2, const_cast< MovingImageType * >( ptr ) );
}

template<class TFixedImage , class TMovingImage , class TDeformationField >

void otb::NCCRegistrationFilter< TFixedImage, TMovingImage, TDeformationField >::SetNCCRadius ( RadiusType radius ) [virtual]

Set the NCC radius

Parameters:

radius

Definition at line 138 of file otbNCCRegistrationFilter.txx.

References itk::FiniteDifferenceFunction< TImageType >::SetRadius().

{

  NCCRegistrationFunctionType *drfp =
    dynamic_cast<NCCRegistrationFunctionType *>
      (this->GetDifferenceFunction().GetPointer());

  if (!drfp)
    {
    itkExceptionMacro(<<
                      "Could not cast difference function to NCCRegistrationFunction");
    }

  drfp->SetRadius(radius);

  this->Modified();
}

void itk::ProcessObject::SetNthInput	(	unsigned int	idx,
		DataObject *	input
	)			`[protected, virtual, inherited]`

Protected methods for setting inputs. Subclasses make use of them for setting input.

Set an Input of this filter. This method does Register()/UnRegister() manually to deal with the fact that smart pointers aren't around to do the reference counting.

Reimplemented in itk::ImageToVectorImageFilter< TInputImage >.

Definition at line 210 of file itkProcessObject.cxx.

{
  // does this change anything?
  if ( idx < m_Inputs.size() && m_Inputs[idx] == input )
    {
    return;
    }
  
  // Expand array if necessary.
  if (idx >= m_Inputs.size())
    {
    this->SetNumberOfInputs(idx + 1);
    }
  
  m_Inputs[idx] = input;

  this->Modified();
}

void itk::ProcessObject::SetNthOutput	(	unsigned int	idx,
		DataObject *	output
	)			`[protected, virtual, inherited]`

Protected methods for setting outputs. Subclasses make use of them for getting output.

Set an output of this filter. This method specifically does not do a Register()/UnRegister() because of the desire to break the reference counting loop.

Definition at line 333 of file itkProcessObject.cxx.

References itk::DataObject::ConnectSource().

Referenced by itk::CoreAtomImageToDistanceMatrixProcess< TSourceImage >::CoreAtomImageToDistanceMatrixProcess(), itk::CoreAtomImageToUnaryCorrespondenceMatrixProcess< TSourceImage >::CoreAtomImageToUnaryCorrespondenceMatrixProcess(), itk::SimplexMeshToTriangleMeshFilter< TInputMesh, TOutputMesh >::CreateTriangles(), otb::FillGapsFilter::FillGapsFilter(), itk::ImageRegistrationMethod< TFixedImage, TMovingImage >::ImageRegistrationMethod(), itk::ImageToSpatialObjectRegistrationMethod< TFixedImage, TMovingSpatialObject >::ImageToSpatialObjectRegistrationMethod(), itk::SimplexMeshAdaptTopologyFilter< TInputMesh, TOutputMesh >::InsertNewCells(), itk::MedialNodePairCorrespondenceProcess< TSourceImage >::MedialNodePairCorrespondenceProcess(), itk::MedialNodeTripletCorrespondenceProcess< TSourceImage >::MedialNodeTripletCorrespondenceProcess(), itk::MultiResolutionImageRegistrationMethod< TFixedImage, TMovingImage >::MultiResolutionImageRegistrationMethod(), itk::MultiScaleHessianBasedMeasureImageFilter< TInputImage, THessianImage, TOutputImage >::MultiScaleHessianBasedMeasureImageFilter(), itk::PointSetToImageRegistrationMethod< TFixedPointSet, TMovingImage >::PointSetToImageRegistrationMethod(), itk::PointSetToPointSetRegistrationMethod< TFixedPointSet, TMovingPointSet >::PointSetToPointSetRegistrationMethod(), otb::ProspectModel::ProspectModel(), itk::RegularSphereMeshSource< TOutputMesh >::RegularSphereMeshSource(), otb::SailModel::SailModel(), itk::watershed::BoundaryResolver< TPixelType, TDimension >::SetEquivalencyTable(), itk::watershed::Relabeler< ScalarType, itkGetStaticConstMacro(ImageDimension)>::SetOutputImage(), itk::watershed::EquivalenceRelabeler< TScalarType, TImageDimension >::SetOutputImage(), itk::SphereMeshSource< TOutputMesh >::SphereMeshSource(), otb::WaveletFilterBank< TInputImage, TOutputImage, TWaveletOperator, Wavelet::FORWARD >::WaveletFilterBank(), otb::WaveletPacketTransform< TInputImage, TOutputImage, TFilter, Wavelet::FORWARD, TCost >::WaveletPacketTransform(), otb::WaveletTransform< TInputImage, TOutputImage, TFilter, Wavelet::INVERSE >::WaveletTransform(), and otb::WaveletTransform< TInputImage, TOutputImage, TFilter, Wavelet::FORWARD >::WaveletTransform().

{
  // does this change anything?
  if ( idx < m_Outputs.size() && output == m_Outputs[idx])
    {
    return;
    }

  // Expand array if necessary.
  if (idx >= m_Outputs.size())
    {
    this->SetNumberOfOutputs(idx + 1);
    }

  // Keep a handle to the original output and disconnect the old output from
  // the pipeline
  DataObjectPointer oldOutput;
  if ( m_Outputs[idx] )
    {
    oldOutput = m_Outputs[idx];
    m_Outputs[idx]->DisconnectSource(this, idx);
    }

  if (output)
    {
    output->ConnectSource(this, idx);
    }
  // save the current reference (which releases the previous reference)
  m_Outputs[idx] = output;

  // if we are clearing an output, we need to create a new blank output
  // so we are prepared for the next Update(). this copies the requested
  // region ivar
  if (!m_Outputs[idx])
    {
    itkDebugMacro( " creating new output object." );
    DataObjectPointer newOutput = this->MakeOutput(idx);
    this->SetNthOutput(idx, newOutput);

    // If we had an output object before, copy the requested region
    // ivars and release data flag to the the new output
    if (oldOutput)
      {
      newOutput->SetRequestedRegion( oldOutput );
      newOutput->SetReleaseDataFlag( oldOutput->GetReleaseDataFlag() );
      }
    }

  this->Modified();
}

void itk::ProcessObject::SetNumberOfInputs ( unsigned int num ) [protected, inherited]

Called to allocate the input array. Copies old inputs.

Called by constructor to set up input array.

Definition at line 96 of file itkProcessObject.cxx.

Referenced by otb::WaveletFilterBank< TInputImage, TOutputImage, TWaveletOperator, Wavelet::FORWARD >::WaveletFilterBank(), otb::WaveletFilterBank< TInputImage, TOutputImage, TWaveletOperator, Wavelet::INVERSE >::WaveletFilterBank(), otb::WaveletPacketTransform< TInputImage, TOutputImage, TFilter, Wavelet::FORWARD, TCost >::WaveletPacketTransform(), otb::WaveletTransform< TInputImage, TOutputImage, TFilter, Wavelet::INVERSE >::WaveletTransform(), and otb::WaveletTransform< TInputImage, TOutputImage, TFilter, Wavelet::FORWARD >::WaveletTransform().

{
  // in case nothing has changed.
  if (num == m_Inputs.size())
    {
    return;
    }
  m_Inputs.resize(num);
  this->Modified();
}

virtual void itk::FiniteDifferenceImageFilter< TDeformationField , TDeformationField >::SetNumberOfIterations ( unsigned int _arg ) [virtual, inherited]

Set/Get the number of iterations that the filter will run.

void itk::ProcessObject::SetNumberOfOutputs ( unsigned int num ) [protected, inherited]

Called to allocate the output array. Copies old outputs.

Called by constructor to set up output array.

Definition at line 424 of file itkProcessObject.cxx.

Referenced by otb::WaveletFilterBank< TInputImage, TOutputImage, TWaveletOperator, Wavelet::FORWARD >::WaveletFilterBank(), otb::WaveletFilterBank< TInputImage, TOutputImage, TWaveletOperator, Wavelet::INVERSE >::WaveletFilterBank(), otb::WaveletPacketTransform< TInputImage, TOutputImage, TFilter, Wavelet::FORWARD, TCost >::WaveletPacketTransform(), otb::WaveletTransform< TInputImage, TOutputImage, TFilter, Wavelet::INVERSE >::WaveletTransform(), and otb::WaveletTransform< TInputImage, TOutputImage, TFilter, Wavelet::FORWARD >::WaveletTransform().

{

  // in case nothing has changed.
  if (num == m_Outputs.size())
    {
    return;
    }
  m_Outputs.resize(num);
  this->Modified();
}

virtual void itk::ProcessObject::SetNumberOfThreads ( int _arg ) [virtual, inherited]

Get/Set the number of threads to create when executing.

Referenced by itk::CannyEdgeDetectionImageFilter< TInputImage, TOutputImage >::Compute2ndDerivative(), itk::CannyEdgeDetectionImageFilter< TInputImage, TOutputImage >::Compute2ndDerivativePos(), itk::FiniteDifferenceSparseImageFilter< TInputImageType, TSparseOutputImageType >::PrecalculateChange(), and otb::WaveletFilterBank< TInputImage, TOutputImage, TWaveletOperator, Wavelet::INVERSE >::WaveletFilterBank().

virtual void itk::ProcessObject::SetProgress ( float _arg ) [virtual, inherited]

Set the execution progress of a process object. The progress is a floating number in [0,1] with 0 meaning no progress and 1 meaning the filter has completed execution. The ProgressEvent is NOT invoked.

void itk::Object::SetReferenceCount ( int ref ) [virtual, inherited]

Sets the reference count (use with care)

If there is a delete method, invoke it.

Reimplemented from itk::LightObject.

Definition at line 347 of file itkObject.cxx.

{
  itkDebugMacro(<< "Reference Count set to " << ref);

  // ReferenceCount in now unlocked.  We may have a race condition to
  // to delete the object.
  if( ref <= 0 )
    {
    this->InvokeEvent(DeleteEvent());
    }

  Superclass::SetReferenceCount(ref);
}

virtual void itk::ProcessObject::SetReleaseDataBeforeUpdateFlag ( bool _arg ) [virtual, inherited]

Turn on/off the flags to control whether the bulk data belonging to the outputs of this ProcessObject are released/reallocated during an Update(). In limited memory scenarios, a user may want to force the elements of a pipeline to release any bulk data that is going to be regenerated anyway during an Update() in order to control peak memory allocation. Note that this flag is different from the ReleaseDataFlag. ReleaseDataFlag manages the deallocation of a ProcessObject's bulk output data once that data has been consumed by a downstream ProcessObject. The ReleaseDataBeforeUpdateFlag manages the deallocation/reallocation of bulk data during a pipeline update to control peak memory utilization. Default value is on.

void itk::ProcessObject::SetReleaseDataFlag ( bool flag ) [virtual, inherited]

Turn on/off the flags to control whether the bulk data belonging to the outputs of this ProcessObject are released after being used by a downstream ProcessObject. Default value is off. Another options for controlling memory utilization is the ReleaseDataBeforeUpdateFlag.

Definition at line 528 of file itkProcessObject.cxx.

{
  unsigned int idx;
  
  for (idx = 0; idx < m_Outputs.size(); idx++)
    {
    if (m_Outputs[idx])
      {
      m_Outputs[idx]->SetReleaseDataFlag(val);
      }
    }
}

virtual void itk::FiniteDifferenceImageFilter< TDeformationField , TDeformationField >::SetRMSChange ( double _arg ) [virtual, inherited]

Set/Get the root mean squared change of the previous iteration. May not be used by all solvers.

template<class TFixedImage , class TMovingImage , class TDeformationField >

virtual void itk::PDEDeformableRegistrationFilter< TFixedImage, TMovingImage, TDeformationField >::SetSmoothDeformationField ( bool _arg ) [virtual, inherited]

Set/Get whether the deformation field is smoothed (regularized). Smoothing the deformation yields a solution elastic in nature. If SmoothDeformationField is on, then the deformation field is smoothed with a Gaussian whose standard deviations are specified with SetStandardDeviations()

template<class TFixedImage , class TMovingImage , class TDeformationField >

virtual void itk::PDEDeformableRegistrationFilter< TFixedImage, TMovingImage, TDeformationField >::SetSmoothUpdateField ( bool _arg ) [virtual, inherited]

Set/Get whether the update field is smoothed (regularized). Smoothing the update field yields a solution viscous in nature. If SmoothUpdateField is on, then the update field is smoothed with a Gaussian whose standard deviations are specified with SetUpdateFieldStandardDeviations()

template<class TFixedImage , class TMovingImage , class TDeformationField >

void itk::PDEDeformableRegistrationFilter< TFixedImage, TMovingImage, TDeformationField >::SetStandardDeviations ( double value ) [virtual, inherited]

Set the standard deviations.

Definition at line 147 of file itkPDEDeformableRegistrationFilter.txx.

{

  unsigned int j;
  for( j = 0; j < ImageDimension; j++ )
    {
    if( value != m_StandardDeviations[j] )
      {
      break;
      }
    }
  if( j < ImageDimension )
    {
    this->Modified();
    for( j = 0; j < ImageDimension; j++ )
      {
      m_StandardDeviations[j] = value;
      }
    }

}

template<class TFixedImage , class TMovingImage , class TDeformationField >

virtual void itk::PDEDeformableRegistrationFilter< TFixedImage, TMovingImage, TDeformationField >::SetStandardDeviations ( double data[] ) [virtual, inherited]

Set the Gaussian smoothing standard deviations for the deformation field. The values are set with respect to pixel coordinates.

virtual void itk::FiniteDifferenceImageFilter< TDeformationField , TDeformationField >::SetState ( FilterStateType _arg ) [virtual, inherited]

Set/Get the state of the filter.

void itk::FiniteDifferenceImageFilter< TDeformationField , TDeformationField >::SetStateToInitialized ( ) [inline, inherited]

Set the state of the filter to INITIALIZED

Definition at line 198 of file itkFiniteDifferenceImageFilter.h.

    {
    this->SetState(INITIALIZED);
    }

void itk::FiniteDifferenceImageFilter< TDeformationField , TDeformationField >::SetStateToUninitialized ( ) [inline, inherited]

Set the state of the filter to UNINITIALIZED

Definition at line 204 of file itkFiniteDifferenceImageFilter.h.

    {
    this->SetState(UNINITIALIZED);
    }

template<class TFixedImage , class TMovingImage , class TDeformationField >

virtual void itk::PDEDeformableRegistrationFilter< TFixedImage, TMovingImage, TDeformationField >::SetUpdateFieldStandardDeviations ( double data[] ) [virtual, inherited]

Set the Gaussian smoothing standard deviations for the update field. The values are set with respect to pixel coordinates.

virtual void itk::FiniteDifferenceImageFilter< TDeformationField , TDeformationField >::SetUseImageSpacing ( bool _arg ) [virtual, inherited]

Use the image spacing information in calculations. Use this option if you want derivatives in physical space. Default is UseImageSpacingOff.

template<class TFixedImage , class TMovingImage , class TDeformationField >

void itk::PDEDeformableRegistrationFilter< TFixedImage, TMovingImage, TDeformationField >::SmoothDeformationField ( ) [protected, virtual, inherited]

Utility to smooth the deformation field (represented in the Output) using a Guassian operator. The amount of smoothing can be specified by setting the StandardDeviations.

Definition at line 414 of file itkPDEDeformableRegistrationFilter.txx.

{

  DeformationFieldPointer field = this->GetOutput();

  // copy field to TempField
  m_TempField->SetOrigin( field->GetOrigin() );
  m_TempField->SetSpacing( field->GetSpacing() );
  m_TempField->SetDirection( field->GetDirection() );
  m_TempField->SetLargestPossibleRegion( 
    field->GetLargestPossibleRegion() );
  m_TempField->SetRequestedRegion(
    field->GetRequestedRegion() );
  m_TempField->SetBufferedRegion( field->GetBufferedRegion() );
  m_TempField->Allocate();
  
  typedef typename DeformationFieldType::PixelType    VectorType;
  typedef typename VectorType::ValueType              ScalarType;
  typedef GaussianOperator<ScalarType,ImageDimension> OperatorType;
  typedef VectorNeighborhoodOperatorImageFilter<
    DeformationFieldType,
    DeformationFieldType>                             SmootherType;

  OperatorType * oper = new OperatorType;
  typename SmootherType::Pointer smoother = SmootherType::New();

  typedef typename DeformationFieldType::PixelContainerPointer 
    PixelContainerPointer;
  PixelContainerPointer swapPtr;

  // graft the output field onto the mini-pipeline
  smoother->GraftOutput( m_TempField );

  for( unsigned int j = 0; j < ImageDimension; j++ )
    {
    // smooth along this dimension
    oper->SetDirection( j );
    double variance = vnl_math_sqr( m_StandardDeviations[j] );
    oper->SetVariance( variance );
    oper->SetMaximumError( m_MaximumError );
    oper->SetMaximumKernelWidth( m_MaximumKernelWidth );
    oper->CreateDirectional();

    // todo: make sure we only smooth within the buffered region
    smoother->SetOperator( *oper );
    smoother->SetInput( field );
    smoother->Update();

    if ( j < ImageDimension - 1 )
      {
      // swap the containers
      swapPtr = smoother->GetOutput()->GetPixelContainer();
      smoother->GraftOutput( field );
      field->SetPixelContainer( swapPtr );
      smoother->Modified();
      }

    }

  // graft the output back to this filter
  m_TempField->SetPixelContainer( field->GetPixelContainer() );
  this->GraftOutput( smoother->GetOutput() );

  delete oper;

}

template<class TFixedImage , class TMovingImage , class TDeformationField >

void itk::PDEDeformableRegistrationFilter< TFixedImage, TMovingImage, TDeformationField >::SmoothUpdateField ( ) [protected, virtual, inherited]

Utility to smooth the UpdateBuffer using a Gaussian operator. The amount of smoothing can be specified by setting the UpdateFieldStandardDeviations.

Definition at line 487 of file itkPDEDeformableRegistrationFilter.txx.

{
  // The update buffer will be overwritten with new data.
  DeformationFieldPointer field = this->GetUpdateBuffer();

  typedef typename DeformationFieldType::PixelType    VectorType;
  typedef typename VectorType::ValueType              ScalarType;
  typedef GaussianOperator<ScalarType,ImageDimension> OperatorType;
  typedef VectorNeighborhoodOperatorImageFilter<
    DeformationFieldType,
    DeformationFieldType>                             SmootherType;
  
  OperatorType opers[ImageDimension];
  typename SmootherType::Pointer smoothers[ImageDimension];

  for( unsigned int j = 0; j < ImageDimension; j++ )
    {
    // smooth along this dimension
    opers[j].SetDirection( j );
    double variance = vnl_math_sqr( this->GetUpdateFieldStandardDeviations()[j] );
    //double variance = vnl_math_sqr( 1.0 );
    opers[j].SetVariance( variance );
    opers[j].SetMaximumError( this->GetMaximumError() );
    opers[j].SetMaximumKernelWidth( this->GetMaximumKernelWidth() );
    opers[j].CreateDirectional();

    smoothers[j] = SmootherType::New();
    smoothers[j]->SetOperator( opers[j] );
    smoothers[j]->ReleaseDataFlagOn();

    if (j > 0)
      {
      smoothers[j]->SetInput( smoothers[j-1]->GetOutput() );
      }
    }
  smoothers[0]->SetInput( field );
  smoothers[ImageDimension-1]->GetOutput()
    ->SetRequestedRegion( field->GetBufferedRegion() );

  smoothers[ImageDimension-1]->Update();
  
  // field to contain the final smoothed data, do the equivalent of a graft
  field->SetPixelContainer( smoothers[ImageDimension-1]->GetOutput()
                            ->GetPixelContainer() );
  field->SetRequestedRegion( smoothers[ImageDimension-1]->GetOutput()
                             ->GetRequestedRegion() );
  field->SetBufferedRegion( smoothers[ImageDimension-1]->GetOutput()
                            ->GetBufferedRegion() );
  field->SetLargestPossibleRegion( smoothers[ImageDimension-1]->GetOutput()
                                   ->GetLargestPossibleRegion() );
  field->CopyInformation( smoothers[ImageDimension-1]->GetOutput() );
}

virtual int itk::ImageSource< TDeformationField >::SplitRequestedRegion	(	int	i,
		int	num,
		OutputImageRegionType &	splitRegion
	)			`[protected, virtual, inherited]`

Split the output's RequestedRegion into "num" pieces, returning region "i" as "splitRegion". This method is called "num" times. The regions must not overlap. The method returns the number of pieces that the routine is capable of splitting the output RequestedRegion, i.e. return value is less than or equal to "num".

template<class TFixedImage , class TMovingImage , class TDeformationField >

virtual void itk::PDEDeformableRegistrationFilter< TFixedImage, TMovingImage, TDeformationField >::StopRegistration ( void ) [inline, virtual, inherited]

Stop the registration after the current iteration.

Definition at line 187 of file itkPDEDeformableRegistrationFilter.h.

    { m_StopRegistrationFlag = true; }

virtual void itk::DenseFiniteDifferenceImageFilter< TDeformationField , TDeformationField >::ThreadedApplyUpdate	(	TimeStepType	dt,
		const ThreadRegionType &	regionToProcess,
		int	threadId
	)			`[protected, virtual, inherited]`

Does the actual work of updating the output from the UpdateContainer over an output region supplied by the multithreading mechanism.

See also:: ApplyUpdate; ApplyUpdateThreaderCallback

virtual TimeStepType itk::DenseFiniteDifferenceImageFilter< TDeformationField , TDeformationField >::ThreadedCalculateChange	(	const ThreadRegionType &	regionToProcess,
		int	threadId
	)			`[protected, virtual, inherited]`

Does the actual work of calculating change over a region supplied by the multithreading mechanism.

See also:: CalculateChange; CalculateChangeThreaderCallback

virtual void itk::ImageSource< TDeformationField >::ThreadedGenerateData	(	const OutputImageRegionType &	outputRegionForThread,
		int	threadId
	)			`[protected, virtual, inherited]`

If an imaging filter can be implemented as a multithreaded algorithm, the filter will provide an implementation of ThreadedGenerateData(). This superclass will automatically split the output image into a number of pieces, spawn multiple threads, and call ThreadedGenerateData() in each thread. Prior to spawning threads, the BeforeThreadedGenerateData() method is called. After all the threads have completed, the AfterThreadedGenerateData() method is called. If an image processing filter cannot support threading, that filter should provide an implementation of the GenerateData() method instead of providing an implementation of ThreadedGenerateData(). If a filter provides a GenerateData() method as its implementation, then the filter is responsible for allocating the output data. If a filter provides a ThreadedGenerateData() method as its implementation, then the output memory will allocated automatically by this superclass. The ThreadedGenerateData() method should only produce the output specified by "outputThreadRegion" parameter. ThreadedGenerateData() cannot write to any other portion of the output image (as this is responsibility of a different thread).

See also:: GenerateData(), SplitRequestedRegion()

virtual bool itk::FiniteDifferenceImageFilter< TDeformationField , TDeformationField >::ThreadedHalt ( void * ) [inline, protected, virtual, inherited]

This method is similar to Halt(), and its default implementation in this class is simply to call Halt(). However, this method takes as a parameter a void pointer to the MultiThreader::ThreadInfoStruct structure. If you override this method instead of overriding Halt, you will be able to get the current thread ID and handle the Halt method accordingly. This is useful if you are doing a lot of processing in Halt that you don't want parallelized. Notice that ThreadedHalt is only called by the multithreaded filters, so you still should implement Halt, just in case a non-threaded filter is used.

Definition at line 284 of file itkFiniteDifferenceImageFilter.h.

{ return this->Halt(); }

static ITK_THREAD_RETURN_TYPE itk::ImageSource< TDeformationField >::ThreaderCallback ( void * arg ) [static, protected, inherited]

Static function used as a "callback" by the MultiThreader. The threading library will call this routine for each thread, which will delegate the control to ThreadedGenerateData().

void itk::Object::UnRegister ( ) const [virtual, inherited]

Decrease the reference count (release by another object).

If there is a delete method, invoke it.

Reimplemented from itk::LightObject.

Definition at line 324 of file itkObject.cxx.

Referenced by itk::Object::New(), and itk::ObjectFactoryBase::UnRegisterFactory().

{
  // call the parent
  itkDebugMacro(<< "UnRegistered, "
                << "ReferenceCount = " << (m_ReferenceCount-1));

  if ( (m_ReferenceCount-1) <= 0)
    {
    this->InvokeEvent(DeleteEvent());
    }

  Superclass::UnRegister();
}

void itk::ProcessObject::Update ( void ) [virtual, inherited]

Bring this filter up-to-date. Update() checks modified times against last execution times, and re-executes objects if necessary. A side effect of this method is that the whole pipeline may execute in order to bring this filter up-to-date. This method updates the currently prescribed requested region. If no requested region has been set on the output, then the requested region will be set to the largest possible region. Once the requested region is set, Update() will make sure the specified requested region is up-to-date. This is a confusing side effect to users who are just calling Update() on a filter. A first call to Update() will cause the largest possible region to be updated. A second call to Update() will update that same region. If a modification to the upstream pipeline cause a filter to have a different largest possible region, this second call to Update() will not cause the output requested region to be reset to the new largest possible region. Instead, the output requested region will be the same as the last time Update() was called. To have a filter always to produce its largest possible region, users should call UpdateLargestPossibleRegion() instead.

The following methods are used to coordinate the execution of the data processing pipeline.

Definition at line 611 of file itkProcessObject.cxx.

Referenced by otb::ImageSeriesFileReader< VectorImage< TPixel, 2 >, VectorImage< TInternalPixel, 2 > >::GenerateData(), otb::ImageSeriesFileReader< Image< TPixel, 2 >, VectorImage< TInternalPixel, 2 > >::GenerateData(), otb::ImageSeriesFileReader< Image< TPixel, 2 >, Image< TInternalPixel, 2 > >::GenerateData(), itk::DerivativeImageFilter< TInputImage, TOutputImage >::GenerateData(), itk::AutoCropLabelMapFilter< TInputImage >::GenerateOutputInformation(), and itk::LabelMapMaskImageFilter< TInputImage, TOutputImage >::GenerateOutputInformation().

{
  if (this->GetOutput(0))
    {
    this->GetOutput(0)->Update();
    }
}

void itk::ProcessObject::UpdateLargestPossibleRegion ( ) [virtual, inherited]

Like Update(), but sets the output requested region to the largest possible region for the output. This is the method users should call if they want the entire dataset to be processed. If a user wants to update the same output region as a previous call to Update() or a previous call to UpdateLargestPossibleRegion(), then they should call the method Update().

Definition at line 1120 of file itkProcessObject.cxx.

{
  this->UpdateOutputInformation();

  if (this->GetOutput(0))
    {
    this->GetOutput(0)->SetRequestedRegionToLargestPossibleRegion();
    this->GetOutput(0)->Update();
    }
}

void itk::ProcessObject::UpdateOutputData ( DataObject * output ) [virtual, inherited]

Actually generate new output

prevent chasing our tail

Prepare all the outputs. This may deallocate previous bulk data.

Propagate the update call - make sure everything we might rely on is up-to-date Must call PropagateRequestedRegion before UpdateOutputData if multiple inputs since they may lead back to the same data object.

Cache the state of any ReleaseDataFlag's on the inputs. While the filter is executing, we need to set the ReleaseDataFlag's on the inputs to false in case the current filter is implemented using a mini-pipeline (which will try to release the inputs). After the filter finishes, we restore the state of the ReleaseDataFlag's before the call to ReleaseInputs().

Tell all Observers that the filter is starting

GenerateData this object - we have not aborted yet, and our progress before we start to execute is 0.0.

Count the number of required inputs which have been assigned

If we ended due to aborting, push the progress up to 1.0 (since it probably didn't end there)

Notify end event observers

Now we have to mark the data as up to date.

Restore the state of any input ReleaseDataFlags

Release any inputs if marked for release

Reimplemented in otb::StreamingImageFileWriter< TInputImage >, and itk::StreamingImageFilter< TInputImage, TOutputImage >.

Definition at line 908 of file itkProcessObject.cxx.

{
  DataObjectPointerArraySizeType idx;

  if (m_Updating)
    {
    return;
    }

  this->PrepareOutputs();

  m_Updating = true;
  if ( m_Inputs.size() == 1 )
    {
    if (m_Inputs[0])
      {
      m_Inputs[0]->UpdateOutputData();
      }
    }
  else
    {
    for (idx = 0; idx < m_Inputs.size(); ++idx)
      {
      if (m_Inputs[idx])
        {
        m_Inputs[idx]->PropagateRequestedRegion();
        m_Inputs[idx]->UpdateOutputData();
        }
      }
    }

  this->CacheInputReleaseDataFlags();
  
  this->InvokeEvent( StartEvent() );

  m_AbortGenerateData = false;
  m_Progress = 0.0f;

  DataObjectPointerArraySizeType ninputs = this->GetNumberOfValidRequiredInputs();
  if (ninputs < m_NumberOfRequiredInputs)
    {
    itkExceptionMacro(<< "At least " << m_NumberOfRequiredInputs 
                      << " inputs are required but only " << ninputs 
                      << " are specified.");
    }
  else
    {
    try
      {
      this->GenerateData();
      }
    catch( ProcessAborted & excp )
      {
      this->InvokeEvent( AbortEvent() );
      this->ResetPipeline();
      this->RestoreInputReleaseDataFlags();
      throw excp;
      }
    catch( ExceptionObject& excp )
      {
      this->ResetPipeline();
      this->RestoreInputReleaseDataFlags();
      throw excp;
      }
    }

  if ( m_AbortGenerateData )
    {
    this->UpdateProgress(1.0f);
    }

  this->InvokeEvent( EndEvent() );

  for (idx = 0; idx < m_Outputs.size(); ++idx)
    {
    if (m_Outputs[idx])
      {
      m_Outputs[idx]->DataHasBeenGenerated();
      }
    }
  
  this->RestoreInputReleaseDataFlags();
  
  this->ReleaseInputs();
  
  // Mark that we are no longer updating the data in this filter
  m_Updating = false;
}

void itk::ProcessObject::UpdateOutputInformation ( void ) [virtual, inherited]

Update the information decribing the output data. This method transverses up the pipeline gathering modified time information. On the way back down the pipeline, this method calls GenerateOutputInformation() to set any necessary information about the output data objects. For instance, a filter that shrinks an image will need to provide an implementation for GenerateOutputInformation() that changes the spacing of the pixels. Such filters should call their superclass' implementation of GenerateOutputInformation prior to changing the information values they need (i.e. GenerateOutputInformation() should call Superclass::GenerateOutputInformation() prior to changing the information.

Watch out for loops in the pipeline

Since we are in a loop, we will want to update. But if we don't modify this filter, then we will not execute because our OutputInformationMTime will be more recent than the MTime of our output.

We now wish to set the PipelineMTime of each output DataObject to the largest of this ProcessObject's MTime, all input DataObject's PipelineMTime, and all input's MTime. We begin with the MTime of this ProcessObject.

Loop through the inputs

Propagate the UpdateOutputInformation call

What is the PipelineMTime of this input? Compare this against our current computation to find the largest one.

Pipeline MTime of the input does not include the MTime of the data object itself. Factor these mtimes into the next PipelineMTime

Call GenerateOutputInformation for subclass specific information. Since UpdateOutputInformation propagates all the way up the pipeline, we need to be careful here to call GenerateOutputInformation only if necessary. Otherwise, we may cause this source to be modified which will cause it to execute again on the next update.

Keep track of the last time GenerateOutputInformation() was called

Reimplemented in itk::watershed::Segmenter< TInputImage >, itk::VTKImageImport< TOutputImage >, and itk::watershed::Segmenter< InputImageType >.

Definition at line 664 of file itkProcessObject.cxx.

References itk::Object::GetMTime(), itk::DataObject::SetPipelineMTime(), and itk::DataObject::UpdateOutputInformation().

Referenced by itk::VTKImageExportBase::UpdateInformationCallback(), and itk::watershed::Segmenter< TInputImage >::UpdateOutputInformation().

{
  unsigned long t1, t2;
  DataObjectPointerArraySizeType idx;
  DataObject *input;
  DataObject *output;

  if ( m_Updating )
    {
    this->Modified();
    return;
    }

  t1 = this->GetMTime();

  for (idx = 0; idx < m_Inputs.size(); ++idx)
    {
    if (m_Inputs[idx])
      {
      input = m_Inputs[idx];

      m_Updating = true;
      input->UpdateOutputInformation();
      m_Updating = false;
      
      t2 = input->GetPipelineMTime();

      if (t2 > t1)
        {
        t1 = t2;
        }

      t2 = input->GetMTime();
      if (t2 > t1)
        {
        t1 = t2;
        }
      }
    }

  if (t1 > m_OutputInformationMTime.GetMTime())
    {
    for (idx = 0; idx < m_Outputs.size(); ++idx)
      {
      output = this->GetOutput( static_cast<int>( idx ) );
      if (output)
        {
        output->SetPipelineMTime(t1);
        }  
      }
    
    this->GenerateOutputInformation();

    m_OutputInformationMTime.Modified();
    }
}

void itk::ProcessObject::UpdateProgress ( float amount ) [inherited]

Update the progress of the process object.

Sets the Progress ivar to amount and invokes any observers for the ProgressEvent. The parameter amount should be in [0,1] and is the cumulative (not incremental) progress.

Update the progress of the process object. If a ProgressMethod exists, execute it. Then set the Progress ivar to amount. The parameter amount should range between (0,1).

Definition at line 500 of file itkProcessObject.cxx.

Referenced by itk::watershed::Segmenter< TInputImage >::GenerateData(), itk::ProgressReporter::ProgressReporter(), and itk::ProgressReporter::~ProgressReporter().

{
  m_Progress = amount;
  this->InvokeEvent( ProgressEvent() );
}

Member Data Documentation

template<class TFixedImage , class TMovingImage , class TDeformationField >

const unsigned int itk::PDEDeformableRegistrationFilter< TFixedImage, TMovingImage, TDeformationField >::ImageDimension = Superclass::ImageDimension [static, inherited]

Inherit some enums and typedefs from the superclass.

Reimplemented from itk::DenseFiniteDifferenceImageFilter< TDeformationField, TDeformationField >.

Definition at line 116 of file itkPDEDeformableRegistrationFilter.h.

const unsigned int itk::InPlaceImageFilter< TDeformationField , TDeformationField >::InputImageDimension [static, inherited]

ImageDimension constants

Reimplemented from itk::ImageToImageFilter< TDeformationField, TDeformationField >.

Definition at line 90 of file itkInPlaceImageFilter.h.

unsigned int itk::FiniteDifferenceImageFilter< TDeformationField , TDeformationField >::m_ElapsedIterations [protected, inherited]

A counter for keeping track of the number of elapsed iterations during filtering.

Definition at line 330 of file itkFiniteDifferenceImageFilter.h.

bool itk::FiniteDifferenceImageFilter< TDeformationField , TDeformationField >::m_ManualReinitialization [protected, inherited]

Indicates whether the filter automatically resets to UNINITIALIZED state after completing, or whether filter must be manually reset

Definition at line 334 of file itkFiniteDifferenceImageFilter.h.

unsigned int itk::FiniteDifferenceImageFilter< TDeformationField , TDeformationField >::m_NumberOfIterations [protected, inherited]

The maximum number of iterations this filter will run

Definition at line 326 of file itkFiniteDifferenceImageFilter.h.

TimeStamp itk::ProcessObject::m_OutputInformationMTime [protected, inherited]

Time when GenerateOutputInformation was last called.

Definition at line 421 of file itkProcessObject.h.

InternalReferenceCountType itk::LightObject::m_ReferenceCount [mutable, protected, inherited]

Number of uses of this object by other objects.

Definition at line 141 of file itkLightObject.h.

SimpleFastMutexLock itk::LightObject::m_ReferenceCountLock [mutable, protected, inherited]

Mutex lock to protect modification to the reference count

Definition at line 144 of file itkLightObject.h.

bool itk::ProcessObject::m_Updating [protected, inherited]

These ivars are made protected so filters like itkStreamingImageFilter can access them directly. This flag indicates when the pipeline is executing. It prevents infinite recursion when pipelines have loops.

Reimplemented in otb::HistogramAndTransferFunctionWidget< THistogram, TPixel >, and otb::ImageViewerBase< TPixel, TLabel >.

Definition at line 418 of file itkProcessObject.h.

const unsigned int itk::InPlaceImageFilter< TDeformationField , TDeformationField >::OutputImageDimension [static, inherited]

ImageDimension constant

Reimplemented from itk::ImageToImageFilter< TDeformationField, TDeformationField >.

Definition at line 92 of file itkInPlaceImageFilter.h.

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

otb::NCCRegistrationFilter< TFixedImage, TMovingImage, TDeformationField > Class Template Reference [Deformable Registration Methods]

Public Types

Public Member Functions

Static Public Member Functions

Static Public Attributes

Protected Types

Protected Member Functions

Static Protected Member Functions

Protected Attributes

Private Member Functions

Detailed Description

template<class TFixedImage, class TMovingImage, class TDeformationField> class otb::NCCRegistrationFilter< TFixedImage, TMovingImage, TDeformationField >

Member Typedef Documentation

Member Function Documentation

Member Data Documentation

otb::NCCRegistrationFilter< TFixedImage, TMovingImage, TDeformationField > Class Template Reference
[Deformable Registration Methods]

template<class TFixedImage, class TMovingImage, class TDeformationField>
class otb::NCCRegistrationFilter< TFixedImage, TMovingImage, TDeformationField >