[Insight-users] I need help please!! DeformableMesh3DFilter

Waltraud Henrich henrich@ira.uka.de
Thu, 9 Jan 2003 11:31:43 +0100


Hallo,

at 16 December I had written my problem in the mailing list.  Unfortunately
I received still no answer and I did not repair the problem .
I have a newer version of ITK and I  downloaded InsightToolkit-1.0.0 to
28.11.2002.
Where can I look for the errors? 
What are the meaning of the parameters: Stiffness, Scale,Time Step and Step
Threshold in the  ItkDeformableMesh3DFilter?
In the example itkDeformableTest.cxx it does not use these parameters. Why?
Is sufficient in such a way? 
I have tried it also without these parameters and it  does not go. 
 I tried with this parameter: 


	double m_stiff[2];
	m_stiff[0]=0.001;
	m_stiff[1]=0.04;
	
	myIndexType m_center={{SEEDX, SEEDY, SEEDZ}};
	
	Deformable3DFilterType::Pointer m_bfilter =
Deformable3DFilterType::New();
	
	m_bfilter->SetCenter(m_center);
	m_bfilter->SetStiffness(m_stiff);
	m_bfilter->SetTimeStep(0.001);
	m_bfilter->SetStepThreshold(20);
	m_bfilter->SetGradient(gradient->GetOutput());
	m_bfilter->SetInput(m_spheresource->GetOutput()); 
	m_bfilter->SetInput(initmesh);
	m_bfilter->Update();

The mail from 16.12 was: 
"
1) I work with itkDeformableMesh3DFilter and now I get the OutputWindow. 

Debug: In
d:\programme\development\insighttoolkit-1.0.0\code\common\itkpointset.txx,
line 74
Mesh (088D9178): returning Points container of VectorContainer (013B4AA8)
  Reference Count: 2
  Modified Time: 9838
  Debug: Off
  Observers: 
    none



Debug: In
D:\Programme\Development\InsightToolkit-1.0.0\Code\Common\itkDataObject.cxx,
line 239
Mesh (088D9178): disconnecting source  088D7D98, source output index 0

Debug: In
D:\Programme\Development\InsightToolkit-1.0.0\Code\Common\itkObject.cxx,
line 288
Mesh (088D9178): UnRegistered, ReferenceCount = 0

Debug: In
d:\programme\development\insighttoolkit-1.0.0\code\common\itkmesh.txx, line
1061
Mesh (088D9178): Mesh Destructor 

Debug: In
d:\programme\development\insighttoolkit-1.0.0\code\common\itkmesh.txx, line
1078
Mesh (088D9178): Mesh  ReleaseCellsMemory method 

Debug: In
d:\programme\development\insighttoolkit-1.0.0\code\common\itkmesh.txx, line
1104
Mesh (088D9178): m_CellsContainer->GetReferenceCount()= 1

Debug: In
d:\programme\development\insighttoolkit-1.0.0\code\common\itkmesh.txx, line
1149
Mesh (088D9178): CellsAllocatedDynamicallyCellByCell

Debug: In
d:\programme\development\insighttoolkit-1.0.0\code\common\itkmesh.txx, line
1168
Mesh (088D9178): Mesh  ReleaseBoundariesMemory method 

Debug: In
d:\programme\development\insighttoolkit-1.0.0\code\common\itkmesh.txx, line
1195
Mesh (088D9178): Mesh  ReleaseBoundariesMemory(unsigned int) method 

Debug: In
d:\programme\development\insighttoolkit-1.0.0\code\common\itkmesh.txx, line
1286
Mesh (088D9178): Boundaries already deleted.

Debug: In
d:\programme\development\insighttoolkit-1.0.0\code\common\itkmesh.txx, line
1195
Mesh (088D9178): Mesh  ReleaseBoundariesMemory(unsigned int) method 

Debug: In
d:\programme\development\insighttoolkit-1.0.0\code\common\itkmesh.txx, line
1286
Mesh (088D9178): Boundaries already deleted.

Debug: In
d:\programme\development\insighttoolkit-1.0.0\code\common\itkmesh.txx, line
1195
Mesh (088D9178): Mesh  ReleaseBoundariesMemory(unsigned int) method 

Debug: In
d:\programme\development\insighttoolkit-1.0.0\code\common\itkmesh.txx, line
1286
Mesh (088D9178): Boundaries already deleted.

Debug: In
D:\Programme\Development\InsightToolkit-1.0.0\Code\Common\itkObject.cxx,
line 447
Object (088D9178): Destructing!


2) When I save vtkUnstructuredGrid in a file and when I look this file I get
following:

# vtk DataFile Version 3.0
vtk output
ASCII
DATASET UNSTRUCTURED_GRID
POINTS 402 float
-1.#INF -1.#INF 1.#INF -1.#INF -1.#INF 1.#INF -1.#INF -1.#INF 1.#INF 
-1.#INF -1.#INF 1.#INF -1.#INF -1.#INF 1.#INF -1.#INF -1.#INF 1.#INF 
-1.#INF -1.#INF 1.#INF -1.#INF -1.#INF 1.#INF -1.#INF -1.#INF 1.#INF 
-1.#INF -1.#INF 1.#INF -1.#INF -1.#INF 1.#INF -1.#INF -1.#INF 1.#INF 
-1.#INF -1.#INF 1.#INF -1.#INF -1.#INF 1.#IN.......etc

My problem is that in vtk2itkDouble.cxx this coordinates are #INF..

pcoordinates[0] = pp[0];
pcoordinates[1] = pp[1];
pcoordinates[2] = pp[2];

Can you tell me where the problem is?

3) I can disply my SphereSource with vtk2itkDouble. The coordinates are OK
here."

Thanks a lot.

Waltraut Henrich