Changeset 8483 in ntrip for trunk/BNC/src/satObs.cpp

Timestamp:

Sep 17, 2018, 2:53:30 PM (6 years ago)

Author:

stuerze

Message:

SSR parameter clock rate, clock drift and URA are added within RTNET format

File:

• trunk/BNC/src/satObs.cpp (modified) (17 diffs)

trunk/BNC/src/satObs.cpp

@@ -18 +18 @@
 }
 
-// 
+//
 ////////////////////////////////////////////////////////////////////////////
 void t_clkCorr::writeEpoch(ostream* out, const QList<t_clkCorr>& corrList) {
@@ -43 +43 @@
 }
 
-// 
+//
 ////////////////////////////////////////////////////////////////////////////
 void t_clkCorr::readEpoch(const string& epoLine, istream& inStream, QList<t_clkCorr>& corrList) {
@@ -62 +62 @@
     getline(inStream, line);
     istringstream in(line.c_str());
-    
+
     in >> corr._prn >> corr._iod >> corr._dClk >> corr._dotDClk >> corr._dotDotDClk;
     if (corr._prn.system() == 'E') {
@@ -85 +85 @@
 }
 
-// 
+//
 ////////////////////////////////////////////////////////////////////////////
 void t_orbCorr::writeEpoch(ostream* out, const QList<t_orbCorr>& corrList) {
@@ -113 +113 @@
 }
 
-// 
+//
 ////////////////////////////////////////////////////////////////////////////
 void t_orbCorr::readEpoch(const string& epoLine, istream& inStream, QList<t_orbCorr>& corrList) {
@@ -132 +132 @@
     getline(inStream, line);
     istringstream in(line.c_str());
-    
+
     in >> corr._prn      >> corr._iod
-       >> corr._xr[0]    >> corr._xr[1]    >> corr._xr[2]   
+       >> corr._xr[0]    >> corr._xr[1]    >> corr._xr[2]
        >> corr._dotXr[0] >> corr._dotXr[1] >> corr._dotXr[2];
 
@@ -144 +144 @@
 }
 
-// 
+// Constructor
+////////////////////////////////////////////////////////////////////////////
+t_URA::t_URA() {
+  _updateInt  = 0;
+  _iod        = 0;
+  _ura        = 0.0;
+}
+
+//
+////////////////////////////////////////////////////////////////////////////
+void t_URA::writeEpoch(ostream* out, const QList<t_URA>& corrList) {
+  if (!out || corrList.size() == 0) {
+    return;
+  }
+  out->setf(ios::fixed);
+  bncTime epoTime;
+  QListIterator<t_URA> it(corrList);
+  while (it.hasNext()) {
+    const t_URA& corr = it.next();
+    if (!epoTime.valid()) {
+      epoTime = corr._time;
+      *out << "> URA " << epoTime.datestr(' ') << ' ' << epoTime.timestr(1,' ') << " "
+          <<  corr._updateInt <<  " "
+           << corrList.size() << ' ' << corr._staID << endl;
+    }
+    *out << corr._prn.toString() << ' ' << setw(11) << corr._iod << ' '
+         << setw(10) << setprecision(4) << corr._ura << endl;
+  }
+  out->flush();
+}
+
+//
+////////////////////////////////////////////////////////////////////////////
+void t_URA::readEpoch(const string& epoLine, istream& inStream, QList<t_URA>& corrList) {
+  bncTime      epoTime;
+  unsigned int updateInt;
+  int          numCorr;
+  string       staID;
+  if (t_corrSSR::readEpoLine(epoLine, epoTime, updateInt, numCorr, staID) != t_corrSSR::URA) {
+    return;
+  }
+  for (int ii = 0; ii < numCorr; ii++) {
+    t_URA corr;
+    corr._time      = epoTime;
+    corr._updateInt = updateInt;
+    corr._staID     = staID;
+
+    string line;
+    getline(inStream, line);
+    istringstream in(line.c_str());
+
+    in >> corr._prn >> corr._iod >> corr._ura;
+
+    corrList.push_back(corr);
+  }
+}
+
+//
 ////////////////////////////////////////////////////////////////////////////
 void t_satCodeBias::writeEpoch(ostream* out, const QList<t_satCodeBias>& biasList) {
@@ -172 +229 @@
 }
 
-// 
+//
 ////////////////////////////////////////////////////////////////////////////
 void t_satCodeBias::readEpoch(const string& epoLine, istream& inStream, QList<t_satCodeBias>& biasList) {
@@ -191 +248 @@
     getline(inStream, line);
     istringstream in(line.c_str());
-    
+
     int numBias;
     in >> satCodeBias._prn >> numBias;
@@ -207 +264 @@
 }
 
-// 
+//
 ////////////////////////////////////////////////////////////////////////////
 void t_satPhaseBias::writeEpoch(ostream* out, const QList<t_satPhaseBias>& biasList) {
@@ -242 +299 @@
   out->flush();
 }
-  
-// 
+
+//
 ////////////////////////////////////////////////////////////////////////////
 void t_satPhaseBias::readEpoch(const string& epoLine, istream& inStream, QList<t_satPhaseBias>& biasList) {
@@ -290 +347 @@
 }
 
-// 
+//
 ////////////////////////////////////////////////////////////////////////////
 void t_vTec::write(ostream* out, const t_vTec& vTec) {
@@ -306 +363 @@
          << setw(2)  << layer._C.Nrows()-1 << ' '
          << setw(2)  << layer._C.Ncols()-1 << ' '
-         << setw(10) << setprecision(1) << layer._height << endl  
-         << setw(10) << setprecision(4) << layer._C 
+         << setw(10) << setprecision(1) << layer._height << endl
+         << setw(10) << setprecision(4) << layer._C
          << setw(10) << setprecision(4) << layer._S;
   }
@@ -313 +370 @@
 }
 
-// 
+//
 ////////////////////////////////////////////////////////////////////////////
 void t_vTec::read(const string& epoLine, istream& inStream, t_vTec& vTec) {
@@ -357 +414 @@
 }
 
-// 
-////////////////////////////////////////////////////////////////////////////
-t_corrSSR::e_type t_corrSSR::readEpoLine(const string& line, bncTime& epoTime, 
+//
+////////////////////////////////////////////////////////////////////////////
+t_corrSSR::e_type t_corrSSR::readEpoLine(const string& line, bncTime& epoTime,
                                          unsigned int& updateInt, int& numEntries,
                                          string& staID) {
@@ -370 +427 @@
   double sec;
 
-  inLine >> epoChar >> typeString 
+  inLine >> epoChar >> typeString
          >> year >> month >> day >> hour >> min >> sec >> updateInt >> numEntries >> staID;
 
@@ -390 +447 @@
       return vTec;
     }
+    else if (typeString == "URA") {
+      return URA;
+    }
   }