3D Morphable Model pca algorithm

function [face,vert,extra] = read_obj(filename) fid = fopen(filename); if fid<0 error(['cannot open ' filename '.']); end face = []; vert = []; vt = []; vn = []; while true line = get_next_line(fid); if line == -1 break; end if length(line) < 2 continue end [format,sz] = get_format(line); if strcmpi(line(1:2), 'vt') fseek(fid,-length(line),'cof'); data = fscanf(fid,format); vt = [vt;reshape(data,[sz,length(data)/sz])']; fseek(fid,-1,'cof'); elseif strcmpi(line(1:2), 'vn') fseek(fid,-length(line),'cof'); data = fscanf(fid,format); vn = [vn;reshape(data,[sz,length(data)/sz])']; fseek(fid,-1,'cof'); elseif strcmpi(line(1), 'v') fseek(fid,-length(line),'cof'); data = fscanf(fid,format); vert = [vert;reshape(data,[sz,length(data)/sz])']; fseek(fid,-1,'cof'); elseif strcmpi(line(1), 'f') fseek(fid,-length(line),'cof'); data = fscanf(fid,format); face = [face;reshape(data,[sz,length(data)/sz])']; fseek(fid,-1,'cof'); end end fclose(fid); has_vt = ~isempty(vt); has_vn = ~isempty(vn); if has_vt extra.vt = vt; end if has_vn extra.vn = vn; end extra.face = face; extra.vert = vert; if size(face,2) == 6 face = face(:,[1 3 5]); elseif size(face,2) == 9 face = face(:,[1 4 7]); end if size(vert,2) > 3 vert = vert(:,1:3); end function line = get_next_line(fid) % read next line, skip comment, blank line and eof line = fgets(fid); if line == -1 return end while isempty(strtrim(line)) || line(1) == '#' line = fgets(fid); if line == -1 line = ''; return; end end function [format,sz] = get_format(str) % determine the format of input str sn = '[\-+]?(?:\d*\.|)\d+\.?(?:[eE][\-+]?\d+|)'; % match number format = regexprep(str,sn,'%f'); [~,splitstr] = regexp(str,sn,'match'); sz = length(splitstr);