Changeset 9249 in ntrip

Timestamp:

Nov 12, 2020, 11:49:23 AM (3 years ago)

Author:

stuerze

Message:

another check is added, to prevent the usage of not updated nav data sets during ssr upload

Location:

branches/BNC_2.12/src

Files:

• bncephuser.cpp (modified) (4 diffs)
• bncutils.cpp (modified) (2 diffs)
• bncutils.h (modified) (2 diffs)
• upload/bncrtnetuploadcaster.cpp (modified) (4 diffs)

branches/BNC_2.12/src/bncephuser.cpp

@@ -166 +166 @@
 
   if (eph->checkState() != t_eph::bad &&
-       eph->checkState() != t_eph::outdated) {
+      eph->checkState() != t_eph::outdated) {
     deque<t_eph*>& qq = _eph[prn];
     qq.push_back(newEph);
@@ -243 +243 @@
     }
     // maximum update interval: 300 sec
-    else if (eph->type() == t_eph::SBAS    && (dt >     600 || dt <       -300)) {
+    else if (eph->type() == t_eph::SBAS    && (dt >     600 || dt <     -300.0)) {
       eph->setCheckState(t_eph::outdated);
       return;
@@ -253 +253 @@
     }
     // update interval: up to 24 hours
-    else if (eph->type() == t_eph::IRNSS   && fabs(dt > 24*3600)) {
+    else if (eph->type() == t_eph::IRNSS   && fabs(dt > 24*3600.0)) {
       eph->setCheckState(t_eph::outdated);
       return;
@@ -279 +279 @@
     // some lines to allow update of ephemeris data sets after an outage
     // update interval: 2h,
-    if      (eph->type() == t_eph::GPS     && dt >  2*3600) {
+    if      (eph->type() == t_eph::GPS     && dt >  2*3600.0) {
       ephL->setCheckState(t_eph::outdated);
       return;
     }
     // update interval: 3h,
-    else if (eph->type() == t_eph::Galileo && dt >  3*3600) {
+    else if (eph->type() == t_eph::Galileo && dt >  3*3600.0) {
       ephL->setCheckState(t_eph::outdated);
       return;
     }
     // updated every 30 minutes
-    else if (eph->type() == t_eph::GLONASS && dt >    1800) {
+    else if (eph->type() == t_eph::GLONASS && dt >    1800.0) {
       ephL->setCheckState(t_eph::outdated);
       return;
     }
     // updated every ?
-    else if (eph->type() == t_eph::QZSS    && dt >  2*3600) {
+    else if (eph->type() == t_eph::QZSS    && dt >  2*3600.0) {
       ephL->setCheckState(t_eph::outdated);
       return;
     }
     // maximum update interval: 300 sec
-    else if  (eph->type() == t_eph::SBAS   && dt >     300) {
+    else if  (eph->type() == t_eph::SBAS   && dt >     300.0) {
       ephL->setCheckState(t_eph::outdated);
       return;
     }
     // updates 1h (GEO) up to 6 hours non-GEO
-    else if  (eph->type() == t_eph::BDS    && dt >  6*3600) {
+    else if  (eph->type() == t_eph::BDS    && dt >  6*3600.0) {
       ephL->setCheckState(t_eph::outdated);
       return;
     }
     // update interval: up to 24 hours
-    else if  (eph->type() == t_eph::IRNSS  && dt > 24*3600) {
+    else if  (eph->type() == t_eph::IRNSS  && dt > 24*3600.0) {
       ephL->setCheckState(t_eph::outdated);
       return;

branches/BNC_2.12/src/bncutils.cpp

@@ -255 +255 @@
   currentGPSWeeks(week, sec);
   bncTime currTime(week, sec);
+
   if (((currTime - obsTime) < 0.0) ||
       (fabs(currTime - obsTime) > maxDt)) {
@@ -261 +262 @@
   return false;
 }
+
+//
+////////////////////////////////////////////////////////////////////////////
+bool outDatedBcep(const t_eph *eph) {
+  bncTime   toc = eph->TOC();
+  QDateTime now = currentDateAndTimeGPS();
+  bncTime currentTime(now.toString(Qt::ISODate).toStdString());
+  double dt = currentTime - toc;
+
+  // update interval: 2h, data sets are valid for 4 hours
+  if      (eph->type() == t_eph::GPS     && dt > 4*3600.0) {
+    return true;
+  }
+  // update interval: 3h, data sets are valid for 4 hours
+  else if (eph->type() == t_eph::Galileo && dt > 4*3600.0) {
+    return true;
+  }
+  // updated every 30 minutes
+  else if (eph->type() == t_eph::GLONASS && dt > 1*3600.0) {
+    return true;
+  }
+  // orbit parameters are valid for 7200 seconds (minimum)
+  else if (eph->type() == t_eph::QZSS    && dt > 2*3600.0) {
+    return true;
+  }
+  // maximum update interval: 300 sec
+  else if (eph->type() == t_eph::SBAS    && dt >    600.0) {
+    return true;
+  }
+  // updates 1h (GEO) up to 6 hours non-GEO
+  else if (eph->type() == t_eph::BDS     && dt > 6*3600.0) {
+    return true;
+  }
+  // update interval: up to 24 hours
+  else if (eph->type() == t_eph::IRNSS   && dt > 24*3600.0) {
+    return true;
+  }
+  return false;
+}
+
 //
 ////////////////////////////////////////////////////////////////////////////

branches/BNC_2.12/src/bncutils.h

@@ -65 +65 @@
 bool         checkForWrongObsEpoch(bncTime obsEpoch);
 
+bool         outDatedBcep(const t_eph *eph);
+
 QByteArray   ggaString(const QByteArray& latitude, const QByteArray& longitude,
                        const QByteArray& height, const QString& ggaType);
@@ -149 +151 @@
 */
 double       lti2sec(int type, int lti);
-
 
 // CRC24Q checksum calculation function (only full bytes supported).

branches/BNC_2.12/src/upload/bncrtnetuploadcaster.cpp

@@ -22 +22 @@
 #include "bncsp3.h"
 #include "gnss.h"
+#include "bncutils.h"
 
 using namespace std;
@@ -384 +385 @@
   for (int ii = 1; ii < lines.size(); ii++) {
     QString key;  // prn or key VTEC, IND (phase bias indicators)
-    double rtnUra = 0.0;    // [m]
+    double rtnUra = 0.0; // [m]
     ColumnVector rtnAPC; rtnAPC.ReSize(3); rtnAPC = 0.0;          // [m, m, m]
     ColumnVector rtnVel; rtnVel.ReSize(3); rtnVel = 0.0;          // [m/s, m/s, m/s]
@@ -467 +468 @@
     }
 
-    if (eph  &&
-        eph->checkState() != t_eph::bad &&
+    if (eph  &&  !outDatedBcep(eph)           &&  // detected from storage because of no update
+        eph->checkState() != t_eph::bad       &&
         eph->checkState() != t_eph::unhealthy &&
-        eph->checkState() != t_eph::outdated) {
+        eph->checkState() != t_eph::outdated) {  // detected during reception (bncephuser)
       QMap<QString, double> codeBiases;
       QList<phaseBiasSignal> phaseBiasList;
@@ -926 +927 @@
 
   _outBuffer += hlpBufferCo + hlpBufferBias + hlpBufferPhaseBias
-      + hlpBufferVtec + '\0';;
+      + hlpBufferVtec + '\0';
 }