wire2data.m

Go to the documentation of this file.

%>@ingroup ioio
%>@file
%> @brief Merges several "Wire txt" files into a dataset
%>
%> @c data.groupcodes is made from file name. File name is trimmed at the trimdot-th dot counted from right to left. E.g., 
%> This allows one to trim "sample.0.dat" at the penultimate dot (ignodedoctount=2) or "sample.dat" at the last dot
%> (trimdot=1)
%>
%> data.obsnames will contain the file names
%>
%> For reference on parameters, please check @ref mergetool.m
%>
%> @sa mergetool.m

%> @param wild
%> @param trimdot See pirtool
%> @param flag_image
%> @param height
%> @return A dataset
function [data, flag_error] = wire2data(wild, trimdot, flag_image, height)

if ~exist('trimdot', 'var')
    trimdot = 2;
end;

[filenames, groupcodes] = resolve_dir(wild, trimdot, flag_image);

no_files = length(filenames);

if flag_image
    if no_files/height ~= floor(no_files/height)
        irerror('Wire2Data: Invalid image height!');
    end;
end;


idxs = find(wild == '/');
if isempty(idxs)
    path_ = '';
else
    path_ = wild(1:idxs(end));
end;
data = irdata();

flag_first = 1;
ipro = progress2_open('WIRE2DATA', [], 0, no_files);
cnt_error = 0;
errors = {};
ii = 0;
for i = 1:no_files
    filename = fullfile(path_, filenames{i});
    try
        M = dlmread(filename);

        if flag_first
            wns = M(:, 1)';
            nf = numel(wns);
        else
            % Checks against previous size
            n1 = size(M, 1);
            if ~(n1 == nf)
                irerror(sprrintf('Wire2Data: Wrong number of data points: expected=%d; found=%d', nf, n1));
            end;
        end;
        
        x = M(:, 2)';
        
        if flag_first
            % Initializes dataset if first row
            data.X = zeros(no_files, nf);
            data.classes = zeros(no_files, 1);
            data.obsnames = cell(no_files, 1);
            data.groupcodes = cell(no_files, 1);
            data.fea_x = wns;
            data.classlabels = {'Class 0'};
            flag_first = 0;
        end;

        ii = ii+1;
        data.X(ii, :) = x;
        data.classes(ii) = 0;
        data.obsnames{ii} = filenames{i};
        data.groupcodes{ii} = groupcodes{i};
    catch ME
        irverbose(['ERROR: ', ME.message]);
        irverbose(ME.getReport());
        cnt_error = cnt_error+1;
        errors{end+1} = filename;
    end;
    ipro = progress2_change(ipro, [], [], i);
end;
progress2_close(ipro);

flag_error = 0;
if cnt_error > 0
    irverbose(sprintf('NOTICE: only %d/%d files were successfully read. Import failed on following files:', no_files-cnt_error, no_files));
    for i = 1:cnt_error
        irverbose(errors{i});
    end;
    flag_error = cnt_error;
    
    % Trims matrices rows because were overdimensioned
    data.X = data.X(1:ii, :);
    data.classes = data.classes(1:ii, :);
    data.obsnames = data.obsnames(1:ii, :);
    data.groupcodes = data.groupcodes(1:ii, :);
end;


if flag_image
    data.height = height;
end;