Added Daily View AHI/hr Graph for all CPAP machines

Graphs/gLineChart.cpp

@@ -443,3 +443,135 @@ void gLineChart::paint(gGraph & w,int left, int top, int width, int height)
     }
 }
 
+
+AHIChart::AHIChart(const QColor col)
+:Layer("ZZZ"),m_color(col)
+{
+    m_miny=m_maxy=0;
+    addGLBuf(lines=new GLShortBuffer(100000,GL_LINES));
+    lines->setColor(col);
+    lines->setAntiAlias(true);
+    lines->setSize(1.5);
+}
+
+AHIChart::~AHIChart()
+{
+}
+
+void AHIChart::paint(gGraph & w,int left, int top, int width, int height)
+{
+    if (!m_visible)
+        return;
+
+    if (!m_day)
+        return;
+
+    // Draw bounding box
+    GLShortBuffer *outlines=w.lines();
+    QColor blk=Qt::black;
+    outlines->add(left, top, left, top+height, blk);
+    outlines->add(left, top+height, left+width,top+height, blk);
+    outlines->add(left+width,top+height, left+width, top, blk);
+    outlines->add(left+width, top, left,top, blk);
+    width--;
+    height-=2;
+
+    EventDataType miny,maxy;
+    double minx,maxx;
+    miny=w.min_y, maxy=w.max_y;
+
+    maxx=w.max_x, minx=w.min_x;
+
+    // hmmm.. subtract_offset..
+
+    w.roundY(miny,maxy);
+
+    double xx=maxx-minx;
+    double xmult=double(width)/xx;
+
+    EventDataType yy=maxy-miny;
+    EventDataType ymult=EventDataType(height-3)/yy;   // time to pixel conversion multiplier
+
+    bool first=false;
+    double px,py;
+    double lastpx,lastpy;
+    double top1=top+height;
+    for (int i=0;i<m_time.size();i++) {
+        qint64 ti=m_time[i];
+        EventDataType v=m_data[i];
+        if ((ti>=minx) && (ti<maxx)) {
+            px=left+(double(ti-minx)*xmult);
+            py=top1-(double(v-miny)*ymult);
+            if (!first) {
+                first=true;
+            } else {
+                lines->add(px,py,lastpx,lastpy,m_color);
+            }
+        }
+        lastpx=px;
+        lastpy=py;
+    }
+}
+
+void AHIChart::SetDay(Day *d)
+{
+    m_day=d;
+    m_data.clear();
+    m_time.clear();
+    m_maxy=0;
+    m_miny=0;
+
+    if (!d) return;
+    m_miny=9999;
+    QVector<Session *>::iterator s;
+    qint64 first=d->first();
+    qint64 last=d->last();
+    qint64 f;
+
+    qint64 winsize=30000; // 30 second windows
+    for (qint64 ti=first;ti<last;ti+=winsize) {
+        f=ti-3600000L;
+        //if (f<first) f=first;
+        EventList *el[4];
+        EventDataType ahi=0;
+        int cnt=0;
+
+        for (s=d->begin();s!=d->end();s++) {
+            Session *sess=*s;
+            if ((ti<sess->first()) || (f>sess->last())) continue;
+
+            if (sess->eventlist.contains(CPAP_Obstructive))
+                el[0]=sess->eventlist[CPAP_Obstructive][0];
+            else el[0]=NULL;
+            if (sess->eventlist.contains(CPAP_Apnea))
+                el[1]=sess->eventlist[CPAP_Apnea][0];
+            else el[1]=NULL;
+            if (sess->eventlist.contains(CPAP_Hypopnea))
+                el[2]=sess->eventlist[CPAP_Hypopnea][0];
+            else el[2]=NULL;
+            if (sess->eventlist.contains(CPAP_ClearAirway))
+                el[3]=sess->eventlist[CPAP_ClearAirway][0];
+            else el[3]=NULL;
+
+            qint64 t;
+            for (int i=0;i<4;i++) {
+                if (!el[i]) continue;
+                for (int j=0;j<el[i]->count();j++) {
+                    t=el[i]->time(j);
+                    if ((t>=f) && (t<=ti)) {
+                        cnt++;
+                    }
+                }
+            }
+        }
+        double g=double(ti-f)/3600000.0;
+        if (g>0) ahi=cnt/g;
+
+        if (ahi<m_miny) m_miny=ahi;
+        if (ahi>m_maxy) m_maxy=ahi;
+        m_time.append(ti);
+        m_data.append(ahi);
+    }
+    m_minx=first;
+    m_maxx=last;
+}

Graphs/gLineChart.h

@@ -13,6 +13,25 @@
 #include "gGraphView.h"
 //#include "graphlayer.h"
 
+class AHIChart:public Layer
+{
+public:
+    AHIChart(const QColor col=QColor("black"));
+    ~AHIChart();
+    virtual void paint(gGraph & w,int left, int top, int width, int height);
+    virtual void SetDay(Day *d);
+    virtual EventDataType Miny() { return m_miny; }
+    virtual EventDataType Maxy() { return m_maxy; }
+    virtual bool isEmpty() { return m_data.size()==0; }
+
+protected:
+    QVector<EventDataType> m_data;
+    QVector<quint64> m_time;
+    EventDataType m_miny,m_maxy;
+    QColor m_color;
+    GLShortBuffer * lines;
+};
+
 class gLineChart:public Layer
 {
     public:

daily.cpp

@@ -74,6 +74,7 @@ Daily::Daily(QWidget *parent,gGraphView * shared, MainWindow *mw)
 
     SF=new gGraph(GraphView,"Event Flags",default_height);
     FRW=new gGraph(GraphView,"Flow Rate",default_height);
+    AHI=new gGraph(GraphView,"AHI",default_height);
     MP=new gGraph(GraphView,"Mask Pressure",default_height);
     PRD=new gGraph(GraphView,"Pressure",default_height);
     LEAK=new gGraph(GraphView,"Leak",default_height);
@@ -152,7 +153,7 @@ Daily::Daily(QWidget *parent,gGraphView * shared, MainWindow *mw)
     FRW->AddLayer(AddCPAP(los));
 
 
-    gGraph *graphs[]={ PRD, LEAK, SNORE, PTB, MP, RR, MV, TV, FLG, IE, TI, TE, TgMV, SPO2, PLETHY, PULSE,INTPULSE, INTSPO2 };
+    gGraph *graphs[]={ PRD, LEAK, AHI, SNORE, PTB, MP, RR, MV, TV, FLG, IE, TI, TE, TgMV, SPO2, PLETHY, PULSE,INTPULSE, INTSPO2 };
     int ng=sizeof(graphs)/sizeof(gGraph*);
     for (int i=0;i<ng;i++){
         graphs[i]->AddLayer(new gXGrid());
@@ -177,6 +178,7 @@ Daily::Daily(QWidget *parent,gGraphView * shared, MainWindow *mw)
     PRD->AddLayer(AddCPAP(new gLineChart(CPAP_EPAP,Qt::blue,square)));
     PRD->AddLayer(AddCPAP(new gLineChart(CPAP_IPAP,Qt::red,square)));
 
+    AHI->AddLayer(AddCPAP(new AHIChart(Qt::darkYellow)));
     LEAK->AddLayer(AddCPAP(new gLineChart(CPAP_Leak,Qt::darkYellow,square)));
     LEAK->AddLayer(AddCPAP(new gLineChart(CPAP_MaxLeak,Qt::darkRed,square)));
     SNORE->AddLayer(AddCPAP(new gLineChart(CPAP_Snore,Qt::darkGray,true)));

daily.h

@@ -77,7 +77,7 @@ private:
 
     gGraph *PRD,*FRW,*GAHI,*TAP,*LEAK,*SF,*TAP_EAP,*TAP_IAP,*PULSE,*SPO2,
            *SNORE,*RR,*MP,*MV,*TV,*FLG,*PTB,*OF,*INTPULSE,*INTSPO2, *THPR,
-           *PLETHY,*TI,*TE, *RE, *IE, *TgMV;
+           *PLETHY,*TI,*TE, *RE, *IE, *TgMV, *AHI;
 
     QList<Layer *> OXIData;
     QList<Layer *> CPAPData;

docs/index.html

@@ -13,8 +13,11 @@ p,a,td,body { font-size: 14px }
 <p>This software is currently being designed to assist you in reviewing data for your CPAP Machine, Oximeter, and Sleep Stage monitors, as well as help you <i>track</i> general issues related to sleep health.</p>
 <p><i><b>This is a developer preview, features are missing and bugs will be plentyful.</b></i></p>
 <p>Currenly supports the following machines:</p>
+<b>CPAP</b>
 <li>Philips Respironics System One</li>
 <li>ResMed S9 models</li>
+<li>DeVilbiss Intellipap models (*new)</li>
+<b>Oximetry</b>
 <li>Contec CMS50 Oximeters (rather poorly still I'm afraid)</li>
 <p>I don't recommend using this built in "web browser" to do any major surfing in, it will work, but it's mainly meant as a help browser.
 (It doesn't support SSL encryption.)</p>