[Insight-users] Problem with DeformableSimplexMesh3DFilter

[Thomas Boettger]

Hello Gunn Marie,

Unfortunately the latest version of the DeformableSimplexMesh3DFilter 
has a bug and would not produce a useful output. We are currently 
working on this. I have no CVS access to commit the fixes right now. I 
will talk to Luis (and Leila, who also needs the classes) next week 
about this topic and hopefully commit the necessary changes.

Now to your problem which what I think has nothing to do with the bug I 
was talking about. How many points does your mesh have? And is your 
program really crashing? Please try to use for example one iteration 
only and see if it finishes.

How did you create the mesh? You were talking about a binary mesh 
source. Could you please post the whole code. Maybe the problem could be 
related to the creation of the mesh.


    Regards,

       Thomas

Gunn Marie Navestad wrote:

>I have a problem when I try to use the DeformableSimplexMesh3DFilter. When
>I try to run my code, everything halts when I update the filter. Can
>somebody please tell me what is wrong? I have tested both the Mesh and the
>gradient image that the filter uses, both look ok.
>
>Here is my code:
>(binaryMeshSource returns a Mesh of type SimplexMeshType, and
>gradientMapFilter an Image of type gradientImageType)
>
>typedef itk::SimplexMesh<double, 3> SimplexMeshType;
>typedef itk::CovariantVector<double, 3> gradientPixelType;
>typedef itk::Image<gradientPixelType, 3> gradientImageType;
>
>itk::DeformableSimplexMesh3DFilter<SimplexMeshType,
>SimplexMeshType>::Pointer deformableModelFilter;
>deformableModelFilter =
>itk::DeformableSimplexMesh3DFilter<SimplexMeshType,
>SimplexMeshType>::New();
>
>deformableModelFilter->SetGradient(gradientMapFilter->GetOutput());
>deformableModelFilter->SetInput(binaryMeshSource->GetOutput());
>deformableModelFilter->SetAlpha(0.1);
>deformableModelFilter->SetBeta(0.3);
>deformableModelFilter->SetGamma(0.03);
>deformableModelFilter->SetDamping(0.4);
>deformableModelFilter->SetRigidity(0.03);
>deformableModelFilter->SetRange(3);
>deformableModelFilter->SetLocality(3.0);
>deformableModelFilter->SetMaxMagnitude(1.0);
>deformableModelFilter->SetIterations(500);
>
>deformableModelFilter->Update();
>
>
>
>Thanks,
>Gunn Marie
>_______________________________________________
>Insight-users mailing list
>Insight-users at itk.org
>http://www.itk.org/mailman/listinfo/insight-users
>
>  
>


-- 
Dipl.-Inform. Thomas Boettger
Deutsches Krebsforschungszentrum         (German Cancer Research Center)
Div. Medical and Biological Informatics B010    Tel: (+49) 6221-42 2328
Im Neuenheimer Feld 280                          Fax: (+49) 6221-42 2345
D-69120 Heidelberg                            e-mail: t.boettger at dkfz.de
Germany                      http://www.dkfz.de/mbi/people/thomasb.shtml