This function I just wrote seems to do the trick:
function A = AVLread(fnm)
% Read the entire file into one string:
txt = fileread(fnm);
% Find sections starting 'key: value':
[S,E,T] = regexpi(txt, '^\s*([A-Z]*):\s*(\w*?)\s*$', 'start', 'end', 'tokens', 'lineanchors');
T = vertcat(T{:});
S = [S,numel(txt)];
% Get all numeric values in the first section:
X = regexpi(txt(E(1):S(2)), '([A-Z]+)\s*=\s*(\d+(\.\d+)?)', 'tokens');
X = [X{:}];
X(2:2:end) = cellfun(@(s) sscanf(s,'%f'), X(2:2:end), 'UniformOutput',false);
A = struct(X{:});
% Get the title and configuration:
A.(T{1,1}) = T{1,2};
A.Title = regexpi(txt, '(?<=^\s*)[A-Z].+?(?=\s*$)', 'match', 'once', 'lineanchors');
% Get each remaining section:
fun = @(z,s) regexpi(z, ['(?<=',s,'\S+\s*=\s*)\S+(?=\s*$)'], 'match', 'once', 'lineanchors');
for k = 2:numel(E)
N = sprintf('%s_%s',T{k,:});
Z = txt(1+E(k):S(k+1));
% Get surface type:
A.(N).Type = regexpi(Z, '(?<=\s*)\w+', 'match', 'once');
% Get min/max:
A.(N).min2Y = sscanf(fun(Z,'min'),'%f');
A.(N).max2Y = sscanf(fun(Z,'max'),'%f');
% Indices of the matrix:
[s,e] = regexpi(Z, '(^\s*\S+\s+[^\s=]+\s+\S+\s*$)+', 'lineanchors');
% Get matrix header:
A.(N).Header = regexpi(Z(s(1):e(1)), '(?<=\s*)\S+', 'match');
% Get matrix data:
A.(N).Data = reshape(sscanf(Z(s(2):e(end)),'%f'), 3,[]).';
end
end
Just call it with the filename of the data file:
A = AVLread('test.vm')
Note that it returns all of data in one structure. All numeric data are converted to double. If you double-click the structure in your workspace then you can browse through it: the arrangement should be fairly self-explanatory.
