#include <stdio.h>
#include <wfdb/wfdb.h>
#include <wfdb/ecgcodes.h>
main(int argc, char **argv)
{
static char *record, *p;
static double map[30], m, msum;
static int i, j, h, mapsig, minutes, n, nsig, seconds, spm, state, v;
static WFDB_Anninfo ai;
static WFDB_Annotation annot;
static WFDB_Sample *x;
static WFDB_Siginfo *si;
static WFDB_Time t, tt;
if (argc < 2) {
fprintf(stderr, "usage: %s RECORD\n", argv[0]);
exit(1);
}
record = argv[1];
if ((nsig = isigopen(record, NULL, 0)) < 1) exit(1);
si = (WFDB_Siginfo *)malloc(nsig * sizeof(WFDB_Siginfo));
nsig = isigopen(record, si, nsig);
x = (WFDB_Sample *)malloc(nsig * sizeof(WFDB_Sample));
if ((mapsig = findsig("ABP Mean")) < 0) {
fprintf(stderr, "%s: no ABP Mean signal for record %s\n",
argv[0], argv[1]);
exit(1);
}
ai.name = "ahe";
ai.stat = WFDB_WRITE;
if (annopen(record, &ai, 1) < 0) exit(1);
spm = strtim("1:0"); /* samples per minute */
for (p = timstr(0L); *p != ':'; p++)
;
minutes = atoi(p+1);
if (spm == 60) {
for (seconds = atoi(p+4); seconds < 60; seconds++, t++)
getvec(x);
}
while (1) { /* process 1 minute of MAP per iteration, until end of data */
annot.anntyp = NOTQRS;
if (++minutes % 60 == 59) annot.anntyp = ARFCT;
if (++i >= 30) i = 0;
if (map[i] > 60.0) n--; /* oldest MAP is normal */
else if (map[i] > 10.0) h--; /* oldest MAP is low */
else if (t > 30*spm) n--; /* assume invalid sample is normal */
for (j = v = 0, msum = 0; j < spm; j++, t++) {
if (getvec(x) < 0) break;
m = aduphys(mapsig, x[mapsig]);
if (m > 0.0) { msum += m; v++; }
}
if (t < 0) t = -t; /* absolute times returned as < 0 */
if (j < spm) break; /* stop at end of data, discard remaining <60 sec */
if (v > 0) map[i] = msum/v;
else map[i] = 0;
// printf("%6ld %g\n", t, map[i]);
if (map[i] > 60.0) n++; /* newest MAP is normal */
else if (map[i] > 10.0) {
h++; /* newest MAP is low */
annot.anntyp = SVPB;
}
else n++; /* assume invalid sample is normal */
switch (state) {
case 0: /* indeterminate state */
if (h >= 27) { /* hypotensive for >= 90% of last 30 minutes */
annot.anntyp = PVC;
for (tt = t - 30*spm; tt < t; tt += spm) {
annot.time = tt;
putann(0, &annot);
}
state = 1;
}
else if (n == 30) { /* normotensive for 100% of last 30 minutes */
state = 2;
}
break;
case 1: /* hypotensive episode in progress */
if (h < 27) /* episode has ended */
state = 0;
else
annot.anntyp = PVC;
break;
case 2: /* normotensive episode in progress */
if (n < 30) /* episode has ended */
state = 0;
break;
}
if (annot.anntyp) {
annot.time = t;
putann(0, &annot);
}
}
wfdbquiet();
wfdbquit();
}