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

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

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

for i=1:5
    Fname = [path,'/random_dataset/subject',subject,'_session',num2str(session),'/random_',sig_type,'_sample',num2str(i),'.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,'/random_dataset/subject',subject,'_session',num2str(session),'/random_',sig_type,'_sample',num2str(i),'.hea'];         % Name with path.
    head_info=textread(Fname,'%s');
    idx = find( strcmp( head_info , ['random_',sig_type,'_sample',num2str(i),'.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{1,i}=data_tmp;
end