[ < ] | [ > ] | [ << ] | [ Up ] | [ >> ] | [Top] | [Contents] | [Index] | [ ? ] |
This program reads an annotation file, determines the intervals between beat annotations (assumed to be the R-R intervals), and accumulates a histogram of them.
1 #include <stdio.h> 2 #include <wfdb/wfdb.h> 3 #include <wfdb/ecgmap.h> 4 5 main(argc, argv) 6 int argc; 7 char *argv[]; 8 { 9 int rr, *rrhist, rrmax; 10 long t; 11 WFDB_Anninfo a; 12 WFDB_Annotation annot; 13 void *calloc(); 14 15 if (argc < 3) { 16 fprintf(stderr, "usage: %s annotator record\n", argv[0]); 17 exit(1); 18 } 19 a.name = argv[1]; a.stat = WFDB_READ; 20 if (annopen(argv[2], &a, 1) < 0) exit(2); 21 if ((rrmax = (int)(3*sampfreq(argv[2]))) <= 0) exit(3); 22 if ((rrhist = (int *)calloc(rrmax+1, sizeof(int))) == NULL) { 23 fprintf(stderr, "%s: insufficient memory\n", argv[0]); 24 exit(4); 25 } 26 while (getann(0, &annot) == 0 && !isqrs(annot.anntyp)) 27 ; 28 t = annot.time; 29 while (getann(0, &annot) == 0) 30 if (isqrs(annot.anntyp)) { 31 if ((rr = annot.time - t) > rrmax) rr = rrmax; 32 rrhist[rr]++; 33 t = annot.time; 34 } 35 for (rr = 1; rr < rrmax; rr++) 36 printf("%4d %s\n", rrhist[rr], mstimstr((long)rr)); 37 printf("%4d %s (or longer)\n", rrhist[rr], mstimstr((long)rr)); 38 exit(0); 39 } |
(See http://physionet.org/physiotools/wfdb/examples/example4.c for a copy of this program.)
Notes:
Here we allocate storage for the histogram. The value returned by
sampfreq
, if positive, specifies the number of sample intervals
per second; we will allocate 3 seconds’ worth of bins, initialized to
zero. See K&R, page 167, for a description of
calloc
.
This code sets t
to the time of the first annotated beat in the
record.
Here we read the remainder of the annotations, skipping any non-beat
annotations. The difference between the values of annot.time
for
consecutive beat annotations defines an R-R interval (rr
). Each
possible value of rr
up to rrmax
is assigned a bin in
rrhist
. Intervals longer than 3 seconds (rrmax
) are
counted in the bin corresponding to rr
= rrmax
.
The histogram is printed as a two-column table, with the number of
intervals in the first column and the length of the interval (with
millisecond resolution) in the second column. (What happens if
rr
starts at 0 rather than 1 in line 35?)
[ < ] | [ > ] | [ << ] | [ Up ] | [ >> ] |
PhysioNet (wfdb@physionet.org)