Hi,<br><br>I am trying to take the inverse FFT of a complex image. However, I couldn't make it work. I searched the mailing list, and found the following mail (<a href="http://public.kitware.com/pipermail/insight-users/2005-June/013559.html">
http://public.kitware.com/pipermail/insight-users/2005-June/013559.html</a>). <br><br>I use the latest ITK version. The mail is quite old, but is it true?<br><br>Thanks in advance...<br><br>haris<br><br><pre>><i> Martin,
<br></i>><i><br></i>><i> I believe that the current FFT implementations will only work with<br></i>><i> taking forward FFT of float and double data types (i.e. not complex,<br></i>><i> integer, RGB, vector, or any other type of data). The current
<br></i>><i> implementation also takes advantage of the complex conjugate symmetry<br></i>><i> when working with real data.<br></i>><i><br></i>><i> Hans<br></i>><i><br></i>><i><br></i>><i><br></i>>
<i><br></i>><i> On Jun 09, 2005 07:41 AM, Martin Kavec <<a href="http://www.itk.org/mailman/listinfo/insight-users">kavec at messi.uku.fi</a>> wrote:<br></i>><i><br></i>><i> > Hi,<br></i>><i> ><br>
</i>><i> > following code does not compile:<br></i>><i> ><br></i>><i> > typedef std::complex<float> ComplexPixelType;<br></i>><i> > const unsigned int Dimension = 3;<br></i>><i>
><br></i>><i> > typedef itk::Image< ComplexPixelType, Dimension > ComplexImageType;<br></i>><i> ><br></i>><i> > typedef itk::ImageFileReader< ComplexImageType > ReaderType;<br></i>>
<i> > ReaderType::Pointer reader = ReaderType::New();<br></i>><i> > reader->SetFileName( "foo.hdr" );<br></i>><i> ><br></i>><i> ><br></i>><i> > typedef itk::VnlFFTComplexConjugateToRealImageFilter< ComplexPixelType
<br></i>><i> > ,<br></i>><i> > Dimension > FFTFilterType;<br></i>><i> > FFTFilterType::Pointer fftFilter = FFTFilterType::New();<br></i>><i> ><br></i>><i> > fftFilter->SetInput( reader->GetOutput() );
<br></i>><i> ><br></i>><i> > So I am in doubts, if Fourier transform in ITK is actually a complex<br></i>><i> > FFT?<br></i>><i> ></i></pre><br>