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

Timestamp:

Dec 2, 2021, 5:18:27 PM (2 years ago)

Author:

stuerze

Message:

update regarding PPP

File:

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

trunk/BNC/src/PPP/pppClient.cpp

@@ -243 +243 @@
   return pseudoObsIono;
 }
+
+//
+//////////////////////////////////////////////////////////////////////////////
+void t_pppClient::useObsWithCodeBiasesOnly(std::vector<t_pppSatObs*>& obsVector) {
+
+  vector<t_pppSatObs*>::iterator it = obsVector.begin();
+  while (it != obsVector.end()) {
+    t_pppSatObs* pppSatObs = *it;
+    char sys = pppSatObs->prn().system();
+    bool codeBiasesAvailable = false;
+    t_frequency::type fType1 = t_lc::toFreq(sys,t_lc::c1);
+    t_frequency::type fType2 = t_lc::toFreq(sys,t_lc::c2);
+    if (pppSatObs->getCodeBias(fType1) &&
+        pppSatObs->getCodeBias(fType2)) {
+      codeBiasesAvailable = true;
+    }
+    if (codeBiasesAvailable) {
+      ++it;
+    }
+    else {
+      it = obsVector.erase(it);
+      delete pppSatObs;
+    }
+  }
+
+}
+
 
 // Compute the Bancroft position, check for blunders
@@ -637 +664 @@
       }
 
-      if (int(_obsRover.size()) < _opt->_minObs) {
-        LOG << "t_pppClient::processEpoch not enough observations" << endl;
-        return finish(failure,5);
-      }
-
       if (_opt->_refSatRequired) {
         if (handleRefSatellites(_obsRover) != success) {
@@ -651 +673 @@
           continue;
         }
+      }
+
+      // Check if obs have code biases
+      // -----------------------------
+      if (OPT->_obsModelType == OPT->DCMcodeBias ||
+          OPT->_obsModelType == OPT->DCMphaseBias) {
+        useObsWithCodeBiasesOnly(_obsRover);
+      }
+
+      if (int(_obsRover.size()) < _opt->_minObs) {
+        LOG << "t_pppClient::processEpoch not enough observations" << endl;
+        return finish(failure,5);
       }