itkTransformToSpatialJacobianSource.h

Go to the documentation of this file.

/*======================================================================

This file is part of the elastix software.

Copyright (c) University Medical Center Utrecht. All rights reserved.
See src/CopyrightElastix.txt or http://elastix.isi.uu.nl/legal.php for
details.

This software is distributed WITHOUT ANY WARRANTY; without even
the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR
PURPOSE. See the above copyright notices for more information.

======================================================================*/
#ifndef __itkTransformToSpatialJacobianSource_h
#define __itkTransformToSpatialJacobianSource_h

#include "itkAdvancedTransform.h"
#include "itkImageSource.h"

namespace itk
{

template< class TOutputImage,
class TTransformPrecisionType = double >
class TransformToSpatialJacobianSource :
  public ImageSource< TOutputImage >
{
public:

  typedef TransformToSpatialJacobianSource Self;
  typedef ImageSource< TOutputImage >      Superclass;
  typedef SmartPointer< Self >             Pointer;
  typedef SmartPointer< const Self >       ConstPointer;

  typedef TOutputImage                           OutputImageType;
  typedef typename OutputImageType::Pointer      OutputImagePointer;
  typedef typename OutputImageType::ConstPointer OutputImageConstPointer;
  typedef typename OutputImageType::RegionType   OutputImageRegionType;

  itkNewMacro( Self );

  itkTypeMacro( TransformToSpatialJacobianSource, ImageSource );

  itkStaticConstMacro( ImageDimension, unsigned int,
    TOutputImage::ImageDimension );

  typedef AdvancedTransform< TTransformPrecisionType,
    itkGetStaticConstMacro( ImageDimension ),
    itkGetStaticConstMacro( ImageDimension ) >     TransformType;
  typedef typename TransformType::ConstPointer        TransformPointerType;
  typedef typename TransformType::SpatialJacobianType SpatialJacobianType;

  typedef typename OutputImageType::PixelType PixelType;
  //typedef typename PixelType::ValueType           PixelValueType;
  typedef typename OutputImageType::RegionType    RegionType;
  typedef typename RegionType::SizeType           SizeType;
  typedef typename OutputImageType::IndexType     IndexType;
  typedef typename OutputImageType::PointType     PointType;
  typedef typename OutputImageType::SpacingType   SpacingType;
  typedef typename OutputImageType::PointType     OriginType;
  typedef typename OutputImageType::DirectionType DirectionType;

  typedef ImageBase< itkGetStaticConstMacro( ImageDimension ) > ImageBaseType;

  itkSetConstObjectMacro( Transform, TransformType );

  itkGetConstObjectMacro( Transform, TransformType );

  virtual void SetOutputSize( const SizeType & size );

  virtual const SizeType & GetOutputSize();

  virtual void SetOutputIndex( const IndexType & index );

  virtual const IndexType & GetOutputIndex();

  itkSetMacro( OutputRegion, OutputImageRegionType );

  itkGetConstReferenceMacro( OutputRegion, OutputImageRegionType );

  itkSetMacro( OutputSpacing, SpacingType );
  virtual void SetOutputSpacing( const double * values );

  itkGetConstReferenceMacro( OutputSpacing, SpacingType );

  itkSetMacro( OutputOrigin, OriginType );
  virtual void SetOutputOrigin( const double * values );

  itkGetConstReferenceMacro( OutputOrigin, OriginType );

  itkSetMacro( OutputDirection, DirectionType );
  itkGetConstReferenceMacro( OutputDirection, DirectionType );

  void SetOutputParametersFromImage( const ImageBaseType * image );

  virtual void GenerateOutputInformation( void );

  virtual void BeforeThreadedGenerateData( void );

  unsigned long GetMTime( void ) const;

protected:

  TransformToSpatialJacobianSource();
  ~TransformToSpatialJacobianSource() {}

  void PrintSelf( std::ostream & os, Indent indent ) const;

  void ThreadedGenerateData(
    const OutputImageRegionType & outputRegionForThread,
    ThreadIdType threadId );

  void NonlinearThreadedGenerateData(
    const OutputImageRegionType & outputRegionForThread,
    ThreadIdType threadId );

  void LinearGenerateData( void );

private:

  TransformToSpatialJacobianSource( const Self & ); // purposely not implemented
  void operator=( const Self & );                   // purposely not implemented

  RegionType           m_OutputRegion;         // region of the output image
  TransformPointerType m_Transform;            // Coordinate transform to use
  SpacingType          m_OutputSpacing;        // output image spacing
  OriginType           m_OutputOrigin;         // output image origin
  DirectionType        m_OutputDirection;      // output image direction cosines

};

} // end namespace itk

#ifndef ITK_MANUAL_INSTANTIATION
#include "itkTransformToSpatialJacobianSource.hxx"
#endif

#endif // end #ifndef __itkTransformToSpatialJacobianSource_h