Changeset 8961 in ntrip for trunk/BNC/src/pppMain.cpp

Timestamp:

Jun 29, 2020, 10:09:16 AM (4 years ago)

Author:

stuerze

Message:

PPP update: pseudo obs tropo added

File:

• trunk/BNC/src/pppMain.cpp (modified) (14 diffs)

trunk/BNC/src/pppMain.cpp

@@ -185 +185 @@
     opt->_obsModelType   = t_pppOptions::IF;
     opt->_pseudoObsIono  = false;
+    opt->_pseudoObsTropo = false;
     opt->_refSatRequired = false;
 #ifdef USE_PPP
     // Pseudo Observations
     if      (settings.value("PPP/pseudoObs").toString() == "Ionosphere") {
-      opt->_pseudoObsIono = true;
+      opt->_pseudoObsIono  = true;
+      opt->_pseudoObsTropo = false;
+    }
+    else if (settings.value("PPP/pseudoObs").toString() == "Iono+Tropo") {
+      opt->_pseudoObsIono  = true;
+      opt->_pseudoObsTropo = true;
     }
     else if (settings.value("PPP/pseudoObs").toString() == "no") {
-      opt->_pseudoObsIono = false;
+      opt->_pseudoObsIono  = false;
+      opt->_pseudoObsTropo = false;
     }
     // Observation Model
     if      (settings.value("PPP/modelObs").toString() == "Ionosphere-free PPP") {
       opt->_obsModelType = t_pppOptions::IF;
-      opt->_pseudoObsIono = false;
+      opt->_pseudoObsIono  = false;
+      opt->_pseudoObsTropo = false;
     }
     else if (settings.value("PPP/modelObs").toString() == "PPP-RTK") {
       opt->_obsModelType = t_pppOptions::PPPRTK;
-      opt->_pseudoObsIono = false;
+      opt->_pseudoObsIono  = false;
+      opt->_pseudoObsTropo = false;
     }
     else if (settings.value("PPP/modelObs").toString() == "Uncombined PPP") {
@@ -229 +238 @@
           opt->_LCsGPS.push_back(t_lc::GIM);
         }
+        if (opt->_pseudoObsTropo) {
+          opt->_LCsGPS.push_back(t_lc::Tz0);
+        }
       }
     }
@@ -240 +252 @@
         if (opt->_pseudoObsIono) {
           opt->_LCsGPS.push_back(t_lc::GIM);
+        }
+        if (opt->_pseudoObsTropo) {
+          opt->_LCsGPS.push_back(t_lc::Tz0);
         }
       }
@@ -256 +271 @@
           opt->_LCsGPS.push_back(t_lc::GIM);
         }
+        if (opt->_pseudoObsTropo) {
+          opt->_LCsGPS.push_back(t_lc::Tz0);
+        }
       }
     }
@@ -268 +286 @@
         if (opt->_pseudoObsIono) {
           opt->_LCsGLONASS.push_back(t_lc::GIM);
+        }
+        if (opt->_pseudoObsTropo) {
+          opt->_LCsGPS.push_back(t_lc::Tz0);
         }
       }
@@ -281 +302 @@
           opt->_LCsGLONASS.push_back(t_lc::GIM);
         }
+        if (opt->_pseudoObsTropo) {
+          opt->_LCsGPS.push_back(t_lc::Tz0);
+        }
       }
     }
@@ -296 +320 @@
           opt->_LCsGLONASS.push_back(t_lc::GIM);
         }
+        if (opt->_pseudoObsTropo) {
+          opt->_LCsGPS.push_back(t_lc::Tz0);
+        }
       }
     }
@@ -309 +336 @@
           opt->_LCsGalileo.push_back(t_lc::GIM);
         }
+        if (opt->_pseudoObsTropo) {
+          opt->_LCsGPS.push_back(t_lc::Tz0);
+        }
       }
     }
@@ -320 +350 @@
         if (opt->_pseudoObsIono) {
           opt->_LCsGalileo.push_back(t_lc::GIM);
+        }
+        if (opt->_pseudoObsTropo) {
+          opt->_LCsGPS.push_back(t_lc::Tz0);
         }
       }
@@ -336 +369 @@
           opt->_LCsGalileo.push_back(t_lc::GIM);
         }
+        if (opt->_pseudoObsTropo) {
+          opt->_LCsGPS.push_back(t_lc::Tz0);
+        }
       }
     }
@@ -349 +385 @@
           opt->_LCsBDS.push_back(t_lc::GIM);
         }
+        if (opt->_pseudoObsTropo) {
+          opt->_LCsGPS.push_back(t_lc::Tz0);
+        }
       }
     }
@@ -360 +399 @@
         if (opt->_pseudoObsIono) {
           opt->_LCsBDS.push_back(t_lc::GIM);
+        }
+        if (opt->_pseudoObsTropo) {
+          opt->_LCsGPS.push_back(t_lc::Tz0);
         }
       }
@@ -375 +417 @@
         if (opt->_pseudoObsIono) {
           opt->_LCsBDS.push_back(t_lc::GIM);
+        }
+        if (opt->_pseudoObsTropo) {
+          opt->_LCsGPS.push_back(t_lc::Tz0);
         }
       }
@@ -417 +462 @@
     opt->_aprSigPhaseBias = 1000.0;
     // TODO: Find realistic values!!!!!!
-    opt->_noiseIon        = 10.0;
+    opt->_noiseIon        = 1.00;
     opt->_noiseCodeBias   = 1.00;
     opt->_noisePhaseBias  = 5.00;
-    opt->_sigmaGIMdiff    = 2.00; //pseudo observation GIM: STEC(ref_sat) - STEC(sat)
+    // pseudo observations
+    opt->_sigmaGIMdiff    = 2.00; // Todo: add to bncWindow
+    opt->_sigmaTz0        = 0.10;
 
     _options << opt;