Changeset 10328 in ntrip

Timestamp:

Feb 9, 2024, 12:00:35 PM (10 months ago)

Author:

stuerze

Message:

hDOP and vDOP values are now computed from the variances of the NEU components of the receiver position estimates

File:

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

trunk/BNC/src/PPP/pppFilter.cpp

@@ -465 +465 @@
     const unsigned numPar = 4;
     Matrix AA(obsVector.size(), numPar);
+    t_pppParam* parX = 0;
+    t_pppParam* parY = 0;
+    t_pppParam* parZ = 0;
     _numSat = 0;
     for (unsigned ii = 0; ii < obsVector.size(); ii++) {
@@ -471 +474 @@
         ++_numSat;
         for (unsigned iPar = 0; iPar < numPar; iPar++) {
-          const t_pppParam* par = _parlist->params()[iPar];
+          t_pppParam* par = _parlist->params()[iPar];
           AA[_numSat - 1][iPar] = par->partial(_epoTime, obs, t_lc::c1);
+          if      (par->type() == t_pppParam::crdX) {
+            parX = par;
+          }
+          else if (par->type() == t_pppParam::crdY) {
+            parY = par;
+          }
+          else if (par->type() == t_pppParam::crdZ) {
+            parZ = par;
+          }
         }
       }
@@ -482 +494 @@
     SymmetricMatrix NN; NN << AA.t() * AA;
     SymmetricMatrix QQ = NN.i();
-
-    _dop.H = sqrt(QQ(1, 1) + QQ(2, 2));
-    _dop.V = sqrt(QQ(3, 3));
+    SymmetricMatrix QQxyz = QQ.SymSubMatrix(1,3);
+
+    ColumnVector xyz(3), neu(3);
+    SymmetricMatrix QQneu(3);
+    const t_pppStation *sta = PPP_CLIENT->staRover();
+    xyz[0] = _xFlt[parX->indexNew()];
+    xyz[1] = _xFlt[parY->indexNew()];
+    xyz[2] = _xFlt[parZ->indexNew()];
+    xyz2neu(sta->ellApr().data(), xyz.data(), neu.data());
+    covariXYZ_NEU(QQxyz, sta->ellApr().data(), QQneu);
+
+    _dop.H = sqrt(QQneu(1, 1) + QQneu(2, 2));
+    _dop.V = sqrt(QQneu(3, 3));
     _dop.P = sqrt(QQ(1, 1) + QQ(2, 2) + QQ(3, 3));
     _dop.T = sqrt(QQ(4, 4));