<div dir="ltr">Hi Luis,<br>Thanks for your detailed answer.<br>According to your answer, the Hessian filter itself is not streamable, However, the main cause for memory insufficient is this filter,<br>since it produces images that are 6-times (in 3D) larger than the original file.<br>
If I'm right, streaming the EigenAnalysisFilter and the vesselness will not solve the problem of large files. ?<br>Moti<br><br><div class="gmail_quote">On Mon, Oct 6, 2008 at 3:34 PM, Luis Ibanez <span dir="ltr"><<a href="mailto:luis.ibanez@kitware.com">luis.ibanez@kitware.com</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;"><br>
Hi Moti,<br>
<br>
Looking at the code of this filter, the piece that looks suspicious is<br>
the internal use of the Symmetric Eigen Analysis filter.<br>
<br>
Both of these filters should be able to stream data, but, in the way<br>
they are nested, the EigenAnalysys filter may be trying to update the<br>
entire image unnecessarily.<br>
<br>
Please try the following:<br>
<br>
1) Remove the EigenAnalysisFilter from inside the<br>
itkHessianSmoothed3DToVesselnessMeasureImageFilter<br>
<br>
2) create a sequential pipeline where you connect:<br>
<br>
Hessian --> Symm.EigenAnalysis --> VesselnessFilter<br>
<br>
The Hessian filter itself is not streamable, but the<br>
Symm.EigenAnalysis and the Vesselness ones, are both<br>
pixel-wise filter, and they should be streamable.<br>
<br>
<br>
Please let us know if you encounter any problem.<br>
<br>
BTW: Once we sort this out, it will be great if you<br>
could add the conclusions of your observations<br>
as a review to this paper, and if you found it<br>
useful, we could try moving this contribution<br>
to the Review directory, in order to include it<br>
with the upcoming release of ITK 3.10.<br>
<br>
<br>
Regards,<br>
<br>
<br>
Luis<br>
<br>
<br>
---------------------<br>
Moti Freiman wrote:<br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;"><div><div></div><div class="Wj3C7c">
Hi,<br>
I want to compute a vesselness measure on a large CT volume (512*512*512, 16bit per voxel) using the itkHessianSmoothed3DToVesselnessMeasureImageFilter<br>
published in the InsightJournal (<a href="http://www.insight-journal.org/browse/publication/163" target="_blank">http://www.insight-journal.org/browse/publication/163</a>).<br>
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.<br>
However it seems that the streaming filter do not divide the computation using the itkHessianSmoothed3DToVesselnessMeasureImageFilter.<br>
Any reason, idea?<br>
Thanks,<br>
Moti<br>
-- <br>
__<br>
Moti Freiman, Ph.D Student.<br>
Medical Image Processing and Computer-Assisted Surgery Laboratory.<br>
School of Computer Science and Engineering.<br>
The Hebrew University of Jerusalem Givat Ram, Jerusalem 91904, Israel<br>
Phone: +(972)-2-658-5371 (laboratory)<br>
WWW site: <a href="http://www.cs.huji.ac.il/%7Efreiman" target="_blank">http://www.cs.huji.ac.il/~freiman</a><br>
<br>
<br></div></div>
------------------------------------------------------------------------<br>
<br>
_______________________________________________<br>
Insight-users mailing list<br>
<a href="mailto:Insight-users@itk.org" target="_blank">Insight-users@itk.org</a><br>
<a href="http://www.itk.org/mailman/listinfo/insight-users" target="_blank">http://www.itk.org/mailman/listinfo/insight-users</a><br>
</blockquote>
</blockquote></div><br><br clear="all"><br>-- <br>__<br>Moti Freiman, Ph.D Student.<br>Medical Image Processing and Computer-Assisted Surgery Laboratory.<br>School of Computer Science and Engineering.<br>The Hebrew University of Jerusalem Givat Ram, Jerusalem 91904, Israel<br>
Phone: +(972)-2-658-5371 (laboratory)<br>WWW site: <a href="http://www.cs.huji.ac.il/~freiman">http://www.cs.huji.ac.il/~freiman</a><br>
</div>