Changeset 10021 in ntrip

Timestamp:

Apr 3, 2023, 11:02:58 AM (13 months ago)

Author:

stuerze

Message:

reset ambiguitis as soon as an outier appears

File:

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

trunk/BNC/src/PPP/pppFilter.cpp

@@ -100 +100 @@
       }
     }
-    if (num < OPT->_minObs) {
-      setNeuNoiseToZero = true;
-    }
+//    if (num < OPT->_minObs) {
+//      setNeuNoiseToZero = true;
+//    }
   }
   setStateVectorAndVarCovMatrix(xFltOld, QFltOld, setNeuNoiseToZero);
@@ -275 +275 @@
     // Check number of observations
     // ----------------------------
-    if (iObs == -1) {
-      LOG << " number of observations == " << iObs + 1 << "\n";
+    if ((iObs +1) < OPT->_minObs) {
+      LOG << "t_pppFilter::processSystem not enough observations: " << iObs + 1 << "\n";
       if (preProcessing) {
         _obsPool->setRefSatChangeRequired(sys, true);
@@ -371 +371 @@
           obs->setOutlier();
         }
-      } else {    // fin-processing
-        LOG << epoTimeStr << " Outlier " << t_lc::toString(maxOutlierLC) << ' '
-            << obs->prn().toString() << ' ' << setw(8) << setprecision(4)
-            << maxOutlier << endl;
-        if (par) {
-          if (par->ambResetCandidate() ||
-              OPT->_obsModelType == OPT->DCMcodeBias ||
-              OPT->_obsModelType == OPT->DCMphaseBias) {
-            resetAmb(par->prn(), obsVector, maxOutlierLC, &QSav, &xSav);
-          }
-          else {
-            par->setAmbResetCandidate();
+      }
+      else {    // fin-processing
+        if (obs->prn() != refPrn) {
+          LOG << epoTimeStr << " Outlier " << t_lc::toString(maxOutlierLC)
+              << ' ' << obs->prn().toString() << ' ' << setw(8)
+              << setprecision(4) << maxOutlier << endl;
+          if (par) {
+              resetAmb(par->prn(), obsVector, maxOutlierLC, &QSav, &xSav);
+            } else {
             obs->setOutlier();
           }
-        }
-        else {
-          obs->setOutlier();
         }
       }