Currently I am using the MedianImageFilter and the SmoothingRecursiveGaussianImageFilter or HessianRecursiveGaussianImageFilter.<br>We also do some other "filtering" using code that is not in the form of an itk filter (doesn't use any ITK yet)<br>
<br>This same project requires some visualization of this image - our plan was to create a 2D projection and create a sliding box that the is used to define the region that is rendered using VTK.<br>For this task I also need to just grab a portion of the image at one time for creating the projection and for creating the volume.<br>
<br>Isaac<br><br><br><div class="gmail_quote">On Thu, Jan 28, 2010 at 3:48 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;">
Hi Isaac,<br>
<br>
As you have already noticed, not all filters are<br>
streamable.<br>
<br>
However, for those that are, it should be possible<br>
to process your image by pieces.<br>
<br>
Brad Lowekamp has recently revamped the streaming<br>
pipeline and verified that is possible to pass the entire<br>
Visible Human dataset through an ITK pipeline.<br>
<br>
<br>
What specific ITK filters are you using ?<br>
<br>
<br>
If you are having problems with the pipeline, we should<br>
identify a reproducible minimal example, in case there<br>
is anything that needs to be fixed.<br>
<br>
<br>
Please let us know,<br>
<br>
<br>
Thanks<br>
<br>
<br>
<br>
Luis<br>
<br>
<br>
-----<br>
PS. I'm copying the ITK users-list, since this<br>
topics will be of interest to many users.<br>
<br>
--------------------------------------------------------------<br>
<div><div></div><div class="h5">On Thu, Jan 28, 2010 at 3:43 PM, Isaac Abbott <<a href="mailto:isaac.abbott@gmail.com">isaac.abbott@gmail.com</a>> wrote:<br>
> Luis,<br>
><br>
> I'm trying to process some very large 3D images. I cannot load these images<br>
> using ImageFileReader because I run out of RAM.<br>
> I would like to apply various filters to only subregions of this image, and<br>
> write the output to the same output image file.<br>
><br>
> This means my input and output images might be something like 10,000 x<br>
> 10,000 x 200 voxels, but I only want 1000x1000x100 in memory.<br>
> I'd load up a region, do the processing, then save the result to file, and<br>
> iterate, saving the results to the same file each time.<br>
><br>
> I have been following the discussion about applying filters to ROIs of an<br>
> itk::Image.<br>
> I've also seen some discussion about the usage of buffered region vs largest<br>
> possible region vs requested region.<br>
> And, I've looked at the StreamingImageFilter - which only seems to work for<br>
> a limited subset of the filters.<br>
><br>
> None of these discussions seem to meet my need. According to the<br>
> documentation I should be able to do something with the buffered region,<br>
> because this defines only the region of the image that is in memory - but I<br>
> haven't seen this work in practice.<br>
><br>
> Thanks for any help you can give,<br>
> Isaac (Roysam lab, RPI)<br>
><br>
><br>
</div></div></blockquote></div><br>