Changeset 7546 in ntrip for trunk/BNC/src/rinex


Ignore:
Timestamp:
Oct 22, 2015, 1:00:03 PM (9 years ago)
Author:
stuerze
Message:

observation types per system are added in header of the QC output

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/BNC/src/rinex/reqcanalyze.cpp

    r6867 r7546  
    7878  connect(this, SIGNAL(dspSkyPlot(const QString&, const QString&, QVector<t_polarPoint*>*,
    7979                                  const QString&, QVector<t_polarPoint*>*,
    80                                   const QByteArray&, double)), 
     80                                  const QByteArray&, double)),
    8181          this, SLOT(slotDspSkyPlot(const QString&, const QString&, QVector<t_polarPoint*>*,
    8282                                    const QString&, QVector<t_polarPoint*>*,
     
    196196  // --------------------
    197197  try {
    198     QMap<QString, bncTime> lastObsTime; 
     198    QMap<QString, bncTime> lastObsTime;
    199199    bool firstEpo = true;
    200200    while ( (_currEpo = obsFile->nextEpoch()) != 0) {
     
    329329//
    330330////////////////////////////////////////////////////////////////////////////
    331 void t_reqcAnalyze::setQcObs(const bncTime& epoTime, const ColumnVector& xyzSta, 
     331void t_reqcAnalyze::setQcObs(const bncTime& epoTime, const ColumnVector& xyzSta,
    332332                             const t_satObs& satObs, QMap<QString, bncTime>& lastObsTime,
    333333                             t_qcSat& qcSat) {
     
    458458      // Loop over all Chunks of Data
    459459      // ----------------------------
    460       for (bncTime chunkStart = _qcFile._startTime; 
     460      for (bncTime chunkStart = _qcFile._startTime;
    461461           chunkStart < _qcFile._endTime; chunkStart += chunkStep) {
    462462
     
    465465        QVector<t_qcFrq*> frqVec;
    466466        QVector<double>   MP;
    467    
     467
    468468        // Loop over all Epochs within one Chunk of Data
    469469        // ---------------------------------------------
     
    494494          }
    495495          meanMP /= MP.size();
    496        
     496
    497497          bool slipMP = false;
    498        
     498
    499499          double stdMP = 0.0;
    500500          for (int ii = 0; ii < MP.size(); ii++) {
     
    506506            stdMP += diff * diff;
    507507          }
    508        
     508
    509509          if (slipMP) {
    510510            stdMP = 0.0;
     
    517517            qcFrqSum._sumMP += stdMP;
    518518          }
    519        
     519
    520520          for (int ii = 0; ii < frqVec.size(); ii++) {
    521521            t_qcFrq* qcFrq = frqVec[ii];
     
    659659    QVector<QWidget*> plots;
    660660    if (data1) {
    661       QwtText title(title1); 
     661      QwtText title(title1);
    662662      QFont font = title.font(); font.setPointSize(font.pointSize()-1); title.setFont(font);
    663663      t_polarPlot* plot1 = new t_polarPlot(title, scaleInterval, BNC_CORE->mainWindow());
     
    793793  // Summary
    794794  // -------
    795   *_log << "Observation File  : " << obsFileName                                   << endl
    796         << "RINEX Version     : " << QString("%1").arg(obsFile->version(),4,'f',2) << endl
    797         << "Marker Name       : " << _qcFile._markerName                           << endl
    798         << "Marker Number     : " << obsFile->markerNumber()                       << endl
    799         << "Receiver          : " << _qcFile._receiverType                         << endl
    800         << "Antenna           : " << _qcFile._antennaName                          << endl
    801         << "Position XYZ      : " << QString("%1 %2 %3").arg(obsFile->xyz()(1), 14, 'f', 4)
     795  *_log << "Observation File   : " << obsFileName                                   << endl
     796        << "RINEX Version      : " << QString("%1").arg(obsFile->version(),4,'f',2) << endl
     797        << "Marker Name        : " << _qcFile._markerName                           << endl
     798        << "Marker Number      : " << obsFile->markerNumber()                       << endl
     799        << "Receiver           : " << _qcFile._receiverType                         << endl
     800        << "Antenna            : " << _qcFile._antennaName                          << endl
     801        << "Position XYZ       : " << QString("%1 %2 %3").arg(obsFile->xyz()(1), 14, 'f', 4)
    802802                                                        .arg(obsFile->xyz()(2), 14, 'f', 4)
    803803                                                        .arg(obsFile->xyz()(3), 14, 'f', 4) << endl
    804         << "Antenna dH/dE/dN  : " << QString("%1 %2 %3").arg(obsFile->antNEU()(3), 8, 'f', 4)
     804        << "Antenna dH/dE/dN   : " << QString("%1 %2 %3").arg(obsFile->antNEU()(3), 8, 'f', 4)
    805805                                                        .arg(obsFile->antNEU()(2), 8, 'f', 4)
    806806                                                        .arg(obsFile->antNEU()(1), 8, 'f', 4) << endl
    807         << "Start Time        : " << _qcFile._startTime.datestr().c_str()         << ' '
     807        << "Start Time         : " << _qcFile._startTime.datestr().c_str()         << ' '
    808808                                  << _qcFile._startTime.timestr(1,'.').c_str()    << endl
    809         << "End Time          : " << _qcFile._endTime.datestr().c_str()           << ' '
     809        << "End Time           : " << _qcFile._endTime.datestr().c_str()           << ' '
    810810                                  << _qcFile._endTime.timestr(1,'.').c_str()      << endl
    811         << "Interval          : " << _qcFile._interval                            << endl;
     811        << "Interval           : " << _qcFile._interval                            << endl;
    812812
    813813  // Number of systems
    814814  // -----------------
    815   QMap<QChar, QVector<const t_qcSatSum*> > systemMap; 
     815  QMap<QChar, QVector<const t_qcSatSum*> > systemMap;
    816816  QMapIterator<t_prn, t_qcSatSum> itSat(_qcFile._qcSatSum);
    817817  while (itSat.hasNext()) {
     
    821821    systemMap[prn.system()].push_back(&qcSatSum);
    822822  }
    823   *_log << "Navigation Systems: " << systemMap.size() << "   ";
     823  *_log << "Navigation Systems : " << systemMap.size() << "   ";
    824824
    825825  QMapIterator<QChar, QVector<const t_qcSatSum*> > itSys(systemMap);
     
    829829  }
    830830  *_log << endl;
     831
     832  // Observation types per system
     833  // -----------------------------
     834  for (int iSys = 0; iSys < obsFile->numSys(); iSys++) {
     835    char sys = obsFile->system(iSys);
     836    if (sys != ' ') {
     837      *_log << "Observation Types " << sys << ":";
     838      for (int iType = 0; iType < obsFile->nTypes(sys); iType++) {
     839        QString type = obsFile->obsType(sys, iType);
     840        *_log << " " << type;
     841      }
     842      *_log << endl;
     843    }
     844  }
    831845
    832846  itSys.toFront();
     
    836850    const QVector<const t_qcSatSum*>& qcSatVec = itSys.value();
    837851    QString prefixSys = QString("  ") + sys + QString(": ");
    838     QMap<QString, QVector<const t_qcFrqSum*> > frqMap; 
     852    QMap<QString, QVector<const t_qcFrqSum*> > frqMap;
    839853    for (int ii = 0; ii < qcSatVec.size(); ii++) {
    840854      const t_qcSatSum* qcSatSum = qcSatVec[ii];
Note: See TracChangeset for help on using the changeset viewer.