Changeset 6755 in ntrip

Timestamp:

Apr 1, 2015, 3:03:06 PM (10 years ago)

Author:

stuerze

Message:

different BDS ephemeris sources ( RTCM stream / RINEX file ) are considered now for RINEX file writing

Location:

trunk/BNC/src

Files:

• ephemeris.cpp (modified) (12 diffs)
• ephemeris.h (modified) (2 diffs)

trunk/BNC/src/ephemeris.cpp

@@ -1336 +1336 @@
   // ------------
   int fieldLen = 19;
+  _URAI = -1; // undefined in RINEX format
 
   int pos[4];
@@ -1342 +1343 @@
   pos[2] = pos[1] + fieldLen;
   pos[3] = pos[2] + fieldLen;
-
-  double TOTs;
-  double TOEs;
-  double TOEw;
 
   // Read eight lines
@@ -1407 +1404 @@
 
     else if ( iLine == 3 ) {
-      if ( readDbl(line, pos[0], fieldLen, TOEs   )  ||
+      if ( readDbl(line, pos[0], fieldLen, _TOEs   )  ||
            readDbl(line, pos[1], fieldLen, _Cic   )  ||
            readDbl(line, pos[2], fieldLen, _OMEGA0)  ||
@@ -1428 +1425 @@
     else if ( iLine == 5 ) {
       if ( readDbl(line, pos[0], fieldLen, _IDOT    ) ||
-           readDbl(line, pos[2], fieldLen, TOEw)    ) {
+           readDbl(line, pos[2], fieldLen, _TOEw)    ) {
         _checkState = bad;
         return;
@@ -1436 +1433 @@
     else if ( iLine == 6 ) {
       double SatH1;
-      if ( readDbl(line, pos[1], fieldLen, SatH1) ||
+      if ( readDbl(line, pos[0], fieldLen, _URA ) ||
+           readDbl(line, pos[1], fieldLen, SatH1) ||
            readDbl(line, pos[2], fieldLen, _TGD1) ||
            readDbl(line, pos[3], fieldLen, _TGD2) ) {
@@ -1447 +1445 @@
     else if ( iLine == 7 ) {
       double aodc;
-      if ( readDbl(line, pos[0], fieldLen, TOTs) ||
+      if ( readDbl(line, pos[0], fieldLen, _TOTs) ||
            readDbl(line, pos[1], fieldLen, aodc) ) {
         _checkState = bad;
         return;
       }
-      if (TOTs == 0.9999e9) {  // 0.9999e9 means not known (RINEX standard)
-        TOTs = TOEs;
+      if (_TOTs == 0.9999e9) {  // 0.9999e9 means not known (RINEX standard)
+        _TOTs = _TOEs;
       }
       _AODC = int(aodc);
@@ -1459 +1457 @@
   }
 
-  TOEw += 1356;  // BDT -> GPS week number
-  _TOE_bdt.set(int(TOEw), TOEs);
+  _TOEw += 1356;  // BDT -> GPS week number
+  _TOE_bdt.set(int(_TOEw), _TOEs);
 
   // GPS->BDT
@@ -1475 +1473 @@
 ////////////////////////////////////////////////////////////////////////////
 void t_ephBDS::set(const bdsephemeris* ee) {
+
+  // RINEX File entries
+  // --------------------
+  _TOTs = 0.0;
+  _TOEs = 0.0;
+  _TOEw = 0.0;
 
   _receptDateTime = currentDateAndTimeGPS();
@@ -1677 +1681 @@
     .arg(_sqrt_A, 19, 'e', 12);
 
-  out << QString(fmt)
-    .arg(_TOE_bdt.gpssec(), 19, 'e', 12)
-    .arg(_Cic,              19, 'e', 12)
-    .arg(_OMEGA0,           19, 'e', 12)
-    .arg(_Cis,              19, 'e', 12);
+  double toes = _TOEs;
+  if (!toes) { // RTCM stream input
+    toes = _TOE_bdt.gpssec();
+  }
+  out << QString(fmt)
+    .arg(toes,   19, 'e', 12)
+    .arg(_Cic,    19, 'e', 12)
+    .arg(_OMEGA0, 19, 'e', 12)
+    .arg(_Cis,    19, 'e', 12);
 
   out << QString(fmt)
@@ -1689 +1697 @@
     .arg(_OMEGADOT, 19, 'e', 12);
 
-  out << QString(fmt)
-    .arg(_IDOT,                            19, 'e', 12)
-    .arg(0.0,                              19, 'e', 12)
-    .arg(double(_TOE_bdt.gpsw() - 1356.0), 19, 'e', 12)
-    .arg(0.0,                              19, 'e', 12);
-
-  double ura = 0.0;
+  double toew = _TOEw;
+  if (!toew) { // RTCM stream input
+    toew = double(_TOE_bdt.gpsw() - 1356.0);
+  }
+  out << QString(fmt)
+    .arg(_IDOT, 19, 'e', 12)
+    .arg(0.0,   19, 'e', 12)
+    .arg(toew,  19, 'e', 12)
+    .arg(0.0,   19, 'e', 12);
+
+  double ura = _URA; // RINEX file input
   if ((_URAI <  6) && (_URAI >= 0)) {
     ura = ceil(10.0 * pow(2.0, ((double)_URAI/2.0) + 1.0)) / 10.0;
@@ -1702 +1714 @@
     ura = ceil(10.0 * pow(2.0, ((double)_URAI/2.0)      )) / 10.0;
   }
-
   out << QString(fmt)
     .arg(ura,            19, 'e', 12)
@@ -1709 +1720 @@
     .arg(_TGD2,          19, 'e', 12);
 
-  out << QString(fmt)
-    .arg(_TOE_bdt.gpssec(), 19, 'e', 12)
-    .arg(double(_AODC),     19, 'e', 12)
-    .arg("",                19, QChar(' '))
-    .arg("",                19, QChar(' '));
+  double tots = _TOTs;
+  if (!tots) { // RTCM stream input
+    tots = _TOE_bdt.gpssec();
+  }
+  out << QString(fmt)
+    .arg(tots,         19, 'e', 12)
+    .arg(double(_AODC), 19, 'e', 12)
+    .arg("",            19, QChar(' '))
+    .arg("",            19, QChar(' '));
   return rnxStr;
 }

trunk/BNC/src/ephemeris.h

@@ -260 +260 @@
   int     _AODE;
   int     _AODC;
-  int     _URAI;             //  [0..15]
+  int     _URAI;             //  [0..15] from RTCM stream
+  double  _URA;              //  [m]     from RINEX file
   double  _clock_bias;       //  [s]    
   double  _clock_drift;      //  [s/s]  
@@ -282 +283 @@
   double  _TGD2;             //  [s]    
   int     _SatH1;            // 
+  double  _TOTs;             //  [s] of BDT week; RINEX file entry
+  double  _TOEs;             //  [s] of BDT week; RINEX file entry
+  double  _TOEw;             //  [-] BDT week; RINEX file entry
 };