Changeset 8372 in ntrip

Timestamp:

Jun 13, 2018, 2:37:07 PM (6 years ago)

Author:

stuerze

Message:

minor changes to allow 10 Hz observation data processing (latency checker) and resampling (RINEX files, feed engine)

Location:

trunk/BNC/src

Files:

• GPSDecoder.cpp (modified) (5 diffs)
• bnccaster.cpp (modified) (1 diff)
• bnctime.cpp (modified) (28 diffs)
• latencychecker.cpp (modified) (1 diff)
• rinex/reqcedit.cpp (modified) (1 diff)

trunk/BNC/src/GPSDecoder.cpp

@@ -35 +35 @@
  * Created:    16-Dec-2011
  *
- * Changes:    
+ * Changes:
  *
  * -----------------------------------------------------------------------*/
@@ -45 +45 @@
 #include "bncsettings.h"
 #include "bncrinex.h"
+#include "bncutils.h"
 
 using namespace std;
@@ -63 +64 @@
 //////////////////////////////////////////////////////////////////////////////
 void GPSDecoder::initRinex(const QByteArray& staID, const QUrl& mountPoint,
-                           const QByteArray& latitude, 
+                           const QByteArray& latitude,
                            const QByteArray& longitude, const QByteArray& nmea,
                            const QByteArray& ntripVersion) {
@@ -70 +71 @@
   }
   bncSettings settings;
-  if ( !settings.value("rnxPath").toString().isEmpty() ) { 
-    _rnx = new bncRinex(staID, mountPoint, latitude, longitude, 
+  if ( !settings.value("rnxPath").toString().isEmpty() ) {
+    _rnx = new bncRinex(staID, mountPoint, latitude, longitude,
                         nmea, ntripVersion);
   }
-} 
+}
 
 // Write RINEX Epoch
@@ -80 +81 @@
 void GPSDecoder::dumpRinexEpoch(const t_satObs& obs, const QByteArray& format) {
   if (_rnx) {
-    long iSec    = long(floor(obs._time.gpssec()+0.5));
-    if (_rnx->samplingRate() == 0 || iSec % _rnx->samplingRate() == 0) {
+    int sec = int(nint(obs._time.gpssec()*10));
+    if (_rnx->samplingRate() == 0 || sec % (_rnx->samplingRate()*10) == 0) {
       _rnx->deepCopy(obs);
     }
     _rnx->dumpEpoch(format, obs._time);
   }
-} 
+}
 
 // Set RINEX Reconnect Flag

trunk/BNC/src/bnccaster.cpp

@@ -303 +303 @@
     if (epoTime <= maxTime) {
       const QList<t_satObs>& allObs = itEpo.value();
-      int sec = int(nint(epoTime.gpssec()));
-      if ( (_out || _sockets) && (_samplingRate == 0 || sec % _samplingRate == 0) ) {
-
+      int sec = int(nint(epoTime.gpssec()*10));
+      if ( (_out || _sockets) && (_samplingRate == 0 || sec % (_samplingRate*10) == 0) ) {
         QListIterator<t_satObs> it(allObs);
         bool firstObs = true;

trunk/BNC/src/bnctime.cpp

@@ -17 +17 @@
   this->set(gpsw, gpssec);
 }
-  
+
 // Constructor (from ISO String yyyy-mm-ddThh:mm:ss)
 //////////////////////////////////////////////////////////////////////////////
@@ -24 +24 @@
     QDateTime dt = QDateTime::fromString(isoString.c_str(), Qt::ISODate);
     this->set(dt.date().year(), dt.date().month(), dt.date().day(),
-              dt.time().hour(), dt.time().minute(), 
+              dt.time().hour(), dt.time().minute(),
               dt.time().second() + dt.time().msec()/1000.0);
   }
@@ -32 +32 @@
 }
 
-// 
+//
 //////////////////////////////////////////////////////////////////////////////
 bncTime& bncTime::set(int gpsw, double gpssec) {
   int  deltad;
   int  dow = 0;
-  while ( gpssec >= 86400 ) {
-    gpssec-=86400;
+  while ( gpssec >= 86400.0 ) {
+    gpssec-=86400.0;
     dow++;
   }
-  while ( gpssec <  0 ) {
-    gpssec+=86400;
+  while ( gpssec <  0.0 ) {
+    gpssec+=86400.0;
     dow--;
   }
@@ -59 +59 @@
   gpssec += 14.0;
   gpsw   += 1356.0;
-  while ( gpssec >= 86400 ) {
-    gpssec-=86400;
+  while ( gpssec >= 86400.0 ) {
+    gpssec-=86400.0;
     dow++;
   }
-  while ( gpssec <  0 ) {
-    gpssec+=86400;
+  while ( gpssec <  0.0 ) {
+    gpssec+=86400.0;
     dow--;
   }
@@ -73 +73 @@
 }
 
-// 
+//
 //////////////////////////////////////////////////////////////////////////////
 bncTime &bncTime::set(int msec) {
@@ -85 +85 @@
 }
 
-// 
+//
 //////////////////////////////////////////////////////////////////////////////
 bncTime &bncTime::setTOD(int msec) {
@@ -108 +108 @@
 }
 
-// 
+//
 //////////////////////////////////////////////////////////////////////////////
 bncTime &bncTime::setTk(int msec) {
@@ -122 +122 @@
 }
 
-// 
+//
 //////////////////////////////////////////////////////////////////////////////
 bncTime &bncTime::setBDS(int msec) {
@@ -137 +137 @@
 }
 
-// 
+//
 //////////////////////////////////////////////////////////////////////////////
 bncTime& bncTime::setmjd(double daysec, int mjd) {
   _sec = daysec;
   _mjd = mjd;
-  while ( _sec >= 86400 ) {
-    _sec-=86400;
+  while ( _sec >= 86400.0 ) {
+    _sec-=86400.0;
     _mjd++;
   }
-  while ( _sec <  0 ) {
-    _sec+=86400;
+  while ( _sec <  0.0 ) {
+    _sec+=86400.0;
     _mjd--;
   }
@@ -153 +153 @@
 }
 
-// 
+//
 //////////////////////////////////////////////////////////////////////////////
 bncTime& bncTime::setmjd(double mjddec) {
@@ -161 +161 @@
 }
 
-// 
+//
 //////////////////////////////////////////////////////////////////////////////
 unsigned int bncTime::mjd() const {
   return _mjd;
 }
- 
+
 //
 //////////////////////////////////////////////////////////////////////////////
@@ -182 +182 @@
 }
 
-// 
+//
 //////////////////////////////////////////////////////////////////////////////
 double bncTime::gpssec() const {
@@ -202 +202 @@
 }
 
-// 
+//
 //////////////////////////////////////////////////////////////////////////////
 double bncTime::bdssec() const {
@@ -215 +215 @@
 }
 
-// 
+//
 //////////////////////////////////////////////////////////////////////////////
 bool bncTime::operator!=(const bncTime &time1) const {
@@ -226 +226 @@
 }
 
-// 
+//
 //////////////////////////////////////////////////////////////////////////////
 bool bncTime::operator==(const bncTime &time1) const {
@@ -237 +237 @@
 }
 
-// 
+//
 //////////////////////////////////////////////////////////////////////////////
 bool bncTime::operator>(const bncTime &time1) const {
@@ -248 +248 @@
 }
 
-// 
+//
 //////////////////////////////////////////////////////////////////////////////
 bool bncTime::operator>=(const bncTime &time1) const {
@@ -259 +259 @@
 }
 
-// 
+//
 //////////////////////////////////////////////////////////////////////////////
 bool bncTime::operator<(const bncTime &time1) const {
@@ -270 +270 @@
 }
 
-// 
+//
 //////////////////////////////////////////////////////////////////////////////
 bool bncTime::operator<=(const bncTime &time1) const {
@@ -281 +281 @@
 }
 
-// 
+//
 //////////////////////////////////////////////////////////////////////////////
 bncTime bncTime::operator+(double sec) const {
@@ -290 +290 @@
 }
 
-// 
+//
 //////////////////////////////////////////////////////////////////////////////
 bncTime bncTime::operator-(double sec) const {
@@ -296 +296 @@
 }
 
-// 
+//
 //////////////////////////////////////////////////////////////////////////////
 double bncTime::operator-(const bncTime &time1) const {
@@ -323 +323 @@
 }
 
-// 
+//
 //////////////////////////////////////////////////////////////////////////////
 void bncTime::civil_date (unsigned int& year, unsigned int& month,
@@ -335 +335 @@
 }
 
-// 
-//////////////////////////////////////////////////////////////////////////////
-void bncTime::civil_time(unsigned int &hour, unsigned int &min, 
+//
+//////////////////////////////////////////////////////////////////////////////
+void bncTime::civil_time(unsigned int &hour, unsigned int &min,
                           double &sec) const {
   hour = static_cast<unsigned int>(_sec/3600.0);
@@ -352 +352 @@
 }
 
-// 
+//
 //////////////////////////////////////////////////////////////////////////////
 string bncTime::timestr(unsigned numdec, char sep) const {
@@ -386 +386 @@
 }
 
-// 
+//
 //////////////////////////////////////////////////////////////////////////////
 string bncTime::datestr(char sep) const {
@@ -402 +402 @@
 }
 
-// 
+//
 //////////////////////////////////////////////////////////////////////////////
 bncTime::operator std::string() const {
@@ -408 +408 @@
 }
 
-// 
-//////////////////////////////////////////////////////////////////////////////
-bncTime& bncTime::set(int year, int month, int day, 
+//
+//////////////////////////////////////////////////////////////////////////////
+bncTime& bncTime::set(int year, int month, int day,
                       int hour, int min, double sec) {
   return set(year, month, day, hour*3600 + min*60 + sec);
 }
 
-// 
-//////////////////////////////////////////////////////////////////////////////
-bncTime& bncTime::setBDS(int year, int month, int day, 
+//
+//////////////////////////////////////////////////////////////////////////////
+bncTime& bncTime::setBDS(int year, int month, int day,
                       int hour, int min, double sec) {
   return set(year, month, day, hour*3600 + min*60 + sec+14.0);
 }
 
-// 
+//
 //////////////////////////////////////////////////////////////////////////////
 bncTime& bncTime::set(int year, int month, int day, double daysec) {
   _sec = daysec;
-  
+
   _mjd = (unsigned int)djul(year, month, day);
-  
+
   while ( _sec >= 86400 ) {
     _sec-=86400;

trunk/BNC/src/latencychecker.cpp

@@ -327 +327 @@
       const t_satObs& obs = it.next();
       bool wrongObservationEpoch = checkForWrongObsEpoch(obs._time);
-      l._newSec = static_cast<int>(obs._time.gpssec());
+      l._newSec = static_cast<int>(nint(obs._time.gpssec()*10));
       if (l._newSec > l._oldSec && !wrongObservationEpoch) {
-        if (l._newSec % _miscIntr < l._oldSec % _miscIntr) {
+        if (l._newSec % _miscIntr < l._oldSec % (_miscIntr * 10)) {
           if (l._numLat > 0) {
             if (l._meanDiff > 0.0) {

trunk/BNC/src/rinex/reqcedit.cpp

@@ -337 +337 @@
         }
 
-        if (_samplingRate == 0 ||
-            fmod(round(epo->tt.gpssec()), _samplingRate) == 0) {
+        int sec = int(nint(epo->tt.gpssec()*10));
+        if (_samplingRate == 0 || sec % (_samplingRate*10) == 0) {
           applyLLI(obsFile, epo);
           outObsFile.writeEpoch(epo);