Mapped Hexahedral Meshing

Objective:

• Develop a tool for mapped meshing that will warp a template mesh onto a surface for a new subject
• Determine the limits that the warping algorithm will succeed and fail
• Apply these techniques to a sample of data to determine the reliability of the algorithm with respect to mesh quality

Progress:

• Warped meshing framework has been developed using VTK and ITK.
  • VTK
    • Support the reading and writing of meshes (vtkUnstructuredGrid) and surfaces (vtkPolyData)
    • Distance measures between the deforming mesh and the subject surface
    • Initial registration - Iterative Closest Point or Procrustes

• Meshing module is now built by default with Slicer3

To Do:

• We should remove the need for the surface extraction when Slicer3 supports the loading of unstructured grids.

Key Investigators:

• Iowa: Nicole Grosland, Vincent Magnotta, Ritesh Bafna

Links:

• Musculoskeletal Imaging, Modelling and Experimentation (MIMX)

Description: This program will convert an image into a hexahedral mesh using the voxels as the hexahedral elements. This is is typically done for only a portion of the image and thus a binary image from a segmentation may be used to define the region for meshing. Since these meshes may be fairly large the user is allowed to downsample the image by an integer value. Currently this must be odd.

This tool was developed as part of an NIH NA-MIC Collaboration Grant EB005973.

Program Usage:

VoxelMeshingModule 

 --imagefilename      Input image File Name
 --maskfilename       Input mask File Name
 --abaqusfilename     Output Abaqus Filename
 --vtkfilename        Output mesh VTK filename
 --vtksurfacefilename Output mesh surface VTK filename (Hack required for slicer3)
 --resamplesize       The image resample downsampling factor - Must be odd
 --meshindexoffset    The mesh index offset for nodes and elements
 --poissonratio       Poisson Ratio for the mesh
 --abaqustitle        Abaqus title for generated mesh
 --numberofbins       The number of histogram bins in the mesh
 --histogamlowerbound The lower bound for the mesh histogram
 --histogamupperbound The upper bound for the mesh histogram
 --histogambinfile    Mesh material properties histogram bin filename

Figures:

Surface of proximal phalanx generated via Slicer3

Voxel mesh generated via Slicer3 execution Model

Slicer3 Voxel Mesh Module