Changeset 5157 in ntrip


Ignore:
Timestamp:
May 26, 2013, 11:24:09 AM (9 years ago)
Author:
mervart
Message:
 
Location:
trunk/BNC/src
Files:
4 edited

Legend:

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

    r5070 r5157  
    110110  }
    111111
     112  // Glonass Offset
     113  // --------------
     114  else if (type == GLONASS_OFFSET) {
     115    if (satData->prn[0] == 'R') {
     116      return 1.0;
     117    }
     118    else {
     119      return 0.0;
     120    }
     121  }
     122
    112123  // Galileo Offset
    113124  // --------------
     
    228239    _params.push_back(new bncParam(bncParam::TROPO, ++nextPar, ""));
    229240  }
     241  if (_opt->useGlonass) {
     242    _params.push_back(new bncParam(bncParam::GLONASS_OFFSET, ++nextPar, ""));
     243  }
    230244  if (_opt->useGalileo) {
    231245    _params.push_back(new bncParam(bncParam::GALILEO_OFFSET, ++nextPar, ""));
     
    246260      _QQ(iPar,iPar) = _opt->sigTrp0 * _opt->sigTrp0;
    247261      pp->xx = lastTrp;
     262    }
     263    else if (pp->type == bncParam::GLONASS_OFFSET) {
     264      _QQ(iPar,iPar) = _opt->sigGlonassOffset0 * _opt->sigGlonassOffset0;
    248265    }
    249266    else if (pp->type == bncParam::GALILEO_OFFSET) {
     
    334351
    335352  double offset = 0.0;
    336   if (satData->prn[0] == 'E') {
     353  if      (satData->prn[0] == 'R') {
     354    offset = Glonass_offset();
     355  }
     356  else if (satData->prn[0] == 'E') {
    337357    offset = Galileo_offset();
    338358  }
     
    483503      }
    484504   
     505      // Glonass Offset
     506      // --------------
     507      else if (pp->type == bncParam::GLONASS_OFFSET) {
     508        _QQ(iPar,iPar) += _opt->sigGlonassOffsetP * _opt->sigGlonassOffsetP;
     509      }
     510
    485511      // Galileo Offset
    486512      // --------------
     
    621647      newPos->xnt[6] = aprTrp + par->xx;
    622648    }
     649    else if (par->type == bncParam::GLONASS_OFFSET) {
     650      strB << "\n    offGlo  = " << setw(10) << setprecision(3) << par->xx
     651           << " +- " << setw(6) << setprecision(3)
     652           << sqrt(_QQ(par->index,par->index));
     653    }
    623654    else if (par->type == bncParam::GALILEO_OFFSET) {
    624       strB << "\n    offset  = " << setw(10) << setprecision(3) << par->xx
     655      strB << "\n    offGal  = " << setw(10) << setprecision(3) << par->xx
    625656           << " +- " << setw(6) << setprecision(3)
    626657           << sqrt(_QQ(par->index,par->index));
  • trunk/BNC/src/bncmodel.h

    r4278 r5157  
    4141class bncParam {
    4242 public:
    43   enum parType {CRD_X, CRD_Y, CRD_Z, RECCLK, TROPO, AMB_L3, GALILEO_OFFSET};
     43  enum parType {CRD_X, CRD_Y, CRD_Z, RECCLK, TROPO, AMB_L3,
     44                GLONASS_OFFSET, GALILEO_OFFSET};
    4445  bncParam(parType typeIn, int indexIn, const QString& prn);
    4546  ~bncParam();
     
    7071      bncParam* pp = _params[ii];
    7172      if (pp->type == bncParam::TROPO) {
     73        return pp->xx;
     74      }
     75    }
     76    return 0.0;
     77  }
     78  double Glonass_offset() const {
     79    for (int ii = 0; ii < _params.size(); ++ii) {
     80      bncParam* pp = _params[ii];
     81      if (pp->type == bncParam::GLONASS_OFFSET) {
    7282        return pp->xx;
    7383      }
  • trunk/BNC/src/pppopt.cpp

    r4278 r5157  
    107107  if (!pppMode)     corrSync   = 0.0;
    108108
     109  sigGlonassOffset0 = 1000.0;
     110  sigGlonassOffsetP =    0.0;
    109111  sigGalileoOffset0 = 1000.0;
    110112  sigGalileoOffsetP =    0.0;
  • trunk/BNC/src/pppopt.h

    r4278 r5157  
    6363  bool    useGlonass;
    6464  bool    useGalileo;
     65  double  sigGlonassOffset0;
     66  double  sigGlonassOffsetP;
    6567  double  sigGalileoOffset0;
    6668  double  sigGalileoOffsetP;
Note: See TracChangeset for help on using the changeset viewer.