[ITK] [ITK-users] 1D Complex To Complex FFT of a 2D Image

[DVigneault]

All--

Sorry for the multiple e-mails--I've changed approach from my last post, and
wanted to send updated code before anyone spent time looking through the
last one.  I'm using the itkFFTWComplexToComplexImageFilter header files as
a starting point.  I copied them into a local directory and am editing them
directly.  To summarize:

1.  Use line iterator to copy one row into a std::complex<double> array.
2.  Convert to FFTW's input form using FFTW's
reinterpret_cast<fftw_complex*>() function
3.  Create the 1d plan
4.  Execute the plan
5.  Use linear iterator to copy into the output image
6.  Destroy the plan
7.  Repeat for the next row

This is giving me a Segmentation Fault on step 4 (when I execute the plan). 
Any ideas where I might be going wrong?

Best, and thanks,

--Davis

itkFFTWComplexToComplexImageFilter.h
<http://itk-users.7.n7.nabble.com/file/n33871/itkFFTWComplexToComplexImageFilter.h>  
itkFFTWComplexToComplexImageFilter.hxx
<http://itk-users.7.n7.nabble.com/file/n33871/itkFFTWComplexToComplexImageFilter.hxx>  





--
View this message in context: http://itk-users.7.n7.nabble.com/ITK-users-1D-Complex-To-Complex-FFT-of-a-2D-Image-tp33868p33871.html
Sent from the ITK - Users mailing list archive at Nabble.com.
_____________________________________
Powered by www.kitware.com

Visit other Kitware open-source projects at
http://www.kitware.com/opensource/opensource.html

Kitware offers ITK Training Courses, for more information visit:
http://www.kitware.com/products/protraining.php

Please keep messages on-topic and check the ITK FAQ at:
http://www.itk.org/Wiki/ITK_FAQ

Follow this link to subscribe/unsubscribe:
http://www.itk.org/mailman/listinfo/insight-users