Headmodel with EEG questions

tmedani · May 20, 2021, 6:13pm

could you change on your computers (both windows and mac) this line

       [node,elem] = surf2mesh(newnode, newelem, min(newnode), max(newnode),...
            OPTIONS.KeepRatio, factor_bst .* OPTIONS.MaxVol, regions, [], [], 'tetgen1.5');

to

[node,elem] = surf2mesh(newnode, newelem, min(newnode), max(newnode),...
OPTIONS.KeepRatio, factor_bst .* OPTIONS.MaxVol, regions, , );

just remove the ''tetgen1.5'' arguments?

you can find it in the function process_fem_mesh at line 423

and run it again, then let us know what you get

github.com

brainstorm-tools/brainstorm3/blob/1a6df44494c4cf91c43666a17509c02bd7facc91/toolbox/process/functions/process_fem_mesh.m#L432


      
          iFace = iFace(end-nBem+1:end);
          dist = dist(end-nBem+1:end);
          % Define region seeds for all the BEM regions: head center, then half-way between each layer
          dist = dist(:);
          distSeed = [0; (dist(1:end-1) + dist(2:end)) .* 0.5];
          regions = repmat(orig, nBem, 1) + distSeed * v0;
          
          % Create tetrahedral mesh
          bst_progress('text', 'Creating 3D mesh (Iso2mesh/surf2mesh)...');
          factor_bst = 1.e-6;
          [node,elem] = surf2mesh(newnode, newelem, min(newnode), max(newnode),...
              OPTIONS.KeepRatio, factor_bst .* OPTIONS.MaxVol, regions, [], [], 'tetgen1.5');
          
          % Removing the label 0 (Tetgen 1.4) or higher than number of layers (Tetgen 1.5)
          bst_progress('text', 'Fixing 3D mesh...');
          iOther = find((elem(:,5) == 0) & (elem(:,5) > nBem));
          if ~isempty(iOther) && (length(iOther) < 0.1 * length(elem))
              elem(iOther,:) = [];
          end
          % Check labelling from 1 to nBem
          allLabels = unique(elem(:,5));