Changeset 5857 in ntrip for trunk/BNC/src/PPP
- Timestamp:
- Aug 7, 2014, 1:58:21 PM (10 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/BNC/src/PPP/pppClient.cpp
r5831 r5857 280 280 281 281 if (OPT->useGlonass()) { 282 while (true) { 282 283 while (obsVector.size() > 0) { 283 284 offGG = 0.0; 284 bool outlierFound = false; 285 unsigned nObs = 0; 285 double maxRes = 0.0; 286 int maxResIndex = -1; 287 unsigned nObs = 0; 286 288 for (unsigned ii = 0; ii < obsVector.size(); ii++) { 287 289 t_pppSatObs* satObs = obsVector.at(ii); 288 if ( !satObs->outlier() &&satObs->isValid() && satObs->prn().system() == 'R' &&290 if ( satObs->isValid() && satObs->prn().system() == 'R' && 289 291 (!satObs->modelSet() || satObs->eleSat() >= OPT->_minEle) ) { 292 double ll = satObs->obsValue(tLC) - satObs->cmpValue(tLC); 290 293 ++nObs; 291 double ll = satObs->obsValue(tLC) - satObs->cmpValue(tLC); 292 if (fabs(ll) > 1000.0) { 293 satObs->setOutlier(); 294 outlierFound = true; 295 LOG << "t_pppClient::cmpOffGG outlier " << satObs->prn().toString() 296 << " " << ll << endl; 294 offGG += ll; 295 if (fabs(ll) > maxRes) { 296 maxRes = fabs(ll); 297 maxResIndex = ii; 297 298 } 298 offGG += ll;299 299 } 300 300 } 301 301 302 if (nObs > 0) { 302 303 offGG = offGG / nObs; … … 305 306 offGG = 0.0; 306 307 } 307 if (!outlierFound) { 308 309 if (maxRes > 1000.0) { 310 obsVector.erase(obsVector.begin() + maxResIndex); 311 } 312 else { 308 313 break; 309 314 }
Note:
See TracChangeset
for help on using the changeset viewer.