[Insight-developers] changes in level set routines since ITK 3.4 breaks itk-snap

[Luca Antiga]

Hey,
  I found out that by using the serial level set class instead of the  
parallel one, i.e. by replacing
     typedef ParallelSparseFieldLevelSetImageFilter<
with
     typedef SparseFieldLevelSetImageFilter<

in function
SNAPLevelSetDriver<VDimension>
::DoCreateLevelSetFilter()

in SNAPLevelSetDriver.txx line 143

the thing works as expected, both on initialization and after evolution.

So the culprit of the problem seems to be  
ParallelSparseFieldLevelSetImageFilter.
I'll now look at the modifications I made to the ParallelSparseField  
solver when I introduced anisotropic spacings.


Luca

--
Luca Antiga, PhD
Head, Medical Imaging Unit,
Bioengineering Department,
Mario Negri Institute
email: antiga at marionegri.it
web: http://villacamozzi.marionegri.it/~luca
mail: Villa Camozzi, 24020, Ranica (BG), Italy
phone: +39 035 4535-381


On Aug 6, 2008, at 1:43 AM, Casey Goodlett wrote:

> Hi Luis,
>
> I submitted a bug report in mantis with the requested information  
> here http://www.itk.org/Bug/view.php?id=7454
>
> Thanks
>
> Casey
>
> On Tue, Aug 5, 2008 at 5:21 PM, Luis Ibanez  
> <luis.ibanez at kitware.com> wrote:
>
> Hi Casey,
>
> Could you please give us a hand by providing a description
> of the circumstances in which this error becomes visible ?
>
> E.g. a step-by-step set of instructions will be ideal.
>
> We should then be able to introduce in ITK a new Nightly
> test that exercise this same case, and then we can proceed
> to fix the problem.
>
>
>  Thanks for any advice,
>
>
>      Luis
>
>
> ---------------------
>
> Bill Lorensen wrote:
> Casey,
>
> Snap is an important ITK customer. We'll work with you to figure  
> out the issues.
>
> Bill
>
> On Tue, Aug 5, 2008 at 5:15 PM, Casey Goodlett  
> <gcasey at sci.utah.edu> wrote:
>
> It's a functionality change.
>
> On Tue, Aug 5, 2008 at 3:14 PM, Bill Lorensen  
> <bill.lorensen at gmail.com>
> wrote:
>
> Is this a compilation error or functionality change?
>
> On Tue, Aug 5, 2008 at 5:06 PM, Casey Goodlett <gcasey at sci.utah.edu>
> wrote:
>
> Hi,
>
> I have recently compiled itk-snap 1.6.0.1 against the latest  
> version of
> ITK,
> and the automatic segmentation routines in SNAP no longer behave
> correctly.
> When initializing the segmentation routine with a bubble, the bubble
> does
> not get correctly converted into a label map using newer versions of
> ITK,
> and the level set evolution does not behave correctly.  Going back to
> ITK
> 3.4 fixes the problem.
>
> I've noticed that there were changes made to the level set routines in
> ITK
> with respect to bug 0006156 (
> http://public.kitware.com/Bug/view.php?id=6156
> ).  I've checked out CVS of ITK before and after these changes and can
> confirm that the problem came about in ITK between 2008-03-01 and
> 2008-03-06.
>
> I tried toggling the ITK_DEPRECATED_LEVELSET_INTERPOLATION flag in the
> cmake
> config but there was no effect on the problem.
>
> I'm unsure how to solve the problem so that snap can be used  
> against the
> current version of ITK.  Can someone explain how the behavior of the
> level
> set filters changed?
>
> Thank you
>
> --
> Casey Goodlett
>
> _______________________________________________
> Insight-developers mailing list
> Insight-developers at itk.org
> http://www.itk.org/mailman/listinfo/insight-developers
>
>
>
>
>
>
> _______________________________________________
> Insight-developers mailing list
> Insight-developers at itk.org
> http://www.itk.org/mailman/listinfo/insight-developers