Changeset 6841 in ntrip for trunk/BNC/src/rinex/rnxobsfile.cpp

Timestamp:

May 19, 2015, 11:24:19 AM (9 years ago)

Author:

stuerze

Message:

consideration of obs header etries for phase shifts, GLONASS slots and GLONASS biases during merging of RINEX files

File:

• trunk/BNC/src/rinex/rnxobsfile.cpp (modified) (7 diffs)

trunk/BNC/src/rinex/rnxobsfile.cpp

@@ -244 +244 @@
         in >> type >> value;
         if (type.size())
-          _gloPhaseBiases[type] = value;
+          _gloBiases[type] = value;
       }
     }
@@ -345 +345 @@
 ////////////////////////////////////////////////////////////////////////////
 void t_rnxObsHeader::set(const t_rnxObsHeader& header, int version,
-                         const QStringList* useObsTypes) {
+                         const QStringList* useObsTypes,
+                         const QStringList* phaseShifts,
+                         const QStringList* gloBiases,
+                         const QStringList* gloSlots) {
 
   if (version <= 2) {
@@ -370 +373 @@
   _comments        = header._comments;
   _usedSystems     = header._usedSystems;
-  if (_version >= 3.0) {
-    _phaseShifts    = header._phaseShifts;
-    _gloPhaseBiases = header._gloPhaseBiases;
-    _gloSlots       = header._gloSlots;
-  }
+
   for (unsigned iPrn = 1; iPrn <= t_prn::MAXPRN_GPS; iPrn++) {
     _wlFactorsL1[iPrn] =  header._wlFactorsL1[iPrn]; 
@@ -456 +455 @@
     }
   }
+
+  if (_version >= 3.0) {
+    // set phase shifts
+    if (!phaseShifts ||  phaseShifts->empty()) {
+      _phaseShifts = header._phaseShifts;
+    }
+    else {
+      foreach (const QString &str, *phaseShifts) {
+        QStringList hlp  = str.split("_", QString::SkipEmptyParts);
+        QStringList hlp1 = hlp.last().split(":", QString::SkipEmptyParts);
+        QString type = hlp.first();
+        double shift = hlp1.first().toDouble();
+        hlp1.removeFirst();
+        QStringList &satList = hlp1;
+        QMap<QString, QPair<double, QStringList> >::iterator it = _phaseShifts.find(type);
+        if ( it != _phaseShifts.end()) {
+          it.value().second.append(satList);
+          it.value().second.removeDuplicates();
+        }
+        else {
+          _phaseShifts.insert(type, QPair<double, QStringList>(shift, satList));
+        }
+      }
+    }
+    // set GLONASS biases
+    if (!gloBiases || gloBiases->empty()) {
+      _gloBiases = header._gloBiases;
+    }
+    else {
+      foreach (const QString &str, *gloBiases) {
+        QStringList hlp = str.split(":", QString::SkipEmptyParts);
+        QString type = hlp.first();;
+        double  value = hlp.last().toDouble();
+        if (type.size())
+          _gloBiases[type] = value;
+      }
+    }
+    // set GLONASS slots
+    if (!gloSlots  || gloSlots->empty()) {
+      _gloSlots = header._gloSlots;
+    }
+    else {
+      foreach (const QString &str, *gloSlots) {
+        QStringList hlp = str.split(":", QString::SkipEmptyParts);
+        QString sat = hlp.first();
+        int    slot = hlp.last().toInt();
+        t_prn prn;
+        prn.set(sat.toStdString());
+        if(sat.size())
+          _gloSlots[prn] = slot;
+      }
+    }
+  }
 }
 
@@ -646 +698 @@
   if (_version >= 3.0) {
     QString hlp = "";
-    QMap<QString, double>::const_iterator it = _gloPhaseBiases.begin();
-    while (it != _gloPhaseBiases.end()){
+    QMap<QString, double>::const_iterator it = _gloBiases.begin();
+    while (it != _gloBiases.end()){
       hlp += QString("%1%2").arg(it.key(), 4).arg(it.value(), 9, 'f', 3);
       it++;
@@ -723 +775 @@
 }
 
+// Number of GLONASS biases
+////////////////////////////////////////////////////////////////////////////
+int t_rnxObsHeader::numGloBiases() const {
+  return _gloBiases.size();
+}
+
+// Number of GLONASS biases
+////////////////////////////////////////////////////////////////////////////
+int t_rnxObsHeader::numGloSlots() const {
+  return _gloSlots.size();
+}
+
 // Observation Type (satellite-system specific)
 ////////////////////////////////////////////////////////////////////////////
@@ -743 +807 @@
   }
   return "";
+}
+
+//
+////////////////////////////////////////////////////////////////////////////
+QStringList t_rnxObsHeader::phaseShifts() const {
+  QStringList strList;
+  QMap<QString, QPair<double, QStringList> >::const_iterator it =  _phaseShifts.begin();
+  while (it != _phaseShifts.end()) {
+    strList.append(QString("%1_%2:%3").arg(it.key(), 3).arg(it.value().first, 9, 'f', 3).arg(it.value().second.join("")));
+    it++;
+  }
+  return strList;
+}
+
+//
+////////////////////////////////////////////////////////////////////////////
+QStringList t_rnxObsHeader::gloBiases() const {
+  QStringList strList;
+  QMap<QString, double>::const_iterator it = _gloBiases.begin();
+  while (it != _gloBiases.end()) {
+    strList.append(QString("%1:%2").arg(it.key(), 3).arg(it.value(), 9, 'f', 3));
+    it++;
+  }
+  return strList;
+}
+
+//
+////////////////////////////////////////////////////////////////////////////
+QStringList t_rnxObsHeader::gloSlots() const {
+  QStringList strList;
+  QMap<t_prn, int>::const_iterator it = _gloSlots.begin();
+  while (it != _gloSlots.end()){
+    QString prn(it.key().toString().c_str());
+    strList.append(QString("%1:%2").arg(prn, 3).arg(it.value()));
+    it++;
+  }
+  return strList;
 }