[Insight-users] Operations on large 3D 8 bits datas

[Dan Mueller]

Hi Yannick,

Indeed, the itkGradientMagnitudeRecursiveGaussianImageFilter uses 
floating point intermediate images: the great thing about open source is 
that you can see for yourself!
   
http://www.itk.org/cgi-bin/viewcvs.cgi/Code/BasicFilters/itkGradientMagnitudeRecursiveGaussianImageFilter.h?root=Insight&view=markup
around line 71.

I too am working with large 3D datasets and regularly received similar 
memory allocation exceptions -- the solution for me was to move to a 
64-bit machine. Are you on a 32-bit machine? How much physical memory 
does your machine have?

Cheers, Dan

Yannick Pannier wrote:
> Hi everobody,
>
> I have some difficulties to adapt  
> GradientMagnitudeRecursiveGaussianImageFilter to be used on 
> 1024*1024*512 8-bits data.
> I have modified the  "InputPixelType" and "OutputPixelType" to be 
> "unsigned char", also :
> "typedef itk::Image< InputPixelType,  3 >   InputImageType;
>  typedef itk::Image< OutputPixelType, 3 >   OutputImageType; "
> and :
> "typedef unsigned char WritePixelType;
>  typedef itk::Image< WritePixelType, 3 > WriteImageType;"
> and when executed on 1024*1024*512 8 bits data, it returns :
> "terminate called after throwing an instance of 'itk::ExceptionObject'
>  what():  .../src/Insight/Code/Common/itkImportImageContainer.txx:188:
> Failed to allocate memory for image."
> I think there stills some hidden float operations due to the gaussian 
> filtering because when I launch the GradientMagnitudeImageFilter (with 
> the same modifications in the code) with the same data, it ends with 
> success.
> Can you tell me what I've missed ?
>
> Thanks
>
> Yannick


-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://public.kitware.com/pipermail/insight-users/attachments/20070728/20f03dd1/attachment.html