#include <fstream>
int main( int argc, char * argv[] )
{
if( argc < 3 )
{
std::cerr << "Missing Parameters " << std::endl;
std::cerr << "Usage: " << argv[0];
std::cerr << " landmarksFile fixedImage ";
std::cerr << "movingImage deformedMovingImage" << std::endl;
return EXIT_FAILURE;
}
const unsigned int Dimension = 2;
typedef float VectorComponentType;
typedef unsigned char PixelType;
fixedReader->SetFileName( argv[2] );
try
{
fixedReader->Update();
}
{
std::cerr << "Exception thrown " << std::endl;
std::cerr << excp << std::endl;
return EXIT_FAILURE;
}
movingReader->SetFileName( argv[3] );
movingWriter->SetFileName( argv[4] );
DisplacementFieldType
> DisplacementSourceType;
deformer->SetOutputSpacing( fixedImage->GetSpacing() );
deformer->SetOutputOrigin( fixedImage->GetOrigin() );
deformer->SetOutputRegion( fixedImage->GetLargestPossibleRegion() );
deformer->SetOutputDirection( fixedImage->GetDirection() );
typedef DisplacementSourceType::LandmarkContainer LandmarkContainerType;
typedef DisplacementSourceType::LandmarkPointType LandmarkPointType;
LandmarkPointType sourcePoint;
LandmarkPointType targetPoint;
std::ifstream pointsFile;
pointsFile.open( argv[1] );
unsigned int pointId = 0;
pointsFile >> sourcePoint;
pointsFile >> targetPoint;
while( !pointsFile.fail() )
{
sourceLandmarks->InsertElement( pointId, sourcePoint );
targetLandmarks->InsertElement( pointId, targetPoint );
pointId++;
pointsFile >> sourcePoint;
pointsFile >> targetPoint;
}
pointsFile.close();
deformer->SetSourceLandmarks( sourceLandmarks.GetPointer() );
deformer->SetTargetLandmarks( targetLandmarks.GetPointer() );
try
{
deformer->UpdateLargestPossibleRegion();
}
{
std::cerr << "Exception thrown " << std::endl;
std::cerr << excp << std::endl;
return EXIT_FAILURE;
}
MovingImageType,
DisplacementFieldType > FilterType;
MovingImageType, double > InterpolatorType;
warper->SetInterpolator( interpolator );
warper->SetOutputSpacing( displacementField->GetSpacing() );
warper->SetOutputOrigin( displacementField->GetOrigin() );
warper->SetDisplacementField( displacementField );
warper->SetInput( movingReader->GetOutput() );
movingWriter->SetInput( warper->GetOutput() );
try
{
movingWriter->Update();
}
{
std::cerr << "Exception thrown " << std::endl;
std::cerr << excp << std::endl;
return EXIT_FAILURE;
}
return EXIT_SUCCESS;
}