Changeset 9395 in ntrip for trunk/BNC/src/PPP/pppFilter.cpp

Timestamp:

Mar 30, 2021, 9:59:17 AM (3 years ago)

Author:

stuerze

Message:

update regarding PPP

File:

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

trunk/BNC/src/PPP/pppFilter.cpp

@@ -138 +138 @@
         }
       }
+      if (!obsVector.size()) {continue;}
       if (processSystem(OPT->LCs(sys), obsVector, _refPrn,
                         epoch->pseudoObsIono(), preProcessing) != success) {
@@ -174 +175 @@
       }
     }
+    if (!num) {continue;}
     LOG << epoTimeStr << " SATNUM " << system << ' ' << right << setw(2) << num << endl;
     if (processSystem(OPT->LCs(system), obsVector, _refPrn,
@@ -266 +268 @@
           if (tLC == t_lc::GIM) {continue;}
           if (tLC == t_lc::Tz0) {continue;}
+          ++iObs;
+          usedObs.push_back(obs);
+          usedTypes.push_back(tLC);
+          for (unsigned iPar = 0; iPar < params.size(); iPar++) {
+            const t_pppParam* par = params[iPar];
+            AA[iObs][iPar] = par->partial(_epoTime, obs, tLC, refPrn);
+          }
+          ll[iObs] = obs->obsValue(tLC) - obs->cmpValue(tLC) - DotProduct(_x0, AA.Row(iObs+1));
+          PP[iObs] = 1.0 / (obs->sigma(tLC) * obs->sigma(tLC));
+        }
+      }
+    }
+
+    // pseudo Obs Tropo
+    // ================
+    if (OPT->_pseudoObsTropo) {
+      for (unsigned ii = 0; ii < obsVector.size(); ii++) {
+        t_pppSatObs* obs = obsVector[ii];
+        if (!obs->isReference()) {continue;}
+        for (unsigned jj = 0; jj < usedLCs; jj++) {
+          const t_lc::type tLC = LCs[jj];
+          if (tLC != t_lc::Tz0) {continue;}
           ++iObs;
           usedObs.push_back(obs);
@@ -309 +333 @@
       }
     }
-    // pseudo Obs Tropo
-    // ================
-    if (OPT->_pseudoObsTropo) {
-      for (unsigned ii = 0; ii < obsVector.size(); ii++) {
-        t_pppSatObs* obs = obsVector[ii];
-        if (!obs->isReference()) {continue;}
-        for (unsigned jj = 0; jj < usedLCs; jj++) {
-          const t_lc::type tLC = LCs[jj];
-          if (tLC != t_lc::Tz0) {continue;}
-          ++iObs;
-          usedObs.push_back(obs);
-          usedTypes.push_back(tLC);
-          for (unsigned iPar = 0; iPar < params.size(); iPar++) {
-            const t_pppParam* par = params[iPar];
-            AA[iObs][iPar] = par->partial(_epoTime, obs, tLC, refPrn);
-          }
-          ll[iObs] = obs->obsValue(tLC) - obs->cmpValue(tLC) - DotProduct(_x0, AA.Row(iObs+1));
-          PP[iObs] = 1.0 / (obs->sigma(tLC) * obs->sigma(tLC));
-        }
-      }
-    }
 
     // Check number of observations, truncate matrices
@@ -561 +564 @@
 }
 
-// Add innfinite noise to iono
+// Add infinite noise to iono
 ////////////////////////////////////////////////////////////////////////////
 t_irc t_pppFilter::addNoiseToIono(char sys) {
@@ -745 +748 @@
       }
     }
+    // pseudo Obs Tropo
+    // ================
+    if (OPT->_pseudoObsTropo) {
+      for (unsigned ii = 0; ii < obsVector.size(); ii++) {
+        t_pppSatObs* obs = obsVector[ii];
+        if (!obs->isReference()) {continue;}
+        for (unsigned jj = 0; jj < usedLCs; jj++) {
+          const t_lc::type tLC = LCs[jj];
+          if (tLC != t_lc::Tz0) {continue;}
+          ++iObs;
+          for (unsigned iPar = 0; iPar < params.size(); iPar++) {
+            const t_pppParam* par = params[iPar];
+            AA[iObs][iPar] = par->partial(_epoTime, obs, tLC, refPrn);
+          }
+        }
+      }
+    }
     _datumTrafo->updateIndices(iObs+1);
     _datumTrafo->prepareAA(AA.SubMatrix(1, iObs+1 , 1, _parlist->nPar()), 2);
@@ -791 +811 @@
 ////////////////////////////////////////////////////////////////////////////
 void t_pppFilter::initDatumTransformation(const std::vector<t_pppSatObs*>& allObs) {
-  unsigned realObs = 0;
+  unsigned trafoObs = 0;
   for (unsigned iSys = 0; iSys < OPT->systems().size(); iSys++) {
     char system = OPT->systems()[iSys];
@@ -809 +829 @@
       realUsedLCs -= 1;
     }
-    realObs += satNum * realUsedLCs;
-  }
-  _datumTrafo->setObsNum(realObs);
+    trafoObs += satNum * realUsedLCs;
+
+    if (OPT->_pseudoObsTropo) {
+      trafoObs += 1;
+    }
+
+  }
+  _datumTrafo->setObsNum(trafoObs);
   _datumTrafo->setParNum(_parlist->nPar());
   _datumTrafo->initAA();