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];
typedef float FloatPixelType;
ReaderType::Pointer reader = ReaderType::New();
reader->SetFileName( inputFileName );
typedef unsigned char UnsignedCharPixelType;
PadFilterType::Pointer padFilter = PadFilterType::New();
padFilter->SetInput( reader->GetOutput() );
padding[0] = 0;
padding[1] = 2;
padding[2] = 6;
padFilter->SetPadUpperBound( padding );
FFTType::Pointer fftFilter = FFTType::New();
fftFilter->SetInput( padFilter->GetOutput() );
typedef FFTType::OutputImageType FFTOutputImageType;
RealFilterType::Pointer realFilter = RealFilterType::New();
realFilter->SetInput(fftFilter->GetOutput());
RescaleFilterType::Pointer realRescaleFilter = RescaleFilterType::New();
realRescaleFilter->SetInput(realFilter->GetOutput());
WriterType::Pointer realWriter = WriterType::New();
realWriter->SetFileName( realFileName );
realWriter->SetInput( realRescaleFilter->GetOutput() );
try
{
realWriter->Update();
}
{
std::cerr << "Error: " << error << std::endl;
return EXIT_FAILURE;
}
ImaginaryFilterType::Pointer imaginaryFilter = ImaginaryFilterType::New();
imaginaryFilter->SetInput(fftFilter->GetOutput());
RescaleFilterType::Pointer imaginaryRescaleFilter = RescaleFilterType::New();
imaginaryRescaleFilter->SetInput(imaginaryFilter->GetOutput());
WriterType::Pointer complexWriter = WriterType::New();
complexWriter->SetFileName( imaginaryFileName );
complexWriter->SetInput( imaginaryRescaleFilter->GetOutput() );
try
{
complexWriter->Update();
}
{
std::cerr << "Error: " << error << std::endl;
return EXIT_FAILURE;
}
ModulusFilterType::Pointer modulusFilter = ModulusFilterType::New();
modulusFilter->SetInput(fftFilter->GetOutput());
RescaleFilterType::Pointer magnitudeRescaleFilter = RescaleFilterType::New();
magnitudeRescaleFilter->SetInput(modulusFilter->GetOutput());
WriterType::Pointer magnitudeWriter = WriterType::New();
magnitudeWriter->SetFileName( modulusFileName );
magnitudeWriter->SetInput( magnitudeRescaleFilter->GetOutput() );
try
{
magnitudeWriter->Update();
}
{
std::cerr << "Error: " << error << std::endl;
return EXIT_FAILURE;
}
return EXIT_SUCCESS;
}