Changeset 9635 in ntrip for trunk/BNC/src/combination/bnccomb.h

Timestamp:

Feb 24, 2022, 5:08:42 PM (2 years ago)

Author:

stuerze

Message:

consideration of incoming Code Biases during clock combination

File:

• trunk/BNC/src/combination/bnccomb.h (modified) (5 diffs)

trunk/BNC/src/combination/bnccomb.h

@@ -4 +4 @@
 
 #include <fstream>
+#include <map>
 #include <newmat.h>
 #include "bncephuser.h"
 #include "satObs.h"
+#include "bncconst.h"
 #include "../RTCM3/clock_and_orbit/clock_orbit_rtcm.h"
 #include "../RTCM3/clock_and_orbit/clock_orbit_igs.h"
@@ -80 +82 @@
       _iod        = 0;
       _dClkResult = 0.0;
+      _codeBiasIF   = 0.0;
     }
     ~cmbCorr() {}
@@ -90 +93 @@
     t_satCodeBias _satCodeBias;
     QString       _acName;
+    double        _codeBiasIF;
     double        _dClkResult;
     ColumnVector  _diffRao;
@@ -106 +110 @@
     QVector<cmbCorr*> corrs;
   };
+
+  class t_lcRefSig {
+   public:
+    enum type {dummy = 0, c1, c2, cIF};
+
+    static t_frequency::type toFreq(char sys, type tt) {
+      switch (tt) {
+      case c1:
+        if      (sys == 'G') return t_frequency::G1;
+        else if (sys == 'R') return t_frequency::R1;
+        else if (sys == 'E') return t_frequency::E1;
+        else if (sys == 'C') return t_frequency::C2;
+        else if (sys == 'J') return t_frequency::J1;
+        else if (sys == 'S') return t_frequency::S1;
+        else                 return t_frequency::dummy;
+      case c2:
+        if      (sys == 'G') return t_frequency::G2;
+        else if (sys == 'R') return t_frequency::R2;
+        else if (sys == 'E') return t_frequency::E5;
+        else if (sys == 'C') return t_frequency::C6;
+        else if (sys == 'J') return t_frequency::J2;
+        else if (sys == 'S') return t_frequency::S5;
+        else                 return t_frequency::dummy;
+      case dummy:
+      case cIF:
+        return t_frequency::dummy;
+      }
+      return t_frequency::dummy;
+    }
+
+    static char toAttrib(char sys, type LC) {
+      switch (LC) {
+        case c1:
+          if      (sys == 'G') return 'W';
+          else if (sys == 'R') return 'P';
+          else if (sys == 'E') return 'C';
+          else if (sys == 'C') return 'I';
+          else if (sys == 'J') return 'C';
+          else if (sys == 'S') return 'C';
+          break;
+        case c2:
+          if      (sys == 'G') return 'W';
+          else if (sys == 'R') return 'P';
+          else if (sys == 'E') return 'Q';
+          else if (sys == 'C') return 'I';
+          else if (sys == 'J') return 'L';
+          else if (sys == 'S') return 'Q';
+          break;
+      case dummy:
+      case cIF:
+        return '_';
+        break;
+      }
+      return '_';
+    }
+
+      static void coeff(char sys, type tLC, double channel, std::map<t_frequency::type, double>& codeCoeff)  {
+      codeCoeff.clear();
+      t_frequency::type fType1 = toFreq(sys, c1);
+      t_frequency::type fType2 = toFreq(sys, c2);
+      double f1 = t_CST::freq(fType1, channel);
+      double f2 = t_CST::freq(fType2, channel);
+      switch (tLC) {
+        case c1:
+          codeCoeff[fType1] = 1.0;
+          return;
+        case c2:
+          codeCoeff[fType2] = 1.0;
+          return;
+        case cIF:
+          codeCoeff[fType1] =  f1 * f1 / (f1 * f1 - f2 * f2);
+          codeCoeff[fType2] = -f2 * f2 / (f1 * f1 - f2 * f2);
+          return;
+        case t_lcRefSig::dummy:
+          return;
+      }
+      return;
+    }
+  };
+
 
   void  processEpoch(char sys);
@@ -134 +218 @@
   e_method                                   _method;
   int                                        _cmbSampl;
+  QString                                    _cmbRefAttributes;
   QMap<QString, QMap<t_prn, t_orbCorr> >     _orbCorrections;
   QMap<QString, QMap<t_prn, t_satCodeBias> > _satCodeBiases;