Changeset 7917 in ntrip

Timestamp:

May 23, 2016, 6:18:50 PM (8 years ago)

Author:

stuerze

Message:

rinex qc: satellite visibility is considered now if expected obs is computed

Location:

trunk/BNC/src/rinex

Files:

• reqcanalyze.cpp (modified) (8 diffs)
• reqcanalyze.h (modified) (2 diffs)

trunk/BNC/src/rinex/reqcanalyze.cpp

@@ -217 +217 @@
       for (unsigned iObs = 0; iObs < _currEpo->rnxSat.size(); iObs++) {
         const t_rnxObsFile::t_rnxSat& rnxSat = _currEpo->rnxSat[iObs];
+        if (_navFileNames.size() &&
+            _numExpObs.find(rnxSat.prn) == _numExpObs.end()) {
+          _numExpObs[rnxSat.prn] = 0;
+        }
         if (_signalTypes.find(rnxSat.prn.system()) == _signalTypes.end()) {
           continue;
@@ -230 +234 @@
 
     analyzeMultipath();
+
+    if (_navFileNames.size()) {
+      setExpectedObs(_qcFile._startTime, _qcFile._endTime, _qcFile._interval, xyzSta);
+    }
 
     preparePlotData(obsFile);
@@ -812 +820 @@
         << "Interval           : " << _qcFile._interval                            << endl;
 
-  int numPossibleObs = int((_qcFile._endTime - _qcFile._startTime + _qcFile._interval)/_qcFile._interval);
-
   // Number of systems
   // -----------------
@@ -847 +853 @@
   }
 
+  // System specific summary
+  // -----------------------
   itSys.toFront();
   while (itSys.hasNext()) {
@@ -852 +860 @@
     const QChar&                      sys      = itSys.key();
     const QVector<const t_qcSatSum*>& qcSatVec = itSys.value();
+    int numExpectedObs = 0;
+    for(QMap<t_prn, int>::iterator it = _numExpObs.begin();
+        it != _numExpObs.end(); it++) {
+      if (sys == it.key().system()) {
+        numExpectedObs += it.value();
+      }
+    }
     QString prefixSys = QString("  ") + sys + QString(": ");
     QMap<QString, QVector<const t_qcFrqSum*> > frqMap;
@@ -881 +896 @@
       const QVector<const t_qcFrqSum*> qcFrqVec = itFrq.value();
       QString prefixFrq = QString("  ") + frqType + QString(": ");
-
       int    numObs          = 0;
       int    numSlipsFlagged = 0;
@@ -907 +921 @@
         sumMP /= numMP;
       }
-      double ratio = double(numObs) / ((double(numPossibleObs) * double(qcSatVec.size())));
+
+      double ratio = (double(numObs) / double(numExpectedObs)) * 100.0;
+
       *_log << endl
-            << prefixSys2 << prefixFrq << "Observations      : "
-            << QString("%1 (%2) %3 \%\n").arg(numObs,           6).arg(numPossibleObs*qcSatVec.size(),           8).arg(ratio*100.0, 8, 'f', 2)
-            << prefixSys2 << prefixFrq << "Slips (file+found): " << QString("%1 +").arg(numSlipsFlagged,  8)
+            << prefixSys2 << prefixFrq << "Observations      : ";
+      if(_navFileNames.isEmpty() || _navFileIncomplete.contains(sys.toLatin1())) {
+        *_log << QString("%1\n").arg(numObs,           6);
+      }
+      else {
+        *_log << QString("%1 (%2) %3 \%\n").arg(numObs,           6).arg(numExpectedObs,           8).arg(ratio, 8, 'f', 2);
+      }
+      *_log << prefixSys2 << prefixFrq << "Slips (file+found): " << QString("%1 +").arg(numSlipsFlagged,  8)
                                                                  << QString("%1\n").arg(numSlipsFound,    8)
             << prefixSys2 << prefixFrq << "Gaps              : " << QString("%1\n").arg(numGaps,          8)
@@ -1060 +1081 @@
   }
 }
+
+void t_reqcAnalyze::setExpectedObs(const bncTime& startTime, const bncTime& endTime,
+                                   double interval, const ColumnVector& xyzSta) {
+
+  for(QMap<t_prn, int>::iterator it = _numExpObs.begin();
+      it != _numExpObs.end(); it++) {
+    t_eph* eph = 0;
+    for (int ie = 0; ie < _ephs.size(); ie++) {
+      if (_ephs[ie]->prn().system() == it.key().system() &&
+          _ephs[ie]->prn().number() == it.key().number()) {
+        eph = _ephs[ie];
+        break;
+      }
+    }
+    if (eph) {
+      int numExpObs = 0;
+      bncTime epoTime;
+      for (epoTime = startTime - interval; epoTime < endTime;
+           epoTime = epoTime + interval) {
+        ColumnVector xc(4);
+        ColumnVector vv(3);
+        if ( xyzSta.size() == 3 && (xyzSta[0] != 0.0 || xyzSta[1] != 0.0 || xyzSta[2] != 0.0) &&
+             eph->getCrd(epoTime, xc, vv, false) == success) {
+          double rho, eleSat, azSat;
+          topos(xyzSta(1), xyzSta(2), xyzSta(3), xc(1), xc(2), xc(3), rho, eleSat, azSat);
+          if ((eleSat * 180.0/M_PI) > 0.0) {
+            numExpObs++;
+          }
+        }
+      }
+      it.value() = numExpObs;
+    }
+    else {
+      if (!_navFileIncomplete.contains(it.key().system())) {
+        qDebug() <<  it.key().system() << it.key().number();
+        _navFileIncomplete.append(it.key().system());
+      }
+    }
+  }
+}

trunk/BNC/src/rinex/reqcanalyze.h

@@ -171 +171 @@
   void   analyzeFile(t_rnxObsFile* obsFile);
 
-  void   updateQcSat(const t_qcSat& qcSat, t_qcSatSum& qcSatSum);
-
-  void   setQcObs(const bncTime& epoTime, const ColumnVector& xyzSta,
-                  const t_satObs& satObs, QMap<QString, bncTime>& lastObsTime, t_qcSat& qcSat);
+  void  updateQcSat(const t_qcSat& qcSat, t_qcSatSum& qcSatSum);
+
+  void  setQcObs(const bncTime& epoTime, const ColumnVector& xyzSta,
+                 const t_satObs& satObs, QMap<QString, bncTime>& lastObsTime, t_qcSat& qcSat);
+
+  void   setExpectedObs(const bncTime& startTime, const bncTime& endTime,
+                        double interval, const ColumnVector& xyzSta);
 
   void   analyzeMultipath();
@@ -192 +195 @@
   QString                       _reqcPlotSignals;
   QMap<char, QVector<QString> > _signalTypes;
+  QMap<t_prn, int>              _numExpObs;
+  QVector<char>                 _navFileIncomplete;
   QStringList                   _defaultSignalTypes;
   QVector<t_eph*>               _ephs;