-
Notifications
You must be signed in to change notification settings - Fork 6
/
load_claris.m
132 lines (102 loc) · 3.9 KB
/
load_claris.m
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% BRAVEHEART - Open source software for electrocardiographic and vectorcardiographic analysis
% load_claris.m -- Load Abbott WorkMate Claris ASCII format ECGs
% Copyright 2016-2024 Hans F. Stabenau and Jonathan W. Waks
%
% Source code/executables: https://github.com/BIVectors/BRAVEHEART
% Contact: [email protected]
%
% BRAVEHEART is free software: you can redistribute it and/or modify it under the terms of the GNU
% General Public License as published by the Free Software Foundation, either version 3 of the License,
% or (at your option) any later version.
%
% BRAVEHEART is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY;
% without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.
% See the GNU General Public License for more details.
%
% You should have received a copy of the GNU General Public License along with this program.
% If not, see <https://www.gnu.org/licenses/>.
%
% This software is for research purposes only and is not intended to diagnose or treat any disease.
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function [hz, I, II, III, avR, avF, avL, V1, V2, V3, V4, V5, V6] = load_claris(filename)
% Parse filename to get structure needed to find all of the corresponding
% text files since the data is stored in 1 lead per file
% It does not matter which lead is chosen - this function will parse all
% leads with a smilar structure:
% e.g. II.Session 18 - Page 1.2.TXT
% Will load all 12 <L>.Session 18 - Page 1.<N>.TXT where <L> is the lead and
% <N> is the order in which the signals appear in the Claris page which has
% no bearing on loading the signals.
% This load module will ONLY work if the file names are in the above format.
[d,f,~] = fileparts(filename);
common_file = strsplit(f,'.');
common_file = common_file{2};
% common_file now contains the structure of the data:
%<L>.common_file.<N>.txt
L =struct();
leads = fieldnames(ECG12);
for i = 3:14 % ECG12 fieldnames
for N = 1:40 % Max number of leads from Claris page (may need to change)
F = strcat(leads{i},'.',common_file,'.',num2str(N),'.txt');
if isfile(fullfile(d,F))
L.(leads{i}) = parse_claris_txt(fullfile(d,F));
break;
end
end
end
% Need to have at least 8 leads
assert(numel(fieldnames(L)) >= 8, "Not enough files to make a 12-lead ECG");
% See if an information file exists
info_file_name = strsplit(common_file,' - ');
info_file_name = info_file_name{1};
info_file_name = strcat(info_file_name,' Information.TXT');
if isfile(fullfile(d,info_file_name))
info_file = textread(fullfile(d,info_file_name),'%s','delimiter','\n');
freq_idx = find(contains(info_file,'Sample Rate='));
freq_cell = strsplit(info_file{freq_idx},'=');
hz = str2double(extract(string(freq_cell{2}), digitsPattern));
gain_idx = find(contains(info_file,'Channel Resolution='));
gain_cell = strsplit(info_file{gain_idx},'=');
gain_val = str2double(extract(string(gain_cell{2}), digitsPattern));
gain_units = strsplit(gain_cell{2},' ');
gain_units = gain_units{2};
switch gain_units
case 'nV/LSB'
gain = gain_val / 1000000;
otherwise
error('gain units not in nV')
end
else
gain = 78/1000000;
hz = 2000;
end
% Form row vectors gained properly
I = L.I' * gain;
II = L.II' * gain;
if isfield(L,'III')
III = L.III' * gain;
else
III = -I + II;
end
if isfield(L,'avR')
avR = L.avR' * gain;
else
avR = -0.5*I - 0.5*II;
end
if isfield(L,'avF')
avF = L.avF' * gain;
else
avF = II - 0.5*I;
end
if isfield(L,'avL')
avL = L.avL' * gain;
else
avL = I - 0.5*II;
end
V1 = L.V1' * gain;
V2 = L.V2' * gain;
V3 = L.V3' * gain;
V4 = L.V4' * gain;
V5 = L.V5' * gain;
V6 = L.V6' * gain;