Internal:Adaptive Meshing For APBNRR

Problem Description

Reducing the error for Non-Rigid Registration (NRR) applications depends among other factors on the Equi-distribution of the registration points within the mesh. The goal of this project is to explore metric adaptive methods from CFD and apply them if possible in the area of NRR.

Mesh Requirements

In a set of application specific mesh requirements for NRR were presented.

• (R1) Equi-distribution of the registration points w.r.t. mesh vertex cells  :

  small number of registration points (e.g., less than 10, but greater than 0) within the vertex cell complex makes the formulation more sensitive to outliers and introduces additional displacement error

• (R2) Minimization of the approximation error at registration points :

  error of the displacements recovered at registration points can be reduced locally by using smaller mesh elements.

• (R3) Prevention of tetrahedron inversion during mesh deformation:

  while the interpolation error shown does not depend on the tetrahedron shape, inversion or collapse of a tetrahedron will result in an unrealistic deformation field, e.g., points inside the different tetrahedra can map to the same image location. We can attempt to remedy this problem by adjusting element size locally according to the expected deformation.

Alongside the application-specific requirements we would also like to adhere to the following generic requirements which are related with FEM calculations and image-based mesh generatio in general

• (R4) The mesh boundary (triangulation) should be close to the surface of the image.
• (R5) Mesh size should be minimized to reduce the computational costs.
• (R6) Mesh elements should not have small angles

some of these rules work against each other. For example we require small elements (R2) to reduce the error but we also need a small size (R5) of a mesh.

Methods

Method 1 :Isotropic Sizing Functions based on knn

Already in two attempts were made.

1. Creating new splitting rules in a BCC-based mesh generator and
2. Creating a sizing function for TetGen and NETGEN.

The sizing function of the second approach is based on a distance criterion. In particular:

The sizing value at each mesh vertex was initialized with the distance to the kth registration point closest to it, to reflect the density of the registration point distribution. The idea here is that the shape of a perfect mesh vertex cell complex is close to a ball with the radius prescribed by the background mesh. Where vertex cell complex is defined as the union of all the elements attached to a vertex.