Changeset 10251 in ntrip for trunk/BNC/src/PPP/pppSatObs.cpp

Timestamp:

Nov 20, 2023, 10:23:42 PM (11 months ago)

Author:

stuerze

Message:

changes regarding PPP: allow single frequency PPP and allow to select the frequency bands that are used

File:

• trunk/BNC/src/PPP/pppSatObs.cpp (modified) (11 diffs)

trunk/BNC/src/PPP/pppSatObs.cpp

@@ -46 +46 @@
   _stecSat    = 0.0;
   _signalPriorities = QString::fromStdString(OPT->_signalPriorities);
-  if (!_signalPriorities.size()) {
-    _signalPriorities = "G:12&CWPSLX R:12&CP E:1&CBX E:5&QIX C:26&IQX";
-  }
-
   for (unsigned ii = 0; ii < t_frequency::max; ii++) {
     _obs[ii] = 0;
@@ -66 +62 @@
 //
 ////////////////////////////////////////////////////////////////////////////
-void t_pppSatObs::prepareObs(const t_satObs& pppSatObs) {
+void t_pppSatObs::prepareObs(const t_satObs& pppSatObs) {//cout << "SATELLITE: " << _prn.toString() << endl;
 
   _model.reset();
+
+  std::vector<char> bb = OPT->frqBands(_prn.system());
+  char frqNum1 = '0';
+  if (bb.size() >= 1) {
+    frqNum1 = bb[0];
+  }
+  char frqNum2 = '0';
+  if (bb.size() == 2) {
+    frqNum2 = bb[1];
+  }
 
   // Select pseudo-ranges and phase observations
@@ -74 +80 @@
   QStringList priorList = _signalPriorities.split(" ", QString::SkipEmptyParts);
   string preferredAttrib;
-  char obsSys = pppSatObs._prn.system();  //cout << "SATELLITE: " << pppSatObs._prn.toString() << endl;
   for (unsigned iFreq = 1; iFreq < t_frequency::max; iFreq++) {
     t_frequency::type frqType = static_cast<t_frequency::type>(iFreq);
     char frqSys = t_frequency::toString(frqType)[0]; //cout << "frqSys: " << frqSys << endl;
     char frqNum = t_frequency::toString(frqType)[1]; //cout << "frqNum: " << frqNum << endl;
-    if (obsSys != frqSys) {
+    if (frqSys != _prn.system()) {
+      continue;
+    }
+    if (frqNum != frqNum1 &&
+        frqNum != frqNum2 ) {
       continue;
     }
@@ -100 +109 @@
             //cout << "observation2char: " << obs->_rnxType2ch << " vs. " << obsType.toStdString().c_str()<< endl;
             if (obs->_rnxType2ch == obsType.toStdString() &&
-                obs->_codeValid && obs->_code &&
+                obs->_codeValid  && obs->_code &&
                 obs->_phaseValid && obs->_phase) {
               _obs[iFreq] = new t_frqObs(*obs); //cout << "================> newObs: " << obs->_rnxType2ch <<endl;
@@ -112 +121 @@
   // Used frequency types
   // --------------------
-  _fType1 = t_lc::toFreq(_prn.system(),t_lc::l1);
-  _fType2 = t_lc::toFreq(_prn.system(),t_lc::l2);
+  _fType1 = t_frqBand::toFreq(_prn.system(), frqNum1);
+  _fType2 = t_frqBand::toFreq(_prn.system(), frqNum2);
 
   // Check whether all required frequencies available
@@ -524 +533 @@
         }
         const t_frqObs* obs = _obs[iFreq];
-        if (obs && obs->_rnxType2ch == bias._rnxType2ch) {
+        if (obs &&
+            obs->_rnxType2ch == bias._rnxType2ch) {
           _model._codeBias[iFreq]  = bias._value;
         }
@@ -551 +561 @@
         }
         const t_frqObs* obs = _obs[iFreq];
-        if (obs && obs->_rnxType2ch == bias._rnxType2ch) {
+        if (obs &&
+            obs->_rnxType2ch == bias._rnxType2ch) {
           _model._phaseBias[iFreq]  = bias._value;
         }
@@ -606 +617 @@
   _model._set = true;
 
-  //printModel();
+  printModel();
 
   return success;
@@ -619 +630 @@
 
       << "\n======================= " << endl
-     << "PPP STRATEGY  : " <<  OPT->_obsmodelTypeStr.at((int)OPT->_obsModelType).toLocal8Bit().constData()
+      << "PPP "
       <<  ((OPT->_pseudoObsIono) ? " with pseudo-observations for STEC" : "")          << endl
       << "RHO           : " << setw(12) << setprecision(3) << _model._rho              << endl
@@ -701 +712 @@
       t_frequency::type tFreq = it->first;
       dispPart += it->second * (_model._antPCO[tFreq] - _model._codeBias[tFreq]);
-      if (OPT->PPPRTK) {
+      if (OPT->PPP_RTK) {
         dispPart += it->second * (_model._ionoCodeDelay[tFreq]);
       }
@@ -709 +720 @@
       dispPart += it->second * (_model._antPCO[tFreq] - _model._phaseBias[tFreq] +
                                 _model._windUp * t_CST::lambda(tFreq, _channel));
-      if (OPT->PPPRTK) {
+      if (OPT->PPP_RTK) {
         dispPart += it->second * (- _model._ionoCodeDelay[tFreq]);
       }