Changeset 10356 in ntrip
- Timestamp:
- Feb 27, 2024, 10:05:36 AM (6 months ago)
- Location:
- trunk/BNC/src
- Files:
-
- 4 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/BNC/src/PPP/pppClient.cpp
r10333 r10356 255 255 ColumnVector& xyzc, bool print) { 256 256 t_lc::type tLC = t_lc::dummy; 257 char sysBancroft = 'G';258 257 int numBancroft = obsVector.size(); 259 260 _usedSystems['G'] = _usedSystems['R'] = _usedSystems['E'] = _usedSystems['C'] = 0;261 for (unsigned jj = 0; jj < obsVector.size(); jj++) {262 const t_pppSatObs* satObs = obsVector[jj];263 char sys = satObs->prn().system();264 _usedSystems[sys]++;265 }266 267 if ((numBancroft = _usedSystems.value('G')) >= _opt->_minObs) {268 sysBancroft = 'G';269 }270 else if ((numBancroft = _usedSystems.value('E')) >= _opt->_minObs) {271 sysBancroft = 'E';272 }273 else if ((numBancroft = _usedSystems.value('C')) >= _opt->_minObs) {274 sysBancroft = 'C';275 }276 else if ((numBancroft = _usedSystems.value('R')) >= _opt->_minObs) {277 sysBancroft = 'R';278 }279 else {280 LOG << "t_pppClient::cmpBancroft not enough observations: " << endl;281 return failure;282 }283 258 284 259 while (_running) { … … 287 262 for (unsigned ii = 0; ii < obsVector.size(); ii++) { 288 263 const t_pppSatObs* satObs = obsVector.at(ii); 289 if (satObs->prn().system() == sysBancroft) { 290 if (tLC == t_lc::dummy) { 291 if (satObs->isValid(t_lc::cIF)) { 292 tLC = t_lc::cIF; 293 } 294 else if (satObs->isValid(t_lc::c1)) { 295 tLC = t_lc::c1; 296 } 297 else if (satObs->isValid(t_lc::c2)) { 298 tLC = t_lc::c2; 299 } 264 if (tLC == t_lc::dummy) { 265 if (satObs->isValid(t_lc::cIF)) { 266 tLC = t_lc::cIF; 300 267 } 301 if ( satObs->isValid(tLC) && 302 (!satObs->modelSet() || satObs->eleSat() >= _opt->_minEle) ) { 303 ++iObs; 304 BB[iObs][0] = satObs->xc()[0]; 305 BB[iObs][1] = satObs->xc()[1]; 306 BB[iObs][2] = satObs->xc()[2]; 307 BB[iObs][3] = satObs->obsValue(tLC) - satObs->cmpValueForBanc(tLC); 268 else if (satObs->isValid(t_lc::c1)) { 269 tLC = t_lc::c1; 308 270 } 271 else if (satObs->isValid(t_lc::c2)) { 272 tLC = t_lc::c2; 273 } 274 } 275 if ( satObs->isValid(tLC) && 276 (!satObs->modelSet() || satObs->eleSat() >= _opt->_minEle) ) { 277 ++iObs; 278 BB[iObs][0] = satObs->xc()[0]; 279 BB[iObs][1] = satObs->xc()[1]; 280 BB[iObs][2] = satObs->xc()[2]; 281 BB[iObs][3] = satObs->obsValue(tLC) - satObs->cmpValueForBanc(tLC); 309 282 } 310 283 } … … 327 300 for (unsigned ii = 0; ii < obsVector.size(); ii++) { 328 301 const t_pppSatObs* satObs = obsVector.at(ii); 329 char sys = satObs->prn().system(); 330 if (satObs->isValid() && sys == sysBancroft && 331 (!satObs->modelSet() || satObs->eleSat() >= _opt->_minEle) ) { 302 if (satObs->isValid() && (!satObs->modelSet() || satObs->eleSat() >= _opt->_minEle) ) { 332 303 ColumnVector rr = satObs->xc().Rows(1,3) - xyzc.Rows(1,3); 333 304 double res = rr.NormFrobenius() - satObs->obsValue(tLC) … … 345 316 if (!_epoTimeRover.undef()) LOG << string(_epoTimeRover); 346 317 LOG << "\n---------------------------------------------------------------\n"; 347 LOG << string(epoTime) << " BANCROFT " << sysBancroft<< ": "318 LOG << string(epoTime) << " BANCROFT " << ": " 348 319 << setw(14) << setprecision(3) << xyzc[0] << ' ' 349 320 << setw(14) << setprecision(3) << xyzc[1] << ' ' -
trunk/BNC/src/PPP/pppFilter.cpp
r10335 r10356 375 375 376 376 // Check Pre-Fit Residuals 377 / *-----------------------377 // ----------------------- 378 378 else { 379 379 ColumnVector AA(params.size()); … … 402 402 resetAmb(obs->prn(), obsVector, tLC); 403 403 } 404 } */404 } 405 405 } 406 406 } -
trunk/BNC/src/PPP/pppParlist.cpp
r10326 r10356 45 45 _noise = 0.0; 46 46 _ambInfo = 0; 47 double fac = 1.0;48 47 49 48 switch (_type) { … … 68 67 break; 69 68 case amb: 70 if (tLC == t_lc::lIF) {fac = 3.0;}71 69 _epoSpec = false; 72 _sigma0 = fac *OPT->_aprSigAmb;70 _sigma0 = OPT->_aprSigAmb; 73 71 _ambInfo = new t_ambInfo(); 74 72 if (obsVector) { … … 499 497 // GLONASS Clock Offset 500 498 // -------------------- 501 if ( _usedSystems.contains('R') && 502 (_usedSystems.contains('G') || _usedSystems.contains('E') || _usedSystems.contains('C'))) { 499 if (OPT->useSystem('R')) { 503 500 required.push_back(new t_pppParam(t_pppParam::offGlo, t_prn(), t_lc::dummy)); 504 }505 else {506 PPP_CLIENT->resetOffGlo();507 501 } 508 502 509 503 // Galileo Clock Offset 510 504 // -------------------- 511 if ( _usedSystems.contains('E') && _usedSystems.contains('G')) {505 if (OPT->useSystem('E')) { 512 506 required.push_back(new t_pppParam(t_pppParam::offGal, t_prn(), t_lc::dummy)); 513 }514 else {515 PPP_CLIENT->resetOffGal();516 507 } 517 508 518 509 // BDS Clock Offset 519 510 // ---------------- 520 if ( _usedSystems.contains('C') && (_usedSystems.contains('G') || _usedSystems.contains('E'))) {511 if (OPT->useSystem('C')) { 521 512 required.push_back(new t_pppParam(t_pppParam::offBds, t_prn(), t_lc::dummy)); 522 }523 else {524 PPP_CLIENT->resetOffBds();525 513 } 526 514 -
trunk/BNC/src/pppMain.cpp
r10327 r10356 474 474 // Some default values 475 475 // ------------------- 476 opt->_aprSigClk = 3 00000.0;477 opt->_aprSigClkOff = 3 00000.0;478 opt->_aprSigAmb = 1 0000.0;479 opt->_aprSigIon = 1.0e 3;480 opt->_noiseIon = 1.0e 3;476 opt->_aprSigClk = 3.0e5; 477 opt->_aprSigClkOff = 3.0e5; 478 opt->_aprSigAmb = 1.0e4; 479 opt->_aprSigIon = 1.0e8; 480 opt->_noiseIon = 1.0e8; 481 481 opt->_aprSigCodeBias = 10000.0; 482 482 opt->_noiseCodeBias = 10000.0;
Note:
See TracChangeset
for help on using the changeset viewer.