Glottal outline fits to x subdivs, neighbor based outlier detection

Author: Henningson

source/GUI/MainMenuWidget.py

@@ -24,7 +24,7 @@ def __init__(self, viewer_palette, parent=None):
             "Tensor Product M5",
             [
                 ("Z Subdivisions", "field", 10), 
-                ("X Subdivisions", "field", 2),
+                ("X Subdivisions", "field", 4),
                 ("R_0", "field", 0.1), # Again from Scherer, for Silicone: 1.0, 2.5 is good
                 ("T", "field",   0.3),
                 ("psi", "field", 0.0),

source/SurfaceReconstruction.py

@@ -154,7 +154,7 @@ def generateARAPAnchors(vertices, points, nPointsU, glottalOutlinePoints, x_subd
     # Fit glottal out and midline
     if glottalOutlinePoints.size != 0:
         for i in range(nPointsV):
-            for j in range(x_subdivisions + 2):
+            for j in range(x_subdivisions + 3):
                 controlPointIndex = (x_subdivisions + 2) + j + i*nPointsU
                 controlPoint = vertices[controlPointIndex]
 

source/Triangulation.py

@@ -43,6 +43,14 @@ def triangulationMatNew(camera, laser, laser_correspondences, points_2d, minInte
 
         #Filter NaNs and badly estimated 3-D Positions
         worldPositions = worldPositions[distances < 5.0]
+        
+        # Filter points that have a large distance from neighboring points of the same frame
+        diff = np.expand_dims(worldPositions, 1) - np.expand_dims(worldPositions, 0)
+        dist = np.linalg.norm(diff, axis=2)
+        np.fill_diagonal(dist, np.inf)
+        sorted_indices = np.argsort(dist, axis=1)
+        closest_distances = np.take_along_axis(dist, sorted_indices[:, 0:1], axis=1).flatten()
+        worldPositions = worldPositions[closest_distances < 3.0]
 
         points3D.append((aPoints + ((bPoints - aPoints) / 2.0)))