[feat] integrate isosurface with binsurface

Author: fangq

binsurface.m

@@ -11,6 +11,8 @@
 %   nface: nface=3 or ignored - for triangular faces,
 %          nface=4 - square faces
 %          nface=0 - return a boundary mask image via node
+%          nface='iso' - call isosurface() function to return a
+%                    marching-cube based isosurface
 %
 % output:
 %   elem: integer array with dimensions of NE x nface, each row represents
@@ -22,6 +24,12 @@
 % -- this function is part of iso2mesh toolbox (http://iso2mesh.sf.net)
 %
 
+if (nargin > 1 && ischar(nface) && strcmp(nface, 'iso'))
+    [elem, node] = isosurface(img);
+    node = node(:, [2, 1, 3]) - 0.5;
+    return
+end
+
 dim = size(img);
 if (length(dim) < 3)
     dim(3) = 1;

vol2surf.m

@@ -91,8 +91,7 @@
         end
         [isoct, octver] = isoctavemesh;
         if (~isoct || sscanf(octver, '%d') > 5)
-            [levelel, levelno] = isosurface(levelmask, 0.5);
-            levelno = levelno(:, [2 1 3]) - 0.5;
+            [levelno, levelel] = binsurface(levelmask, 'iso');
         else
             [levelno, levelel] = binsurface(levelmask);
         end