deformation_field.cxx

#include <itkVectorContainer.h>
#include "itkImage.h"
#include "itkImageFileReader.h"

const unsigned int                          Dimension = 3;
typedef double CoordinateRepType;
typedef itk::Vector< float, Dimension >      VectorType;
typedef itk::Image< VectorType, Dimension >  DisplacementFieldType;
typedef itk::VectorContainer<double,DisplacementFieldType::Pointer>DFTransformContainerType;

int main( int argc, char *argv[] )
{
       int inputId = 0;
       typedef itk::ImageFileReader< DisplacementFieldType >  fieldReaderType;
       DFTransformContainerType::Pointer DeffieldContainer = DFTransformContainerType::New() ;
       DeffieldContainer->Reserve( 1 );
       fieldReaderType::Pointer deffield = fieldReaderType::New();
       deffield->SetFileName(  argv[1]);
       deffield->Update();

       std::cout<<"Updating container!"<<std::endl;
              try
              {
                     DeffieldContainer->SetElement( inputId++ , deffield->GetOutput() ); //Fails here
              }
              catch (itk::ExceptionObject &excp)
              {
              std::cerr << "Exception thrown while updating container" << std::endl;
              std::cerr << excp << std::endl;
              return EXIT_FAILURE;
              }
       std::cout<<"DONE!"<<std::endl;

}