Changeset 10614 in ntrip for trunk/BNC/src/rinex/reqcedit.cpp

Timestamp:

Feb 28, 2025, 5:53:12 PM (9 hours ago)

Author:

stuerze

Message:

bug fixes egarding RINEX Editing and QC

File:

• trunk/BNC/src/rinex/reqcedit.cpp (modified) (9 diffs)

trunk/BNC/src/rinex/reqcedit.cpp

@@ -75 +75 @@
   _begTime        = bncTime(settings.value("reqcStartDateTime").toString().toLatin1().data());
   _endTime        = bncTime(settings.value("reqcEndDateTime").toString().toLatin1().data());
+
+  _checkEph = false;
 
 }
@@ -171 +173 @@
   while (it.hasNext()) {
     QString fileName = it.next();
-    if (fileName.indexOf('*') != -1 || fileName.indexOf('?') != -1) {
+    if (fileName.indexOf('*') != -1 ||
+        fileName.indexOf('?') != -1) {
       QFileInfo fileInfo(fileName);
       QDir dir = fileInfo.dir();
@@ -316 +319 @@
     t_rnxObsFile* obsFile = _rnxObsFiles[ii];
     if (_log) {
-      *_log << "Processing File: " << obsFile->fileName() << "  start: "
+      *_log << "Input Obs File: " << obsFile->fileName() << "  start: "
             << obsFile->startTime().datestr().c_str() << ' '
             << obsFile->startTime().timestr(0).c_str() << Qt::endl;
@@ -520 +523 @@
 ////////////////////////////////////////////////////////////////////////////
 void t_reqcEdit::readEphemerides(const QStringList& navFileNames,
-                                 QVector<t_eph*>& ephs) {
+                                 QVector<t_eph*>& ephs, QTextStream* log,
+                                 bool checkEph) {
 
   QStringListIterator it(navFileNames);
   while (it.hasNext()) {
     QString fileName = it.next();
-    if (fileName.indexOf('*') != -1 || fileName.indexOf('?') != -1) {
+    if (fileName.indexOf('*') != -1 ||
+        fileName.indexOf('?') != -1) {
       QFileInfo fileInfo(fileName);
       QDir dir = fileInfo.dir();
@@ -532 +537 @@
       while (it.hasNext()) {
         QString filePath = it.next().filePath();
-        appendEphemerides(filePath, ephs);
+        appendEphemerides(filePath, ephs, log, checkEph);
       }
     }
     else {
-      appendEphemerides(fileName, ephs);
+      appendEphemerides(fileName, ephs, log, checkEph);
     }
   }
@@ -565 +570 @@
   while (it.hasNext()) {
     QString fileName = it.next();
+
     t_rnxNavFile rnxNavFile(fileName, t_rnxNavFile::input);
     QStringListIterator itCmnt(rnxNavFile.comments());
@@ -580 +586 @@
   // Read Ephemerides
   // ----------------
-  t_reqcEdit::readEphemerides(_navFileNames, _ephs);
+  t_reqcEdit::readEphemerides(_navFileNames, _ephs, _log, _checkEph);
 
   // Check Satellite Systems
@@ -709 +715 @@
 ////////////////////////////////////////////////////////////////////////////
 void t_reqcEdit::appendEphemerides(const QString& fileName,
-                                   QVector<t_eph*>& ephs) {
+                                   QVector<t_eph*>& ephs, QTextStream* log,
+                                   bool checkEph) {
   t_rnxNavFile rnxNavFile(fileName, t_rnxNavFile::input);
+  unsigned numOK  = 0;
+  unsigned numBad = 0;
+  unsigned numUnhealthy = 0;
+  bncEphUser ephUser(false);
+
+  if (log) {
+    QFileInfo navFi(rnxNavFile.fileName());
+      *log << "Input Nav File     : " << navFi.fileName()      << Qt::endl;
+    if (checkEph) {
+      *log << "RINEX Version      : " << rnxNavFile.version()  << Qt::endl;
+    }
+  }
+
   for (unsigned ii = 0; ii < rnxNavFile.ephs().size(); ii++) {
     t_eph* eph   = rnxNavFile.ephs()[ii];
@@ -724 +744 @@
 
     if (isNew) {
-      if      (eph->system() == t_eph::GPS &&
-               eph->type()   == t_eph::LNAV) {
+
+      if (checkEph) {
+        ephUser.putNewEph(eph, false);
+        if      (eph->checkState() == t_eph::bad) {
+          ++numBad;
+        }
+        else if (eph->checkState() == t_eph::unhealthy) {
+          ++numUnhealthy;
+        }
+        else {
+          ++numOK;
+        }
+        if (eph->checkState() == t_eph::bad) {
+          continue;
+        }
+      }
+
+      if      (eph->system() == t_eph::GPS) {
         ephs.append(new t_ephGPS(*dynamic_cast<t_ephGPS*>(eph)));
       }
-      else if (eph->system() == t_eph::GLONASS &&
-               eph->type()   == t_eph::FDMA_M) {
+      else if (eph->system() == t_eph::GLONASS) {
         ephs.append(new t_ephGlo(*dynamic_cast<t_ephGlo*>(eph)));
       }
-      else if (eph->system() == t_eph::Galileo &&
-               eph->type()   ==  t_eph::INAV) {
+      else if (eph->system() == t_eph::Galileo) {
         ephs.append(new t_ephGal(*dynamic_cast<t_ephGal*>(eph)));
       }
-      else if (eph->system() == t_eph::QZSS &&
-               eph->type()   == t_eph::LNAV) {
+      else if (eph->system() == t_eph::QZSS) {
         ephs.append(new t_ephGPS(*dynamic_cast<t_ephGPS*>(eph)));
       }
-      else if (eph->system() == t_eph::SBAS &&
-               eph->type()   == t_eph::SBASL1) {
+      else if (eph->system() == t_eph::SBAS) {
         ephs.append(new t_ephSBAS(*dynamic_cast<t_ephSBAS*>(eph)));
       }
-      else if (eph->system() == t_eph::BDS &&
-               (eph->type()  == t_eph::D1 ||
-                eph->type()  == t_eph::D2)) {
+      else if (eph->system() == t_eph::BDS) {
         ephs.append(new t_ephBDS(*dynamic_cast<t_ephBDS*>(eph)));
       }
-      else if (eph->system() == t_eph::IRNSS &&
-               eph->type()   == t_eph::LNAV) {
+      else if (eph->system() == t_eph::IRNSS) {
         ephs.append(new t_ephGPS(*dynamic_cast<t_ephGPS*>(eph)));
       }
     }
+  }
+
+  if (log && checkEph) {
+    *log << "Ephemeris Check    : " << numOK        << " OK   "
+                                    << numUnhealthy << " UNHEALTHY   "
+                                    << numBad       << " WRONG\n";
+    if (numBad > 0) {
+      for (unsigned ii = 0; ii < rnxNavFile.ephs().size(); ii++) {
+        t_eph* eph = rnxNavFile.ephs()[ii];
+        QFileInfo navFi(fileName);
+        if (eph->checkState() == t_eph::bad) {
+          if (log) {
+            *log << "                   : "
+                 << QString("WRONG %2:%3\n")
+                         .arg(eph->typeStr(eph->type(), eph->prn(), 99.0))
+                         .arg(eph->rinexDateStr(eph->TOC(), eph->prn(), 99.0)).toLatin1();
+          }
+        }
+      }
+    }
+    *log << Qt::endl;
   }
 }