ITK  5.2.0
Insight Toolkit
SphinxExamples/src/Core/Mesh/WorkingWithPointAndCellData/Code.cxx
/*=========================================================================
*
* Copyright NumFOCUS
*
* Licensed under the Apache License, Version 2.0 (the "License");
* you may not use this file except in compliance with the License.
* You may obtain a copy of the License at
*
* http://www.apache.org/licenses/LICENSE-2.0.txt
*
* Unless required by applicable law or agreed to in writing, software
* distributed under the License is distributed on an "AS IS" BASIS,
* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
* See the License for the specific language governing permissions and
* limitations under the License.
*
*=========================================================================*/
// Include the relevant header files.
#include "itkMesh.h"
// We define the dimension and coordinate type...
constexpr unsigned int Dimension = 3;
using TCoordinate = float;
// ...and then type alias the mesh, sphere, and writer.
using TMeshWriter = itk::MeshFileWriter<TMesh>;
int
main(int, char *[])
{
// 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;
}
itk::RegularSphereMeshSource
Inputs are the center of the mesh, the scale (radius in each dimension) of the mesh and a resolution ...
Definition: itkRegularSphereMeshSource.h:41
itkRegularSphereMeshSource.h
itk::MeshFileWriter
Writes mesh data to a single file.
Definition: itkMeshFileWriter.h:52
itkMesh.h
itk::Mesh
Implements the N-dimensional mesh structure.
Definition: itkMesh.h:114
itkMeshFileWriter.h
itk::GTest::TypedefsAndConstructors::Dimension2::Dimension
constexpr unsigned int Dimension
Definition: itkGTestTypedefsAndConstructors.h:44