[Insight-users] Problem PointSetToImageFilter

[Karthik Krishnan]

Vera Naumburger wrote:

> Hi Karthik,
>
> thanks you very much for your response. Yes, you are right. I´m 
> working on  the 3D Curves Extractor based on the Insight Application 
> 2D Curves  Extractor. Luis already told me, that you are working on 
> the 3D Version of  the application, which is great!
> Your comments helped me a lot. Actually downloading the new version 
> of  itkImageToParametricSpaceFilter.txx solved the problem in one 
> second :-)))  I´m so happy that this problem is fixed !!!
>
Glad you found it useful.

> You advised me to update the Auxiallary/FltkImageViewer sources. But 
> I  could not find new versions in the CVS?
>
It is there... Here's the webcvs:
http://www.itk.org/cgi-bin/viewcvs.cgi/Auxiliary/FltkImageViewer/?root=InsightApplications

> Thanks,
> Vera
>
thanks
karthik

>
>> Vera
>>
>> Are you using the Curves2DExtractor application or the 
>> Curves3DExtractor  application ? The sources for Curves3D were 
>> committed yesterday, so  there may be a few issues.
>>
>> Did you update InsightApplications/Curves3DExtractor. We just 
>> updated  the sources to display curves 3D for MRA images.
>> You can find a screenshot at  
>> www.ece.arizona.edu/~karthikk/ITK/Curves3DScreenshot.png
>>
>> Here are some of the issues you may be having, ( We had similar 
>> issues a  few days ago).
>>
>> 1. ImageToParametricSpaceImageFilter was creating new point data 
>> every  time erasing existing point data if present. This was fixed a 
>> few days  ago. Please update itkImageToParametricSpaceImageFilter.txx
>>
>> 2. InteriorExteriorMeshFilter uses a spatial function. However this  
>> wasn't an Input 'per say'. So changes to the spatial function do not  
>> cause the filter pipeline to re-run. We fixed this by using a 
>> dataobject  decorator and passing it as an input. Please update  
>> itkInteriorExteriorMeshFilter.txx .h
>>
>> 3. FLTK image viewer, as of now takes only int, char, short images 
>> as  overlays... the output from your point set is probably double... 
>> Please  cast (see ceExtractorConsoleBase.h)
>> //snip
>>   m_ThresholdImageFilter->UpdateLargestPossibleRegion();
>>   m_Viewer_Extracted_Points.SetOverlayColorIndex( 0 ); // Overlay 
>> with  red points
>>   m_Viewer_Extracted_Points.SetImage( m_Reader->GetOutput() );    
>> m_Viewer_Extracted_Points.SetOverlay(  
>> m_ThresholdImageFilter->GetOutput() );  4. Update all the sources in 
>> Auxiallary/FltkImageViewer
>>
>> 5. Update Code/Common/itkFrustumSpatialFunction.txx .h Please let us  
>> know if you have any issues with the Curves3DExtractor application.  
>> Looks like you are trying to do the same thing.
>>
>> You may also find the following filter useful:  
>> BasicFilters/itkHessian3DToVesselnessMeasureImageFilter   This is 
>> the  line filter from Sato's paper.
>>
>> Ther are still some issues with fltk draw events. If you close some 
>> of  the windows, you might find that the point set is not redrawn.
>>
>> Thanks
>> Regards
>> Karthik
>>
>> Vera Naumburger wrote:
>>
>>> Hi,
>>>
>>> I have a problem when using PointSetToImageFilter.
>>> I want to extract points from an Image. To filter them I have a   
>>> ParametricSpace with a FrustrumSpatialFunction.
>>> Points within the pyramid are extracted using the   
>>> InteriorExteriorMeshFilter.
>>> I get the PointSets by using ParametricSpaceToImageSpaceMeshFilter.
>>>
>>> Now I use PointSetToImageFilter to transform the PointSets into an  
>>> image  in order to display them with an FLTKImageViewer.
>>> When I press the output-button a function is called, that updates   
>>> PointSetToImageFilter:
>>>
>>> void
>>> ceExtractorConsole
>>> ::ShowPoints( void )
>>> {
>>>   m_resampler->UpdateLargestPossibleRegion();
>>>   m_PointSetToImageFilter->UpdateLargestPossibleRegion();
>>>
>>>   m_Viewer_Extracted_Points->SetImage( m_resampler->GetOutput() );
>>>   m_Viewer_Extracted_Points->SetOverlay(  
>>> m_PointSetToImageFilter->GetOutput());
>>>   m_Viewer_Extracted_Points->Show();
>>> }
>>>
>>>
>>> The problem is that the points are shown, when the button is 
>>> pressed  the  first time. If I press it again there are no points 
>>> shown anymore  and the  PointSetImage is not available for further 
>>> processing.
>>>
>>> Thanks for any help!!!
>>>
>>> Kind Regards,
>>> Vera
>>> _______________________________________________
>>> Insight-users mailing list
>>> Insight-users at itk.org
>>> http://www.itk.org/mailman/listinfo/insight-users
>>>
>>
>
>
>