function reportECG(ECG_w, detailLevel, report_mode, win_lengths, report_format, filename)
% This function creates a report of the signals handled by the ECG wrapper
% object ECG_w. The report includes several views of the signals, from a
% wide to a narrow scale. Some aspects of the report can be configured as
% the detail degree, the length of each time scale and the report format.
%
% Arguments:
%
% +ECG_w: [ECGwrapper object] REQUIRED. The signal handler.
%
% +detailLevel: [char] OPTIONAL. The report detail level:
% 'HighDetail', 'MediumDetail', 'LowDetail'. A higher
% detail level means report the whole recording at every
% time resolution defined in "win_lengths". High
% resolution also means larger reports. Default 'LowDetail'.
%
% +report_mode: [char] OPTIONAL. Information from other tasks like QRS
% detection/delineation/classification added to the
% signals in case available mode. Possible values are:
% 'full' 'ECG only' 'QRS detection' 'Wave delineation'
% 'Heartbeat classification'
% Default "ECG only".
%
% +win_lengths: [numeric] OPTIONAL. The amount and size (in seconds) of
% each scale length present in the report. Default
% [60*60 30*60 60 7] it means 1 hour - 30 min - 1 min
% and 7 seconds.
%
% +report_format: [char] OPTIONAL. The report format of the document.
% Default PDF.
%
% +filename: [char] OPTIONAL. The report filename. Default same folder
% of the original recording handled by ECG_w.
%
%
% Author: Mariano Llamedo Soria (llamedom at {electron.frba.utn.edu.ar; unizar.es}
% Version: 0.1 beta
% Last update: 17/6/2014
% Birthdate : 15/8/2012
%% constants
cKnownModes = {'full' 'ECG only' 'QRS detection' 'Wave delineation' 'Heartbeat classification' };
cKnownFormats = {'pdf' 'eps' 'png' 'tiff' 'jpg' 'bmp' };
cDetails = {'HighDetail' 'MediumDetail' 'LowDetail'};
%% argument parsing
if( nargin < 5 || ~any(strcmpi(cKnownFormats, report_format)) )
report_format = 'pdf';
end
heasig = ECG_w.ECG_header;
if( nargin < 3 || isempty(intersect(cKnownModes, report_mode)) )
report_mode = 'ECG only';
end
file_suffix = rowvec(colvec(char(report_mode)'));
prev_length = length(file_suffix);
file_suffix = strrep(file_suffix, ' ', ' ');
while( length(file_suffix) ~= prev_length )
prev_length = length(file_suffix);
file_suffix = strrep(file_suffix, ' ', ' ');
end
file_suffix = ['_' file_suffix];
if( any(strcmpi(report_mode, 'full')) )
report_mode = cKnownModes(2:end);
end
if( nargin < 6 || isempty(filename) )
[~, file_recname ]= fileparts(ECG_w.recording_name);
report_path = ECG_w.output_path;
if( isfield(heasig, 'recname') )
filename = [report_path heasig.recname file_suffix '.' report_format];
else
filename = [report_path file_recname file_suffix '.' report_format];
end
end
if( nargin < 2 || ~any(strcmpi(cDetails, detailLevel)) )
detailLevel = 'LowDetail'; % seconds
end
if( nargin < 4 || isempty(win_lengths) || ~isnumeric(win_lengths) )
win_lengths = [60*60 30*60 60 7]; % seconds
end
if( ~isobject(ECG_w) || ~isa(ECG_w, 'ECGwrapper') )
error('reportECG:BadArg', 'Need an ECGwrapper object to work.')
end
init_ghostscript();
win_lengths = bsxfun( @min, win_lengths, floor(heasig.nsamp/heasig.freq) );
win_lengths = sort(unique(win_lengths), 'descend');
% Activate the progress_struct bar.
pb = progress_bar('Report ECG');
pb.Loops2Do = length(win_lengths);
ii = 1;
fig_hdl = [];
for win_length = win_lengths
pb.start_loop();
if( (win_length * heasig.freq) > heasig.nsamp )
win_length = round(heasig.nsamp / heasig.freq);
end
switch detailLevel
case 'LowDetail'
max_cant_figs = 10;
case 'MediumDetail'
max_cant_figs = 30;
case 'HighDetail'
max_cant_figs = realmax;
end
cant_figs = min(max_cant_figs, ceil(heasig.nsamp / win_length / heasig.freq));
starts = linspace(0, round(heasig.nsamp / heasig.freq) - win_length, cant_figs);
str_aux = ['Reporting for window size: ' Seconds2HMS(win_length) ];
pb1 = pb.AddInerLoop( str_aux );
pb1.Loops2Do = cant_figs;
jj = 1;
for start_time = starts
pb1.start_loop();
fig_hdl = prepare_fig_hdl( fig_hdl );
eval_str = 'plot_ecg_strip(ECG_w, ''PrettyPrint'', true, ''Figure_handle'', fig_hdl, ''Start_time'', start_time, ''End_time'', start_time + win_length ';
if( any(strcmpi(report_mode, 'QRS detection')) )
eval_str = [ eval_str ',''QRS_locations'', ECG_w' ];
end
if( any(strcmpi(report_mode, 'Wave delineation')) )
eval_str = [ eval_str ',''ECG_delineation_single_lead'', ECG_w, ''ECG_delineation_multilead'', ECG_w' ];
end
if( any(strcmpi(report_mode, 'Heartbeat classification')) )
eval_str = [ eval_str ',''Heartbeat_classification'', ECG_w' ];
end
eval_str = [ eval_str ');'];
eval(eval_str);
set(fig_hdl, 'Visible', 'off');
str_aux2 = [ str_aux ' - Start: ' Seconds2HMS(start_time) ];
pb1.checkpoint(str_aux2);
switch( report_format )
case {'pdf' 'tiff' }
if( jj == 1 )
% export_fig(filename, '-nocrop', '-bookmark', Seconds2HMS(win_length), '-append', ['-' report_format], fig_hdl);
export_fig(filename, '-nocrop', '-append', ['-' report_format], fig_hdl);
else
% export_fig(filename, '-nocrop', '-bookmark', [' ' Seconds2HMS(start_time)], '-append', ['-' report_format], fig_hdl);
export_fig(filename, '-nocrop', '-append', ['-' report_format], fig_hdl);
end
otherwise
[aux_path, aux_filename] = fileparts(filename);
aux_filename = [aux_path filesep aux_filename];
aux_filename = sprintf( '%s_%04d.%s', aux_filename, ii, report_format);
export_fig(aux_filename, '-nocrop', ['-' report_format], fig_hdl);
end
pb1.checkpoint(str_aux2);
ii = ii + 1;
jj = jj + 1;
pb1.end_loop();
end
pb.end_loop();
end
clear pb
function fig_hdl = prepare_fig_hdl( fig_hdl )
paper_size = [0 0 29.7 21.0 ];
if( ishandle(fig_hdl) )
close(fig_hdl)
end
fig_hdl = figure('Visible', 'off');
set(fig_hdl, 'NumberTitle', 'off');
set(fig_hdl, 'PaperUnits', 'centimeters');
set(fig_hdl, 'PaperType', 'A4');
set(fig_hdl, 'PaperOrientation', 'landscape');
set(fig_hdl, 'PaperPosition', paper_size);
set(fig_hdl, 'Units', 'centimeters');
set(fig_hdl, 'Position', paper_size);