Changeset 6280 in ntrip


Ignore:
Timestamp:
Oct 31, 2014, 7:24:17 PM (10 years ago)
Author:
mervart
Message:
 
File:
1 edited

Legend:

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

    r6279 r6280  
    419419      bncTime prevTime;
    420420      for (int iEpo = 0; iEpo < _qcFile._qcEpo.size(); iEpo++) {
    421         const t_qcEpo& qcEpo = _qcFile._qcEpo[iEpo];
     421        t_qcEpo& qcEpo = _qcFile._qcEpo[iEpo];
    422422        if (qcEpo._epoTime < chunkStart) {
    423423          continue;
     
    427427        }
    428428     
    429         const t_qcObs& qcObs = qcEpo._qcObs[prn];
     429        t_qcObs& qcObs = qcEpo._qcObs[prn];
    430430     
    431431        // Compute the Azimuth and Zenith Distance
     
    453453                double rho, eleSat, azSat;
    454454                topos(xyzSta(1), xyzSta(2), xyzSta(3), xc(1), xc(2), xc(3), rho, eleSat, azSat);
    455                 qcObsSampled._azDeg  = azSat * 180.0/M_PI;
     455                qcObsSampled._azDeg  = azSat  * 180.0/M_PI;
    456456                qcObsSampled._eleDeg = eleSat * 180.0/M_PI;
     457                qcObs._azDeg         = azSat  * 180.0/M_PI;
     458                qcObs._eleDeg        = eleSat * 180.0/M_PI;
    457459              }
    458460            }
     
    647649void t_reqcAnalyze::slotDspAvailPlot(const QString& fileName, const QByteArray& title) {
    648650
     651  t_plotData              plotData;
    649652  QMap<t_prn, t_plotData> plotDataMap;
    650653
    651654  for (int ii = 0; ii < _qcFile._qcEpo.size(); ii++) {
    652655    const t_qcEpo& qcEpo = _qcFile._qcEpo[ii];
     656    double mjdX24 = qcEpo._epoTime.mjddec() * 24.0;
     657
     658    plotData._mjdX24 << mjdX24;
     659    plotData._PDOP   << qcEpo._PDOP;
     660    plotData._numSat << qcEpo._qcObs.size();
     661
    653662    QMapIterator<t_prn, t_qcObs> it(qcEpo._qcObs);
    654663    while (it.hasNext()) {
     
    656665      const t_prn&   prn   = it.key();
    657666      const t_qcObs& qcObs = it.value();
    658       t_plotData& plotData = plotDataMap[prn];
    659       double mjdX24 = qcEpo._epoTime.mjddec() * 24.0;
    660       if (qcObs._hasL1)  plotData._L1ok   << mjdX24;
    661       if (qcObs._hasL2)  plotData._L2ok   << mjdX24;
    662       if (qcObs._slipL1) plotData._L1slip << mjdX24;
    663       if (qcObs._slipL2) plotData._L2slip << mjdX24;
    664       if (qcObs._gapL1)  plotData._L1gap  << mjdX24;
    665       if (qcObs._gapL2)  plotData._L2gap  << mjdX24;
    666     }
    667   }
    668 
    669   t_plotData plotData;
     667      t_plotData&    data  = plotDataMap[prn];
     668      if (qcObs._hasL1)  data._L1ok   << mjdX24;
     669      if (qcObs._hasL2)  data._L2ok   << mjdX24;
     670      if (qcObs._slipL1) data._L1slip << mjdX24;
     671      if (qcObs._slipL2) data._L2slip << mjdX24;
     672      if (qcObs._gapL1)  data._L1gap  << mjdX24;
     673      if (qcObs._gapL2)  data._L2gap  << mjdX24;
     674      data._eleDeg << qcObs._eleDeg;
     675    }
     676  }
     677
    670678
    671679  if (BNC_CORE->GUIenabled()) {
Note: See TracChangeset for help on using the changeset viewer.