mri_scalp.m

来自「Matlab下的EEG处理程序库」· M 代码 · 共 46 行

function [vertices] = mri_scalp
% MRI_SCALP: Find the scalp surface of an MRI volume
%
%

% $Revision: 1.2 $ $Date: 2003/03/02 03:20:44 $

% Licence:  GNU GPL, no express or implied warranties
% History:  03/2002, [e-mail]
%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

load mri;

D = squeeze(D);

Ds = smooth3(D);
FV = isosurface(Ds,5);

vertices = FV.vertices;


% delaunayn can be used to tesselate these vertices
% dsearchn can be used to identify nearest neighbours

r = 10000;
NFV = reducepatch(FV,r); % reduces the faces of struct fv.


Nvertices = size(NFV.vertices,1);

fprintf('\n...Found %d vertices.\n\n',Nvertices);

plot = 1;
if isequal(plot,1),
    P = patch(NFV,'FaceColor',[1,.75,.65],'EdgeColor','none');
    isonormals(Ds,P);
    view(45,30), axis tight, daspect([1,1,.4]), rotate3D;
    lightangle(45,30); 
    set(gcf,'Renderer','zbuffer'); lighting phong
    set(P,'SpecularColorReflectance',0,'SpecularExponent',50);
end


return