Changeset 7142 in ntrip

Timestamp:

Jul 22, 2015, 4:18:09 PM (9 years ago)

Author:

stuerze

Message:

some more changes to consider the user's decision regarding the observations that have to be used

File:

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

trunk/BNC/src/PPP_SSR_I/pppFilter.cpp

@@ -531 +531 @@
       // ----------
       else if (pp->type == t_pppParam::BDS_OFFSET) {
-        _QQ(iPar,iPar) += 0.1 * 0.1;    //TODO: TEST
+        _QQ(iPar,iPar) = 1000.0 * 1000.0;    //TODO: TEST
       }
     }
@@ -848 +848 @@
 void t_pppFilter::addAmb(t_satData* satData) {
   Tracer tracer("t_pppFilter::addAmb");
+  if (!OPT->ambLCs(satData->system()).size()){
+    return;
+  }
   bool    found = false;
   for (int iPar = 1; iPar <= _params.size(); iPar++) {
@@ -909 +912 @@
   else {
     double sigP3 = 2.98 * OPT->_sigmaC1;
+    if  (satData->system() == 'C') {
+      sigP3 *= BDS_WEIGHT_FACTOR;
+    }
     ll(iObs)      = satData->P3 - cmpValue(satData, false);
     PP(iObs,iObs) = 1.0 / (sigP3 * sigP3) / (ellWgtCoef * ellWgtCoef);
@@ -926 +932 @@
   ostringstream str;
   str.setf(ios::fixed);
+  bool useObs;
         
   QMapIterator<QString, t_satData*> it(satDataMap);
@@ -931 +938 @@
     it.next();
     t_satData* satData = it.value();
-    if (satData->obsIndex != 0) {
+    (iPhase == 0) ? useObs = OPT->codeLCs(satData->system()).size() :
+                    useObs = OPT->ambLCs(satData->system()).size();
+    if (satData->obsIndex != 0 && useObs) {
       str << _time.timestr(1)
           << " RES " << satData->prn.mid(0,3).toAscii().data()
@@ -1018 +1027 @@
       nObs = epoData->sizeAll();
       bool useObs = false;
-      char additionalSys[] ={'R', 'E', 'C'};
-      for (unsigned ii = 0; ii < sizeof(additionalSys); ii++) {
-        const char s = additionalSys[ii];
+      char sys[] ={'G', 'R', 'E', 'C'};
+      for (unsigned ii = 0; ii < sizeof(sys); ii++) {
+        const char s = sys[ii];
         (iPhase == 0) ? useObs = OPT->codeLCs(s).size() : useObs = OPT->ambLCs(s).size();
         if (!useObs) {
@@ -1039 +1048 @@
         t_satData* satData = it.value();
         QString prn = satData->prn;
-        char sys =   satData->system();
-        (iPhase == 0) ? useObs = OPT->codeLCs(sys).size() : useObs = OPT->ambLCs(sys).size();
-        if (sys == 'G' || useObs) {
+        (iPhase == 0) ? useObs = OPT->codeLCs(satData->system()).size() :
+                        useObs = OPT->ambLCs(satData->system()).size();
+        if (useObs) {
           addObs(iPhase, iObs, satData, AA, ll, PP);
         }
@@ -1054 +1063 @@
       // Print Residuals
       // ---------------
-      if (iPhase == 0) {
+      if      (iPhase == 0) {
         strResCode  = printRes(iPhase, vv, epoData->satData);
       }