[ITK-users] why can the itkStatisticsImageFilter not stream?

[Bradley Lowekamp]

Hi,

I was able to easy update a the StatisticsImageFilter to use my streaming base class for reduction algorithms [1]. It runs basically the same expect that you can set the number of stream divisions[2], there is the potential for setting the way the input image is split via the RegionSplitter in a base class [3].

You should be able to just clone the repository into you ITK/Modules/External directory, then enable it in the CMake configuration to use.

Please let me know if it works for you.

This should help :)
Brad

[1] https://github.com/blowekamp/itkStreamingSinc/blob/master/include/itkStreamingStatisticsImageFilter.h
[2] https://github.com/blowekamp/itkStreamingSinc/blob/master/test/itkStreamingStatisticsImageFilterTest2.cxx#L61
[3] https://github.com/blowekamp/itkStreamingSinc/blob/master/include/itkImageSinc.h#L99



On Apr 17, 2015, at 7:14 AM, Bradley Lowekamp <blowekamp at mail.nih.gov> wrote:

> Hello fellow mailing list member,
> 
> Currently in ITK there are two base classes for streaming images in ITK. The ImageToImageFilter provide the framework of streaming one input region to one output region. The other is the StreamingImageFilter which provides the framework to stream multiple image regions and merge them together into a single output image. This StatisticsImageFilter needs a framework which consumes multiple input image regions and reduces the data produced. I happen to have implemented such a framework here [1]. It currently has the LabelStatisticsImageFilter implemented in a streaming fashion.
> 
> Let me give it a try and see how quickly I can get the StatisticsImageFilter in this framework. Help writing tests would be appreciated.
> 
> Thanks,
> Brad
> 
> 
> [1] https://github.com/blowekamp/itkStreamingSinc
> 
> On Apr 17, 2015, at 5:08 AM, Dr. Roman Grothausmann <grothausmann.roman at mh-hannover.de> wrote:
> 
>> Dear mailing list members,
>> 
>> 
>> I'm wondering why the itkStatisticsImageFilter cannot stream?
>> Looking into itkStatisticsImageFilter.hxx it requests the whole input in GenerateInputRequestedRegion and EnlargeOutputRequestedRegion but in ThreadedGenerateData only iterates over the outputRegionForThread.
>> Where in the code is the whole input image needed as one single block?
>> 
>> Thanks for any hints
>> Roman
>> 
>> -- 
>> Dr. Roman Grothausmann
>> 
>> Tomographie und Digitale Bildverarbeitung
>> Tomography and Digital Image Analysis
>> 
>> Institut für Funktionelle und Angewandte Anatomie, OE 4120
>> Medizinische Hochschule Hannover
>> Carl-Neuberg-Str. 1
>> D-30625 Hannover
>> 
>> Tel. +49 511 532-2900
>> _____________________________________
>> Powered by www.kitware.com
>> 
>> Visit other Kitware open-source projects at
>> http://www.kitware.com/opensource/opensource.html
>> 
>> Kitware offers ITK Training Courses, for more information visit:
>> http://www.kitware.com/products/protraining.php
>> 
>> Please keep messages on-topic and check the ITK FAQ at:
>> http://www.itk.org/Wiki/ITK_FAQ
>> 
>> Follow this link to subscribe/unsubscribe:
>> http://public.kitware.com/mailman/listinfo/insight-users
> 
> _____________________________________
> Powered by www.kitware.com
> 
> Visit other Kitware open-source projects at
> http://www.kitware.com/opensource/opensource.html
> 
> Kitware offers ITK Training Courses, for more information visit:
> http://www.kitware.com/products/protraining.php
> 
> Please keep messages on-topic and check the ITK FAQ at:
> http://www.itk.org/Wiki/ITK_FAQ
> 
> Follow this link to subscribe/unsubscribe:
> http://public.kitware.com/mailman/listinfo/insight-users