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


Ignore:
Timestamp:
Nov 1, 2014, 4:11:16 PM (10 years ago)
Author:
mervart
Message:
 
File:
1 edited

Legend:

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

    r6291 r6292  
    313313    // Compute the Multipath Linear Combination
    314314    // ----------------------------------------
    315     if (frqObs->_phaseValid && frqObs->_codeValid) {
     315    if (frqObs->_codeValid) {
    316316      t_frequency::type fA;
    317317      t_frequency::type fB;
     
    348348
    349349      if (fA != t_frequency::dummy && fB != t_frequency::dummy) {
     350        double f_a = t_CST::freq(fA, qcSat._slotNum);
     351        double f_b = t_CST::freq(fB, qcSat._slotNum);
     352        double C_a = frqObs->_code;
     353
     354        bool   foundA = false;
     355        double L_a    = 0.0;
     356        bool   foundB = false;
     357        double L_b    = 0.0;
    350358        for (unsigned jj = 0; jj < satObs._obs.size(); jj++) {
    351           const t_frqObs* frqObsB = satObs._obs[jj];
    352           if (frqObsB->_rnxType2ch[0] == t_frequency::toString(fB)[1] &&
    353               frqObsB->_phaseValid) {
    354 
    355             double f_a = t_CST::freq(fA, qcSat._slotNum);
    356             double f_b = t_CST::freq(fB, qcSat._slotNum);
    357 
    358             double L_a = frqObs->_phase  * t_CST::c / f_a;
    359             double C_a = frqObs->_code;
    360             double L_b = frqObsB->_phase * t_CST::c / f_b;
    361 
    362             qcFrq._setMP = true;
    363             qcFrq._rawMP = C_a - L_a - 2.0*f_b*f_b/(f_a*f_a-f_b*f_b) * (L_a - L_b);
    364           }
     359          const t_frqObs* frqObsHlp = satObs._obs[jj];
     360          if      (frqObsHlp->_rnxType2ch[0] == t_frequency::toString(fA)[1] &&
     361              frqObsHlp->_phaseValid) {
     362            foundA = true;
     363            L_a    = frqObsHlp->_phase * t_CST::c / f_a;
     364          }
     365          else if (frqObsHlp->_rnxType2ch[0] == t_frequency::toString(fB)[1] &&
     366              frqObsHlp->_phaseValid) {
     367            foundB = true;
     368            L_b    = frqObsHlp->_phase * t_CST::c / f_b;
     369          }
     370        }
     371        if (foundA && foundB) {
     372          qcFrq._setMP = true;
     373          qcFrq._rawMP = C_a - L_a - 2.0*f_b*f_b/(f_a*f_a-f_b*f_b) * (L_a - L_b);
    365374        }
    366375      }
Note: See TracChangeset for help on using the changeset viewer.