[Insight-users] Re: Insight-users Digest, Vol 37, Issue 49

[Catherine Peloquin]

Andriy,

Thank you for your assistance.  I am working to integrate your code snippet
with my current program.  I included itkCommand.h, but am getting the error
messages "expected primary-expression before "template"" and "expected `;'
before "template"" at the line "template<class TInputMesh> void".

I did not forget a ';' in the lines prior to the class creation.  I even
tried re-creating the makefile with ccmake before attempting to compile the
code, but get the same error messages.

Do you have any suggestions as to what I may be missing?

Thank you,
Catherine Peloquin


On 5/18/07, Andriy Fedorov <fedorov at bwh.harvard.edu> wrote:
>
> Catherine,
>
> Below is a code snippet, which saves a mesh in RAW format (very
> simple), given ITK surface mesh. Hope this helps you.
>
> Andriy Fedorov
>
>
>
> template<class TInputMesh> void
> TriangularSurfaceMeshWriter<TInputMesh>::WriteRawMesh(){
>       std::ofstream raw_mesh(m_FileName.c_str());
>       raw_mesh << m_InputMesh->GetPoints()->Size() << " "
>         << m_InputMesh->GetCells()->Size() << std::endl;
>
>       typename InputMeshType::PointsContainer::ConstIterator inPointsI;
>       typename InputMeshType::CellsContainer::ConstIterator inCellsI;
>       inPointsI = m_InputMesh->GetPoints()->Begin();
>       while(inPointsI != m_InputMesh->GetPoints()->End()){
>         typename InputMeshType::PointType curPoint;
>         curPoint = inPointsI.Value();
>         raw_mesh << curPoint[0] << " " << curPoint[1] << " " <<
> curPoint[2] << std::endl;
>         inPointsI++;
>       }
>
>       inCellsI = m_InputMesh->GetCells()->Begin();
>       typedef itk::TriangleCell<typename InputMeshType::CellType>
> TriangleType;
>       while(inCellsI != m_InputMesh->GetCells()->End()){
>         TriangleType *curTri;
>         typename TriangleType::PointIdIterator ptI;
>         if(inCellsI.Value()->GetType() !=
> InputMeshType::CellType::TRIANGLE_CELL)
>           itkExceptionMacro(<<"Input surface mesh is not triangular");
>         curTri = dynamic_cast<TriangleType*>(inCellsI.Value());
>         ptI = curTri->PointIdsBegin();
>         raw_mesh << *ptI++ << " " << *ptI++ << " " << *ptI++ << std::endl;
>         inCellsI++;
>       }
>       raw_mesh.close();
> }
>
>
> > Message: 5
> > Date: Fri, 18 May 2007 10:48:03 -0500
> > From: "Catherine Peloquin" <cepeloquin at gmail.com>
> > Subject: [Insight-users] Retrieving coordinates from
> >         BinaryMask3DMeshSource  mesh
> > To: insight-users at itk.org
> > Message-ID:
> >         <3a72f1140705180848h35f387f2h6e418e3aa2a62ba3 at mail.gmail.com>
> > Content-Type: text/plain; charset="iso-8859-1"
> >
> > I am using the BinaryMask3DMeshSource filter example in
> > Insight/Examples/Filters/SurfaceExtraction.cxx to create an isosurface
> of a
> > series of DICOM images.  I need to retrieve the coordinates of the
> output
> > mesh's nodes (or the vertices of the triangles formed by the mesh) so I
> can
> > then draw the surface in Coin3D.  I realize that I can also draw the
> surface
> > in VTK, but I need to integrate the program into an application that
> uses
> > Coin.
> >
> > In a previous post (located at
> > http://public.kitware.com/pipermail/insight-users/2007-May/022211.html)
> , I
> > attempted to follow the methods of retrieving coordinates in the
> PointSet
> > example in Insight/Examples/DataRepresentation/Mesh/PointSet2.cxx, but
> could
> > not find a way to get the mesh points for the points container.
> >
> > How is it possible to get the x, y, and z coordinates of a mesh
> generated by
> > the BinaryMask3DMeshSource filter?
> >
> > Thank you,
> > Catherine Peloquin
> > -------------- next part --------------
> > An HTML attachment was scrubbed...
> > URL:
> http://public.kitware.com/pipermail/insight-users/attachments/20070518/4bc94481/attachment-0001.htm
> >
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://public.kitware.com/pipermail/insight-users/attachments/20070518/c0171bd8/attachment.html