|
|
| (One intermediate revision by one other user not shown) |
| Line 1: |
Line 1: |
| This example demonstrates how to add point and cell data to an itk::Mesh and write it to a .vtk mesh file. Please note that, currently, the only mesh file formats supported by ITK which allow for cell and point data are .vtk and .gii.
| | {{warning|1=The media wiki content on this page is no longer maintained. The examples presented on the https://itk.org/Wiki/* pages likely require ITK version 4.13 or earlier releases. In many cases, the examples on this page no longer conform to the best practices for modern ITK versions.}} |
| | |
| Contributed by: Davis Vigneault
| |
| | |
| ==PointAndCellData.cxx==
| |
| | |
| <source lang="cpp">
| |
| | |
| // Include the relevant header files. | |
| #include "itkMesh.h"
| |
| #include "itkRegularSphereMeshSource.h"
| |
| #include "itkMeshFileWriter.h"
| |
| | |
| // We define the dimension and coordinate type...
| |
| const unsigned int Dimension = 3;
| |
| typedef float TCoordinate;
| |
| | |
| // ...and then typedef the mesh, sphere, and writer. | |
| typedef itk::Mesh< TCoordinate, Dimension > TMesh;
| |
| typedef itk::RegularSphereMeshSource< TMesh > TSphere;
| |
| typedef itk::MeshFileWriter< TMesh > TMeshWriter;
| |
| | |
| int main()
| |
| {
| |
| | |
| // Create the sphere source.
| |
| TSphere::Pointer sphere = TSphere::New();
| |
| sphere->Update();
| |
| | |
| // We now assign it to a mesh pointer.
| |
| TMesh::Pointer mesh = sphere->GetOutput();
| |
|
| |
| // It is necessary to disconnect the mesh from the pipeline;
| |
| // otherwise, the point and cell data will be deallocated | |
| // when we call "Update()" on the writer later in the program.
| |
| mesh->DisconnectPipeline();
| |
| | |
| // Let's assign a value to each of the mesh's points...
| |
| for (unsigned int i = 0; i < mesh->GetNumberOfPoints(); ++i)
| |
| mesh->SetPointData( i, 5.0 );
| |
| | |
| // ...and assign a different value to each of the mesh's cells.
| |
| for (unsigned int i = 0; i < mesh->GetNumberOfCells(); ++i)
| |
| mesh->SetCellData( i, 10.0 );
| |
| | |
| // We'll print out some data about the points...
| |
| std::cout << mesh->GetNumberOfPoints() << std::endl; // 66
| |
| std::cout << mesh->GetPointData()->Size() << std::endl; // 66
| |
| std::cout << mesh->GetPointData()->ElementAt( 0 ) << std::endl << std::endl; // 5.0
| |
|
| |
| // ...and about the cells.
| |
| std::cout << mesh->GetNumberOfCells() << std::endl; // 128
| |
| std::cout << mesh->GetCellData()->Size() << std::endl; // 128
| |
| std::cout << mesh->GetCellData()->ElementAt( 0 ) << std::endl << std::endl; // 10.0
| |
| | |
| // Finally, we'll write the data to file. Note that the only mesh file
| |
| // formats supported by ITK which support cell and point data are .vtk and .gii.
| |
| TMeshWriter::Pointer meshWriter = TMeshWriter::New();
| |
| meshWriter->SetFileName( "mesh.vtk" );
| |
| meshWriter->SetInput( mesh );
| |
| meshWriter->Update();
| |
|
| |
| return EXIT_SUCCESS;
| |
|
| |
| }
| |
| | |
| | |
| </source>
| |
| | |
| {{ITKCMakeLists|{{SUBPAGENAME}}}}
| |
