[Insight-users] itkDeformableMesh3DFilter

[Pieter Vos]

Hi,

I have a 2 questions about the itkDeformableMesh3DFilter.

The algorithm is running ok, thanks to the example 
itkDeformableMeshTest.cxx. 
1) The filter has an input mesh generated by the class
itkSpheresource. When setting up the sphere you have to set the X and Y
resolution. When the deformable filter is being set up you have to set op
the X, Y an Z resolution. In the example X and Y are 10 resp. 30 and Z=1.
My question is: must the resolutions of the sphere and filter be the
same? When I tried to change the resolutions of the deformable filter I
got the error:
% mesh: /data/cemri/Insight/Code/Algorithms/itkDeformableMesh3DFilter.txx:1335: void
% itk::DeformableMesh3DFilter<TInputMesh, TOutputMesh>::ComputeNormals
% () [with TInputMesh = FMesh, TOutputMesh = FMesh]: Assertion `absvec!=0'
% failed.
% Abort

So what is the meaning of resolutions of the filter? Is it
Multiresolution? How does it relate to the spheresource?

question 2:

As an input mesh for the deformable filter, I would like to input a mesh
created by myself. I used the example vtk2itk.cxx where a mesh is filled
with a pointscontainer and a cellcontainer. However when I connected it to
the deformable filter, I got a segmentation fault. With the help of ddd (a
linux debugger) I found out that a vectorcontainer tried to connect to the
datacontainer of the mesh I made. But what data does it need? Or is the
initialisation of a datacontainer next to the points and cells enough?

My question is: how should I set up a self defined mesh (like in
vtk2itk.cxx) and not a spheresource and use it as input for the deformable
filter?

Many thanks for the good help!

-- 
Pieter Vos
UMC Radboud Nijmegen
The Netherlands
mailto:pieterv@radiology.azn.nl