Changeset 6288 in ntrip
 Timestamp:
 Nov 1, 2014, 2:54:28 PM (9 years ago)
 Location:
 trunk/BNC/src/rinex
 Files:

 2 edited
Legend:
 Unmodified
 Added
 Removed

trunk/BNC/src/rinex/reqcanalyze.cpp
r6287 r6288 163 163 t_satObs satObs; 164 164 t_rnxObsFile::setObsFromRnx(obsFile, _currEpo, rnxSat, satObs); 165 t_qc Obs& qcObs = qcEpo._qcObs[satObs._prn];166 setQcObs(qcEpo._epoTime, xyzSta, satObs, qc Obs);167 updateQcSat(qc Obs, _qcFile._qcSatSum[satObs._prn]);165 t_qcSat& qcSat = qcEpo._qcSat[satObs._prn]; 166 setQcObs(qcEpo._epoTime, xyzSta, satObs, qcSat); 167 updateQcSat(qcSat, _qcFile._qcSatSum[satObs._prn]); 168 168 } 169 169 _qcFile._qcEpo.push_back(qcEpo); … … 245 245 // 246 246 //////////////////////////////////////////////////////////////////////////// 247 void t_reqcAnalyze::updateQcSat(const t_qc Obs& qcObs, t_qcSatSum& qcSatSum) {248 249 for (int ii = 0; ii < qc Obs._qcFrq.size(); ii++) {250 const t_qcFrq& qcFrq = qc Obs._qcFrq[ii];247 void t_reqcAnalyze::updateQcSat(const t_qcSat& qcSat, t_qcSatSum& qcSatSum) { 248 249 for (int ii = 0; ii < qcSat._qcFrq.size(); ii++) { 250 const t_qcFrq& qcFrq = qcSat._qcFrq[ii]; 251 251 t_qcFrqSum& qcFrqSum = qcSatSum._qcFrqSum[qcFrq._rnxType2ch]; 252 252 qcFrqSum._numObs += 1; … … 263 263 //////////////////////////////////////////////////////////////////////////// 264 264 void t_reqcAnalyze::setQcObs(const bncTime& epoTime, const ColumnVector& xyzSta, 265 const t_satObs& satObs, t_qc Obs& qcObs) {265 const t_satObs& satObs, t_qcSat& qcSat) { 266 266 267 267 t_eph* eph = 0; … … 278 278 double rho, eleSat, azSat; 279 279 topos(xyzSta(1), xyzSta(2), xyzSta(3), xc(1), xc(2), xc(3), rho, eleSat, azSat); 280 qc Obs._eleSet = true;281 qc Obs._azDeg = azSat * 180.0/M_PI;282 qc Obs._eleDeg = eleSat * 180.0/M_PI;280 qcSat._eleSet = true; 281 qcSat._azDeg = azSat * 180.0/M_PI; 282 qcSat._eleDeg = eleSat * 180.0/M_PI; 283 283 } 284 284 if (satObs._prn.system() == 'R') { 285 qc Obs._slotSet = true;286 qc Obs._slotNum = eph>slotNum();285 qcSat._slotSet = true; 286 qcSat._slotNum = eph>slotNum(); 287 287 } 288 288 } … … 294 294 const t_frqObs* frqObs = satObs._obs[ii]; 295 295 296 qc Obs._qcFrq.push_back(t_qcFrq());297 t_qcFrq& qcFrq = qc Obs._qcFrq.back();296 qcSat._qcFrq.push_back(t_qcFrq()); 297 t_qcFrq& qcFrq = qcSat._qcFrq.back(); 298 298 299 299 qcFrq._rnxType2ch = QString(frqObs>_rnxType2ch.c_str()); … … 353 353 frqObsB>_phaseValid && frqObsB>_codeValid) { 354 354 355 double f_a = t_CST::freq(fA, qc Obs._slotNum);356 double f_b = t_CST::freq(fB, qc Obs._slotNum);355 double f_a = t_CST::freq(fA, qcSat._slotNum); 356 double f_b = t_CST::freq(fB, qcSat._slotNum); 357 357 358 358 double L_a = frqObs>_phase * t_CST::c / f_a; … … 410 410 t_qcEpo& qcEpo = _qcFile._qcEpo[iEpo]; 411 411 if (chunkStart <= qcEpo._epoTime && qcEpo._epoTime < chunkEnd) { 412 if (qcEpo._qc Obs.contains(prn)) {413 t_qc Obs& qcObs = qcEpo._qcObs[prn];414 for (int iFrq = 0; iFrq < qc Obs._qcFrq.size(); iFrq++) {415 t_qcFrq& qcFrq = qc Obs._qcFrq[iFrq];412 if (qcEpo._qcSat.contains(prn)) { 413 t_qcSat& qcSat = qcEpo._qcSat[prn]; 414 for (int iFrq = 0; iFrq < qcSat._qcFrq.size(); iFrq++) { 415 t_qcFrq& qcFrq = qcSat._qcFrq[iFrq]; 416 416 if (qcFrq._rnxType2ch == frqType) { 417 417 frqVec << &qcFrq; … … 503 503 for (int iEpo = 0; iEpo < _qcFile._qcEpo.size(); iEpo++) { 504 504 t_qcEpo& qcEpo = _qcFile._qcEpo[iEpo]; 505 QMapIterator<t_prn, t_qc Obs> it(qcEpo._qcObs);505 QMapIterator<t_prn, t_qcSat> it(qcEpo._qcSat); 506 506 while (it.hasNext()) { 507 507 it.next(); 508 508 const t_prn& prn = it.key(); 509 const t_qc Obs& qcObs= it.value();509 const t_qcSat& qcSat = it.value(); 510 510 if ( (prn.system() == 'G' && plotGPS)  511 511 (prn.system() == 'R' && plotGlo)  512 512 (prn.system() == 'E' && plotGal) ) { 513 513 514 if (qc Obs._eleSet) {514 if (qcSat._eleSet) { 515 515 QString frqType1; 516 516 QString frqType2; 517 for (int iFrq = 0; iFrq < qc Obs._qcFrq.size(); iFrq++) {518 const t_qcFrq& qcFrq = qc Obs._qcFrq[iFrq];517 for (int iFrq = 0; iFrq < qcSat._qcFrq.size(); iFrq++) { 518 const t_qcFrq& qcFrq = qcSat._qcFrq[iFrq]; 519 519 if (qcFrq._rnxType2ch[0] == '1' && frqType1.isEmpty()) { 520 520 frqType1 = qcFrq._rnxType2ch; … … 524 524 } 525 525 if (qcFrq._rnxType2ch == frqType1) { 526 (*dataSNR1) << (new t_polarPoint(qc Obs._azDeg, 90.0  qcObs._eleDeg, qcFrq._SNR));527 (*dataMP1) << (new t_polarPoint(qc Obs._azDeg, 90.0  qcObs._eleDeg, qcFrq._stdMP));526 (*dataSNR1) << (new t_polarPoint(qcSat._azDeg, 90.0  qcSat._eleDeg, qcFrq._SNR)); 527 (*dataMP1) << (new t_polarPoint(qcSat._azDeg, 90.0  qcSat._eleDeg, qcFrq._stdMP)); 528 528 } 529 529 else if (qcFrq._rnxType2ch == frqType2) { 530 (*dataSNR2) << (new t_polarPoint(qc Obs._azDeg, 90.0  qcObs._eleDeg, qcFrq._SNR));531 (*dataMP2) << (new t_polarPoint(qc Obs._azDeg, 90.0  qcObs._eleDeg, qcFrq._stdMP));530 (*dataSNR2) << (new t_polarPoint(qcSat._azDeg, 90.0  qcSat._eleDeg, qcFrq._SNR)); 531 (*dataMP2) << (new t_polarPoint(qcSat._azDeg, 90.0  qcSat._eleDeg, qcFrq._stdMP)); 532 532 } 533 533 } … … 637 637 plotData._mjdX24 << mjdX24; 638 638 plotData._PDOP << qcEpo._PDOP; 639 plotData._numSat << qcEpo._qc Obs.size();640 641 QMapIterator<t_prn, t_qc Obs> it(qcEpo._qcObs);639 plotData._numSat << qcEpo._qcSat.size(); 640 641 QMapIterator<t_prn, t_qcSat> it(qcEpo._qcSat); 642 642 while (it.hasNext()) { 643 643 it.next(); 644 644 const t_prn& prn = it.key(); 645 const t_qc Obs& qcObs= it.value();645 const t_qcSat& qcSat = it.value(); 646 646 t_plotData& data = plotDataMap[prn]; 647 647 648 648 data._mjdX24 << mjdX24; 649 data._eleDeg << qc Obs._eleDeg;649 data._eleDeg << qcSat._eleDeg; 650 650 651 651 QString frqType1; 652 652 QString frqType2; 653 for (int iFrq = 0; iFrq < qc Obs._qcFrq.size(); iFrq++) {654 const t_qcFrq& qcFrq = qc Obs._qcFrq[iFrq];653 for (int iFrq = 0; iFrq < qcSat._qcFrq.size(); iFrq++) { 654 const t_qcFrq& qcFrq = qcSat._qcFrq[iFrq]; 655 655 if (qcFrq._rnxType2ch[0] == '1' && frqType1.isEmpty()) { 656 656 frqType1 = qcFrq._rnxType2ch; … … 684 684 } 685 685 } 686 687 686 688 687 if (BNC_CORE>GUIenabled()) { 
trunk/BNC/src/rinex/reqcanalyze.h
r6285 r6288 86 86 }; 87 87 88 class t_qc Obs{88 class t_qcSat { 89 89 public: 90 t_qc Obs() {90 t_qcSat() { 91 91 _slotSet = false; 92 92 _eleSet = false; … … 107 107 bncTime _epoTime; 108 108 double _PDOP; 109 QMap<t_prn, t_qc Obs> _qcObs;109 QMap<t_prn, t_qcSat> _qcSat; 110 110 }; 111 111 … … 153 153 void analyzeFile(t_rnxObsFile* obsFile); 154 154 155 void updateQcSat(const t_qc Obs& qcObs, t_qcSatSum& qcSatSum);155 void updateQcSat(const t_qcSat& qcSat, t_qcSatSum& qcSatSum); 156 156 157 157 void setQcObs(const bncTime& epoTime, const ColumnVector& xyzSta, 158 const t_satObs& satObs, t_qc Obs& qcObs);158 const t_satObs& satObs, t_qcSat& qcSat); 159 159 160 160 void analyzeMultipath();
Note:
See TracChangeset
for help on using the changeset viewer.