Changeset 6322 in ntrip for trunk/BNC/src/rinex/reqcanalyze.cpp

Timestamp:

Nov 12, 2014, 5:20:13 PM (9 years ago)

Author:

stuerze

Message:

systems J,S,C added, parameter type of plot function changed

File:

• trunk/BNC/src/rinex/reqcanalyze.cpp (modified) (12 diffs)

trunk/BNC/src/rinex/reqcanalyze.cpp

@@ -69 +69 @@
   _navFileNames = settings.value("reqcNavFile").toString().split(",", QString::SkipEmptyParts);
 
-  connect(this, SIGNAL(dspSkyPlot(const QString&, const QByteArray&, QVector<t_polarPoint*>*, 
-                                  const QByteArray&, QVector<t_polarPoint*>*,
+  connect(this, SIGNAL(dspSkyPlot(const QString&, const QString&, QVector<t_polarPoint*>*,
+                                  const QString&, QVector<t_polarPoint*>*,
                                   const QByteArray&, double)), 
-          this, SLOT(slotDspSkyPlot(const QString&, const QByteArray&, QVector<t_polarPoint*>*, 
-                                    const QByteArray&, QVector<t_polarPoint*>*,
+          this, SLOT(slotDspSkyPlot(const QString&, const QString&, QVector<t_polarPoint*>*,
+                                    const QString&, QVector<t_polarPoint*>*,
                                     const QByteArray&, double)));
 
@@ -297 +297 @@
   // ---------------------------
   for (unsigned ii = 0; ii < satObs._obs.size(); ii++) {
-
     const t_frqObs* frqObs = satObs._obs[ii];
 
@@ -308 +307 @@
     qcFrq._phaseValid = frqObs->_phaseValid;
     qcFrq._codeValid  = frqObs->_codeValid;
-
     // Check Gaps
     // ----------
@@ -355 +353 @@
         }
       }
-
+      else if (satObs._prn.system() == 'J') {
+        if      (frqObs->_rnxType2ch[0] == '1') {
+          fA = t_frequency::J1;
+          fB = t_frequency::J2;
+        }
+        else if (frqObs->_rnxType2ch[0] == '2') {
+          fA = t_frequency::J2;
+          fB = t_frequency::J1;
+        }
+      }
+      else if (satObs._prn.system() == 'S') {
+        if      (frqObs->_rnxType2ch[0] == '1') {
+          fA = t_frequency::S1;
+          fB = t_frequency::S5;
+        }
+        else if (frqObs->_rnxType2ch[0] == '5') {
+          fA = t_frequency::S5;
+          fB = t_frequency::S1;
+        }
+      }
+      else if (satObs._prn.system() == 'C') {
+        if      (frqObs->_rnxType2ch[0] == '1') {
+          fA = t_frequency::C1;
+          fB = t_frequency::C7;
+        }
+        else if (frqObs->_rnxType2ch[0] == '7') {
+          fA = t_frequency::C7;
+          fB = t_frequency::C1;
+        }
+      }
       if (fA != t_frequency::dummy && fB != t_frequency::dummy) {
         double f_a = t_CST::freq(fA, qcSat._slotNum);
@@ -409 +436 @@
       const QString& frqType  = itFrq.key();
       t_qcFrqSum&    qcFrqSum = itFrq.value();
-
 
       // Loop over all Chunks of Data
@@ -494 +520 @@
   bncSettings settings;
   QString reqSkyPlotSystems = settings.value("reqcSkyPlotSystems").toString();
-  bool plotGPS = false;
-  bool plotGlo = false;
-  bool plotGal = false;
+  bool plotGPS  = false;
+  bool plotGlo  = false;
+  bool plotGal  = false;
+  bool plotQZSS = false;
+  bool plotSBAS = false;
+  bool plotBDS  = false;
+  QString mp1Title = reqSkyPlotSystems + ": MP";
+  QString mp2Title = reqSkyPlotSystems + ": MP";
+  QString snr1Title = reqSkyPlotSystems + ": SNR";
+  QString snr2Title = reqSkyPlotSystems + ": SNR";
+  char freq1 = '1';
+  char freq2 = '2';
   if      (reqSkyPlotSystems == "GPS") {
     plotGPS = true;
@@ -505 +540 @@
   else if (reqSkyPlotSystems == "Galileo") {
     plotGal = true;
+    freq2 = '5';
+  }
+  else if (reqSkyPlotSystems == "QZSS") {
+    plotQZSS = true;
+  }
+  else if (reqSkyPlotSystems == "SBAS") {
+    plotSBAS = true;
+    freq2 = '5';
+  }
+  else if (reqSkyPlotSystems == "BDS") {
+    plotBDS = true;
+    freq2 = '7';
   }
   else if (reqSkyPlotSystems == "ALL") {
@@ -510 +557 @@
     plotGlo = true;
     plotGal = true;
+    plotQZSS = true;
+    plotSBAS = true;
+    plotBDS = true;
   }
   else {
           return;
   }
+  mp1Title += freq1;
+  mp2Title += freq2;
+  snr1Title += freq1;
+  snr2Title += freq2;
 
   QVector<t_polarPoint*>* dataMP1  = new QVector<t_polarPoint*>;
@@ -531 +585 @@
       if ( (prn.system() == 'G' && plotGPS) ||
            (prn.system() == 'R' && plotGlo) ||
-           (prn.system() == 'E' && plotGal) ) {
+           (prn.system() == 'E' && plotGal) ||
+           (prn.system() == 'J' && plotQZSS) ||
+           (prn.system() == 'S' && plotSBAS) ||
+           (prn.system() == 'C' && plotBDS) ) {
 
         if (qcSat._eleSet) {
@@ -538 +595 @@
           for (int iFrq = 0; iFrq < qcSat._qcFrq.size(); iFrq++) {
             const t_qcFrq& qcFrq = qcSat._qcFrq[iFrq];
-            if (qcFrq._rnxType2ch[0] == '1' && frqType1.isEmpty()) {
+            if (qcFrq._rnxType2ch[0] == freq1 && frqType1.isEmpty()) {
               frqType1 = qcFrq._rnxType2ch;
             }
-            if (qcFrq._rnxType2ch[0] == '2' && frqType2.isEmpty()) {
+            if (qcFrq._rnxType2ch[0] == freq2 && frqType2.isEmpty()) {
               frqType2 = qcFrq._rnxType2ch;
             }
@@ -563 +620 @@
     QFileInfo  fileInfo(obsFile->fileName());
     QByteArray title = fileInfo.fileName().toAscii();
-    emit dspSkyPlot(obsFile->fileName(), "MP1",  dataMP1,  "MP2",  dataMP2,  "Meters",  2.0);
-    emit dspSkyPlot(obsFile->fileName(), "SNR1", dataSNR1, "SNR2", dataSNR2, "dbHz",   54.0);
+    emit dspSkyPlot(obsFile->fileName(), mp1Title,  dataMP1,  mp2Title,  dataMP2,  "Meters",  2.0);
+    double mean = 0.0;
+    for (int ii = 0; ii < dataSNR1->size(); ii++) {
+      const t_polarPoint* point = dataSNR1->at(ii);
+      mean += point->_value;
+    }
+    mean /= dataSNR1->size();
+    double max = (mean > 9.0) ? 54.0 : 9.0;
+    QByteArray str = (mean > 9.0) ? "dbHz" : "";
+    emit dspSkyPlot(obsFile->fileName(), snr1Title, dataSNR1, snr2Title, dataSNR2, str,  max);
     emit dspAvailPlot(obsFile->fileName(), title);
   }
@@ -589 +654 @@
 //
 ////////////////////////////////////////////////////////////////////////////
-void t_reqcAnalyze::slotDspSkyPlot(const QString& fileName, const QByteArray& title1,
-                                   QVector<t_polarPoint*>* data1, const QByteArray& title2,
+void t_reqcAnalyze::slotDspSkyPlot(const QString& fileName, const QString& title1,
+                                   QVector<t_polarPoint*>* data1, const QString& title2,
                                    QVector<t_polarPoint*>* data2, const QByteArray& scaleTitle,
                                    double maxValue) {