Open Access Dataset and Toolbox of High-Density Surface Electromyogram Recordings 1.0.0

function data=load_mvc(path,subject,session,sig_type)

if(strcmp(sig_type,'force'))
    M=5;
else
    M=256;
end

for i=1:10
    if(mod(i,2)==1)
        direction='flexion';
    else
        direction='extension';
    end
    
    Fname = [path,'/mvc_dataset/subject',subject,'_session',num2str(session),'/mvc_',sig_type,'_finger',num2str(ceil(i/2)),'_',direction,'.dat'];         % Name with path.
    fid = fopen(Fname,'r','n');            % Open for reading.
    if fid<0, error(['Failed to open: ' Fname]); end
    data_tmp = fread(fid, [M,inf], 'int16');  % Read.
    fclose(fid);
    data_tmp=data_tmp';
        
    Fname = [path,'/mvc_dataset/subject',subject,'_session',num2str(session),'/mvc_',sig_type,'_finger',num2str(ceil(i/2)),'_',direction,'.hea'];         % Name with path.
    head_info=textread(Fname,'%s');
    idx = find( strcmp( head_info , ['mvc_',sig_type,'_finger',num2str(ceil(i/2)),'_',direction,'.dat'] ));
    
    for u=1:length(idx)
        str_tmp=head_info(idx(u)+2);
        idx2=strfind( str_tmp , '(' );
        gain=str2num(str_tmp{1,1}(1:(idx2{1,1}-1)));
        idx3=strfind( str_tmp , ')' );
        baseline=str2num(str_tmp{1,1}((idx2{1,1}+1):(idx3{1,1}-1)));
        
        data_tmp(:,u)=(data_tmp(:,u)-baseline)/gain;
    end
    data{i,1}=data_tmp;
end