Changeset 7474 in ntrip for trunk/BNC/src/rinex/rnxobsfile.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/rnxobsfile.cpp (modified) (26 diffs)

trunk/BNC/src/rinex/rnxobsfile.cpp

@@ -35 +35 @@
  * Created:    24-Jan-2012
  *
- * Changes:    
+ * Changes:
  *
  * -----------------------------------------------------------------------*/
@@ -310 +310 @@
   _obsTypes.clear();
   if (_version < 3.0) {
-    _obsTypes['G'] << "C1" << "P1" << "L1" << "S1" 
+    _obsTypes['G'] << "C1" << "P1" << "L1" << "S1"
                    << "C2" << "P2" << "L2" << "S2";
     _obsTypes['R'] = _obsTypes['G'];
@@ -319 +319 @@
   }
   else {
-    _obsTypes['G'] << "C1C" << "L1C"  << "S1C" 
+    _obsTypes['G'] << "C1C" << "L1C"  << "S1C"
                    << "C1W" << "L1W"  << "S1W"
                    << "C2X" << "L2X"  << "S2X"
-                   << "C2W" << "L2W"  << "S2W" 
+                   << "C2W" << "L2W"  << "S2W"
                    << "C5X" << "L5X"  << "S5X";
 
     _obsTypes['J'] = _obsTypes['G'];
-    
-    _obsTypes['R'] << "C1C" << "L1C" << "S1C" 
+
+    _obsTypes['R'] << "C1C" << "L1C" << "S1C"
                    << "C2P" << "L2P" << "S2P";
-    
+
     _obsTypes['E'] << "C1X" << "L1X" << "SX1"
                    << "C5X" << "L5X" << "SX5"
                    << "C7X" << "L7X" << "SX7"
                    << "C8X" << "L8X" << "SX8";
-    
+
     _obsTypes['S'] << "C1C" << "L1C" << "S1C"
                    << "C5I" << "L5I" << "S5I"
                    << "C5Q" << "L5Q" << "S5Q";
-    
+
     _obsTypes['C'] << "C2I" << "L2I" << "S2I"
                    << "C6I" << "L6I" << "S6I"
@@ -359 +359 @@
     _version = t_rnxObsHeader::defaultRnxObsVersion3;
   }
-  _interval        = header._interval;    
-  _antennaNumber   = header._antennaNumber; 
-  _antennaName     = header._antennaName; 
-  _markerName      = header._markerName;  
-  _markerNumber    = header._markerNumber;  
+  _interval        = header._interval;
+  _antennaNumber   = header._antennaNumber;
+  _antennaName     = header._antennaName;
+  _markerName      = header._markerName;
+  _markerNumber    = header._markerNumber;
   _markerType      = header._markerType;
-  _antNEU          = header._antNEU;      
-  _antXYZ          = header._antXYZ;      
-  _antBSG          = header._antBSG;      
-  _xyz             = header._xyz;         
-  _observer        = header._observer;       
-  _agency          = header._agency;         
-  _receiverNumber  = header._receiverNumber; 
-  _receiverType    = header._receiverType;   
+  _antNEU          = header._antNEU;
+  _antXYZ          = header._antXYZ;
+  _antBSG          = header._antBSG;
+  _xyz             = header._xyz;
+  _observer        = header._observer;
+  _agency          = header._agency;
+  _receiverNumber  = header._receiverNumber;
+  _receiverType    = header._receiverType;
   _receiverVersion = header._receiverVersion;
-  _startTime       = header._startTime;   
+  _startTime       = header._startTime;
   _comments        = header._comments;
   _usedSystems     = header._usedSystems;
 
   for (unsigned iPrn = 1; iPrn <= t_prn::MAXPRN_GPS; iPrn++) {
-    _wlFactorsL1[iPrn] =  header._wlFactorsL1[iPrn]; 
-    _wlFactorsL2[iPrn] =  header._wlFactorsL2[iPrn]; 
+    _wlFactorsL1[iPrn] =  header._wlFactorsL1[iPrn];
+    _wlFactorsL2[iPrn] =  header._wlFactorsL2[iPrn];
   }
 
@@ -444 +444 @@
         for (int iSys = 0; iSys < _usedSystems.length(); iSys++) {
           char    sys  = _usedSystems[iSys].toAscii();
-          QString type = _version >= 3.0 ? t_rnxObsFile::type2to3(sys, useObsTypes->at(iType)) : 
+          QString type = _version >= 3.0 ? t_rnxObsFile::type2to3(sys, useObsTypes->at(iType)) :
                                            t_rnxObsFile::type3to2(sys, useObsTypes->at(iType));
           if (!_obsTypes[sys].contains(type)) {
@@ -523 +523 @@
 
   if (txtMap) {
-    QMapIterator<QString, QString> it(*txtMap);
+   QMapIterator<QString, QString> it(*txtMap);
     while (it.hasNext()) {
       it.next();
@@ -529 +529 @@
         runBy = it.value();
       }
-      else if (it.key() == "COMMENT") {
-        newComments = it.value().split("\\n", QString::SkipEmptyParts);
+      else if ((it.key().indexOf("COMMENT")) != -1) {
+        newComments += it.value().split("\\n", QString::SkipEmptyParts);
       }
     }
@@ -570 +570 @@
     .leftJustified(60)
            << "MARKER TYPE\n";
-           
+
   *stream << QString("%1%2")
     .arg(_observer, -20)
@@ -766 +766 @@
 }
 
-// 
+//
 ////////////////////////////////////////////////////////////////////////////
 char t_rnxObsHeader::system(int iSys) const {
@@ -781 +781 @@
 }
 
+//
+////////////////////////////////////////////////////////////////////////////
+QString t_rnxObsHeader::usedSystems(void) const {
+  return _usedSystems;
+}
+
+QStringList t_rnxObsHeader::obsTypes(char sys) const {
+  if (_obsTypes.contains(sys)) {
+    return _obsTypes[sys];
+  }
+  else {
+    return QStringList();
+  }
+}
+
 // Number of Observation Types (satellite-system specific)
 ////////////////////////////////////////////////////////////////////////////
@@ -873 +888 @@
     QTextStream(&hlp) << QString("%1").arg(_obsTypes[sys0].size(), 6);
     for (int ii = 0; ii < _obsTypes[sys0].size(); ii++) {
-      QTextStream(&hlp) << QString("%1").arg(_obsTypes[sys0][ii], 6);   
+      QTextStream(&hlp) << QString("%1").arg(_obsTypes[sys0][ii], 6);
       if ((ii+1) % 9 == 0 || ii == _obsTypes[sys0].size()-1) {
         strList.append(hlp.leftJustified(60) + "# / TYPES OF OBSERV\n");
@@ -887 +902 @@
       for (int iType = 0; iType < nTypes(sys); iType++) {
         QString type = obsType(sys, iType);
-        QTextStream(&hlp) << QString(" %1").arg(type, -3);   
+        QTextStream(&hlp) << QString(" %1").arg(type, -3);
         if ((iType+1) % 13 == 0 || iType == nTypes(sys)-1) {
           strList.append(hlp.leftJustified(60) + "SYS / # / OBS TYPES\n");
@@ -985 +1000 @@
 // Handle Special Epoch Flag
 ////////////////////////////////////////////////////////////////////////////
-void t_rnxObsFile::handleEpochFlag(int flag, const QString& line, 
+void t_rnxObsFile::handleEpochFlag(int flag, const QString& line,
                                    bool& headerReRead) {
 
@@ -1003 +1018 @@
 
   // Re-Read Header
-  // -------------- 
+  // --------------
   else if (flag == 3 || flag == 4 || flag == 5) {
     int numLines = 0;
@@ -1077 +1092 @@
     int numSat;
     readInt(line, 32, 3, numSat);
-  
+
     _currEpo.rnxSat.resize(numSat);
 
@@ -1090 +1105 @@
         int pos = 3 + 16*iType;
         double obsValue = 0.0;
-        int    lli      = 0; 
+        int    lli      = 0;
         int    snr      = 0;
         readDbl(line, pos,     14, obsValue);
@@ -1154 +1169 @@
     int numSat;
     readInt(line, 29, 3, numSat);
-  
+
     _currEpo.rnxSat.resize(numSat);
 
@@ -1172 +1187 @@
       int satNum; readInt(line, pos + 1, 2, satNum);
       _currEpo.rnxSat[iSat].prn.set(sys, satNum);
- 
+
       pos += 3;
     }
@@ -1206 +1221 @@
       }
     }
- 
+
     _flgPowerFail = false;
 
     return &_currEpo;
   }
- 
+
   return 0;
 }
@@ -1226 +1241 @@
     const t_rnxSat& rnxSat = epo->rnxSat[ii];
     if (_header._obsTypes[rnxSat.prn.system()].size() > 0) {
-      epoLocal.rnxSat.push_back(rnxSat);    
+      epoLocal.rnxSat.push_back(rnxSat);
     }
   }
@@ -1240 +1255 @@
 // Write Data Epoch (RINEX Version 2)
 ////////////////////////////////////////////////////////////////////////////
-void t_rnxObsFile::writeEpochV2(QTextStream* stream, const t_rnxObsHeader& header, 
+void t_rnxObsFile::writeEpochV2(QTextStream* stream, const t_rnxObsHeader& header,
                                 const t_rnxEpo* epo) {
 
@@ -1278 +1293 @@
       QString typeV2 = header.obsType(sys, iTypeV2);
       bool    found  = false;
-   
+
       QString preferredAttrib = signalPriorities(sys);
       for (int iPref = 0; iPref < preferredAttrib.length(); iPref++) {
@@ -1324 +1339 @@
 // Write Data Epoch (RINEX Version 3)
 ////////////////////////////////////////////////////////////////////////////
-void t_rnxObsFile::writeEpochV3(QTextStream* stream, const t_rnxObsHeader& header, 
+void t_rnxObsFile::writeEpochV3(QTextStream* stream, const t_rnxObsHeader& header,
                                 const t_rnxEpo* epo) {
 
@@ -1431 +1446 @@
 // Set Observations from RINEX File
 ////////////////////////////////////////////////////////////////////////////
-void t_rnxObsFile::setObsFromRnx(const t_rnxObsFile* rnxObsFile, const t_rnxObsFile::t_rnxEpo* epo, 
+void t_rnxObsFile::setObsFromRnx(const t_rnxObsFile* rnxObsFile, const t_rnxObsFile::t_rnxEpo* epo,
                                  const t_rnxObsFile::t_rnxSat& rnxSat, t_satObs& obs) {
   obs._staID = rnxObsFile->markerName().toAscii().constData();
@@ -1461 +1476 @@
       if (rnxObs.value != 0.0) {
         string type2ch(typeV3.mid(1).toAscii().data());
-        
+
         t_frqObs* frqObs = 0;
         for (unsigned iFrq = 0; iFrq < obs._obs.size(); iFrq++) {
@@ -1474 +1489 @@
           obs._obs.push_back(frqObs);
         }
-        
+
         switch( typeV3.toAscii().data()[0] ) {
         case 'C':