tsurf

function tsurf(F,V,vertex_indices,face_indices)  % trisurf wrapper  % F: list of faces #F x 3  % V: vertex positiosn #V x 3 or #V x 2  % vertex_indices: show vertex indices on plot  %                 0 -> off  %                 1 -> text and grey background  %                >1 -> text  % face_indices: show face indices on plot  %                 0 -> off  %                 1 -> text and grey background  %                >1 -> text  %  % See corresponding paper: "Mixed finite elements for variational surface  % modeling" by Alec Jacobson, Elif Tosun, Olga Sorkine, and Denis Zorin, SGP  % 2010  %  % Copyright 2010, Alec Jacobson, NYU  %  if(~exist('vertex_indices'))    % off by default    vertex_indices = 0;  end  if(~exist('face_indices'))    % off by default    face_indices = 0;  end  is_2d = false;  if(size(V,2)==2 || (size(V,2) ==3 && sum(abs(V(:,3))) == 0))    V = [V(:,1) V(:,2) 0*V(:,1)];    is_2d = true;  elseif(size(V,2)>3 || size(V,2)<2 )     error('V must be #V x 3 or #V x 2');    return;  end      %trisurf(Tri,X,Y,Z,C)  displays triangles defined in the m-by-3 face   %matrix Tri as a surface. Each row of Tri defines a single triangular  %face by indexing into the vectors or matrices that contain the X, Y, and  %Z vertices. The color is defined by the vector C.  trisurf(F,V(:,1),V(:,2),V(:,3));  % if 2d then set to view (x,y) plane  if( is_2d )    view(2);  end  %标注面索引  if(face_indices==1)    FC = (V(F(:,1),:)+V(F(:,2),:)+V(F(:,3),:))./3;    text(FC(:,1),FC(:,2),FC(:,3),num2str((1:size(F,1))'),'BackgroundColor',[.8 .8 .8]);  elseif(face_indices)    FC = (V(F(:,1),:)+V(F(:,2),:)+V(F(:,3),:))./3;    text(FC(:,1),FC(:,2),FC(:,3),num2str((1:size(F,1))'));  end  %标注点索引  if(vertex_indices==1)    text(V(:,1),V(:,2),V(:,3),num2str((1:size(V,1))'),'BackgroundColor',[.8 .8 .8]);  elseif(vertex_indices)    text(V(:,1),V(:,2),V(:,3),num2str((1:size(V,1))'));  end  % uncomment these to switch to a better 3d surface viewing mode  %axis equal; axis vis3d;end