[Insight-users] ImageToParametricSpaceFilter Problem

[Luis Ibanez]

Hi Ashley,

Welcome to ITK !,

Your code, using the ImageToParametricSpaceFilter, looks good.

Could you please run the following tests:


a) Connect a FilterWatcher, to the m_ParametricSpace object.
    and rerun your code to see if you get any progress report.

    This will let us know if the filter is actually going
    through the pixels.

    You only need to add:

       #include "itkFilterWatcher.h"

       FilterWatcher  watcher( m_ParametricSpace );



b) Manually call Update() in the filter that are providing
    inputs to the m_ParametricSpace one.

      m_RescaleIntensitySmoothed->Update();
      m_RescaleIntensityMedialness->Update();
      m_RescaleIntensityMaxEigen->Update();
      m_ParametricSpace->Update();


    Then rerun the filter.


    Normally this shouldn't be necessary, but in this case
    since you are not observing any output, this test will
    let us know if the ImageToParametricSpaceFilter is missing
    to trigger the invocation the Update() methods in its
    preceding filters.




Please let us know what you find after experiments (a) and (b).
It may be that you have found a bug in this filter class.


    Thanks


      Luis


------------------------------
Ashley G. Anderson III wrote:
> Hello everyone,
> I am very new to ITK and I am trying to learn it for a class project. 
> Most of it seems to make a lot of sense so for, but I am stuck at one 
> point. I am playing around trying to implement a simple vessel-detection 
> program similar to the Curves2DExtraction application, but I am having 
> trouble with the ImageToParametricSpaceFilter. It seems to not produce 
> any output. Here is the code I am using, which is basically taken 
> straight from the Curves2DExtraction application:
> 
>    ParametricSpaceFilterType::Pointer m_ParametricSpace = 
> ParametricSpaceFilterType::New();
> 
>    m_ParametricSpace->SetInput( 0, 
> m_RescaleIntensitySmoothed->GetOutput() );
>    m_ParametricSpace->SetInput( 1, 
> m_RescaleIntensityMedialness->GetOutput() );
>    m_ParametricSpace->SetInput( 2, 
> m_RescaleIntensityMaxEigen->GetOutput() );
>    m_ParametricSpace->Update();
> 
>    std::cout << m_ParametricSpace->GetOutput()->GetNumberOfPoints();
> 
> The cout statement at the end tells me there are zero points in the 
> output. Is there something obvious I am doing wrong? I have checked the 
> output of all three RescaleIntensityFilter images, and they look fine.
> 
> Thanks in advance,
> Ashley
> 
> _______________________________________________
> Insight-users mailing list
> Insight-users at itk.org
> http://www.itk.org/mailman/listinfo/insight-users
>