Changeset 9797 in ntrip


Ignore:
Timestamp:
Aug 10, 2022, 4:52:43 PM (8 weeks ago)
Author:
stuerze
Message:

minor changes

Location:
trunk/BNC/src
Files:
3 edited

Legend:

Unmodified
Added
Removed
  • trunk/BNC/src/bncutils.cpp

    r9787 r9797  
    780780
    781781  eleSat = acos( sqrt(neu[0]*neu[0] + neu[1]*neu[1]) / rho );
    782   if (neu[2] < 0) {
     782  if (neu[2] < 0.0) {
    783783    eleSat *= -1.0;
    784784  }
  • trunk/BNC/src/rinex/reqcanalyze.cpp

    r9775 r9797  
    120120
    121121  if (_log) {
    122     *_log << "QC Format Version  : " << QString("%1").arg(QC_FORMAT_VERSION,3,'f',1)  << endl << endl;
     122    *_log << "QC Format Version  : " << QString("%1").arg(QC_FORMAT_VERSION,3,'f',1)  << '\n' << '\n';
    123123  }
    124124
     
    203203      for (unsigned iObs = 0; iObs < _currEpo->rnxSat.size(); iObs++) {
    204204        const t_rnxObsFile::t_rnxSat& rnxSat = _currEpo->rnxSat[iObs];
    205         if (_navFileNames.size() &&
    206             _numExpObs.find(rnxSat.prn) == _numExpObs.end()) {
     205        if (_navFileNames.size() && _numExpObs.find(rnxSat.prn) == _numExpObs.end()) {
    207206          _numExpObs[rnxSat.prn] = 0;
    208207        }
     
    231230  catch (QString str) {
    232231    if (_log) {
    233       *_log << "Exception " << str << endl;
     232      *_log << "Exception " << str << '\n';
    234233    }
    235234    else {
     
    342341      double rho, eleSat, azSat;
    343342      topos(xyzSta(1), xyzSta(2), xyzSta(3), xc(1), xc(2), xc(3), rho, eleSat, azSat);
    344       qcSat._eleSet = true;
    345       qcSat._azDeg  = azSat  * 180.0/M_PI;
    346       qcSat._eleDeg = eleSat * 180.0/M_PI;
     343      if (round(eleSat * 180.0/M_PI) >= 0.0) {
     344        qcSat._eleSet = true;
     345        qcSat._azDeg  = azSat  * 180.0/M_PI;
     346        qcSat._eleDeg = eleSat * 180.0/M_PI;
     347      }
    347348    }
    348349    if (satObs._prn.system() == 'R') {
     
    381382      t_frequency::type fA = t_frequency::dummy;
    382383      t_frequency::type fB = t_frequency::dummy;
    383       char sys             = satObs._prn.system();
     384      char sys = satObs._prn.system();
    384385      if (_signalTypes.find(sys) != _signalTypes.end()) {
    385386        for (int iSig = 0; iSig < _signalTypes[sys].size(); iSig++) {
     
    412413          for (unsigned jj = 0; jj < satObs._obs.size(); jj++) {
    413414            const t_frqObs* frqObsHlp = satObs._obs[jj];
    414             if      (frqObsHlp->_rnxType2ch[0] == t_frequency::toString(fA)[1] &&
    415                 frqObsHlp->_phaseValid) {
     415            if      (frqObsHlp->_rnxType2ch[0] == t_frequency::toString(fA)[1] && frqObsHlp->_phaseValid) {
    416416              foundA = true;
    417417              L_a    = frqObsHlp->_phase * t_CST::c / f_a;
    418418            }
    419             else if (frqObsHlp->_rnxType2ch[0] == t_frequency::toString(fB)[1] &&
    420                 frqObsHlp->_phaseValid) {
     419            else if (frqObsHlp->_rnxType2ch[0] == t_frequency::toString(fB)[1] && frqObsHlp->_phaseValid) {
    421420              foundB = true;
    422421              L_b    = frqObsHlp->_phase * t_CST::c / f_b;
    423422            }
    424423          }
    425           if (foundA && foundB) {
     424          if (foundA && foundB && C_a) {
    426425            qcFrq._setMP = true;
    427426            qcFrq._rawMP = C_a - L_a - 2.0*f_b*f_b/(f_a*f_a-f_b*f_b) * (L_a - L_b);
     
    748747  // Summary
    749748  // -------
    750   *_log << "Observation File   : " << obsFileName                                   << endl
    751         << "RINEX Version      : " << QString("%1").arg(obsFile->version(),4,'f',2) << endl
    752         << "Marker Name        : " << _qcFile._markerName                           << endl
    753         << "Marker Number      : " << obsFile->markerNumber()                       << endl
    754         << "Receiver           : " << _qcFile._receiverType                         << endl
    755         << "Antenna            : " << _qcFile._antennaName                          << endl
     749  *_log << "Observation File   : " << obsFileName                                   << '\n'
     750        << "RINEX Version      : " << QString("%1").arg(obsFile->version(),4,'f',2) << '\n'
     751        << "Marker Name        : " << _qcFile._markerName                           << '\n'
     752        << "Marker Number      : " << obsFile->markerNumber()                       << '\n'
     753        << "Receiver           : " << _qcFile._receiverType                         << '\n'
     754        << "Antenna            : " << _qcFile._antennaName                          << '\n'
    756755        << "Position XYZ       : " << QString("%1 %2 %3").arg(obsFile->xyz()(1), 14, 'f', 4)
    757756                                                        .arg(obsFile->xyz()(2), 14, 'f', 4)
    758                                                         .arg(obsFile->xyz()(3), 14, 'f', 4) << endl
     757                                                        .arg(obsFile->xyz()(3), 14, 'f', 4) << '\n'
    759758        << "Antenna dH/dE/dN   : " << QString("%1 %2 %3").arg(obsFile->antNEU()(3), 8, 'f', 4)
    760759                                                        .arg(obsFile->antNEU()(2), 8, 'f', 4)
    761                                                         .arg(obsFile->antNEU()(1), 8, 'f', 4) << endl
     760                                                        .arg(obsFile->antNEU()(1), 8, 'f', 4) << '\n'
    762761        << "Start Time         : " << _qcFile._startTime.datestr().c_str()         << ' '
    763                                    << _qcFile._startTime.timestr(1,'.').c_str()    << endl
     762                                   << _qcFile._startTime.timestr(1,'.').c_str()    << '\n'
    764763        << "End Time           : " << _qcFile._endTime.datestr().c_str()           << ' '
    765                                    << _qcFile._endTime.timestr(1,'.').c_str()      << endl
    766         << "Interval           : " << _qcFile._interval << " sec"                  << endl;
     764                                   << _qcFile._endTime.timestr(1,'.').c_str()      << '\n'
     765        << "Interval           : " << _qcFile._interval << " sec"                  << '\n';
    767766
    768767  // Number of systems
     
    783782    *_log << ' ' << itSys.key();
    784783  }
    785   *_log << endl;
     784  *_log << '\n';
    786785
    787786  // Observation types per system
     
    795794        *_log << " " << type;
    796795      }
    797       *_log << endl;
     796      *_log << '\n';
    798797    }
    799798  }
     
    804803  while (itSys.hasNext()) {
    805804    itSys.next();
    806     const QChar&                      sys      = itSys.key();
     805    const QChar&  sys = itSys.key();
    807806    const QVector<const t_qcSatSum*>& qcSatVec = itSys.value();
    808807    int numExpectedObs = 0;
    809     for(QMap<t_prn, int>::iterator it = _numExpObs.begin();
    810         it != _numExpObs.end(); it++) {
     808    for(QMap<t_prn, int>::iterator it = _numExpObs.begin(); it != _numExpObs.end(); it++) {
    811809      if (sys == it.key().system()) {
    812810        numExpectedObs += it.value();
     
    820818      while (itFrq.hasNext()) {
    821819        itFrq.next();
    822         QString           frqType  = itFrq.key(); if (frqType.length() < 2) frqType += '?';
     820        QString frqType  = itFrq.key(); if (frqType.length() < 2) frqType += '?';
    823821        const t_qcFrqSum& qcFrqSum = itFrq.value();
    824822        frqMap[frqType].push_back(&qcFrqSum);
    825823      }
    826824    }
    827     *_log << endl
    828           << prefixSys << "Satellites: " << qcSatVec.size() << endl
     825    *_log << '\n'
     826          << prefixSys << "Satellites: " << qcSatVec.size() << '\n'
    829827          << prefixSys << "Signals   : " << frqMap.size() << "   ";
    830828    QMapIterator<QString, QVector<const t_qcFrqSum*> > itFrq(frqMap);
     
    834832      *_log << ' ' << frqType;
    835833    }
    836     *_log << endl;
     834    *_log << '\n';
    837835    QString prefixSys2 = "    " + prefixSys;
    838836    itFrq.toFront();
    839837    while (itFrq.hasNext()) {
    840838      itFrq.next();
    841       QString                          frqType  = itFrq.key(); if (frqType.length() < 2) frqType += '?';
     839      QString frqType  = itFrq.key(); if (frqType.length() < 2) frqType += '?';
    842840      const QVector<const t_qcFrqSum*> qcFrqVec = itFrq.value();
    843841      QString prefixFrq = QString("  ") + frqType + QString(": ");
     
    870868      double ratio = (double(numObs) / double(numExpectedObs)) * 100.0;
    871869
    872       *_log << endl
     870      *_log << '\n'
    873871            << prefixSys2 << prefixFrq << "Observations      : ";
    874       if(_navFileNames.isEmpty() || _navFileIncomplete.contains(sys.toLatin1())) {
     872      if(_navFileNames.isEmpty() || numExpectedObs == 0.0  || _navFileIncomplete.contains(sys.toLatin1())) {
    875873        *_log << QString("%1\n").arg(numObs,           6);
    876874      }
     
    892890    return;
    893891  }
    894   *_log << endl;
     892  *_log << '\n';
    895893  for (int iEpo = 0; iEpo < _qcFile._qcEpo.size(); iEpo++) {
    896894    const t_qcEpo& qcEpo = _qcFile._qcEpo[iEpo];
     
    912910    *_log << dateStr << QString(" %1").arg(qcEpo._qcSat.size(), 2)
    913911          << QString(" %1").arg(qcEpo._PDOP, 4, 'f', 1)
    914           << endl;
     912          << '\n';
    915913
    916914    QMapIterator<t_prn, t_qcSat> itSat(qcEpo._qcSat);
     
    964962        }
    965963      }
    966       *_log << endl;
     964      *_log << '\n';
    967965    }
    968966  }
     
    988986  }
    989987  if (_log) {
    990     *_log << "Navigation File(s) : " << navFileName  << endl;
     988    *_log << "Navigation File(s) : " << navFileName  << '\n';
    991989  }
    992990  QStringListIterator it(_navFileNames);
     
    10401038  }
    10411039  if (_log) {
    1042     *_log << endl;
     1040    *_log << '\n';
    10431041  }
    10441042}
     
    10471045                                   double interval, const ColumnVector& xyzSta) {
    10481046
    1049   for(QMap<t_prn, int>::iterator it = _numExpObs.begin();
    1050       it != _numExpObs.end(); it++) {
     1047
     1048  for(QMap<t_prn, int>::iterator it = _numExpObs.begin(); it != _numExpObs.end(); it++) {
    10511049    t_eph* eph = 0;
    10521050    for (int ie = 0; ie < _ephs.size(); ie++) {
     
    10601058      int numExpObs = 0;
    10611059      bncTime epoTime;
    1062       for (epoTime = startTime - interval; epoTime < endTime;
    1063            epoTime = epoTime + interval) {
     1060      for (epoTime = startTime - interval; epoTime < endTime; epoTime = epoTime + interval) {
    10641061        ColumnVector xc(6);
    10651062        ColumnVector vv(3);
     
    10681065          double rho, eleSat, azSat;
    10691066          topos(xyzSta(1), xyzSta(2), xyzSta(3), xc(1), xc(2), xc(3), rho, eleSat, azSat);
    1070           if ((eleSat * 180.0/M_PI) > 0.0) {
     1067          if (round(eleSat * 180.0/M_PI) >= 0.0) {
    10711068            numExpObs++;
    10721069          }
  • trunk/BNC/src/rinex/reqcedit.cpp

    r9770 r9797  
    680680                                   QVector<t_eph*>& ephs) {
    681681  t_rnxNavFile rnxNavFile(fileName, t_rnxNavFile::input);
     682  bncEphUser ephUser(false);
    682683  for (unsigned ii = 0; ii < rnxNavFile.ephs().size(); ii++) {
    683684    t_eph* eph   = rnxNavFile.ephs()[ii];
Note: See TracChangeset for help on using the changeset viewer.