[Insight-users] Writing mesh to Metafile

[Julien Jomier]

Hi Dominique,

We are having the same issue with the SimplexMesh class which derives 
from itk::Mesh.
This is a problem with classes deriving from itk::Mesh and because of 
some cast issues there is no way to write them out at this point using 
the MeshSpatialObject.

So you have two options:

1) If you plan to make the itkIndexedTriangleMesh class part of ITK, I'd 
  be happy to create a new itkIndexedTriangleMeshSpatialObject so you 
can write this type of object out.

2) If you can convert the itkIndexedTriangleMesh to an itkMesh structure 
then you can use the current writer.

I'll investigate if there is a better way to make the writer more 
generic in the future.

let us know if that would work for you,

Julien

Dominique Zosso wrote:
> Hello,
> 
> I'm trying to save an itkIndexedTriangleMesh (author : Tobias Heimann.
> Division Medical and Biological Informatics,  German Cancer Research Center,
> Heidelberg, Germany ) in some meta-file.
> 
> itkIndexedTriangleMesh is a class derived from itkMesh offering additional
> CellLink-functionalities.
> 
> According to the SoftwareGuide, I put the mesh into a MeshSpatialObject and
> try to write this latter into a meta-file. The source code doing this is
> given below.
> 
> The evolveMesh-Output is already calculated when it is assigned to the
> meshSpatialObject. No need to update. (see Meta-Output)
> 
> The program crashes as soon as I try to update the SpatialObjectWriter, with
> a segmentation fault. There is not even an itkException being catched by the
> block. Does anybody have an idea what could be wrong? Or is there another
> way to save a mesh?
> 
> Thank you for every piece of advice!!
> 
> For your information, I'm using ITK 2.2.0, and c++ (GCC) 4.0.2 20050808
> (prerelease)  on Ubuntu 4.0.1-4ubuntu9.
> 
> 
> Thanks a lot, Dominique
> 
> 
> Piece of code producing the crash:
> **************************
> 
>   MeshSpatialObjectType::Pointer meshSpatialObject =
> MeshSpatialObjectType::New();
>   meshSpatialObject->SetMesh( evolveMesh->GetOutput() );
> 
>   MeshWriterType::Pointer meshWriter = MeshWriterType::New();
>   meshWriter->SetInput( meshSpatialObject );
>   meshWriter->SetFileName( "mesh.meta" );
> 
>   std::cout << std::endl << " meshSpatialObject : " << std::endl <<
> meshSpatialObject << std::endl;
>   std::cout << std::endl << " meshWriter : " << std::endl << meshWriter <<
> std::endl;
> 
>   try
>   {
>     meshWriter->Update();
>   }
>   catch ( itk::ExceptionObject & err )
>   {
>     std::cerr << " Could not write mesh meta file. Going down *tilt* " <<
> std::endl;
>     std::cerr << err << std::endl;
>     return EXIT_FAILURE;
>   }
> 
> 
> 
> And the console output produced before crashing:
> **************************************
> 
> meshSpatialObject :
> MeshSpatialObject (0x82e0a60)
>   RTTI typeinfo:   itk::MeshSpatialObject<itk::IndexedTriangleMesh<short,
> unsigned int, double> >
>   Reference Count: 3
>   Modified Time: 26050
>   Debug: Off
>   Observers:
>     none
>   Source: (none)
>   Source output index:  0
>   Release Data: Off
>   Data Released: False
>   Global Release Data: Off
>   PipelineMTime: 0
>   UpdateMTime: 0
> Bounding Box:
>   Object (0x82e08f0)
>   RTTI typeinfo:   itk::BoundingBox<unsigned long, 3, double,
> itk::VectorContainer<unsigned long, itk::Point<double, 3u> > >
>   Reference Count: 2
>   Modified Time: 26004
>   Debug: Off
>   Observers:
>     none
>   Bounding Box: ( 60.3259,140.075 60.3259,140.072 60.2,140.6  )
> 
> Geometric properties:
>   Object to World Transform: ScalableAffineTransform (0x8412150)
>   RTTI typeinfo:   itk::ScalableAffineTransform<double, 3u>
>   Reference Count: 2
>   Modified Time: 26009
>   Debug: Off
>   Observers:
>     none
>   Matrix:
>     1 0 0
>     0 1 0
>     0 0 1
>   Offset: [0, 0, 0]
>   Center: [0, 0, 0]
>   Translation: [0, 0, 0]
>   Inverse:
>     1 0 0
>     0 1 0
>     0 0 1
>   Singular: 0
>   Scale : 1 1 1
>   MatrixScale : 1 1 1
> 
>   Index to World Transform: ScalableAffineTransform (0x83e6930)
>   RTTI typeinfo:   itk::ScalableAffineTransform<double, 3u>
>   Reference Count: 3
>   Modified Time: 26017
>   Debug: Off
>   Observers:
>     none
>   Matrix:
>     1 0 0
>     0 1 0
>     0 0 1
>   Offset: [0, 0, 0]
>   Center: [0, 0, 0]
>   Translation: [0, 0, 0]
>   Inverse:
>     1 0 0
>     0 1 0
>     0 0 1
>   Singular: 0
>   Scale : 1 1 1
>   MatrixScale : 1 1 1
> 
> 
> 
>   Bounding Box Children Depth: 9999999
>   Bounding Box Children Name:
> Object properties:
> SpatialObjectProperty (0x82e0960)
>   RTTI typeinfo:   itk::SpatialObjectProperty<float>
>   Reference Count: 2
>   Name:
>   RGBA: 1 1 1
> 
> Mesh:
> m_IsInsidePrecision: 1
>   IndexedTriangleMesh (0x82e0b80)
>   RTTI typeinfo:   itk::IndexedTriangleMesh<short, unsigned int, double>
>   Reference Count: 3
>   Modified Time: 26050
>   Debug: Off
>   Observers:
>     none
>   Source: (0x83151e0)
>   Source output index: 0
>   Release Data: Off
>   Data Released: False
>   Global Release Data: Off
>   PipelineMTime: 22445
>   UpdateMTime: 25842
>   Number Of Points: 402
>   Requested Number Of Regions: 1
>   Requested Region: 0
>   Buffered Region: 0
>   Maximum Number Of Regions: 1
>   Number Of Points: 402
>   Number Of Cell Links: 402
>   Number Of Cells: 800
>   Size of Cell Data Container: 800
>   Number of explicit cell boundary assignments: 2
>   CellsAllocationMethod: 3
>   Number Of Edges: 1200
>   Number Of Faces: 800
> 
> 
> 
>  meshWriter :
> SpatialObjectWriter (0x8432e38)
>   RTTI typeinfo:   itk::SpatialObjectWriter<3u, short,
> itk::DefaultStaticMeshTraits<short, 3u, 2u, double, double, short> >
>   Reference Count: 2
>   Modified Time: 26056
>   Debug: Off
>   Observers:
>     none
> 
> Erreur de segmentation
> 
> 
> ---
> Dominique Zosso
> cand. ing. el. EPFL
> 
> 
> _______________________________________________
> Insight-users mailing list
> Insight-users at itk.org
> http://www.itk.org/mailman/listinfo/insight-users
> 
>