Changeset 10356 in ntrip for trunk/BNC/src/PPP/pppClient.cpp
- Timestamp:
- Feb 27, 2024, 10:05:36 AM (5 months ago)
- File:
-
- 1 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] << ' '
Note:
See TracChangeset
for help on using the changeset viewer.