#include <stdio.h>
#include <wfdb/wfdb.h>
tprint(double tin)
{
static char tstring[32];
WFDB_Time t = tin;
printf("\t%s", timstr(-t)); /* time and date */
printf(" %s", timstr(t)); /* elapsed time since start of record */
}
main(int argc, char **argv)
{
char buf[40], tstring[20], *p = buf;
static double map[30], m, msum;
double t;
static int i, j, h, n, spm, state, v;
spm = 60 * sampfreq(argv[1]); /* samples per minute */
if (spm <= 0.) exit(1);
printf("Episode\t.-----------Start------------.\t"
".----------- End ------------.\n");
printf("type \t time date elapsed\t"
" time date elapsed\n");
while (1) { /* process 1 minute of MAP per iteration, until end of data */
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 && fgets(buf, sizeof(buf), stdin); j++) {
for (p = buf; *p; p++)
if (*p == '\t') {
m = atof(p+1);
break;
}
if (m > 0.0) { msum += m; v++; }
}
*p = '\0';
t = strtim(buf);
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;
if (map[i] > 60.0) n++; /* newest MAP is normal */
else if (map[i] > 10.0) h++; /* newest MAP is low */
else n++; /* assume invalid sample is normal */
// printf("state = %d, t = %s, map[%d] = %g, h = %d, n = %d, v = %d\n",
// state, timstr(-t), i, map[i], h, n, v);
switch (state) {
case 0: /* indeterminate state */
if (h >= 27) { /* hypotensive for >= 90% of last 30 minutes */
printf("AHE");
tprint(t > 30*spm ? t-30*spm : 1.0);
state = 1;
}
else if (n == 30) { /* normotensive for 100% of last 30 minutes */
printf(">60mmHg");
tprint(t > 30*spm ? t-30*spm : 1.0);
state = 2;
}
break;
case 1: /* hypotensive episode in progress */
// printf("t = %s, h = %d\n", timstr(-t), h);
if (h < 27) { /* episode has ended */
tprint(t-1); /* AHE definition last satisfied at t-1 */
printf("\n");
state = 0;
}
break;
case 2: /* normotensive episode in progress */
if (n < 30) { /* episode has ended */
tprint(t-2); /* first low MAP was at t-1, so last known
non-low was at t-2 */
printf("\n");
state = 0;
}
break;
}
}
if (state == 0)
printf("End");
tprint(t);
printf("\n");
}