int main( int argc, char* argv[] )
{
if( argc != 5 )
{
std::cerr << "Usage: "<< std::endl;
std::cerr << argv[0];
std::cerr << " <InputFileName> <Real filename> <Imaginary filename> <Modulus filename>";
std::cerr << std::endl;
return EXIT_FAILURE;
}
const char * inputFileName = argv[1];
const char * realFileName = argv[2];
const char * imaginaryFileName = argv[3];
const char * modulusFileName = argv[4];
const unsigned int Dimension = 3;
typedef float FloatPixelType;
reader->SetFileName( inputFileName );
typedef unsigned char UnsignedCharPixelType;
padFilter->SetInput( reader->GetOutput() );
PadFilterType::SizeType padding;
padding[0] = 0;
padding[1] = 2;
padding[2] = 6;
padFilter->SetPadUpperBound( padding );
fftFilter->SetInput( padFilter->GetOutput() );
typedef FFTType::OutputImageType FFTOutputImageType;
realFilter->SetInput(fftFilter->GetOutput());
realRescaleFilter->SetInput(realFilter->GetOutput());
realWriter->SetFileName( realFileName );
realWriter->SetInput( realRescaleFilter->GetOutput() );
try
{
realWriter->Update();
}
{
std::cerr << "Error: " << error << std::endl;
return EXIT_FAILURE;
}
imaginaryFilter->SetInput(fftFilter->GetOutput());
imaginaryRescaleFilter->SetInput(imaginaryFilter->GetOutput());
complexWriter->SetFileName( imaginaryFileName );
complexWriter->SetInput( imaginaryRescaleFilter->GetOutput() );
try
{
complexWriter->Update();
}
{
std::cerr << "Error: " << error << std::endl;
return EXIT_FAILURE;
}
modulusFilter->SetInput(fftFilter->GetOutput());
magnitudeRescaleFilter->SetInput(modulusFilter->GetOutput());
magnitudeWriter->SetFileName( modulusFileName );
magnitudeWriter->SetInput( magnitudeRescaleFilter->GetOutput() );
try
{
magnitudeWriter->Update();
}
{
std::cerr << "Error: " << error << std::endl;
return EXIT_FAILURE;
}
return EXIT_SUCCESS;
}