[Insight-developers] vnl FFT filter doesn't build when dim=4

Fri Aug 26 11:42:02 EDT 2011

It actually would be fairly simple to add a method templated over the
dimension, that would subsume the 1,2, & 3D cases. I have not done this
because I'm no longer funded for any ITK work that isn't directly
applicable to the work I am funded for.  In all cases, it's simply a
matter of using the image size array to set the dimensions for the vnl
FFT.

As to whether the output would be 'correct' I'm confident that it would be
reversible -- the output of IFFT(FFT(Image)) would be be close or
identical to Image.  I'm not sure there's an application crying out for a
N-Dimensional FFT for N > 3, though.

On 8/26/11 8:44 AM, "Bradley Lowekamp" <blowekamp at mail.nih.gov> wrote:

>I would agree with Kent. If it didn't produce correct output, then it's
>better not to be able to build it.
>It's a shame that we are uncertain if the output is correct :)
>
>Brad
>
>On Aug 26, 2011, at 9:37 AM, Gaëtan Lehmann wrote:
>
>
>
>I can't say if it was producing anything useful, but at least it was
>used to build.
>I'm not sure who should look at that, but it would be nice to have it
>buildable again!
>
>Gaëtan
>
>
>Le 19 août 11 à 18:16, Williams, Norman K a écrit :
>
>
>Though in fact, the code never worked, as far as I know for dim > 3.
>
>
>
>
>
>On 8/19/11 10:41 AM, "Cory Quammen" <cquammen at cs.unc.edu> wrote:
>
>
>
>Gaetan,
>
>
>
>I think this may have been caused by Kent's recent change:
>
>http://www.itk.org/gitweb?p=ITK.git;a=commit;h=abb9e14a
>
>
>
>Cory
>
>
>
>2011/8/19 Gaëtan Lehmann <gaetan.lehmann at jouy.inra.fr>:
>
>
>
>Hi Dirk,
>
>
>
>The vnl FFT filters are not building anymore when the image
>
>dimension
>
>is 4.
>
>
>
>http://www.cdash.org/CDash/viewBuildError.php?buildid=1444835
>
>
>
>Do you know how to fix that problem?
>
>
>
>Regards,
>
>
>
>Gaëtan
>
>
>
>
>
>/home/glehmann/src/tests/ITK/Modules/Filtering/FFT/include/
>
>itkVnlFFTCompl
>
>exConjugateToRealImageFilter.hxx:
>
>In member function 'void
>
>itk::VnlFFTComplexConjugateToRealImageFilter<TInputImage,
>
>TOutputImage>::GenerateData() [with TInputImage =
>
>itk::Image<std::complex<float>, 4u>, TOutputImage =
>
>itk::Image<float,
>
>4u>]':
>
>
>
>/home/glehmann/src/tests/WrapITK-build/Wrapping/WrapITK/Libraries/
>
>ITKFFT/
>
>itkVnlFFTComplexConjugateToRealImageFilterExplicit.cxx:6:
>
> instantiated from here
>
>
>
>/home/glehmann/src/tests/ITK/Modules/Filtering/FFT/include/
>
>itkVnlFFTCompl
>
>exConjugateToRealImageFilter.hxx:134:
>
>error: no matching function for call to
>
>
>
>'itk
>
>::VnlFFTComplexConjugateToRealImageFilter
>
><itk::Image<std::complex<flo
>
>at>,
>
>4u>, itk::Image<float, 4u>
>
>::FFTND_transform(vnl_vector<std::complex<float>
>
>&, const itk::Size<4u>&,
>
>
>
>itk
>
>::VnlFFTComplexConjugateToRealImageFilter
>
><itk::Image<std::complex<floa
>
>t>,
>
>4u>, itk::Image<float, 4u> >::DimDiscriminator<4u>*)'
>
>
>
>/home/glehmann/src/tests/ITK/Modules/Filtering/FFT/include/
>
>itkVnlFFTCompl
>
>exConjugateToRealImageFilter.hxx:61:
>
>note: candidates are: void
>
>itk::VnlFFTComplexConjugateToRealImageFilter<TInputImage,
>
>TOutputImage>::FFTND_transform(vnl_vector<typename
>
>TInputImage::PixelType>&,
>
>const typename TOutputImage::SizeType&,
>
>itk::VnlFFTComplexConjugateToRealImageFilter<TInputImage,
>
>TOutputImage>::DimDiscriminator<1u>*) [with TInputImage =
>
>itk::Image<std::complex<float>, 4u>, TOutputImage =
>
>itk::Image<float,
>
>4u>]
>
>
>
>/home/glehmann/src/tests/ITK/Modules/Filtering/FFT/include/
>
>itkVnlFFTCompl
>
>exConjugateToRealImageFilter.hxx:70:
>
>note:                 void
>
>itk::VnlFFTComplexConjugateToRealImageFilter<TInputImage,
>
>TOutputImage>::FFTND_transform(vnl_vector<typename
>
>TInputImage::PixelType>&,
>
>const typename TOutputImage::SizeType&,
>
>itk::VnlFFTComplexConjugateToRealImageFilter<TInputImage,
>
>TOutputImage>::DimDiscriminator<2u>*) [with TInputImage =
>
>itk::Image<std::complex<float>, 4u>, TOutputImage =
>
>itk::Image<float,
>
>4u>]
>
>
>
>/home/glehmann/src/tests/ITK/Modules/Filtering/FFT/include/
>
>itkVnlFFTCompl
>
>exConjugateToRealImageFilter.hxx:79:
>
>note:                 void
>
>itk::VnlFFTComplexConjugateToRealImageFilter<TInputImage,
>
>TOutputImage>::FFTND_transform(vnl_vector<typename
>
>TInputImage::PixelType>&,
>
>const typename TOutputImage::SizeType&,
>
>itk::VnlFFTComplexConjugateToRealImageFilter<TInputImage,
>
>TOutputImage>::DimDiscriminator<3u>*) [with TInputImage =
>
>itk::Image<std::complex<float>, 4u>, TOutputImage =
>
>itk::Image<float,
>
>4u>]
>
>
>
>
>
>--
>
>Gaëtan Lehmann
>
>Biologie du Développement et de la Reproduction
>
>INRA de Jouy-en-Josas (France)
>
>tel: +33 1 34 65 29 66    fax: 01 34 65 29 09
>
>http://mima2.jouy.inra.fr  http://www.itk.org
>
>http://www.bepo.fr
>
>
>
>
>
>
>
>
>
>
>
>--
>
>Cory Quammen
>
>Research Associate
>
>Department of Computer Science
>
>The University of North Carolina at Chapel Hill
>
>_______________________________________________
>
>Powered by www.kitware.com <http://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://kitware.com/products/protraining.html
>
>
>
>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-developers
>
>
>
>
>
>
>
>________________________________
>
>Notice: This UI Health Care e-mail (including attachments) is
>
>covered by the Electronic Communications Privacy Act, 18 U.S.C.
>
>2510-2521, is confidential and may be legally privileged.  If you
>
>are not the intended recipient, you are hereby notified that any
>
>retention, dissemination, distribution, or copying of this
>
>communication is strictly prohibited.  Please reply to the sender
>
>that you have received the message in error, then delete it.  Thank
>
>you.
>
>________________________________
>
>
>
>--
>Gaëtan Lehmann
>Biologie du Développement et de la Reproduction
>INRA de Jouy-en-Josas (France)
>tel: +33 1 34 65 29 66    fax: 01 34 65 29 09
>http://mima2.jouy.inra.fr  http://www.itk.org
>http://www.bepo.fr
>
><PGP.sig>_______________________________________________
>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://kitware.com/products/protraining.html
>
>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-developers
>
>
>
>
>
>========================================================
>Bradley Lowekamp
>Lockheed Martin Contractor for
>Office of High Performance Computing and Communications
>National Library of Medicine
>blowekamp at mail.nih.gov
>
>
>

________________________________
Notice: This UI Health Care e-mail (including attachments) is covered by the Electronic Communications Privacy Act, 18 U.S.C. 2510-2521, is confidential and may be legally privileged.  If you are not the intended recipient, you are hereby notified that any retention, dissemination, distribution, or copying of this communication is strictly prohibited.  Please reply to the sender that you have received the message in error, then delete it.  Thank you.
________________________________