Changeset 6288 in ntrip for trunk/BNC

Timestamp:

Nov 1, 2014, 2:54:28 PM (10 years ago)

Author:

mervart

Message:

 

Location:

trunk/BNC/src/rinex

Files:

• reqcanalyze.cpp (modified) (11 diffs)
• reqcanalyze.h (modified) (3 diffs)

trunk/BNC/src/rinex/reqcanalyze.cpp

@@ -163 +163 @@
         t_satObs satObs;
         t_rnxObsFile::setObsFromRnx(obsFile, _currEpo, rnxSat, satObs);
-        t_qcObs& qcObs = qcEpo._qcObs[satObs._prn];
-        setQcObs(qcEpo._epoTime, xyzSta, satObs, qcObs);
-        updateQcSat(qcObs, _qcFile._qcSatSum[satObs._prn]);
+        t_qcSat& qcSat = qcEpo._qcSat[satObs._prn];
+        setQcObs(qcEpo._epoTime, xyzSta, satObs, qcSat);
+        updateQcSat(qcSat, _qcFile._qcSatSum[satObs._prn]);
       }
       _qcFile._qcEpo.push_back(qcEpo);
@@ -245 +245 @@
 //
 ////////////////////////////////////////////////////////////////////////////
-void t_reqcAnalyze::updateQcSat(const t_qcObs& qcObs, t_qcSatSum& qcSatSum) {
-
-  for (int ii = 0; ii < qcObs._qcFrq.size(); ii++) {
-    const t_qcFrq& qcFrq    = qcObs._qcFrq[ii];
+void t_reqcAnalyze::updateQcSat(const t_qcSat& qcSat, t_qcSatSum& qcSatSum) {
+
+  for (int ii = 0; ii < qcSat._qcFrq.size(); ii++) {
+    const t_qcFrq& qcFrq    = qcSat._qcFrq[ii];
     t_qcFrqSum&    qcFrqSum = qcSatSum._qcFrqSum[qcFrq._rnxType2ch];
     qcFrqSum._numObs += 1;
@@ -263 +263 @@
 ////////////////////////////////////////////////////////////////////////////
 void t_reqcAnalyze::setQcObs(const bncTime& epoTime, const ColumnVector& xyzSta, 
-                             const t_satObs& satObs, t_qcObs& qcObs) {
+                             const t_satObs& satObs, t_qcSat& qcSat) {
 
   t_eph* eph = 0;
@@ -278 +278 @@
       double rho, eleSat, azSat;
       topos(xyzSta(1), xyzSta(2), xyzSta(3), xc(1), xc(2), xc(3), rho, eleSat, azSat);
-      qcObs._eleSet = true;
-      qcObs._azDeg  = azSat  * 180.0/M_PI;
-      qcObs._eleDeg = eleSat * 180.0/M_PI;
+      qcSat._eleSet = true;
+      qcSat._azDeg  = azSat  * 180.0/M_PI;
+      qcSat._eleDeg = eleSat * 180.0/M_PI;
     }
     if (satObs._prn.system() == 'R') {
-      qcObs._slotSet = true;
-      qcObs._slotNum = eph->slotNum();
+      qcSat._slotSet = true;
+      qcSat._slotNum = eph->slotNum();
     }
   }
@@ -294 +294 @@
     const t_frqObs* frqObs = satObs._obs[ii];
 
-    qcObs._qcFrq.push_back(t_qcFrq());
-    t_qcFrq& qcFrq = qcObs._qcFrq.back();
+    qcSat._qcFrq.push_back(t_qcFrq());
+    t_qcFrq& qcFrq = qcSat._qcFrq.back();
 
     qcFrq._rnxType2ch = QString(frqObs->_rnxType2ch.c_str());
@@ -353 +353 @@
               frqObsB->_phaseValid && frqObsB->_codeValid) {
 
-            double f_a = t_CST::freq(fA, qcObs._slotNum);
-            double f_b = t_CST::freq(fB, qcObs._slotNum);
+            double f_a = t_CST::freq(fA, qcSat._slotNum);
+            double f_b = t_CST::freq(fB, qcSat._slotNum);
 
             double L_a = frqObs->_phase  * t_CST::c / f_a;
@@ -410 +410 @@
           t_qcEpo& qcEpo = _qcFile._qcEpo[iEpo];
           if (chunkStart <= qcEpo._epoTime && qcEpo._epoTime < chunkEnd) {
-            if (qcEpo._qcObs.contains(prn)) {
-              t_qcObs& qcObs = qcEpo._qcObs[prn];
-              for (int iFrq = 0; iFrq < qcObs._qcFrq.size(); iFrq++) {
-                t_qcFrq& qcFrq = qcObs._qcFrq[iFrq];
+            if (qcEpo._qcSat.contains(prn)) {
+              t_qcSat& qcSat = qcEpo._qcSat[prn];
+              for (int iFrq = 0; iFrq < qcSat._qcFrq.size(); iFrq++) {
+                t_qcFrq& qcFrq = qcSat._qcFrq[iFrq];
                 if (qcFrq._rnxType2ch == frqType) {
                   frqVec << &qcFrq;
@@ -503 +503 @@
   for (int iEpo = 0; iEpo < _qcFile._qcEpo.size(); iEpo++) {
     t_qcEpo& qcEpo = _qcFile._qcEpo[iEpo];
-    QMapIterator<t_prn, t_qcObs> it(qcEpo._qcObs);
+    QMapIterator<t_prn, t_qcSat> it(qcEpo._qcSat);
     while (it.hasNext()) {
       it.next();
       const t_prn&   prn   = it.key();
-      const t_qcObs& qcObs = it.value();
+      const t_qcSat& qcSat = it.value();
       if ( (prn.system() == 'G' && plotGPS) ||
            (prn.system() == 'R' && plotGlo) ||
            (prn.system() == 'E' && plotGal) ) {
 
-        if (qcObs._eleSet) {
+        if (qcSat._eleSet) {
           QString frqType1;
           QString frqType2;
-          for (int iFrq = 0; iFrq < qcObs._qcFrq.size(); iFrq++) {
-            const t_qcFrq& qcFrq = qcObs._qcFrq[iFrq];
+          for (int iFrq = 0; iFrq < qcSat._qcFrq.size(); iFrq++) {
+            const t_qcFrq& qcFrq = qcSat._qcFrq[iFrq];
             if (qcFrq._rnxType2ch[0] == '1' && frqType1.isEmpty()) {
               frqType1 = qcFrq._rnxType2ch;
@@ -524 +524 @@
             }
             if      (qcFrq._rnxType2ch == frqType1) {
-              (*dataSNR1) << (new t_polarPoint(qcObs._azDeg, 90.0 - qcObs._eleDeg, qcFrq._SNR));
-              (*dataMP1)  << (new t_polarPoint(qcObs._azDeg, 90.0 - qcObs._eleDeg, qcFrq._stdMP));
+              (*dataSNR1) << (new t_polarPoint(qcSat._azDeg, 90.0 - qcSat._eleDeg, qcFrq._SNR));
+              (*dataMP1)  << (new t_polarPoint(qcSat._azDeg, 90.0 - qcSat._eleDeg, qcFrq._stdMP));
             }
             else if (qcFrq._rnxType2ch == frqType2) {
-              (*dataSNR2) << (new t_polarPoint(qcObs._azDeg, 90.0 - qcObs._eleDeg, qcFrq._SNR));
-              (*dataMP2)  << (new t_polarPoint(qcObs._azDeg, 90.0 - qcObs._eleDeg, qcFrq._stdMP));
+              (*dataSNR2) << (new t_polarPoint(qcSat._azDeg, 90.0 - qcSat._eleDeg, qcFrq._SNR));
+              (*dataMP2)  << (new t_polarPoint(qcSat._azDeg, 90.0 - qcSat._eleDeg, qcFrq._stdMP));
             }
           }
@@ -637 +637 @@
     plotData._mjdX24 << mjdX24;
     plotData._PDOP   << qcEpo._PDOP;
-    plotData._numSat << qcEpo._qcObs.size();
-
-    QMapIterator<t_prn, t_qcObs> it(qcEpo._qcObs);
+    plotData._numSat << qcEpo._qcSat.size();
+
+    QMapIterator<t_prn, t_qcSat> it(qcEpo._qcSat);
     while (it.hasNext()) {
       it.next();
       const t_prn&   prn   = it.key();
-      const t_qcObs& qcObs = it.value();
+      const t_qcSat& qcSat = it.value();
       t_plotData&    data  = plotDataMap[prn];
 
       data._mjdX24 << mjdX24;
-      data._eleDeg << qcObs._eleDeg;
+      data._eleDeg << qcSat._eleDeg;
 
       QString frqType1;
       QString frqType2;
-      for (int iFrq = 0; iFrq < qcObs._qcFrq.size(); iFrq++) {
-        const t_qcFrq& qcFrq = qcObs._qcFrq[iFrq];
+      for (int iFrq = 0; iFrq < qcSat._qcFrq.size(); iFrq++) {
+        const t_qcFrq& qcFrq = qcSat._qcFrq[iFrq];
         if (qcFrq._rnxType2ch[0] == '1' && frqType1.isEmpty()) {
           frqType1 = qcFrq._rnxType2ch;
@@ -684 +684 @@
     }
   }
-
 
   if (BNC_CORE->GUIenabled()) {

trunk/BNC/src/rinex/reqcanalyze.h

@@ -86 +86 @@
   };
 
-  class t_qcObs {
+  class t_qcSat {
    public:
-    t_qcObs() {
+    t_qcSat() {
       _slotSet = false;
       _eleSet  = false;
@@ -107 +107 @@
     bncTime              _epoTime;
     double               _PDOP;
-    QMap<t_prn, t_qcObs> _qcObs;
+    QMap<t_prn, t_qcSat> _qcSat;
   };
   
@@ -153 +153 @@
   void   analyzeFile(t_rnxObsFile* obsFile);
 
-  void   updateQcSat(const t_qcObs& qcObs, t_qcSatSum& qcSatSum);
+  void   updateQcSat(const t_qcSat& qcSat, t_qcSatSum& qcSatSum);
 
   void   setQcObs(const bncTime& epoTime, const ColumnVector& xyzSta, 
-                  const t_satObs& satObs, t_qcObs& qcObs);
+                  const t_satObs& satObs, t_qcSat& qcSat);
 
   void   analyzeMultipath();