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

Timestamp:

Jun 18, 2026, 4:02:38 PM (5 days ago)

Author:

stuerze

Message:

bugfix regarding ionospheric constraints

File:

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

trunk/BNC/src/PPP/pppClient.cpp

@@ -230 +230 @@
 
   if (_opt->_pseudoObsIono) {
-    vector<t_pppSatObs*>::iterator it = obsVector.begin();
-    while (it != obsVector.end()) {
-      t_pppSatObs* satObs = *it;
-      pseudoObsIono = satObs->setPseudoObsIono(t_frequency::G1);
-      it++;
+
+    // Select reference satellite per system: keep the current one if still
+    // visible; only switch to highest elevation when the current one disappears
+    // -------------------------------------------------------------------------
+    map<char, t_pppSatObs*> refSatMap;
+    map<char, t_pppSatObs*> bestEleMap;
+    for (t_pppSatObs* satObs : obsVector) {
+      satObs->resetReference();
+      double stec = satObs->getIonoCodeDelay(t_frequency::G1);
+      if (stec == 0.0) continue;
+      char sys = satObs->prn().system();
+      if (!bestEleMap.count(sys) || satObs->eleSat() > bestEleMap[sys]->eleSat()) {
+        bestEleMap[sys] = satObs;
+      }
+      if (_refPrnMap.count(sys) && satObs->prn() == _refPrnMap[sys]) {
+        refSatMap[sys] = satObs;
+      }
+    }
+    for (auto& kv : bestEleMap) {
+      if (!refSatMap.count(kv.first)) {
+        refSatMap[kv.first] = kv.second;
+      }
+    }
+    for (auto& kv : refSatMap) {
+      _refPrnMap[kv.first] = kv.second->prn();
+      kv.second->setAsReference();
+    }
+
+    if (bestEleMap.empty()) {
+      LOG << "GIM pseudo-obs unavailable: no valid STEC (vTec data missing?)" << endl;
+    }
+
+    // Set STEC pseudo-observations (single-differenced vs reference satellite)
+    // -------------------------------------------------------------------------
+    int nGIM = 0;
+    for (t_pppSatObs* satObs : obsVector) {
+      char sys = satObs->prn().system();
+      double stecRefSat = refSatMap.count(sys)
+                          ? refSatMap[sys]->getIonoCodeDelay(t_frequency::G1)
+                          : 0.0;
+      if (satObs->setPseudoObsIono(t_frequency::G1, stecRefSat)) {
+        pseudoObsIono = true;
+        nGIM++;
+      }
+    }
+    if (nGIM > 0) {
+      LOG << "GIM pseudo-obs: " << nGIM << " satellites, ref:";
+      for (auto& kv : refSatMap) {
+        LOG << ' ' << kv.second->prn().toString();
+      }
+      LOG << endl;
     }
   }
@@ -598 +644 @@
 
   LOG << "pppClient: reset" << endl;
+  _refPrnMap.clear();
   // to delete old orbit and clock corrections
   delete _ephPool;