Changeset 1243 in ntrip for trunk/BNS/bns.cpp
- Timestamp:
- Nov 28, 2008, 4:19:06 PM (15 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/BNS/bns.cpp
r1217 r1243 491 491 void t_bns::crdTrafo(int GPSWeek, ColumnVector& xyz) { 492 492 493 // Current epoch minus 2000.0 in years 494 // ------------------------------------ 495 double dt = (GPSWeek - (1042.0+6.0/7.0)) / 365.2422 * 7.0; 496 493 497 ColumnVector dx(3); 494 dx(1) = 0.054; 495 dx(2) = 0.051; 496 dx(3) = -0.048; 498 dx(1) = 0.0541 - dt * 0.0002; 499 dx(2) = 0.0502 + dt * 0.0001; 500 dx(3) = -0.0538 - dt * 0.0018; 501 497 502 static const double arcSec = 180.0 * 3600.0 / M_PI; 498 static const double ox = 0.000081 / arcSec; 499 static const double oy = 0.000490 / arcSec; 500 static const double oz = -0.000792 / arcSec; 503 504 static const double ox = ( 0.000891 + dt * 0.000081) / arcSec; 505 static const double oy = ( 0.005390 + dt * 0.000490) / arcSec; 506 static const double oz = (-0.008712 - dt * 0.000792) / arcSec; 501 507 502 508 Matrix rMat(3,3); rMat = 0.0; … … 508 514 rMat(3,2) = ox; 509 515 510 // Current epoch minus 1989.0 in years 511 // ------------------------------------ 512 double dt = (GPSWeek - 469.0) / 365.2422 * 7.0; 513 514 xyz += dx + dt * rMat * xyz; 515 } 516 517 xyz += dx + rMat * xyz; 518 }
Note:
See TracChangeset
for help on using the changeset viewer.