[Insight-users] Problem with DeformableSimplexMesh3DFilter

[Gunn Marie Navestad]

On Fri, 22 Apr 2005, Thomas Boettger wrote:

> 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.
>

The mesh has 29776 points, and it is created from a binary 3D image
through the itkBinaryMeshSource. The program is also crashing when I try
to use one iteration only. I have also tried to change the snake locality
as Leila suggested.

>
>     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
>
>