[Insight-users] itkStreamingImageFilter with itkHessianSmoothed3DToVesselnessMeasureImageFilter

[Moti Freiman]

Hi Luis,
Thanks for your detailed answer.
According to your answer, the Hessian filter itself is not streamable,
However, the main cause for memory insufficient is this filter,
since it produces images that are 6-times (in 3D) larger than the original
file.
If I'm right, streaming the  EigenAnalysisFilter and the vesselness will not
solve the problem of large files. ?
Moti

On Mon, Oct 6, 2008 at 3:34 PM, Luis Ibanez <luis.ibanez at kitware.com> wrote:

>
> Hi Moti,
>
> Looking at the code of this filter, the piece that looks suspicious is
> the internal use of the Symmetric Eigen Analysis filter.
>
> Both of these filters should be able to stream data, but, in the way
> they are nested, the EigenAnalysys filter may be trying to update the
> entire image unnecessarily.
>
> Please try the following:
>
> 1) Remove the EigenAnalysisFilter from inside the
>   itkHessianSmoothed3DToVesselnessMeasureImageFilter
>
> 2) create a sequential pipeline where you connect:
>
>   Hessian --> Symm.EigenAnalysis --> VesselnessFilter
>
> The Hessian filter itself is not streamable, but the
> Symm.EigenAnalysis and the Vesselness ones, are both
> pixel-wise filter, and they should be streamable.
>
>
> Please let us know if you encounter any problem.
>
> BTW: Once we sort this out, it will be great if you
>     could add the conclusions of your observations
>     as a review to this paper, and if you found it
>     useful, we could try moving this contribution
>     to the Review directory, in order to include it
>     with the upcoming release of ITK 3.10.
>
>
>  Regards,
>
>
>     Luis
>
>
> ---------------------
> Moti Freiman wrote:
>
>> Hi,
>> I want to compute a vesselness measure on a large CT volume (512*512*512,
>> 16bit per voxel) using the
>> itkHessianSmoothed3DToVesselnessMeasureImageFilter
>> published in the InsightJournal (
>> http://www.insight-journal.org/browse/publication/163).
>> Since my machine memory is limited I tried to use the
>> itkStreamingImageFilter  to divide the computation for several small parts
>> which do not require so much memory, and then merge them.
>> However it seems that the streaming filter do not divide the computation
>> using the itkHessianSmoothed3DToVesselnessMeasureImageFilter.
>> Any reason, idea?
>> Thanks,
>> Moti
>> --
>> __
>> Moti Freiman, Ph.D Student.
>> Medical Image Processing and Computer-Assisted Surgery Laboratory.
>> School of Computer Science and Engineering.
>> The Hebrew University of Jerusalem Givat Ram, Jerusalem 91904, Israel
>> Phone: +(972)-2-658-5371 (laboratory)
>> WWW site: http://www.cs.huji.ac.il/~freiman<http://www.cs.huji.ac.il/%7Efreiman>
>>
>>
>> ------------------------------------------------------------------------
>>
>> _______________________________________________
>> Insight-users mailing list
>> Insight-users at itk.org
>> http://www.itk.org/mailman/listinfo/insight-users
>>
>


-- 
__
Moti Freiman, Ph.D Student.
Medical Image Processing and Computer-Assisted Surgery Laboratory.
School of Computer Science and Engineering.
The Hebrew University of Jerusalem Givat Ram, Jerusalem 91904, Israel
Phone: +(972)-2-658-5371 (laboratory)
WWW site: http://www.cs.huji.ac.il/~freiman
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://www.itk.org/pipermail/insight-users/attachments/20081012/ed8677f4/attachment.htm>