[ITK-users] [ITK] Geodesic level set error

[Nicolas Gallego]

Hi Sindhoor,

It's hard to tell. Could you provide a simple code that reproduces your
error?
A priori I would say you could print the image information on each
connection of your pipeline, pay special attention to filters that may
change it, as the ROI filter...
Refer also to the concepts of regions: there are
LargestPossibleRegion
BufferedRegion
RequestedRegion

LargestPossibleRegion refers to the whole dataset. Refer to section 13.1 of
the software guide for more details.
I hope that helps,

Nicolás Gallego-Ortiz
Université catholique de Louvain, Belgium


2014-06-17 10:36 GMT+02:00 Sindhoor Bhat <bhatsindhoor at ymail.com>:

> Hi nicgallego,
> Thanks for reply.
> I've read that section. Here is some detail on steps i followed. This may
> help you to give some hints.
>
> 1.read dicom series
> 2.define region of interest (ROI filter)
> 3. compute gradient & speed image.
> 4.Set seed using fast marching filter.
> 5.Give this inital level set image---- to geodesic segmentation filter
>
>   geodesicActiveContour->SetInput(  fastMarching->GetOutput() );
>   geodesicActiveContour->SetFeatureImage( sigmoidFilter->GetOutput() );
>
>
> fastMarching->SetOutputSize(roi_filter->GetOutput()->GetBufferedRegion().GetSize()
> ); //roi=region of interest
>
>  Even after reading from same dicom series, why is this error? it should
> have same origin / spacing right?
>  can you pls help me ?
>
>
>
>   On Tuesday, 17 June 2014 12:59 PM, Nicolas Gallego <nicgallego at gmail.com>
> wrote:
>
>
> Good morning Sindhoor,
>
> In general itk pipelines process images in the physical space, so that
> pixelwise operations are between samples on the same physical position.
> Section 4.1.4 of the Software guide [1] and its correspondent example
> (Image4.cxx) may help you,
>
> [1] http://www.itk.org/ItkSoftwareGuide.pdf
>
> Nicolás Gallego-Ortiz
> Université catholique de Louvain, Belgium
>
>
> 2014-06-17 7:28 GMT+02:00 Sindhoor Bhat <bhatsindhoor at ymail.com>:
>
> Dear Matt,
>
> I've got error in Geodesic active contour level set segmentation :  saying
>  itk::ERROR: GeodesicActiveContourLevelSetImageFilter(0000000004AC7160):
> Inputs do not occupy the same physical space!
>
> I searched a lot on solving this error. but none have replied. Can you pls
> tell how to solve this error?
> (I'm trying to segment 3D dicom images)
>
> thank you.
>
> _____________________________________
> 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
>
> _______________________________________________
> Community mailing list
> Community at itk.org
> http://public.kitware.com/mailman/listinfo/community
>
>
>
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://public.kitware.com/pipermail/insight-users/attachments/20140617/af628586/attachment.html>