[Insight-users] Convert a set of 2D images to 3D

[Luis Ibanez]

Hi Luis,

It is a bit more complicated than that.

You will need to put the filter in a for loop,
where you iterate over all the slices.

Inside the loop you take the output of the
PasteImageFilter and reconnect it as a new
input.


In Pseudo code it will look something like:

       ImageType::Pointer wholeImage;

       for-loop
       {
       pasteFilter->SetInput1( wholeImage );
       pasteFilter->SetInput2( slice_i );
       pasteFilter->Update();
       wholeImage = pasterFilter->GetOutupt();
       wholeImage->DisconnectPipeline();
       }



    Regards,


       Luis



--------------------------
Luis Roberto P. Paula wrote:
> After that, all I need to do is to call the SetInput for each slice and, 
> at the end, call the Update method?
> 
> Thanks,
> Luis
> 
> On Tue, Apr 21, 2009 at 10:06 AM, Luis Ibanez <luis.ibanez at kitware.com 
> <mailto:luis.ibanez at kitware.com>> wrote:
> 
> 
>     Hi Luis,
> 
>     You should instantiate the Paste filter in 3D.
> 
>     Since you are going to paste 2D slices into a 3D image,
>     you need to set up the filter in such a way that a 2D
>     image is seen as a 3D image of a single slice.
> 
> 
> 
>       Regards,
> 
> 
>           Luis
> 
> 
>     -----------------------------
>     Luis Roberto P. Paula wrote:
> 
>         Hi Luis,
> 
>         I'm a little confused on how to use the itkPasteImageFilter.
> 
>         At first, I created the typedef:
> 
>         typedef itk::PasteImageFilter< ImageType2D > PasteFilter2D;
> 
> 
>         Then:
> 
>           PasteFilter2D::Pointer paste2D = PasteFilter2D::New();
> 
>           // I'm lost in this part
> 
>           for (int i; i < VECTOR_SIZE; i++) {
>              paste2D->SetInput( IMAGE2D_VECTOR[i] );
>           }
> 
>         How do I make the output of the paste2D be a 3D image?
> 
>         Do I need to set to the destination image a 3D blank image?
> 
>         The destination index affects the output?
> 
>         I thought it would work like the add method of the std:vector
>         class with an output as a 3D image, but it seems that it can
>         return an image with the same dimension.
> 
>         Thanks,
>         Luis
> 
>         On Fri, Apr 17, 2009 at 3:36 PM, Luis Ibanez
>         <luis.ibanez at kitware.com <mailto:luis.ibanez at kitware.com>
>         <mailto:luis.ibanez at kitware.com
>         <mailto:luis.ibanez at kitware.com>>> wrote:
> 
> 
>            Hi Luis,
> 
>            Yes,
>            the same should work for pasting 3D images into to 4D images.
> 
> 
> 
>              Regards,
> 
> 
>                 Luis
> 
> 
>            ---------------------------
>            Luis Roberto P. Paula wrote:
> 
>                Thanks Luis!! This is a much better approach
>         (itkPasteImageFilter).
> 
>                I guess it works the same way from 3D to 4D, right?
> 
>                Regards,
>                Luis
> 
>                On Thu, Apr 16, 2009 at 1:59 PM, Luis Ibanez
>                <luis.ibanez at kitware.com <mailto:luis.ibanez at kitware.com>
>         <mailto:luis.ibanez at kitware.com <mailto:luis.ibanez at kitware.com>>
>                <mailto:luis.ibanez at kitware.com
>         <mailto:luis.ibanez at kitware.com>
>                <mailto:luis.ibanez at kitware.com
>         <mailto:luis.ibanez at kitware.com>>>> wrote:
> 
> 
>                   Hi Luis,
> 
>                   If your images are in individual 2D files, then
>                   you could simply use the ImageSeriesReader.
> 
>                   If your images are in an array of 2D images in
>                   memory, then you may want to use the Paste
>                   filter:
> 
>                      http://itk.org/Doxygen/html/itkPasteImageFilter_8h.html
> 
>                   You will have to call it multiple times...
>                   once per input slice.
> 
> 
>                      Regards,
> 
> 
>                          Luis
> 
> 
>                   ---------------------------
> 
>                   Luis Roberto P. Paula wrote:
> 
>                       Hi All,
> 
>                       I have a set of 2D images in axial orientation and
>         I want to
>                       create a 3D volume with this images.
> 
>                       It would be de reverse operation of the class
>                ExtractImageFilter.
> 
>                       Does anybody knows how to do it?
> 
>                       Thanks in advance,
>                       Luis
> 
> 
>                            
>         ------------------------------------------------------------------------
> 
>                       _____________________________________
>                       Powered by www.kitware.com
>         <http://www.kitware.com> <http://www.kitware.com>
>                <http://www.kitware.com>
> 
> 
>                       Visit other Kitware open-source projects at
>                       http://www.kitware.com/opensource/opensource.html
> 
>                       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://www.itk.org/mailman/listinfo/insight-users
> 
> 
> 
>