[Insight-users] 2D Slice to 3D Volume

[Dário Oliveira]

Luis,

Just after I sent the email, I managed to run it with no errors, but I can't
have the right resulting image anyway. Here is the code:

 *typedef* itk::CastImageFilter < 2D ImageType, 3DImageType > FilterTypeCA;
FilterTypeCA::Pointer filterCast = FilterTypeCA::New();

*typedef* itk::PasteImageFilter< 3D ImageType, 3DImageType, 3DImageType >
FilterTypeP;
FilterTypeP::Pointer filterP = FilterTypeP::New();

3DImageType::RegionType outputRegion = image->GetLargestPossibleRegion();
3DImageType::SizeType size2 = outputRegion.GetSize();
3DImageType::IndexType start = outputRegion.GetIndex();
size2[2] = 1;
outputRegion.SetSize ( size2 );

filterP->SetDestinationImage(image);
filterP->SetInput(image);

*for* (*int* i=0; i<numSlices; i++)
{
 {some processing using the ith slice}
 start[2] = i;
 filter2D->*Update*();
 filterCast-> SetInput(filter2D->GetOutput());
 filterCast->* Update*();
 outputRegion.SetIndex( start );
 outputRegion = filterCast->GetOutput()->GetRequestedRegion();
 filterP->SetSourceImage(filterCast->GetOutput ());
 filterP->SetDestinationIndex(start);
 filterP->SetSourceRegion (outputRegion);
 filterP->*Update*();
}

When I see the filterP->getOutput() the result is exactly the same as the
filterP Input and Destination (which are the same image). And the processing
in the slice certainly change the image.

Thank you again,

Dário Oliveira.

2007/9/8, Dário Oliveira <dariodisk at gmail.com>:
>
> Hi Luis,
>
> thank you for your reply.
>
> Yes, I'm using the PasteImageFilter and I've tried two different
> approaches:
>
> 1. To define the filter with an 2D image as source and a 3D image as the
> other input, just like that:
> *
> typedef* itk::PasteImageFilter< 3DImageType , 2DImageType, 3DImageType >
> FilterTypeP ;
> FilterTypeP::Pointer filterP = FilterTypeP ::New();
>
> filterP->
> SetDestinationImage(image3D);
> filterP->SetInput(image3D);
>
> *for* (*int* i=0; i<numSlices; i++)
> {
>  {some processing using the ith slice}
>  filterP->SetSourceImage(filter2D->GetOutput());
>  filterP-> SetDestinationIndex(start);
>  filterP->SetSourceRegion(filter2D-> GetOutput
> ()->GetLargestPossibleRegion());
>  filterP->*Update*();
> }
> In that approach I got a compile error.
>
> 2. So I changed my approach using a cast image filter, like that:
>
> *typedef* itk::CastImageFilter < 2D ImageType, 3DImageType > FilterTypeCA;
> FilterTypeCA::Pointer filterCast = FilterTypeCA::New();
>
> *typedef* itk::PasteImageFilter< 3D ImageType, 3DImageType, 3DImageType >
> FilterTypeP;
> FilterTypeP::Pointer filterP = FilterTypeP::New();
>
> 3DImageType::RegionType outputRegion = image->GetLargestPossibleRegion();
> 3DImageType::SizeType size2 = outputRegion.GetSize();
> 3DImageType::IndexType start = outputRegion.GetIndex();
> size2[2] = 1;
> outputRegion.SetSize ( size2 );
>
> filterP->
> SetDestinationImage(image);
> filterP->SetInput(image);
>
> *for* (*int* i=0; i<numSlices; i++)
> {
>  {some processing using the ith slice}
>  start[2] = i;
>  filter2D->*Update*();
>  filterCast-> SetInput(filter2D->GetOutput());
>  filterCast->* Update*();
>  outputRegion = filterCast->GetOutput()->GetLargestPossibleRegion();
>   outputRegion.SetIndex( start );
>  filterP->SetSourceImage(filterCast->GetOutput ());
>  filterP->SetDestinationIndex(start);
>  filterP->SetSourceRegion (outputRegion);
>  filterP->*Update*();
> }
>
> In this case I got a runtime error in "filterP->*Update*();"
>
> I've just started to work with ITK, so I'm probably making some silly
> mistake.
>
> Thank you!
>
> Dário Olivera.
> 2007/9/8, Luis Ibanez <luis.ibanez at kitware.com>:
> >
> >
> >
> > Hi Dario,
> >
> > 1) Are you using the "PasteImageFilter" ?
> >
> > http://www.itk.org/Insight/Doxygen/html/classitk_1_1PasteImageFilter.html
> >
> >
> > 2) When you say that "it doesn't seem to work" ?
> >    Do you mean that
> >
> >    a) It doesn't compile ?
> >    b) It doesn't link ?
> >    c) It produces a segmentation fault at run time ?
> >    d) It throws an exception at run time ?
> >    e) It produces an image that when you visualize
> >       doesn't have the content that you expect ?
> >
> >
> > 3) Could you post to the list the minimal example of your code ?
> >
> >
> >
> >    Thanks
> >
> >
> >       Luis
> >
> >
> >
> > -----------------------
> > Dário Oliveira wrote:
> > > Hi friends,
> > >
> > > I'm having trouble to paste (actually replace) a processed 2D slice
> > into
> > > its original volume. I've read some previous itk-users email in the
> > > mailing list, but I still couldn't manage to suceed. Is there any
> > simple
> > > way to perfom this task?
> > >
> > > I'm trying to cast the processed 2D image to a 3D one, and get its 3D
> > > region to paste into the original 3D volume, but it doesn't seem to
> > work.
> > >
> > > Any help will be very apreciated,
> > >
> > > Regards
> > > --
> > > Dário Oliveira
> > >
> > >
> > >
> > ------------------------------------------------------------------------
> > >
> > > _______________________________________________
> > > Insight-users mailing list
> > > Insight-users at itk.org
> > > http://www.itk.org/mailman/listinfo/insight-users
> >
>
>
>
> --
> Dário Oliveira




-- 
Dário Oliveira
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://public.kitware.com/pipermail/insight-users/attachments/20070908/d187f4cc/attachment.html