loadPly.m
fb70e1933436feat_muscle_subdivisions

loadPly.m
View Options

			function [vertices,faces] = loadPly(filename,face_index_start)
			%LOADPLY Load the ply files and output verticies and faces
			% Adapted by Alexandre Terrier from
			% https://ch.mathworks.com/matlabcentral/fileexchange/55968-ply-file-utilities

			%% Parse the Header
			fid = fopen(filename,'r');
			readFile = true;
			ii = 0;
			while readFile == true
			fLine = fgetl(fid);
			ii = ii + 1;
			if regexpi(fLine,'element vertex')
			fLine = strsplit(fLine,' ');
			nVertex = str2double(fLine{end});
			continue
			end

			if regexpi(fLine,'element face')
			fLine = strsplit(fLine,' ');
			nFace = str2double(fLine{end});
			continue
			end

			if regexpi(fLine,'end_header')
			readFile = false;
			endHeader = ii;
			continue
			end
			end
			fclose(fid);
			%% Textscan the vertex info
			delimiter = ' ';
			startRow = endHeader+1;
			endRow = startRow + nVertex - 1;
			formatSpec = '%f%f%f%*s%[^\n\r]';
			fid = fopen(filename,'r');
			dataArray = textscan(fid, formatSpec, nVertex, 'Delimiter', delimiter, 'MultipleDelimsAsOne', true, 'EmptyValue' ,NaN,'HeaderLines', startRow-1, 'ReturnOnError', false);
			fclose(fid);
			vertices = [dataArray{1:end-1}];

			%% Textscan the face connectivity
			delimiter = ' ';
			startRow = endRow + 1;
			formatSpec = '%q%f%f%f%s%[^\n\r]';
			fid = fopen(filename,'r');
			dataArray = textscan(fid, formatSpec, 'Delimiter', delimiter, 'MultipleDelimsAsOne', true, 'EmptyValue' ,NaN,'HeaderLines' ,startRow-1, 'ReturnOnError', false);
			fclose(fid);
			faces = [dataArray{1:end-1}];
			faces = faces + face_index_start;


			end