Changeset 1776 in ntrip
- Timestamp:
- Mar 26, 2009, 8:24:32 AM (16 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/BNS/bns.cpp
r1772 r1776 505 505 506 506 if (trafo == "ETRF2000") { 507 _dx = 0.0541;508 _dy = 0.0502;509 _dz = -0.0538;510 _dxr = -0.0002;511 _dyr = 0.0001;512 _dzr = -0.0018;513 _ox = 0.000891;514 _oy = 0.005390;515 _oz = -0.008712;516 _oxr = 0.000081;517 _oyr = 0.000490;518 _ozr = -0.000792;519 _sc = 0.40;520 _scr = 0.08;521 _t0 = 2000.0;507 _dx = 0.0541; 508 _dy = 0.0502; 509 _dz = -0.0538; 510 _dxr = -0.0002; 511 _dyr = 0.0001; 512 _dzr = -0.0018; 513 _ox = 0.000891; 514 _oy = 0.005390; 515 _oz = -0.008712; 516 _oxr = 0.000081; 517 _oyr = 0.000490; 518 _ozr = -0.000792; 519 _sc = 0.40; 520 _scr = 0.08; 521 _t0 = 2000.0; 522 522 } 523 523 else if (trafo == "Custom") { 524 _dx = settings.value("trafo_dx").toInt();525 _dy = settings.value("trafo_dy").toInt();526 _dz = settings.value("trafo_dz").toInt();527 _dxr = settings.value("trafo_dxr").toInt();528 _dyr = settings.value("trafo_dyr").toInt();529 _dzr = settings.value("trafo_dzr").toInt();530 _ox = settings.value("trafo_ox").toInt();531 _oy = settings.value("trafo_oy").toInt();532 _oz = settings.value("trafo_oz").toInt();533 _oxr = settings.value("trafo_oxr").toInt();534 _oyr = settings.value("trafo_oyr").toInt();535 _ozr = settings.value("trafo_ozr").toInt();536 _sc = settings.value("trafo_sc").toInt();537 _scr = settings.value("trafo_scr").toInt();538 _t0 = settings.value("trafo_t0").toInt();524 _dx = settings.value("trafo_dx").toDouble(); 525 _dy = settings.value("trafo_dy").toDouble(); 526 _dz = settings.value("trafo_dz").toDouble(); 527 _dxr = settings.value("trafo_dxr").toDouble(); 528 _dyr = settings.value("trafo_dyr").toDouble(); 529 _dzr = settings.value("trafo_dzr").toDouble(); 530 _ox = settings.value("trafo_ox").toDouble(); 531 _oy = settings.value("trafo_oy").toDouble(); 532 _oz = settings.value("trafo_oz").toDouble(); 533 _oxr = settings.value("trafo_oxr").toDouble(); 534 _oyr = settings.value("trafo_oyr").toDouble(); 535 _ozr = settings.value("trafo_ozr").toDouble(); 536 _sc = settings.value("trafo_sc").toDouble(); 537 _scr = settings.value("trafo_scr").toDouble(); 538 _t0 = settings.value("trafo_t0").toDouble(); 539 539 } 540 540 541 541 // Current epoch minus 2000.0 in years 542 542 // ------------------------------------ 543 //double dt = (GPSWeek - (1042.0+6.0/7.0)) / 365.2422 * 7.0;544 543 double dt = (GPSWeek - (1042.0+6.0/7.0)) / 365.2422 * 7.0 + 2000.0 - _t0; 545 544 546 545 ColumnVector dx(3); 547 546 548 //dx(1) = 0.0541 - dt * 0.0002;549 //dx(2) = 0.0502 + dt * 0.0001;550 //dx(3) = -0.0538 - dt * 0.0018;551 547 dx(1) = _dx + dt * _dxr; 552 548 dx(2) = _dy + dt * _dyr; … … 555 551 static const double arcSec = 180.0 * 3600.0 / M_PI; 556 552 557 //double ox = ( 0.000891 + dt * 0.000081) / arcSec;558 //double oy = ( 0.005390 + dt * 0.000490) / arcSec;559 //double oz = (-0.008712 - dt * 0.000792) / arcSec;560 553 double ox = (_ox + dt * _oxr) / arcSec; 561 554 double oy = (_oy + dt * _oyr) / arcSec; 562 555 double oz = (_oz + dt * _ozr) / arcSec; 563 556 564 //double sc = 1.0 + 0.4e-9 + dt * 0.08e-9;565 557 double sc = 1.0 + _sc * 1e-9 + dt * _scr * 1e-9; 566 558
Note:
See TracChangeset
for help on using the changeset viewer.