Changeset 6981 in ntrip for trunk

Timestamp:

Jun 29, 2015, 2:35:47 PM (9 years ago)

Author:

stuerze

Message:

changes to prevent the upload of ephemeris, which are detected as erroneous; further ephemris checks regarding their age are added

Location:

trunk/BNC/src

Files:

• bncephuser.cpp (modified) (3 diffs)
• bncephuser.h (modified) (1 diff)

trunk/BNC/src/bncephuser.cpp

@@ -154 +154 @@
   const t_eph* ephOld = ephLast(prn);
 
-  if (ephOld == 0 || newEph->isNewerThan(ephOld)) {
+  if (ephOld && ephOld->checkState() == t_eph::bad) {
+    ephOld = 0;
+  }
+  
+  if ((ephOld == 0 || newEph->isNewerThan(ephOld)) &&
+      eph->checkState() != t_eph::bad) {
     deque<t_eph*>& qq = _eph[prn];
     qq.push_back(newEph);
@@ -195 +200 @@
     return;
   }
+
+  // Check whether the epoch is too far away the current time
+  // --------------------------------------------------------
+  bncTime   toc = eph->TOC();
+  QDateTime now = currentDateAndTimeGPS();
+  bncTime currentTime(now.toString(Qt::ISODate).toStdString());
+  double timeDiff = fabs(toc - currentTime);
+
+  if      (eph->type() == t_eph::GPS  || t_eph::Galileo) {
+    if (timeDiff > 4*3600) { // data sets are valid for 4 hours
+      eph->setCheckState(t_eph::bad);
+      return;
+    }
+  }
+  else if (eph->type() == t_eph::GLONASS) {
+    if (timeDiff > 1*3600) { // updated every 30 minutes
+      eph->setCheckState(t_eph::bad);
+      return;
+    }
+  }
+  else if (eph->type() == t_eph::QZSS) {
+    if (timeDiff > 1*3600) {
+      eph->setCheckState(t_eph::bad);
+      return;
+    }
+  }
+  else if (eph->type() == t_eph::SBAS) {
+    if (timeDiff > 300) { // maximum update interval
+      eph->setCheckState(t_eph::bad);
+      return;
+    }
+  }
+  else if (eph->type() == t_eph::BDS) {
+    if (timeDiff > 6*3600) { // updates 1 (GEO) up to 6 hours
+      eph->setCheckState(t_eph::bad);
+      return;
+    }
+  }
+
 
   // Check consistency with older ephemerides
@@ -209 +253 @@
     double diff = (xc.Rows(1,3) - xcL.Rows(1,3)).norm_Frobenius();
 
+    // some lines to allow update of ephemeris data sets after outage
+    if      ((eph->type() == t_eph::GPS ||
+              eph->type() == t_eph::Galileo)  && dt > 4*3600) {
+      ephL->setCheckState(t_eph::bad);
+      return;
+    }
+    else if ((eph->type() == t_eph::GLONASS ||
+              eph->type() == t_eph::QZSS)     && dt > 2*3600) {
+      ephL->setCheckState(t_eph::bad);
+      return;
+    }
+    else if  (eph->type() == t_eph::SBAS      && dt > 1*3600)   {
+      ephL->setCheckState(t_eph::bad);
+      return;
+    }
+    else if  (eph->type() == t_eph::BDS       && dt > 6*3600)   {
+      ephL->setCheckState(t_eph::bad);
+      return;
+    }
+
     if (diff < MAXDIFF) {
-      if(dt != 0.0) {
+      if (dt != 0.0) {
         eph->setCheckState(t_eph::ok);
         ephL->setCheckState(t_eph::ok);

trunk/BNC/src/bncephuser.h

@@ -32 +32 @@
 #include "bncconst.h"
 #include "bnctime.h"
+#include "bncutils.h"
 #include "ephemeris.h"