[ITK-users] Splitting and repasting a volume with overlap

[Pol Monsó Purtí]

Hello everyone,

We need to split a volume into pieces, overlaped by a certain amount, do
some computation and then repaste them.

The computation requires the pieces' region indexes to be set to zero. The
original volume might not have index at zero, but be a subregion itself.

What is the best way to achieve this?

I've thought of the following pipeline:

1. imageRegionSplitterSlowDimension
2. Store the index and size of each region
3. manually increase the size of each region (and modify index)
4. crop with the largestPossibleRegion of the original volume (to prevent
requesting more than available)
iterate over regions:
      5. extractImageFilter
      6. disconnect pipeline for each region
      7. reset indexes
      8. perform computation
      9. crop with the regions in step 1 to remove the added overlap
      10. restore indexes
11. pasteImageFilter

I'll have to think of a smart way to deal with the index restore/cropping
of steps 9 and 10.

I somehow feel somebody will have already dealt with this... And there's
certainly a better way to do it.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://public.kitware.com/pipermail/insight-users/attachments/20150316/2edf79a6/attachment.html>