Changeset 4700 in ntrip for trunk/BNC


Ignore:
Timestamp:
Sep 9, 2012, 6:33:18 PM (12 years ago)
Author:
mervart
Message:
 
File:
1 edited

Legend:

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

    r4699 r4700  
    5757using namespace std;
    5858
    59 const double SLIPTRESH = 5.0; // cycle-slip threshold (meters)
     59const double SLIPTRESH = 10.0; // cycle-slip threshold (meters)
    6060
    6161// Constructor
     
    545545    // Compute the Multipath
    546546    // ---------------------
    547     bool slipMP = false;
    548     meanMP1 /= numEpo;
    549     meanMP2 /= numEpo;
    550     double MP1 = 0.0;
    551     double MP2 = 0.0;
    552     for (int ii = 0; ii < numEpo; ii++) {
    553       int iEpo = chunkStart + ii;
    554       const t_oneObs* oneObs = allObs._oneObsVec[iEpo];
    555       double diff1 = oneObs->_MP1 - meanMP1;
    556       double diff2 = oneObs->_MP2 - meanMP2;
    557 
    558       // Check Slip Threshold
    559       // --------------------
    560       if (fabs(diff1) > SLIPTRESH || fabs(diff2) > SLIPTRESH) {
    561         slipMP = true;
    562         break;
    563       }
    564 
    565       MP1 += diff1 * diff1;
    566       MP2 += diff2 * diff2;
    567     }
    568     if (slipMP) {
    569       slipL1 = true;
    570       slipL2 = true;
    571     }
    572     else {
    573       MP1 = sqrt(MP1 / (numEpo-1));
    574       MP2 = sqrt(MP2 / (numEpo-1));
    575       (*dataMP1)  << (new t_polarPoint(aziDeg, zenDeg, MP1));
    576       (*dataMP2)  << (new t_polarPoint(aziDeg, zenDeg, MP2));
     547    if (prn[0] != 'R') { // TODO
     548      bool slipMP = false;
     549      meanMP1 /= numEpo;
     550      meanMP2 /= numEpo;
     551      double MP1 = 0.0;
     552      double MP2 = 0.0;
     553      for (int ii = 0; ii < numEpo; ii++) {
     554        int iEpo = chunkStart + ii;
     555        const t_oneObs* oneObs = allObs._oneObsVec[iEpo];
     556        double diff1 = oneObs->_MP1 - meanMP1;
     557        double diff2 = oneObs->_MP2 - meanMP2;
     558     
     559        // Check Slip Threshold
     560        // --------------------
     561        if (fabs(diff1) > SLIPTRESH || fabs(diff2) > SLIPTRESH) {
     562          slipMP = true;
     563          break;
     564        }
     565     
     566        MP1 += diff1 * diff1;
     567        MP2 += diff2 * diff2;
     568      }
     569      if (slipMP) {
     570        slipL1 = true;
     571        slipL2 = true;
     572        _obsStat._prnStat[prn]._numSlips += 1;
     573      }
     574      else {
     575        MP1 = sqrt(MP1 / (numEpo-1));
     576        MP2 = sqrt(MP2 / (numEpo-1));
     577        (*dataMP1)  << (new t_polarPoint(aziDeg, zenDeg, MP1));
     578        (*dataMP2)  << (new t_polarPoint(aziDeg, zenDeg, MP2));
     579      }
    577580    }
    578581
Note: See TracChangeset for help on using the changeset viewer.