freImageTypes.txx

Go to the documentation of this file.

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

  Program:   F.R.E.E. - flexible registration evaluation engine
  Version:   v.1.0.0
  Date:      $Date: 2006/09/01 12:00:00 $
  Module:    $RCSfile: freImageTypes.txx,v $
  Language:  C++



  Copyright (c) 2007 Ralf o Floca (Department of Medical Informatics,
  Institute for Medical Biometry and Informatics, University of Heidelberg,
  Germany). All rights reserved.
  See FREECopyright.txt or http://www.mi.med.uni-hd.de/free/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 __freImageTypes_txx
#define __freImageTypes_txx

#include "freImageTypes.h"

#include "itkImageRegionIterator.h"

namespace FREE
{

template < class TImageType >
typename DefaultImages<TImageType>::ImagePointer
DefaultImages<TImageType>::
CreateGridImage(const RegionType& imageRegion, const SpacingType& imageSpacing, const PixelType& gridValue, const PixelType& backValue, const unsigned int& iGridSize)
{
        typedef itk::ImageRegionIterator< ImageType > IteratorType;
        ImagePointer gridImage = ImageType::New();

        //Generate imagegrid
        gridImage->SetRegions(imageRegion);
        gridImage->SetSpacing(imageSpacing);
        gridImage->Allocate();
        gridImage->FillBuffer(backValue);

  IteratorType iterator(gridImage,gridImage->GetLargestPossibleRegion());
        for (iterator.GoToBegin(); !(iterator.IsAtEnd()); ++iterator)
        {
                unsigned int iDimInGrid = 0;
                typename ImageType::IndexType index = iterator.GetIndex();

                for (int iIndex=0; iIndex<gridImage->GetImageDimension(); iIndex++)
                {
                        if (index[iIndex]%iGridSize==0) iDimInGrid++;
                };

                if (iDimInGrid>=(gridImage->GetImageDimension()-1)) iterator.Set(gridValue); 
                //at least in n-1 dimensions the index is on a gridline
        };
        return gridImage;
}

}//End of Namespace free

#endif

---