Changeset 7050 in ntrip for trunk/BNC/src/PPP_SSR_I/pppFilter.cpp

Timestamp:

Jul 14, 2015, 4:13:56 PM (9 years ago)

Author:

stuerze

Message:

small changes which allow to the user to decide whether he wants to use code and/or phase observations of GLONASS, Galileo and BDS in SSR I PPP mode

File:

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

trunk/BNC/src/PPP_SSR_I/pppFilter.cpp

@@ -1007 +1007 @@
 
     for (int iPhase = 0; iPhase <= (usePhase ? 1 : 0); iPhase++) {
-    
+
       // Status Prediction
       // -----------------
@@ -1017 +1017 @@
       unsigned nObs = 0;
       nObs = epoData->sizeAll();
+      bool useObs = false;
+      char additionalSys[] ={'R', 'E', 'C'};
+      for (unsigned ii = 0; ii < sizeof(additionalSys); ii++) {
+        const char s = additionalSys[ii];
+        (iPhase == 0) ? useObs = OPT->codeLCs(s).size() : useObs = OPT->ambLCs(s).size();
+        if (!useObs) {
+          nObs -= epoData->sizeSys(s);
+        }
+      }
       
       // Prepare first-design Matrix, vector observed-computed
@@ -1030 +1039 @@
         t_satData* satData = it.value();
         QString prn = satData->prn;
-        addObs(iPhase, iObs, satData, AA, ll, PP);
+        char sys =   satData->system();
+        (iPhase == 0) ? useObs = OPT->codeLCs(sys).size() : useObs = OPT->ambLCs(sys).size();
+        if (sys == 'G' || useObs) {
+          addObs(iPhase, iObs, satData, AA, ll, PP);
+        }
       }