Skip to main content

πŸ•ΈοΈ Mesh Issues and Fixes

There are a variety of issues that may be present in a mesh:

  • small components
  • disconnected vertices
  • degenerate faces
  • multiple edges
  • self-intersections
  • holes

⚠️ Non-manifoldness is already fixed on load or construct, so it cannot be an issue here.

πŸ”Ή Small components​

# remove small components first
area = mesh.area()
bigComps = mm.MeshComponents.getLargeByAreaComponents(mesh, area * 0.002, None)
mesh.deleteFaces(mesh.topology.getValidFaces() - bigComps)

πŸ”Ή Disconnected vertices​

(from test_python/test_meshBuilder.py, line 77 in 09b543e)

mrmesh.uniteCloseVertices(mesh, 0)

πŸ”Ή Degenerate faces​

tolerance = mesh.computeBoundingBox().diagonal() * 1e-2

degFaces = mm.findDegenerateFaces(mesh, 1e3)
shortEdges = mm.findShortEdges(mesh, tolerance * 1e-2)
degFaces |= mm.getIncidentFaces(mesh.topology, shortEdges)
mm.expand(mesh.topology, degFaces, 3)

dSettings = mm.DecimateSettings()
dSettings.strategy = mm.DecimateStrategy.ShortestEdgeFirst
dSettings.maxError = tolerance * 0.1
dSettings.criticalTriAspectRatio = 1e3
dSettings.tinyEdgeLength = tolerance * 1e-2
dSettings.stabilizer = mm.ResolveMeshDegenSettings().stabilizer
dSettings.optimizeVertexPos = False
dSettings.region = degFaces
dSettings.maxAngleChange = mm.ResolveMeshDegenSettings().maxAngleChange

mm.decimateMesh(mesh, dSettings)

πŸ”Ή Multiple edges​

# just find them
multipleEdges = mm.findMultipleEdges(mesh.topology)

# fix possible multiple edges
mm.fixMultipleEdges(mesh)

πŸ”Ή Self-intersections​

# remove self-intersections
selfIntersSettings = mm.SelfIntersections.Settings()
selfIntersSettings.method = mm.SelfIntersections.Settings.Method.CutAndFill
selfIntersSettings.maxExpand = 2
selfIntersSettings.relaxIterations = 2

mm.SelfIntersections.fix(mesh, selfIntersSettings)

πŸ”Ή Holes​

# fill all holes if needed
for e in mesh.topology.findHoleRepresentiveEdges():
    mm.fillHole(mesh, e)

πŸ”Ή One-shot fix (voxel rebuild)​

There is also a single function that fixes holes, self-intersections, multiple edges, and most degenerations (though it may leave some small components):

params = mm.RebuildMeshSettings()
params.voxelSize = mm.suggestVoxelSize(mesh, 5e6)
params.fwn = mc.FastWindingNumber(mesh)  # Enables CUDA (mc = `from meshlib import mrcudapy as mc`)
mesh = mm.rebuildMesh(mesh, params)