[feat] accelerate vol2surf by switching from binsurface to isosurface

Author: fangq

tools/PoissonRecon

@@ -82,13 +82,14 @@
         maxlevel = length(isovalues);
     end
 
+    % compute the region seed for each region
     for i = 1:maxlevel
         if (i < maxlevel)
             levelmask = int8(newimg >= isovalues(i) & newimg < isovalues(i + 1));
         else
             levelmask = int8(newimg >= isovalues(i));
         end
-        [levelno, levelel] = binsurface(levelmask);
+        [levelel, levelno] = isosurface(levelmask, 0.5);
         if (~isempty(levelel))
             if (isstruct(opt) && isfield(opt, 'autoregion'))
                 if (opt.autoregion)
@@ -195,7 +196,6 @@
             end
             [v0, f0] = vol2restrictedtri(newimg, isovalues(i) - perturb, regions(i, :), ...
                                          sum(newdim .* newdim) * 2, 30, radbound, distbound, maxsurfnode);
-
             if (~isempty(surfside))
                 newimg = newimg0;
                 clear newimg0;