Changeset 7474 in ntrip for trunk/BNC/src/rinex/reqcedit.cpp

Timestamp:

Sep 29, 2015, 9:29:30 AM (9 years ago)

Author:

stuerze

Message:

in case of RINEX3 to RINEX 2 conversion now some comment lines are written into the RINEX header

File:

• trunk/BNC/src/rinex/reqcedit.cpp (modified) (21 diffs)

trunk/BNC/src/rinex/reqcedit.cpp

@@ -35 +35 @@
  * Created:    11-Apr-2012
  *
- * Changes:    
+ * Changes:
  *
  * -----------------------------------------------------------------------*/
@@ -88 +88 @@
 }
 
-//  
+//
 ////////////////////////////////////////////////////////////////////////////
 void t_reqcEdit::run() {
- 
+
   // Open Log File
   // -------------
@@ -120 +120 @@
           << _samplingRate << endl;
     *_log << QByteArray("Start time").leftJustified(15) << ": "
-          << _begTime.datestr().c_str() << ' ' 
+          << _begTime.datestr().c_str() << ' '
           << _begTime.timestr(0).c_str() << endl;
     *_log << QByteArray("End time").leftJustified(15) << ": "
-          << _endTime.datestr().c_str() << ' ' 
+          << _endTime.datestr().c_str() << ' '
           << _endTime.timestr(0).c_str() << endl;
     *_log << QByteArray("Input Obs Files").leftJustified(15) << ": "
@@ -159 +159 @@
 // Initialize input observation files, sort them according to start time
 ////////////////////////////////////////////////////////////////////////////
-void t_reqcEdit::initRnxObsFiles(const QStringList& obsFileNames, 
+void t_reqcEdit::initRnxObsFiles(const QStringList& obsFileNames,
                                  QVector<t_rnxObsFile*>& rnxObsFiles,
                                  QTextStream* log) {
@@ -172 +172 @@
       QListIterator<QFileInfo> it(dir.entryInfoList(filters));
       while (it.hasNext()) {
-        QString filePath = it.next().filePath(); 
+        QString filePath = it.next().filePath();
         t_rnxObsFile* rnxObsFile = 0;
         try {
@@ -199 +199 @@
     }
   }
-  qStableSort(rnxObsFiles.begin(), rnxObsFiles.end(), 
+  qStableSort(rnxObsFiles.begin(), rnxObsFiles.end(),
               t_rnxObsFile::earlierStartTime);
 }
 
-//  
+//
 ////////////////////////////////////////////////////////////////////////////
 void t_reqcEdit::editObservations() {
@@ -218 +218 @@
   // ----------------------------------
   t_rnxObsFile outObsFile(_outObsFileName, t_rnxObsFile::output);
-  
+
   // Select observation types
   // ------------------------
@@ -296 +296 @@
     t_rnxObsFile* obsFile = _rnxObsFiles[ii];
     if (_log) {
-      *_log << "Processing File: " << obsFile->fileName() << "  start: " 
-            << obsFile->startTime().datestr().c_str() << ' ' 
+      *_log << "Processing File: " << obsFile->fileName() << "  start: "
+            << obsFile->startTime().datestr().c_str() << ' '
             << obsFile->startTime().timestr(0).c_str() << endl;
     }
@@ -320 +320 @@
         txtMap["COMMENT"]  = comment;
       }
+      if (int(_rnxVersion) < int(obsFile->header().version())) {
+        addRnxConversionDetails(obsFile, txtMap);
+      }
       outObsFile.header().write(outObsFile.stream(), &txtMap);
     }
@@ -331 +334 @@
           break;
         }
-      
-        if (_samplingRate == 0 || 
+
+        if (_samplingRate == 0 ||
             fmod(round(epo->tt.gpssec()), _samplingRate) == 0) {
           applyLLI(obsFile, epo);
@@ -347 +350 @@
       }
       else {
-        qDebug() << str;    
+        qDebug() << str;
       }
       return;
@@ -363 +366 @@
 }
 
-// Change RINEX Header Content  
+// Change RINEX Header Content
 ////////////////////////////////////////////////////////////////////////////
 void t_reqcEdit::editRnxObsHeader(t_rnxObsFile& obsFile) {
@@ -372 +375 @@
   QString newMarkerName   = settings.value("reqcNewMarkerName").toString();
   if (!newMarkerName.isEmpty()) {
-    if (oldMarkerName.isEmpty() || 
+    if (oldMarkerName.isEmpty() ||
         QRegExp(oldMarkerName).exactMatch(obsFile.markerName())) {
       obsFile.setMarkerName(newMarkerName);
@@ -381 +384 @@
   QString newAntennaName  = settings.value("reqcNewAntennaName").toString();
   if (!newAntennaName.isEmpty()) {
-    if (oldAntennaName.isEmpty() || 
+    if (oldAntennaName.isEmpty() ||
         QRegExp(oldAntennaName).exactMatch(obsFile.antennaName())) {
       obsFile.setAntennaName(newAntennaName);
@@ -426 +429 @@
   QString newReceiverType = settings.value("reqcNewReceiverName").toString();
   if (!newReceiverType.isEmpty()) {
-    if (oldReceiverType.isEmpty() || 
+    if (oldReceiverType.isEmpty() ||
         QRegExp(oldReceiverType).exactMatch(obsFile.receiverType())) {
       obsFile.setReceiverType(newReceiverType);
@@ -442 +445 @@
 }
 
-// 
-////////////////////////////////////////////////////////////////////////////
-void t_reqcEdit::rememberLLI(const t_rnxObsFile* obsFile, 
+//
+////////////////////////////////////////////////////////////////////////////
+void t_reqcEdit::rememberLLI(const t_rnxObsFile* obsFile,
                              const t_rnxObsFile::t_rnxEpo* epo) {
 
@@ -467 +470 @@
   }
 }
-  
-// 
-////////////////////////////////////////////////////////////////////////////
-void t_reqcEdit::applyLLI(const t_rnxObsFile* obsFile, 
+
+//
+////////////////////////////////////////////////////////////////////////////
+void t_reqcEdit::applyLLI(const t_rnxObsFile* obsFile,
                           t_rnxObsFile::t_rnxEpo* epo) {
- 
+
  if (_samplingRate == 0) {
     return;
@@ -509 +512 @@
       QListIterator<QFileInfo> it(dir.entryInfoList(filters));
       while (it.hasNext()) {
-        QString filePath = it.next().filePath(); 
+        QString filePath = it.next().filePath();
         appendEphemerides(filePath, ephs);
       }
@@ -520 +523 @@
 }
 
-//  
+//
 ////////////////////////////////////////////////////////////////////////////
 void t_reqcEdit::editEphemerides() {
@@ -594 +597 @@
 }
 
-//  
+//
 ////////////////////////////////////////////////////////////////////////////
 void t_reqcEdit::appendEphemerides(const QString& fileName,
@@ -633 +636 @@
   }
 }
+
+void t_reqcEdit::addRnxConversionDetails(const t_rnxObsFile* obsFile,
+                                          QMap<QString, QString>& txtMap) {
+
+  int key = 0;
+  QString systems = obsFile->header().usedSystems();
+  QString comment = QString("RINEX 3 => 2 CONVERSION DETAILS:");
+  QString commentKey = QString("COMMENT %1").arg(key, 3, 10, QChar('0'));
+  txtMap.insert(commentKey, comment);
+
+  for(int ii = 0; ii < obsFile->numSys(); ii++) {
+    key++;
+    char sys = systems[ii].toAscii();
+    QString preferredAttrib = obsFile->signalPriorities(sys);
+    comment = QString("%1: Signal priority = %2").arg(sys).arg(preferredAttrib);
+    commentKey = QString("COMMENT %1").arg(key, 3, 10, QChar('0'));
+    txtMap.insert(commentKey, comment);
+    QStringList types = obsFile->header().obsTypes(sys);
+    for (int jj = 0; jj < types.size(); jj++) {
+      key++;
+      QString inType = types[jj];
+      for (int iPref = 0; iPref < preferredAttrib.length(); iPref++) {
+        if (preferredAttrib[iPref] == '?'                             ||
+            (inType.length() == 2 && preferredAttrib[iPref] == '_'    ) ||
+            (inType.length() == 3 && preferredAttrib[iPref] == inType[2]) ) {
+          QString outType = t_rnxObsFile::type3to2(sys, inType);
+          comment = QString("%1: %2 => %3").arg(sys).arg(inType).arg(outType);
+          commentKey = QString("COMMENT %1").arg(key, 3, 10, QChar('0'));
+          txtMap.insert(commentKey, comment);
+          break;
+        }
+      }
+    }
+  }
+}