[Insight-users] Problem with DeformableSimplexMesh3DFilter

[Leila baghdadi]

Hi,

As Thomas pointed out we are working on fixing this class, My changes
are almost done and I am in the process of testing them, unfortunately I
am involved with other projects so things are a bit slow, I will however
commit my changes soon and e-mail the list,

a few things that I can point out,

you need to make sure you fully understand how the deformable models
with simplex mesh and specifically the parameters work (there will be
references mentioned in the documentation)

for instance the locality parameter ( 0 <= locality <= 1), this is
exclusively for switching between global forces ( model-based
registration) versus free form (local deformations), so it simply does
not make sense to set it to 3.0.

We also have an Application (it is part of

InsightApplications/DeformableModelSimplexMesh )

which is also being updated and modified as we speak, you might want to
try at least looking at the help and the sample data before you proceed,

once again these have *not* been committed to cvs but will be soon


Leila

On Fri, 2005-04-22 at 13:05, 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.
> 
> 
>     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
> >
> >  
> >
>