chris/OSCAR-code
1
0
Fork 0
mirror of https://gitlab.com/pholy/OSCAR-code.git synced 2025-04-05 18:50:44 +00:00
Code Issues Packages Projects Releases Wiki Activity

Initial TidalVolume & MinuteVent calculations for PRS1 machines. (Need to reimport.) Also added overview counterparts

Browse Source
This commit is contained in:
Mark Watkins 2011-12-01 01:19:31 +10:00
parent 9ebcdf7b1a
commit 69dde95b81
5 changed files with 95 additions and 22 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

84
SleepLib/calcs.cpp
View File

@ -29,21 +29,31 @@ int CalcRespRate::calculate(Session *session)
if (!session->eventlist.contains(CPAP_FlowRate)) return 0; //need flow waveform
EventList *flow, *rr;
EventList *flow, *rr, *tv=NULL, *mv=NULL;
if (!session->eventlist.contains(CPAP_TidalVolume)) {
tv=new EventList(EVL_Event);
session->eventlist[CPAP_TidalVolume].push_back(tv);
}
if (!session->eventlist.contains(CPAP_MinuteVent)) {
mv=new EventList(EVL_Event);
session->eventlist[CPAP_MinuteVent].push_back(mv);
}
int cnt=0;
for (int ws=0; ws < session->eventlist[CPAP_FlowRate].size(); ws++) {
flow=session->eventlist[CPAP_FlowRate][ws];
if (flow->count() > 5) {
rr=new EventList(EVL_Event);
session->eventlist[CPAP_RespRate].push_back(rr);
cnt+=filterFlow(flow,rr,flow->rate());
cnt+=filterFlow(flow,rr,tv,mv,flow->rate());
}
}
return cnt;
}
int CalcRespRate::filterFlow(EventList *in, EventList *out, double rate)
int CalcRespRate::filterFlow(EventList *in, EventList *out, EventList *tv, EventList *mv, double rate)
{
int size=in->count();
EventDataType *stage1=new EventDataType [size];
EventDataType *stage2=new EventDataType [size];
@ -51,7 +61,7 @@ int CalcRespRate::filterFlow(EventList *in, EventList *out, double rate)
QVector<EventDataType> med;
med.reserve(8);
EventDataType r;
EventDataType r,tmp;
int cnt;
EventDataType c;
@ -115,9 +125,12 @@ int CalcRespRate::filterFlow(EventList *in, EventList *out, double rate)
qint64 time=in->first();
qint64 u1=0,u2=0,len,l1=0,l2=0;
int z1=0,z2=0;
EventDataType lastc=0,thresh=0;
QVector<int> breaths;
QVector<EventDataType> TV;
QVector<qint64> breaths_start;
int lasti=0;
for (i=0;i<size;i++) {
c=stage1[i];
@ -126,11 +139,18 @@ int CalcRespRate::filterFlow(EventList *in, EventList *out, double rate)
u2=u1;
u1=time;
if (u2>0) {
z2=i;
len=abs(u2-u1);
//if (len>1500) {
breaths_start.push_back(time);
breaths.push_back(len);
//}
if (tv) { // && z1>0) { // Tidal Volume Calculations
EventDataType t=0;
for (int g=z1;g<z2;g++) {
tmp=-stage1[g];
t+=tmp;
}
TV.push_back(t);
}
breaths_start.push_back(time);
breaths.push_back(len);
}
}
} else {
@ -138,11 +158,20 @@ int CalcRespRate::filterFlow(EventList *in, EventList *out, double rate)
l2=l1;
l1=time;
if (l2>0) {
z1=i;
len=abs(l2-l1);
//if (len>1500) {
// breaths2_start.push_back(time);
// breaths2.push_back(len);
//}
if (tv) {
// Average the other half of the breath to increase accuracy.
EventDataType t=0;
for (int g=z2;g<z1;g++) {
tmp=stage1[g];
t+=tmp;
}
int ts=TV.size()-2;
if (ts>=0) {
// TV[ts]=(TV[ts]+t)/2.0;
}
}
}
}
@ -162,11 +191,17 @@ int CalcRespRate::filterFlow(EventList *in, EventList *out, double rate)
QVector<int> breaths2;
QVector<qint64> breaths2_start;
int fir=0;
QVector<EventDataType> TV2;
QVector<qint64> TV2_start;
int fir=0,fir2=0;
EventDataType T,L;
do {
br=0;
bool first=true;
bool cont=false;
T=0;
L=0;
for (int i=fir;i<breaths.size();i++) {
t=breaths_start[i];
l=breaths[i];
@ -194,10 +229,15 @@ int CalcRespRate::filterFlow(EventList *in, EventList *out, double rate)
continue;
} else
br+=1.0;
T+=TV[i]/2.0;
L+=l/1000.0;
}
if (cont) continue;
breaths2.push_back(br);
breaths2_start.push_back(t1+window/2);
//TV2_start.push_back(t2);
TV2.push_back(T);
//out->AddEvent(t,br);
//stage2[z++]=br;
@ -205,16 +245,26 @@ int CalcRespRate::filterFlow(EventList *in, EventList *out, double rate)
t2+=window/2.0;
} while (t2<in->last());
for (int i=1;i<breaths2.size()-2;i++) {
for (int i=2;i<breaths2.size()-2;i++) {
t=breaths2_start[i];
med.clear();
for (int j=0;j<4;j++) {
med.push_back(breaths2[i+j-1]);
for (int j=0;j<5;j++) {
med.push_back(breaths2[i+j-2]);
}
qSort(med);
br=med[2];
out->AddEvent(t,br);
//t=TV2_start[i];
med.clear();
for (int j=0;j<5;j++) {
med.push_back(TV2[i+j-2]);
}
qSort(med);
tmp=med[3];
tv->AddEvent(t,tmp);
mv->AddEvent(t,(tmp*br)/1000.0);
}
delete [] stage2;

2
SleepLib/calcs.h
View File

@ -25,7 +25,7 @@ public:
CalcRespRate(ChannelID id=CPAP_RespRate);
virtual int calculate(Session *session);
protected:
int filterFlow(EventList *in, EventList *out,double rate);
int filterFlow(EventList *in, EventList *out,EventList *tv, EventList *mv, double rate);
};
class CalcAHIGraph:public Calculation

3
daily.cpp
View File

@ -201,7 +201,8 @@ Daily::Daily(QWidget *parent,gGraphView * shared, MainWindow *mw)
MP->AddLayer(AddCPAP(new gLineChart(CPAP_MaskPressure,Qt::blue,false)));
RR->AddLayer(AddCPAP(new gLineChart(CPAP_RespRate,Qt::darkMagenta,square)));
MV->AddLayer(AddCPAP(new gLineChart(CPAP_MinuteVent,Qt::darkCyan,square)));
TV->AddLayer(AddCPAP(new gLineChart(CPAP_TidalVolume,Qt::magenta,square)));
TV->AddLayer(AddCPAP(new gLineChart(CPAP_TidalVolume,Qt::yellow,square)));
TV->AddLayer(AddCPAP(new gLineChart("TidalVolume2",Qt::magenta,square)));
FLG->AddLayer(AddCPAP(new gLineChart(CPAP_FLG,Qt::darkBlue,true)));
//RE->AddLayer(AddCPAP(new gLineChart(CPAP_RespiratoryEvent,Qt::magenta,true)));
IE->AddLayer(AddCPAP(new gLineChart(CPAP_IE,Qt::darkRed,square)));

24
overview.cpp
View File

@ -98,6 +98,9 @@ Overview::Overview(QWidget *parent,gGraphView * shared) :
SES=createGraph("Sessions");
NPB=createGraph("% in PB");
RR=createGraph("Resp. Rate");
TV=createGraph("Tidal Volume");
MV=createGraph("Minute Vent.");
PTB=createGraph("Pat. Trig. Br.");
uc=new SummaryChart("Hours",GT_BAR);
uc->addSlice("",QColor("green"),ST_HOURS);
@ -128,12 +131,31 @@ Overview::Overview(QWidget *parent,gGraphView * shared) :
SET->setRecMaxY(5);
SET->AddLayer(set);
rr=new SummaryChart("bpm",GT_LINE);
rr=new SummaryChart("breaths/min",GT_LINE);
rr->addSlice(CPAP_RespRate,QColor("light blue"),ST_MIN);
rr->addSlice(CPAP_RespRate,QColor("light green"),ST_90P);
rr->addSlice(CPAP_RespRate,QColor("blue"),ST_WAVG);
RR->AddLayer(rr);
tv=new SummaryChart("L/b",GT_LINE);
tv->addSlice(CPAP_TidalVolume,QColor("light blue"),ST_MIN);
tv->addSlice(CPAP_TidalVolume,QColor("light green"),ST_90P);
tv->addSlice(CPAP_TidalVolume,QColor("blue"),ST_WAVG);
TV->AddLayer(tv);
mv=new SummaryChart("L/m",GT_LINE);
mv->addSlice(CPAP_MinuteVent,QColor("light blue"),ST_MIN);
mv->addSlice(CPAP_MinuteVent,QColor("light green"),ST_90P);
mv->addSlice(CPAP_MinuteVent,QColor("blue"),ST_WAVG);
MV->AddLayer(mv);
ptb=new SummaryChart("%PTB",GT_LINE);
ptb->addSlice(CPAP_PTB,QColor("yellow"),ST_MIN);
ptb->addSlice(CPAP_PTB,QColor("light gray"),ST_90P);
ptb->addSlice(CPAP_PTB,QColor("orange"),ST_WAVG);
PTB->AddLayer(ptb);
pr=new SummaryChart("cmH2O",GT_LINE);
//PR->setRecMinY(4.0);
//PR->setRecMaxY(12.0);

4
overview.h
View File

@ -36,8 +36,8 @@ public:
gGraph * createGraph(QString name);
gGraph *AHI,*UC, *US, *PR,*LK,*NPB,*SET,*SES,*RR;
SummaryChart *bc,*uc, *us, *pr,*lk,*npb,*set,*ses,*rr;
gGraph *AHI,*UC, *US, *PR,*LK,*NPB,*SET,*SES,*RR,*MV,*TV,*PTB;
SummaryChart *bc,*uc, *us, *pr,*lk,*npb,*set,*ses,*rr,*mv,*tv,*ptb;
QVector<SummaryChart *> OverviewCharts;
public slots: