[Insight-users] Is FFT in ITK a complex FFT

[Martin Kavec]

Hi,

following code does not compile:

  typedef std::complex<float> ComplexPixelType;
  const   unsigned int        Dimension = 3;

  typedef itk::Image< ComplexPixelType, Dimension > ComplexImageType;

  typedef itk::ImageFileReader< ComplexImageType >  ReaderType;
  ReaderType::Pointer reader = ReaderType::New();
  reader->SetFileName( "foo.hdr" );


  typedef itk::VnlFFTComplexConjugateToRealImageFilter< ComplexPixelType , 
Dimension >  FFTFilterType;
  FFTFilterType::Pointer fftFilter = FFTFilterType::New();

  fftFilter->SetInput( reader->GetOutput() );

So I am in doubts, if Fourier transform in ITK is actually a complex FFT?

Thanks for explanation.

Martin