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

Use rolling average when computing TV from flow graph

Browse Source 
Also, start with a guessed-at reasonable value to try to
minimized spike artifacts resulting from poor calculation
of initial breath volume.
This commit is contained in:
Guy Scharf 2021-07-17 16:22:40 -07:00
parent 0e3f669210
commit 04d7945241
1 changed files with 8 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

10
oscar/SleepLib/calcs.cpp
View File

@ -480,7 +480,7 @@ void FlowParser::calc(bool calcResp, bool calcTv, bool calcTi, bool calcTe, bool
quint32 *tv_tptr = nullptr; quint32 *tv_tptr = nullptr;
EventStoreType *tv_dptr = nullptr; EventStoreType *tv_dptr = nullptr;
int tv_count = 0; int tv_count = 0;
double tvlast, tvlast2, tvlast3; double tvlast = 0, tvlast2 = 0, tvlast3 = 0;
if (calcTv) { if (calcTv) {
TV = m_session->AddEventList(CPAP_TidalVolume, EVL_Event); TV = m_session->AddEventList(CPAP_TidalVolume, EVL_Event);
@ -600,8 +600,14 @@ void FlowParser::calc(bool calcResp, bool calcTv, bool calcTi, bool calcTe, bool
//val2=x; //val2=x;
// Average TV over last three data points // Average TV over last three data points
if (tv_count == 0) if (tv_count == 0) {
if (tv > 800) // Very much a Q&D patch to handle the first TV value not being calculated well
tv = 300; // If unreasonable, just set to a "reasonable" number.
tvlast = tvlast2 = tvlast3 = tv; tvlast = tvlast2 = tvlast3 = tv;
}
if (tv_count < 4) {
qDebug() << "tv" << tv << tvlast << tvlast2 << tvlast3 << "avg" << (tvlast + tvlast2 + tvlast3 + tv*2)/5;
}
tv = (tvlast + tvlast2 + tvlast3 + tv*2)/5; tv = (tvlast + tvlast2 + tvlast3 + tv*2)/5;
tvlast3 = tvlast2; tvlast3 = tvlast2;
tvlast2 = tvlast; tvlast2 = tvlast;