Changeset 3521 in ntrip
- Timestamp:
- Dec 14, 2011, 12:32:41 PM (13 years ago)
- Location:
- trunk/BNC
- Files:
-
- 5 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/BNC/bncephuser.cpp
r3518 r3521 165 165 if (!relevantMessageType(messageType)) { 166 166 return failure; 167 }168 169 if (messageType < 0) {170 messageType = -messageType;171 xyzCorr = true; // correction in xyz instead of rao172 167 } 173 168 -
trunk/BNC/bncephuser.h
r3508 r3521 44 44 eph = 0; 45 45 hrClk = 0.0; 46 xyzCorr = false;47 46 } 48 47 bool ready() {return raoSet && dClkSet;} 49 48 50 49 static bool relevantMessageType(int msgType) { 51 if (msgType < 0) {52 msgType = -msgType;53 }54 50 return ( msgType == COTYPE_GPSCOMBINED || 55 51 msgType == COTYPE_GLONASSCOMBINED || … … 76 72 bool raoSet; 77 73 bool dClkSet; 78 bool xyzCorr;79 74 const t_eph* eph; 80 75 }; -
trunk/BNC/bncpppclient.cpp
r3519 r3521 376 376 if (eLast && eLast->IOD() == cc->iod) { 377 377 eLast->position(tt.gpsw(), tt.gpssec(), xc.data(), vv.data()); 378 return applyCorr(tt, cc, xc, vv , eLast);378 return applyCorr(tt, cc, xc, vv); 379 379 } 380 380 else if (ePrev && ePrev->IOD() == cc->iod) { 381 381 ePrev->position(tt.gpsw(), tt.gpssec(), xc.data(), vv.data()); 382 return applyCorr(tt, cc, xc, vv , ePrev);382 return applyCorr(tt, cc, xc, vv); 383 383 } 384 384 } … … 399 399 //////////////////////////////////////////////////////////////////////////// 400 400 t_irc bncPPPclient::applyCorr(const bncTime& tt, const t_corr* cc, 401 ColumnVector& xc, ColumnVector& vv, 402 const t_eph* eph) { 401 ColumnVector& xc, ColumnVector& vv) { 403 402 404 403 double dt = tt - cc->tt; … … 409 408 } 410 409 411 if (cc->xyzCorr) { 412 xc[0] += raoHlp[0]; 413 xc[1] += raoHlp[1]; 414 xc[2] += raoHlp[2]; 415 ColumnVector xcHlp(4); 416 ColumnVector vvHlp(3); 417 eph->position(cc->tt.gpsw(), cc->tt.gpssec(), xcHlp.data(), vvHlp.data()); 418 xc[3] += cc->dClk - xcHlp[3]; 419 } 420 else { 421 ColumnVector dx(3); 422 RSW_to_XYZ(xc.Rows(1,3), vv, raoHlp, dx); 423 xc[0] -= dx[0]; 424 xc[1] -= dx[1]; 425 xc[2] -= dx[2]; 426 xc[3] += cc->dClk + cc->dotDClk * dt + cc->dotDotDClk * dt * dt 427 + cc->hrClk; 428 } 410 ColumnVector dx(3); 411 RSW_to_XYZ(xc.Rows(1,3), vv, raoHlp, dx); 412 xc[0] -= dx[0]; 413 xc[1] -= dx[1]; 414 xc[2] -= dx[2]; 415 xc[3] += cc->dClk + cc->dotDClk * dt + cc->dotDotDClk * dt * dt 416 + cc->hrClk; 429 417 430 418 return success; -
trunk/BNC/bncpppclient.h
r3519 r3521 125 125 void putNewObs(const t_obs& pp); 126 126 static t_irc applyCorr(const bncTime& tt, const t_corr* cc, ColumnVector& xc, 127 ColumnVector& vv , const t_eph* eph);127 ColumnVector& vv); 128 128 129 129 public slots: -
trunk/BNC/combination/bnccomb.cpp
r3519 r3521 658 658 corr->eph->position(time12.gpsw(), time12.gpssec(), 659 659 xc.data(), vv.data()); 660 bncPPPclient::applyCorr(time12, corr, xc, vv , corr->eph);660 bncPPPclient::applyCorr(time12, corr, xc, vv); 661 661 662 662 // Relativistic Correction
Note:
See TracChangeset
for help on using the changeset viewer.