Changeset 7546 in ntrip for trunk/BNC/src/rinex

Timestamp:

Oct 22, 2015, 1:00:03 PM (9 years ago)

Author:

stuerze

Message:

observation types per system are added in header of the QC output

File:

• trunk/BNC/src/rinex/reqcanalyze.cpp (modified) (13 diffs)

trunk/BNC/src/rinex/reqcanalyze.cpp

@@ -78 +78 @@
   connect(this, SIGNAL(dspSkyPlot(const QString&, const QString&, QVector<t_polarPoint*>*,
                                   const QString&, QVector<t_polarPoint*>*,
-                                  const QByteArray&, double)), 
+                                  const QByteArray&, double)),
           this, SLOT(slotDspSkyPlot(const QString&, const QString&, QVector<t_polarPoint*>*,
                                     const QString&, QVector<t_polarPoint*>*,
@@ -196 +196 @@
   // --------------------
   try {
-    QMap<QString, bncTime> lastObsTime; 
+    QMap<QString, bncTime> lastObsTime;
     bool firstEpo = true;
     while ( (_currEpo = obsFile->nextEpoch()) != 0) {
@@ -329 +329 @@
 //
 ////////////////////////////////////////////////////////////////////////////
-void t_reqcAnalyze::setQcObs(const bncTime& epoTime, const ColumnVector& xyzSta, 
+void t_reqcAnalyze::setQcObs(const bncTime& epoTime, const ColumnVector& xyzSta,
                              const t_satObs& satObs, QMap<QString, bncTime>& lastObsTime,
                              t_qcSat& qcSat) {
@@ -458 +458 @@
       // Loop over all Chunks of Data
       // ----------------------------
-      for (bncTime chunkStart = _qcFile._startTime; 
+      for (bncTime chunkStart = _qcFile._startTime;
            chunkStart < _qcFile._endTime; chunkStart += chunkStep) {
 
@@ -465 +465 @@
         QVector<t_qcFrq*> frqVec;
         QVector<double>   MP;
-    
+
         // Loop over all Epochs within one Chunk of Data
         // ---------------------------------------------
@@ -494 +494 @@
           }
           meanMP /= MP.size();
-        
+
           bool slipMP = false;
-        
+
           double stdMP = 0.0;
           for (int ii = 0; ii < MP.size(); ii++) {
@@ -506 +506 @@
             stdMP += diff * diff;
           }
-        
+
           if (slipMP) {
             stdMP = 0.0;
@@ -517 +517 @@
             qcFrqSum._sumMP += stdMP;
           }
-        
+
           for (int ii = 0; ii < frqVec.size(); ii++) {
             t_qcFrq* qcFrq = frqVec[ii];
@@ -659 +659 @@
     QVector<QWidget*> plots;
     if (data1) {
-      QwtText title(title1); 
+      QwtText title(title1);
       QFont font = title.font(); font.setPointSize(font.pointSize()-1); title.setFont(font);
       t_polarPlot* plot1 = new t_polarPlot(title, scaleInterval, BNC_CORE->mainWindow());
@@ -793 +793 @@
   // Summary
   // -------
-  *_log << "Observation File  : " << obsFileName                                   << endl
-        << "RINEX Version     : " << QString("%1").arg(obsFile->version(),4,'f',2) << endl
-        << "Marker Name       : " << _qcFile._markerName                           << endl
-        << "Marker Number     : " << obsFile->markerNumber()                       << endl
-        << "Receiver          : " << _qcFile._receiverType                         << endl
-        << "Antenna           : " << _qcFile._antennaName                          << endl
-        << "Position XYZ      : " << QString("%1 %2 %3").arg(obsFile->xyz()(1), 14, 'f', 4)
+  *_log << "Observation File   : " << obsFileName                                   << endl
+        << "RINEX Version      : " << QString("%1").arg(obsFile->version(),4,'f',2) << endl
+        << "Marker Name        : " << _qcFile._markerName                           << endl
+        << "Marker Number      : " << obsFile->markerNumber()                       << endl
+        << "Receiver           : " << _qcFile._receiverType                         << endl
+        << "Antenna            : " << _qcFile._antennaName                          << endl
+        << "Position XYZ       : " << QString("%1 %2 %3").arg(obsFile->xyz()(1), 14, 'f', 4)
                                                         .arg(obsFile->xyz()(2), 14, 'f', 4)
                                                         .arg(obsFile->xyz()(3), 14, 'f', 4) << endl
-        << "Antenna dH/dE/dN  : " << QString("%1 %2 %3").arg(obsFile->antNEU()(3), 8, 'f', 4)
+        << "Antenna dH/dE/dN   : " << QString("%1 %2 %3").arg(obsFile->antNEU()(3), 8, 'f', 4)
                                                         .arg(obsFile->antNEU()(2), 8, 'f', 4)
                                                         .arg(obsFile->antNEU()(1), 8, 'f', 4) << endl
-        << "Start Time        : " << _qcFile._startTime.datestr().c_str()         << ' '
+        << "Start Time         : " << _qcFile._startTime.datestr().c_str()         << ' '
                                   << _qcFile._startTime.timestr(1,'.').c_str()    << endl
-        << "End Time          : " << _qcFile._endTime.datestr().c_str()           << ' '
+        << "End Time           : " << _qcFile._endTime.datestr().c_str()           << ' '
                                   << _qcFile._endTime.timestr(1,'.').c_str()      << endl
-        << "Interval          : " << _qcFile._interval                            << endl;
+        << "Interval           : " << _qcFile._interval                            << endl;
 
   // Number of systems
   // -----------------
-  QMap<QChar, QVector<const t_qcSatSum*> > systemMap; 
+  QMap<QChar, QVector<const t_qcSatSum*> > systemMap;
   QMapIterator<t_prn, t_qcSatSum> itSat(_qcFile._qcSatSum);
   while (itSat.hasNext()) {
@@ -821 +821 @@
     systemMap[prn.system()].push_back(&qcSatSum);
   }
-  *_log << "Navigation Systems: " << systemMap.size() << "   ";
+  *_log << "Navigation Systems : " << systemMap.size() << "   ";
 
   QMapIterator<QChar, QVector<const t_qcSatSum*> > itSys(systemMap);
@@ -829 +829 @@
   }
   *_log << endl;
+
+  // Observation types per system
+  // -----------------------------
+  for (int iSys = 0; iSys < obsFile->numSys(); iSys++) {
+    char sys = obsFile->system(iSys);
+    if (sys != ' ') {
+      *_log << "Observation Types " << sys << ":";
+      for (int iType = 0; iType < obsFile->nTypes(sys); iType++) {
+        QString type = obsFile->obsType(sys, iType);
+        *_log << " " << type;
+      }
+      *_log << endl;
+    }
+  }
 
   itSys.toFront();
@@ -836 +850 @@
     const QVector<const t_qcSatSum*>& qcSatVec = itSys.value();
     QString prefixSys = QString("  ") + sys + QString(": ");
-    QMap<QString, QVector<const t_qcFrqSum*> > frqMap; 
+    QMap<QString, QVector<const t_qcFrqSum*> > frqMap;
     for (int ii = 0; ii < qcSatVec.size(); ii++) {
       const t_qcSatSum* qcSatSum = qcSatVec[ii];