Changeset 849 in ntrip

Timestamp:

Apr 25, 2008, 11:44:33 AM (16 years ago)

Author:

mervart

Message:

* empty log message *

Location:

trunk/BNS

Files:

• bns.cpp (modified) (1 diff)
• bnssp3.cpp (modified) (3 diffs)
• bnssp3.h (modified) (1 diff)

trunk/BNS/bns.cpp

@@ -116 +116 @@
     delete it.value();
   }
+  delete _rnx;
+  delete _sp3;
 }
 

trunk/BNS/bnssp3.cpp

@@ -19 +19 @@
 
 #include "bnssp3.h"
+#include "bnsutils.h"
 
 using namespace std;
@@ -25 +26 @@
 ////////////////////////////////////////////////////////////////////////////
 bnsSP3::bnsSP3() {
+  QSettings settings;
+
+  _headerWritten = false;
+  _ID4           = "BNS_";
+  _ext           = ".SP3";
+  _samplingRate  = settings.value("sp3Sampl").toInt();
+  _intr          = settings.value("rnxIntr").toString();
+  _path          = settings.value("rnxPath").toString();
+  expandEnvVar(_path);
+  if ( _path.length() > 0 && _path[_path.length()-1] != QDir::separator() ) {
+    _path += QDir::separator();
+  }
 }
 
@@ -30 +43 @@
 ////////////////////////////////////////////////////////////////////////////
 bnsSP3::~bnsSP3() {
+  _out.close();
 }
 
-// 
+// Close the Old RINEX File
+////////////////////////////////////////////////////////////////////////////
+void bnsSP3::closeFile() {
+  _out.close();
+}
+
+// Next File Epoch (static)
+////////////////////////////////////////////////////////////////////////////
+QString bnsSP3::nextEpochStr(const QDateTime& datTim, 
+                             const QString& intStr, QDateTime* nextEpoch) {
+
+  QString epoStr;
+
+  QTime nextTime;
+  QDate nextDate;
+
+  int indHlp = intStr.indexOf("min");
+
+  if ( indHlp != -1) {
+    int step = intStr.left(indHlp-1).toInt();
+    char ch = 'A' + datTim.time().hour();
+    epoStr = ch;
+    if (datTim.time().minute() >= 60-step) {
+      epoStr += QString("%1").arg(60-step, 2, 10, QChar('0'));
+      if (datTim.time().hour() < 23) {
+        nextTime.setHMS(datTim.time().hour() + 1 , 0, 0);
+        nextDate = datTim.date();
+      }
+      else {
+        nextTime.setHMS(0, 0, 0);
+        nextDate = datTim.date().addDays(1);
+      }
+    }
+    else {
+      for (int limit = step; limit <= 60-step; limit += step) {
+        if (datTim.time().minute() < limit) {
+          epoStr += QString("%1").arg(limit-step, 2, 10, QChar('0'));
+          nextTime.setHMS(datTim.time().hour(), limit, 0);
+          nextDate = datTim.date();
+          break;
+        }
+      }
+    }
+  }
+  else if (intStr == "1 hour") {
+    char ch = 'A' + datTim.time().hour();
+    epoStr = ch;
+    if (datTim.time().hour() < 23) {
+      nextTime.setHMS(datTim.time().hour() + 1 , 0, 0);
+      nextDate = datTim.date();
+    }
+    else {
+      nextTime.setHMS(0, 0, 0);
+      nextDate = datTim.date().addDays(1);
+    }
+  }
+  else {
+    epoStr = "0";
+    nextTime.setHMS(0, 0, 0);
+    nextDate = datTim.date().addDays(1);
+  }
+
+  if (nextEpoch) {
+   *nextEpoch = QDateTime(nextDate, nextTime);
+  }
+
+  return epoStr;
+}
+
+// File Name according to RINEX Standards
+////////////////////////////////////////////////////////////////////////////
+void bnsSP3::resolveFileName(const QDateTime& datTim) {
+
+  QString hlpStr = nextEpochStr(datTim, _intr, &_nextCloseEpoch);
+
+  _fName = (_ID4
+            + QString("%1").arg(datTim.date().dayOfYear(), 3, 10, QChar('0'))
+            + hlpStr 
+            + _ext).toAscii();
+}
+
+// Write Header
+////////////////////////////////////////////////////////////////////////////
+void bnsSP3::writeHeader(const QDateTime& datTim) {
+
+  // Open the Output File
+  // --------------------
+  resolveFileName(datTim);
+
+  _out.open(_fName.data());
+  _out.setf(ios::showpoint | ios::fixed);
+
+  _out << "THIS IS A DUMMY HEADER" << endl;
+
+  _headerWritten = true;
+}
+
+// Write One Epoch
 ////////////////////////////////////////////////////////////////////////////
 void bnsSP3::write(int GPSweek, double GPSweeks, const QString& prn, 
                    const ColumnVector& xx) {
+
+  QDateTime datTim = dateAndTimeFromGPSweek(GPSweek, GPSweeks);
+
+  // Close the file
+  // --------------
+  if (_nextCloseEpoch.isValid() && datTim >= _nextCloseEpoch) {
+    closeFile();
+    _headerWritten = false;
+  }
+
+  // Write Header
+  // ------------
+  if (!_headerWritten) {
+    writeHeader(datTim);
+  }
+
+  int year, month, day, hour, min;
+  double sec;
+
+  _out << "*  " << setw(4) << year 
+       << setw(3) << month 
+       << setw(3) << day
+       << setw(3) << hour 
+       << setw(3) << min
+       << setw(12) << setprecision(8) << sec << endl; 
+  _out << "P" << prn.toAscii().data()
+       << setw(14) << setprecision(6) << xx(1) / 1000.0
+       << setw(14) << setprecision(6) << xx(2) / 1000.0
+       << setw(14) << setprecision(6) << xx(3) / 1000.0
+       << " 999999.999999" << endl;
 }

trunk/BNS/bnssp3.h

@@ -8 +8 @@
 class bnsSP3 {
  public:
-   bnsSP3();
-   ~bnsSP3();
-   void write(int GPSweek, double GPSweeks, const QString& prn, 
-              const ColumnVector& xx);
+  bnsSP3();
+  ~bnsSP3();
+  void write(int GPSweek, double GPSweeks, const QString& prn, 
+             const ColumnVector& xx);
+
+  static QString nextEpochStr(const QDateTime& datTim,
+                              const QString& intStr, 
+                              QDateTime* nextEpoch = 0);
+
  private:
+  void resolveFileName(const QDateTime& datTim);
+  void writeHeader(const QDateTime& datTim);
+  void closeFile();
+
+  int           _samplingRate;
+  bool          _headerWritten;
+  QDateTime     _nextCloseEpoch;
+  std::ofstream _out;
+  QString       _path;
+  QString       _intr;
+  QString       _ext;
+  QString       _ID4;
+  QByteArray    _fName;
 };