Bugfix: Surface processing: No reducepatch if Nvertices too low

brainstorm-tools · Jan 8, 2022 · 184e8cb · 184e8cb
1 parent bf5e7f4
commit 184e8cb
Show file tree

Hide file tree

Showing 2 changed files with 9 additions and 5 deletions.
diff --git a/toolbox/anatomy/tess_fillholes.m b/toolbox/anatomy/tess_fillholes.m
@@ -84,8 +84,10 @@
 VertConn = tess_vertconn(Vertices, Faces);
 Vertices = tess_smooth(Vertices, 1, 10, VertConn, 0);
 % Downsampling isosurface
-bst_progress('text', 'Fill: Downsampling surface...');
-[Faces, Vertices] = reducepatch(Faces, Vertices, nVertices ./ length(Vertices));
+if (length(Vertices) > nVertices)
+    bst_progress('text', 'Fill: Downsampling surface...');
+    [Faces, Vertices] = reducepatch(Faces, Vertices, nVertices ./ length(Vertices));
+end
 % Remove small components
 bst_progress('text', 'Fill: Removing small components...');
 [Vertices, Faces] = tess_remove_small(Vertices, Faces);

diff --git a/toolbox/anatomy/tess_isohead.m b/toolbox/anatomy/tess_isohead.m
@@ -134,9 +134,11 @@
 bst_progress('inc', 10);
 
 % Downsampling isosurface
-bst_progress('text', 'Downsampling surface...');
-[sHead.Faces, sHead.Vertices] = reducepatch(sHead.Faces, sHead.Vertices, nVertices./length(sHead.Vertices));
-bst_progress('inc', 10);
+if (length(sHead.Vertices) > nVertices)
+    bst_progress('text', 'Downsampling surface...');
+    [sHead.Faces, sHead.Vertices] = reducepatch(sHead.Faces, sHead.Vertices, nVertices./length(sHead.Vertices));
+    bst_progress('inc', 10);
+end
 % Convert to millimeters
 sHead.Vertices = sHead.Vertices(:,[2,1,3]);
 sHead.Faces    = sHead.Faces(:,[2,1,3]);