[Insight-users] vessel enhancing diffusion filter (release 200)

[Luis Ibanez]

Hi Laura,

This is annoying, but normal.

This code computes Hessians of images, which in 3D requires
the allocation of sizeof(double) x 6 bytes per pixel.

That is, you will need 48 bytes per pixel of your image
in order to store the resulting Hessian alone. There will be
of course additional intermediate allocations.

If you need to process large images you may need a 64bits
machine with a larger memory...

What is the actual size (in pixels) of the image that
you are processing ?

Can you process selected regions of the image ?

Usually there is a lot of empty (or at least, non interesting)
space in medical images. You could use the RegionOfInterest
filter to reduce the vessel enhancing processing to smaller
section of the image.


Please let us know,


     Thanks


        Luis





---------------------------------
Laura Fernandez de Manuel wrote:
> Hi all,
>  
> We have been checking the implementation of the "Vessel Enhancing 
> Diffusion Filter" depicted here:
>  
> http://insight-journal.org/midas/handle.php?handle=1926/558
>  
> although it works properly with the example 3D images provided (ranging 
> from around 30 to 250 KB) we didn't succeed to make it work in images 
> any larger (5MB images failed already for instance). We work in a system 
> with 4GB RAM so we don't know which can be the source of the "Failed to 
> allocate memory for image" errors that we get. Here, I attach the error 
> message we get:
>  
> ------------------------------------------------------------------------------------------------------
> ./itkAnisotropicDiffusionVesselEnhancementImageFilterTest.exe 
> image00_2.mhd image_2Enhanced.mhd
> Reading input image : image00_2.mhd
> Enhancing vessels.........: image00_2.mhd
> Iteration:      0
> Computing vesselness for scale with sigma= 0.2
> Exception caught:
> itk::ExceptionObject (0138FB20)
> Location: "class itk::SymmetricSecondRankTensor<double,3> *__thiscall 
> itk::ImportImageContainer<unsigned long,class 
> itk::SymmetricSecondRankTensor<double,3> >:: AllocateElements (unsigned 
> long) const"
> File: itk3.6.0\code\common\itkImportImageContainer.txx
> Line: 193
> Description: Failed to allocate memory for image.
> ------------------------------------------------------------------------------------------------------
>  
> Thanks a lot!
>  
> carlos & laura
> 
>  
>  
> 
> 
> ------------------------------------------------------------------------
> 
> _______________________________________________
> Insight-users mailing list
> Insight-users at itk.org
> http://www.itk.org/mailman/listinfo/insight-users