Changeset 9288 in ntrip for trunk/BNC/src/PPP/pppClient.cpp

Timestamp:

Nov 27, 2020, 10:38:03 AM (3 years ago)

Author:

stuerze

Message:

updates regarding inter system biases

File:

• trunk/BNC/src/PPP/pppClient.cpp (modified) (9 diffs)

trunk/BNC/src/PPP/pppClient.cpp

@@ -73 +73 @@
     }
   }
-  _offGG = 0.0;
-  _offGB = 0.0;
+  _offGR = 0.0;
+  _offGE = 0.0;
+  _offGC = 0.0;
   CLIENTS.setLocalData(this);  // CLIENTS takes ownership over "this"
 }
@@ -332 +333 @@
 // Compute A Priori GPS-Glonass Offset
 //////////////////////////////////////////////////////////////////////////////
-double t_pppClient::cmpOffGG(vector<t_pppSatObs*>& obsVector) {
+double t_pppClient::cmpOffGR(vector<t_pppSatObs*>& obsVector) {
 
   t_lc::type tLC   = t_lc::dummy;
-  double     offGG = 0.0;
+  double     offGR = 0.0;
 
   if (_opt->useSystem('R')) {
     while (obsVector.size() > 0) {
-      offGG = 0.0;
+      offGR = 0.0;
       double   maxRes      = 0.0;
       int      maxResIndex = -1;
@@ -353 +354 @@
             double ll = satObs->obsValue(tLC) - satObs->cmpValue(tLC);
             ++nObs;
-            offGG += ll;
+            offGR += ll;
             if (fabs(ll) > fabs(maxRes)) {
               maxRes      = ll;
@@ -364 +365 @@
 
       if (nObs > 0) {
-        offGG = offGG / nObs;
+        offGR = offGR / nObs;
       }
       else {
-        offGG = 0.0;
+        offGR = 0.0;
       }
 
       if (fabs(maxRes) > 1000.0) {
-        LOG << "t_pppClient::cmpOffGG outlier " << maxResPrn.toString() << " " << maxRes << endl;
+        LOG << "t_pppClient::cmpOffGR outlier " << maxResPrn.toString() << " " << maxRes << endl;
         delete obsVector.at(maxResIndex);
         obsVector.erase(obsVector.begin() + maxResIndex);
@@ -380 +381 @@
     }
   }
-  return offGG;
+  return offGR;
+}
+
+// Compute A Priori GPS-Galileo Offset
+//////////////////////////////////////////////////////////////////////////////
+double t_pppClient::cmpOffGE(vector<t_pppSatObs*>& obsVector) {
+
+  t_lc::type tLC   = t_lc::dummy;
+  double     offGE = 0.0;
+
+  if (_opt->useSystem('E')) {
+    while (obsVector.size() > 0) {
+      offGE = 0.0;
+      double   maxRes      = 0.0;
+      int      maxResIndex = -1;
+      t_prn    maxResPrn;
+      unsigned nObs        = 0;
+      for (unsigned ii = 0; ii < obsVector.size(); ii++) {
+        const t_pppSatObs* satObs = obsVector.at(ii);
+        if (satObs->prn().system() == 'E') {
+          if (tLC == t_lc::dummy) {
+            tLC = satObs->isValid(t_lc::cIF) ? t_lc::cIF : t_lc::c1;
+          }
+          if (satObs->isValid(tLC) && (!satObs->modelSet() || satObs->eleSat() >= _opt->_minEle)) {
+            double ll = satObs->obsValue(tLC) - satObs->cmpValue(tLC);
+            ++nObs;
+            offGE += ll;
+            if (fabs(ll) > fabs(maxRes)) {
+              maxRes      = ll;
+              maxResIndex = ii;
+              maxResPrn   = satObs->prn();
+            }
+          }
+        }
+      }
+
+      if (nObs > 0) {
+        offGE = offGE / nObs;
+      }
+      else {
+        offGE = 0.0;
+      }
+
+      if (fabs(maxRes) > 1000.0) {
+        LOG << "t_pppClient::cmpOffGE outlier " << maxResPrn.toString() << " " << maxRes << endl;
+        delete obsVector.at(maxResIndex);
+        obsVector.erase(obsVector.begin() + maxResIndex);
+      }
+      else {
+        break;
+      }
+    }
+  }
+  return offGE;
 }
 
 // Compute A Priori GPS-BDS Offset
 //////////////////////////////////////////////////////////////////////////////
-double t_pppClient::cmpOffGB(vector<t_pppSatObs*>& obsVector) {
+double t_pppClient::cmpOffGC(vector<t_pppSatObs*>& obsVector) {
 
   t_lc::type tLC   = t_lc::dummy;
-  double     offGB = 0.0;
+  double     offGC = 0.0;
 
   if (_opt->useSystem('C')) {
     while (obsVector.size() > 0) {
-      offGB = 0.0;
+      offGC = 0.0;
       double   maxRes      = 0.0;
       int      maxResIndex = -1;
@@ -406 +460 @@
             double ll = satObs->obsValue(tLC) - satObs->cmpValue(tLC);
             ++nObs;
-            offGB += ll;
+            offGC += ll;
             if (fabs(ll) > fabs(maxRes)) {
               maxRes      = ll;
@@ -417 +471 @@
 
       if (nObs > 0) {
-        offGB = offGB / nObs;
+        offGC = offGC / nObs;
       }
       else {
-        offGB = 0.0;
+        offGC = 0.0;
       }
 
       if (fabs(maxRes) > 1000.0) {
-        LOG << "t_pppClient::cmpOffGB outlier " << maxResPrn.toString() << " " << maxRes << endl;
+        LOG << "t_pppClient::cmpOffGC outlier " << maxResPrn.toString() << " " << maxRes << endl;
         delete obsVector.at(maxResIndex);
         obsVector.erase(obsVector.begin() + maxResIndex);
@@ -433 +487 @@
     }
   }
-  return offGB;
+  return offGC;
 }
 
@@ -570 +624 @@
       }
 
-      _offGG = cmpOffGG(_obsRover);
-      _offGB = cmpOffGB(_obsRover);
+      _offGR = cmpOffGR(_obsRover);
+      _offGE = cmpOffGE(_obsRover);
+      _offGC = cmpOffGC(_obsRover);
 
       if (_opt->_refSatRequired) {