itkImage.h

Go to the documentation of this file.

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

  Program:   Insight Segmentation & Registration Toolkit
  Module:    $RCSfile: itkImage.h,v $
  Language:  C++
  Date:      $Date: 2002/09/11 19:57:12 $
  Version:   $Revision: 1.104 $

  Copyright (c) 2002 Insight Consortium. All rights reserved.
  See ITKCopyright.txt or http://www.itk.org/HTML/Copyright.htm 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 __itkImage_h
#define __itkImage_h

#include "itkImageBase.h"
#include "itkImageRegion.h"
#include "itkImportImageContainer.h"
#include "itkDefaultPixelAccessor.h"
#include "itkPoint.h"
#include "itkContinuousIndex.h"
#include "itkTransform.h"

namespace itk
{

template <class TPixel, unsigned int VImageDimension=2>
class ITK_EXPORT Image : public ImageBase<VImageDimension>
{
public:
  typedef Image               Self;
  typedef ImageBase<VImageDimension>  Superclass;
  typedef SmartPointer<Self>  Pointer;
  typedef SmartPointer<const Self>  ConstPointer;
  
  itkNewMacro(Self);  

  itkTypeMacro(Image, ImageBase);

  typedef TPixel PixelType;

  typedef TPixel ValueType ;

  typedef TPixel InternalPixelType;

  typedef DefaultPixelAccessor< PixelType > AccessorType;

  itkStaticConstMacro(ImageDimension, unsigned int, VImageDimension);
  
  typedef ImportImageContainer<unsigned long, PixelType> PixelContainer;

  typedef Index<VImageDimension>  IndexType;

  typedef Offset<VImageDimension>  OffsetType;

  typedef Size<VImageDimension>  SizeType;

  typedef ImageRegion<VImageDimension>  RegionType;
  
  typedef typename PixelContainer::Pointer PixelContainerPointer;

   typedef Transform< double, VImageDimension, VImageDimension > TransformType;
   typedef typename TransformType::Pointer                       TransformPointer;


  void Allocate();

  void SetRegions(RegionType region)
    {
    this->SetLargestPossibleRegion(region);
    this->SetBufferedRegion(region);
    this->SetRequestedRegion(region);
    };

  void SetRegions(SizeType size)
    {
    RegionType region; region.SetSize(size);
    this->SetLargestPossibleRegion(region);
    this->SetBufferedRegion(region);
    this->SetRequestedRegion(region);
    };

  virtual void Initialize();

  void FillBuffer (const TPixel& value);
  
  void SetPixel(const IndexType &index, const TPixel& value)
    {
    typename Superclass::OffsetValueType offset = this->ComputeOffset(index);
    (*m_Buffer)[offset] = value;
    }
  
  const TPixel& GetPixel(const IndexType &index) const
  {
    typename Superclass::OffsetValueType offset = this->ComputeOffset(index);
    return ( (*m_Buffer)[offset] );
  }

  TPixel& GetPixel(const IndexType &index)
    {
    typename Superclass::OffsetValueType offset = this->ComputeOffset(index);
    return ( (*m_Buffer)[offset] );
    }
    
  TPixel & operator[](const IndexType &index)
     { return this->GetPixel(index); }
  
  const TPixel& operator[](const IndexType &index) const
     { return this->GetPixel(index); }

  TPixel *GetBufferPointer()
    { return m_Buffer ? m_Buffer->GetBufferPointer() : 0; }
  const TPixel *GetBufferPointer() const
    { return m_Buffer ? m_Buffer->GetBufferPointer() : 0; }
  
  PixelContainer* GetPixelContainer()
    { return m_Buffer.GetPointer(); }

  void SetPixelContainer( PixelContainer *container );
  
  AccessorType GetPixelAccessor( void ) 
    { return AccessorType(); }
    
  const AccessorType GetPixelAccessor( void ) const
    { return AccessorType(); }
    
  virtual void SetSpacing( const double spacing[VImageDimension] );
  virtual void SetSpacing( const float spacing[VImageDimension] );
  
  virtual void SetOrigin( const double origin[VImageDimension] );
  virtual void SetOrigin( const float origin[VImageDimension] );
  
  TransformPointer GetIndexToPhysicalTransform(void) const;

  TransformPointer GetPhysicalToIndexTransform(void) const;

  itkSetObjectMacro(IndexToPhysicalTransform, TransformType );
  itkSetObjectMacro(PhysicalToIndexTransform, TransformType );
    
  virtual void RebuildTransforms(void) throw ( std::exception );  

  template<class TCoordRep> 
  bool TransformPhysicalPointToContinuousIndex(
              const Point<TCoordRep, VImageDimension>& point, 
              ContinuousIndex<TCoordRep, VImageDimension>& index   ) const
    {
    // If no current transforms exist, throw an exception
    if ( !m_PhysicalToIndexTransform ) 
      { 
      itkExceptionMacro("The Image lacks a PhysicalToIndexTransform");
      }

    typename TransformType::InputPointType inputPoint =
                          m_PhysicalToIndexTransform->TransformPoint(point) ;

    // Update the output index
    for (unsigned int i = 0 ; i < VImageDimension ; i++)
      { 
      index[i] = static_cast<TCoordRep>(inputPoint[i]);
      }
    
    // Now, check to see if the index is within allowed bounds
    const bool isInside = 
      this->GetLargestPossibleRegion().IsInside( index );

    return isInside;
    }

  template<class TCoordRep> 
  bool TransformPhysicalPointToIndex(
            const Point<TCoordRep, VImageDimension>& point, 
            IndexType & index                                ) const
    {
    // If no current transforms exist, throw an exception
    if ( !m_PhysicalToIndexTransform ) 
      { 
      itkExceptionMacro("The Image lacks a PhysicalToIndexTransform");
      }

    typename TransformType::InputPointType inputPoint =
                          m_PhysicalToIndexTransform->TransformPoint(point) ;

    typedef typename IndexType::IndexValueType IndexValueType;

    // Update the output index
    for (unsigned int i = 0 ; i < VImageDimension ; i++)
      { 
      index[i] = static_cast<IndexValueType>(inputPoint[i]);
      }
    
    // Now, check to see if the index is within allowed bounds
    const bool isInside = 
      this->GetLargestPossibleRegion().IsInside( index );

    return isInside;
    }

  template<class TCoordRep> 
  void TransformContinuousIndexToPhysicalPoint( 
            const ContinuousIndex<TCoordRep, VImageDimension>& index, 
            Point<TCoordRep, VImageDimension>& point        ) const
    {
    // If no current transforms exist, throw an exception
    if ( !m_IndexToPhysicalTransform ) 
      { 
      itkExceptionMacro("The Image lacks a IndexToPhysicalTransform");
      }

    typename TransformType::InputPointType inputPoint;

    // Update the input index
    for (unsigned int i = 0 ; i < VImageDimension ; i++)
      { inputPoint[i] = index[i]; }

    // Transform the point
    typename TransformType::OutputPointType outputPoint =
      m_IndexToPhysicalTransform->TransformPoint(inputPoint) ;
  
    // Update the output point
    point = outputPoint;
    }

  template<class TCoordRep> 
  void TransformIndexToPhysicalPoint(  
                      const IndexType & index, 
                      Point<TCoordRep, VImageDimension>& point ) const 
    {
    if ( !m_IndexToPhysicalTransform ) 
      { 
      itkExceptionMacro("The Image lacks a IndexToPhysicalTransform");
      }

    typename TransformType::InputPointType inputPoint;

    // Update the input index
    for (unsigned int i = 0 ; i < VImageDimension ; i++)
      { 
      inputPoint[i] = index[i]; 
      }

    // Transform the point
    typename TransformType::OutputPointType outputPoint =
      m_IndexToPhysicalTransform->TransformPoint(inputPoint) ;
  
    // Update the output point
    for (unsigned int i = 0 ; i < VImageDimension ; ++i)
      { 
      point[i] = outputPoint[i]; 
      }
    }

  virtual void CopyInformation(const DataObject *data);

protected:
  Image();
  virtual ~Image();
  void PrintSelf(std::ostream& os, Indent indent) const;

private:
  Image(const Self&); //purposely not implemented
  void operator=(const Self&); //purposely not implemented

  PixelContainerPointer m_Buffer;

  TransformPointer  m_IndexToPhysicalTransform;
  TransformPointer  m_PhysicalToIndexTransform;
};

} // end namespace itk
  
#ifndef ITK_MANUAL_INSTANTIATION
#include "itkImage.txx"
#endif

#endif

---