[ITK] [ITK-users] it.IsAtEnd() is true immediately after setting it.GoToBegin()

Matt McCormick matt.mccormick at kitware.com
Tue Apr 15 23:34:36 EDT 2014


Hi David,

Update() has to be called on fourierImage in order for its contents to
be available.

Hope this helps,
Matt

PS. Tip: ForwardFFTImageFilter can be used instead of
VnlForwardFFTImageFilter, and ForwardFFTImageFilter will automatically
use VnlForwardFFTImageFilter if FFTW is not available.

On Tue, Apr 15, 2014 at 10:28 AM, DVigneault <davis.vigneault at gmail.com> wrote:
> All--
>
> In connection with my previous  phase unwrapping question
> <http://itk-users.7.n7.nabble.com/Is-there-a-phase-unwrapping-filter-in-ITK-td33731.html>
> , I'm trying to write a filter with the following structure:
>
> --Take the vnl forward Fourier transform
> --Alter the pixel values using itkImageRegionIteratorWithIndex
> --Take the vnl inverse Fourier transform
> --Multiply by a constant
>
> However, in the second step, it.IsAtEnd() is true immediately after setting
> it.GoToBegin(), suggesting that I've done something wrong in assigning the
> output of the vnl forward FFT to the pointer that I'm trying to iterate
> over.
>
> Does anyone see where I've gone wrong?
>
> Best,
>
> --Davis
>
> *** Code ***
>
> #ifndef DV_FORWARDLAPLACIANIMAGEFILTER_H
> #define DV_FORWARDLAPLACIANIMAGEFILTER_H
>
> #include "itkImageToImageFilter.h"
> #include "itkVnlForwardFFTImageFilter.h"
> #include "itkVnlInverseFFTImageFilter.h"
> #include "itkMultiplyImageFilter.h"
> #include "itkImageRegionIteratorWithIndex.h"
> #include "vnl/vnl_math.h"
>
> namespace itk {
>
> template < typename TWorkImage >
> class ITK_EXPORT ForwardLaplacianImageFilter :
> public ImageToImageFilter< TWorkImage, TWorkImage >
> {
> public:
>
> //  Standard declarations
> //  Used for object creation with the object factory:
>
>   typedef ForwardLaplacianImageFilter                  Self;
>   typedef ImageToImageFilter< TWorkImage, TWorkImage > Superclass;
>   typedef SmartPointer<Self>                           Pointer;
>   typedef SmartPointer<const Self>                     ConstPointer;
>
>   // Typedefs
>         typedef TWorkImage                        WorkImageType;
>         typedef std::complex< double >            ComplexPixelType;
>         typedef itk::Image< ComplexPixelType, 2 > ComplexImageType;
>
>   /** Method for creation through object factory */
>   itkNewMacro(Self);
>
>   /** Run-time type information */
>   itkTypeMacro(ForwardLaplacianImageFilter, ImageToImageFilter);
>
>   /** Display */
>   void PrintSelf( std::ostream& os, Indent indent ) const;
>
> protected:
>
>   ForwardLaplacianImageFilter();
>
> //  Declare the component filter types:
>
> protected:
>
>         typedef VnlForwardFFTImageFilter< WorkImageType, ComplexImageType >
> ForwardFFTType;
>         typedef VnlInverseFFTImageFilter< ComplexImageType, WorkImageType >
> InverseFFTType;
>         typedef MultiplyImageFilter< WorkImageType, WorkImageType, WorkImageType >
> MultiplyType;
>         typedef ImageRegionIteratorWithIndex< ComplexImageType > ItType;
>
>     void GenerateData();
>
> private:
>
>   ForwardLaplacianImageFilter(Self&);   // intentionally not implemented
>   void operator=(const Self&);          // intentionally not implemented
>
> //  Component filters are declared as data members, all using the smart
> //  pointer types.
>
>         typename ForwardFFTType::Pointer m_FFT_Forward;
>         typename InverseFFTType::Pointer m_FFT_Inverse;
>         typename MultiplyType::Pointer   m_Multiply;
>
> };
>
> }
>
> //  The constructor sets up the pipeline, which involves creating the
> //  stages, connecting them together, and setting default parameters.
>
> namespace itk
> {
>
> template < typename WorkImageType >
> ForwardLaplacianImageFilter< WorkImageType >::ForwardLaplacianImageFilter()
> {
>
>         m_FFT_Forward = ForwardFFTType::New();
>         m_FFT_Inverse = InverseFFTType::New();
>         m_Multiply =    MultiplyType::New();
>
> }
>
> template < typename WorkImageType >
> void
> ForwardLaplacianImageFilter< WorkImageType >::GenerateData()
> {
>
>         // Description:
>         //              Take the FORWARD Fourier transform
>         //              Modify the values
>         //              Take the INVERSE Fourier transform
>         //              Multiply by a constant
>
>         // Set the input to the filter to the forward FFT
>         m_FFT_Forward->SetInput( this->GetInput() );
>
>         // Assign the output of this filter to a pointer
>         ComplexImageType::Pointer fourierImage = m_FFT_Forward->GetOutput();
>
>         // Iterate over this image to modify the values based on the indices
>         ItType fourierIt(fourierImage, fourierImage->GetLargestPossibleRegion() );
>     fourierIt.GoToBegin();
>     std::complex< double > a( 4, 0 );
>
>     // IsAtEnd() seems to be true, even though we're at the beginning
>     // Does this mean that the FFT output wasn't really assigned?
>     std::cout << fourierIt.IsAtEnd() << std::endl; // Prints "1"
>
>     // The while loop is not executed, as fourierIt.IsAtEnd() is true from
> the start
>     while ( !fourierIt.IsAtEnd() ) {
>         fourierIt.Set( a );
>         std::cout << fourierIt.Get() << std::endl; // Doesn't print
>         ++fourierIt;
>     }
>
>         // Assign that pointer to the inverse FFT filter
>         m_FFT_Inverse->SetInput( fourierImage );
>
>         // Multiply by a constant
>         m_Multiply->SetInput1( m_FFT_Inverse->GetOutput() );
>         m_Multiply->SetConstant( 2*pow(vnl_math::pi, 2)/pow(160, 2) );
>
>         // Assign the output of the multiplication filter to the output of this
> filter
>         m_Multiply->GraftOutput( this->GetOutput() );
>         m_Multiply->Update();
>         this->GraftOutput( m_Multiply->GetOutput() );
>
> }
>
> //  PrintSelf method prints parameters
>
> template < typename WorkImageType >
> void
> ForwardLaplacianImageFilter< WorkImageType >::PrintSelf( std::ostream& os,
> Indent indent ) const
> {
>   Superclass::PrintSelf(os,indent);
>
>   os
>     << indent << "No parameters."
>     << std::endl;
> }
>
> } /* end namespace itk */
>
> #endif
>
> *** End Code ***
>
>
>
> --
> View this message in context: http://itk-users.7.n7.nabble.com/it-IsAtEnd-is-true-immediately-after-setting-it-GoToBegin-tp33751.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
> _______________________________________________
> Community mailing list
> Community at itk.org
> http://public.kitware.com/cgi-bin/mailman/listinfo/community
_____________________________________
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



More information about the Community mailing list