[Insight-users] itkVTKPolyDataWriter doesn't write edges?

[Bill Lorensen]

I would think you would have to create a line to see the lines in paraview.

On Thu, Sep 17, 2009 at 7:02 PM, David Doria <daviddoria+itk at gmail.com> wrote:
> It seems as though edges are not converted to lines or whatever it is
> necessary to do in a vtk file to view a "wireframe" in Paraview.
>
> An example would be something like this:
>
> #include "itkQuadEdgeMesh.h"
> #include "itkLineCell.h"
> #include "/home/doriad/src/Insight/Code/Review/itkVTKPolyDataWriter.h"
>
> #include <string>
>
> const unsigned int Dimension = 3;
> typedef itk::QuadEdgeMesh< float, Dimension >   MeshType;
>
>
> int main(int argc, char * argv[])
> {
>     if(argc != 2)
>     {
>         std::cout << "Required arguments: Filename" << std::endl;
>         exit(-1);
>     }
>
>     std::string Filename = argv[1];
>
>     MeshType::Pointer  mesh = MeshType::New();
>
>     //create points
>     MeshType::PointType p0,p1,p2,p3;
>
>     p0[0] = -1.0; p0[1]= -1.0; p0[2]= 0.0; // first  point ( -1, -1, 0 )
>     p1[0] =  1.0; p1[1]= -1.0; p1[2]= 0.0; // second point (  1, -1, 0 )
>     p2[0] =  1.0; p2[1]=  1.0; p2[2]= 0.0; // third  point (  1,  1, 0 )
>     p3[0] =  1.0; p3[1]=  1.0; p3[2]= 1.0; // third  point (  1,  1, 1 )
>
>     mesh->SetPoint( 0, p0 );
>     mesh->SetPoint( 1, p1 );
>     mesh->SetPoint( 2, p2 );
>     mesh->SetPoint( 3, p3 );
>
>     std::cout << "NumPoints = " << mesh->GetNumberOfPoints() << std::endl;
>
>     //create 2 edges
>     mesh->AddEdge(0, 1);
>     mesh->AddEdge(1, 2);
>
>     std::cout << "NumEdges = " << mesh->GetNumberOfEdges() << std::endl;
>
>     typedef itk::VTKPolyDataWriter< MeshType >  WriterType;
>     WriterType::Pointer SubmeshWriter = WriterType::New();
>
>     SubmeshWriter->SetFileName(Filename);
>     SubmeshWriter->SetInput(mesh);
>     SubmeshWriter->Update();
>     SubmeshWriter->Write();
>
>     return 0;
> }
>
> I would expect to be able to see these two edges in Paraview. Anyone
> agree/disagree?
>
> Thanks,
>
> David
>
> _____________________________________
> Powered by www.kitware.com
>
> Visit other Kitware open-source projects at
> http://www.kitware.com/opensource/opensource.html
>
> Please keep messages on-topic and check the ITK FAQ at:
> http://www.itk.org/Wiki/ITK_FAQ
>
> Follow this link to subscribe/unsubscribe:
> http://www.itk.org/mailman/listinfo/insight-users
>
>