Homec4science
Diffusion Instrument Control (confint-fix) Browse

rINSTCONTROL/@MyMetadata/MyMetadata.m
f762f246ff0econfint-fix

Last Change
Raw File

MyMetadata.m
View Options

classdef MyMetadata < handle
properties (Access = public)
% Header sections are separated by [hdr_spec, title, hdr_spec]
title = ''
hdr_spec = '=='
% Columns are separated by this symbol (space-tab by default)
column_sep = ' \t'
% Comments start from this symbol
comment_sep = '%'
line_sep = '\r\n'
% Limit for column padding. Variables which take more space than
% this limit are ignored when calculating the padding length.
pad_lim = 15
end
properties (GetAccess = public, SetAccess = protected)
ParamList = struct()
end
methods (Access = public)
function this = MyMetadata(varargin)
P = MyClassParser(this);
processInputs(P, this, varargin{:});
end
%Adds a new metadata parameter.
function addParam(this, param_name, value, varargin)
assert(isvarname(param_name), ['Parameter name must be a ' ...
'valid variable name.']);
p = inputParser();
% Format specifier for printing the value
addParameter(p, 'fmt_spec', '', @ischar);
% Comment to be added to the line
addParameter(p, 'comment', '', @ischar);
addParameter(p, 'SubStruct', struct('type',{},'subs',{}),...
@isstruct)
parse(p,varargin{:});
% Make sure that the comment does not contain new line or
% carriage return characters, which would mess up formating
% when saving the metadata
[comment, is_mod] = toSingleLine(p.Results.comment);
this.ParamList.(param_name).comment = comment;
if is_mod
warning(['Comment string for ''%s'' has been ' ...
'converted to single line.'], param_name);
end
S = p.Results.SubStruct;
if isempty(S)
% Assign value directly
this.ParamList.(param_name).value = value;
else
% Assign using subref structure
tmp = feval([class(value),'.empty']);
this.ParamList.(param_name).value = subsasgn(tmp,S,value);
end
this.ParamList.(param_name).fmt_spec = p.Results.fmt_spec;
end
% The function below is useful to ensure the correspondence between
% metadata parameter names and object property names. It spares
% some lines of code.
function addObjProp(this, Obj, tag, varargin)
addParam(this, tag, Obj.(tag), varargin{:});
end
% Print metadata in a readable form
function str = print(this)
% Compose the list of parameter names expanded over subscripts
% except for those which are already character arrays
par_names = fieldnames(this.ParamList);
% Expand parameters over subscripts, except for the character
% arrays
exp_par_names = cell(1, length(par_names));
maxnmarr = zeros(1, length(par_names));
for i=1:length(par_names)
tmpval = this.ParamList.(par_names{i}).value;
exp_par_names{i} = printSubs(tmpval, ...
'own_name', par_names{i}, ...
'expansion_test',@(y) ~ischar(y));
% Max name length for this parameter including subscripts
maxnmarr(i)=max(cellfun(@(x) length(x), exp_par_names{i}));
end
% Calculate width of the name column
name_pad_length = min(max(maxnmarr), this.pad_lim);
% Compose list of parameter values converted to char strings
par_strs = cell(1, length(par_names));
% Width of the values column will be the maximum parameter
% string width
val_pad_length = 0;
for i=1:length(par_names)
TmpPar = this.ParamList.(par_names{i});
for j=1:length(exp_par_names{i})
tmpnm = exp_par_names{i}{j};
TmpS = str2substruct(tmpnm);
if isempty(TmpS)
tmpval = TmpPar.value;
else
tmpval = subsref(TmpPar.value, TmpS);
end
%Do the check to detect unsupported data type
if ischar(tmpval)&&~isvector(tmpval)&&~isempty(tmpval)
warning(['Argument ''%s'' is a multi-dimensional ',...
'character array. It will be converted to ',...
'single string during saving. Use cell',...
'arrays to save data as a set of separate ',...
'strings.'],tmpnm)
% Flatten
tmpval = tmpval(:);
end
%Check for new line symbols in strings
if (ischar(tmpval)||isstring(tmpval)) && ...
any(ismember({newline,sprintf('\r')},tmpval))
warning(['String value must not contain ',...
'''\\n'' and ''\\r'' symbols, replacing ',...
'them with '' ''.']);
tmpval=replace(tmpval,{newline,sprintf('\r')},' ');
end
if isempty(TmpPar.fmt_spec)
% Convert to string with format specifier
% extracted from the varaible calss
par_strs{i}{j} = var2str(tmpval);
else
par_strs{i}{j} = sprintf(TmpPar.fmt_spec, tmpval);
end
% Find maximum length to determine the colum width,
% but, for beauty, do not account for variables with
% excessively long value strings
tmplen = length(par_strs{i}{j});
if (val_pad_length<tmplen) && (tmplen<=this.pad_lim)
val_pad_length = tmplen;
end
end
end
cs = this.column_sep;
ls = this.line_sep;
% Make the output string. Start by printing the title.
str = sprintf([this.hdr_spec, this.title, this.hdr_spec, ls]);
par_fmt_spec = [sprintf('%%-%is', name_pad_length),...
cs, sprintf('%%-%is', val_pad_length)];
% Print parameters
for i=1:length(par_names)
% Capitalize the first letter of comment
comment = this.ParamList.(par_names{i}).comment;
if ~isempty(comment)
fmt_comment = [this.comment_sep, ' '...
upper(comment(1)), comment(2:end)];
else
fmt_comment = '';
end
% Iterate over the parameter subscripts
for j=1:length(exp_par_names{i})
if j==1
% Add the comment to first line
str = [str, ...
sprintf([par_fmt_spec, cs, '%s', ls], ...
exp_par_names{i}{j}, par_strs{i}{j}, ...
fmt_comment)]; %#ok<AGROW>
else
str = [str, sprintf([par_fmt_spec, ls],...
exp_par_names{i}{j}, par_strs{i}{j})]; %#ok<AGROW>
end
end
end
% Prints an extra line separator at the end
str = [str, sprintf(ls)];
end
% Save metadata to a file
function save(this, filename)
fileID = fopen(filename,'a');
fprintf(fileID, print(this));
fclose(fileID);
end
end
methods (Access = public, Static = true)
% Create metadata indicating the present moment of time
function TimeMdt = time(title)
if nargin()>0
assert(ischar(title)&&isvector(title),...
'Time field name must be a character vector')
else
title = 'Time';
end
TimeMdt = MyMetadata();
TimeMdt.title = title;
dv = datevec(datetime('now'));
addParam(TimeMdt, 'Year', dv(1), 'fmt_spec','%i');
addParam(TimeMdt, 'Month', dv(2), 'fmt_spec','%i');
addParam(TimeMdt, 'Day', dv(3), 'fmt_spec','%i');
addParam(TimeMdt, 'Hour', dv(4), 'fmt_spec','%i');
addParam(TimeMdt, 'Minute', dv(5), 'fmt_spec','%i');
addParam(TimeMdt, 'Second', floor(dv(6)), 'fmt_spec','%i');
addParam(TimeMdt, 'Millisecond',...
round(1000*(dv(6)-floor(dv(6)))),'fmt_spec','%i');
end
% Load metadata from file. Return all the entries found and
% the number of the last line read.
function [MdtArr, n_end_line] = load(filename, varargin)
fileID = fopen(filename,'r');
TmpMdt = MyMetadata(varargin{:});
title_exp = [TmpMdt.hdr_spec, '(\w.*)', TmpMdt.hdr_spec];
%Loop initialization
line_no=0;
curr_title='';
%Loop continues until we reach the next header or we reach
%the end of the file
while ~feof(fileID)
line_no=line_no+1;
% Grabs the current line
curr_line = fgetl(fileID);
% Gives an error if the file is empty, i.e. if fgetl
% returns -1
if curr_line == -1
error('Tried to read empty file. Aborting.')
end
% Skips if current line is empty
if isempty(curr_line)
continue
end
title_token = regexp(curr_line, title_exp,'once','tokens');
% If we find a title, first check if it is the specified
% end header. Then change the title if a title was found,
% then if no title was found, put the data under the current
% current title.
if ismember(this.end_header, title_token)
break
elseif ~isempty(title_token)
% Apply genvarname for sefety in case the title string
% is not a proper variable name
curr_title=genvarname(title_token{1});
addField(this, curr_title);
%This runs if there was no match for the header regular
%expression, i.e. the current line is not a filed
%separator, and the current line is not empty. We then
%add this line to the current field (curr_title), possibly
%iterating over the parameter subscripts.
elseif ~isempty(curr_title)
% First separate the comment if present
tmp = regexp(curr_line, this.comment_sep, ...
'split', 'once');
if length(tmp)>1
comment_str = tmp{2}; % There is a comment
else
comment_str = ''; % There is no comment
end
% Then process name-value pair. Regard everything after
% the first column separator as value.
tmp = regexp(tmp{1}, this.column_sep, 'split', 'once');
if length(tmp)<2
% Ignore the line if a name-value pair is not found
continue
else
% Attempt convertion of value to number
val = str2doubleHedged(strtrim(tmp{2}));
end
% Infer the variable name and subscript reference
try
[S, name] = str2substruct(strtrim(tmp{1}));
catch
name = '';
end
if isempty(name)
% Ignore the line if variable name is not missing
continue
elseif ismember(name, fieldnames(this.(curr_title)))
% If the variable name already presents among
% parameters, add new subscript value
this.(curr_title).(name).value= ...
subsasgn(this.(curr_title).(name).value,S,val);
else
% Add new parameter with comment
addParam(this, curr_title, name, val,...
'SubStruct', S, 'comment', comment_str);
end
end
end
fclose(fileID);
if isempty(this.field_names)
warning('No metadata found, continuing without metadata.')
n_end_header=1;
else
n_end_header=line_no;
end
end
end
methods (Access = protected)
% Parse string and determine the type of string
function S = parseLine(this, str)
S = struct( ...
'type', '', ... % title, paramval, other
'match', {}); % parsed output
% Check if the line contains a parameter - value pair.
% First separate the comment if present
pv_token = regexp(str, this.comment_sep, 'split', 'once');
if length(pv_token)>1
comment = pv_token{2}; % There is a comment
else
comment = ''; % There is no comment
end
% Then process name-value pair. Regard everything after
% the first column separator as value.
pv_token = regexp(pv_token{1}, this.column_sep, 'split', ...
'once');
if length(pv_token)>=2 && isvarname(pv_token{1})
% A parameter-value pair is found
S.type = 'paramval';
S.match = {pv_token{1}, pv_token{2}, comment};
return
end
% Check if the line contains a title
title_exp = [this.hdr_spec, '(\w.*)', this.hdr_spec];
title_token = regexp(str, title_exp, 'once', 'tokens');
if ~isempty(title_token)
% Title expression found
S.type = 'title';
S.match = title_token{1};
return
end
% No match found
S.type = 'other';
S.match = {};
end
end
end
c4science · Help