Changeset 4992 in ntrip
- Timestamp:
- Mar 20, 2013, 11:30:12 AM (12 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/BNC/src/combination/bnccomb.cpp
r4978 r4992 678 678 t_corr* corr = it.value(); 679 679 680 double dT = 60.0;681 682 680 ColumnVector xc(4); 683 681 ColumnVector vv(3); … … 686 684 bncPPPclient::applyCorr(_resTime, corr, xc, vv); 687 685 688 // Relativistic Correction689 // -----------------------690 double relCorr = - 2.0 * DotProduct(xc.Rows(1,3),vv) / t_CST::c / t_CST::c;691 xc(4) -= relCorr;692 693 // Code Biases694 // -----------695 double dcbP1C1 = 0.0;696 double dcbP1P2 = 0.0;697 698 686 // Correction Phase Center --> CoM 699 687 // ------------------------------- … … 701 689 if (_antex) { 702 690 double Mjd = _resTime.mjd() + _resTime.daysec()/86400.0; 703 if (_antex->satCoMcorrection(corr->prn, Mjd, xc.Rows(1,3), dx) == success) { 704 xc(1) -= dx(1); 705 xc(2) -= dx(2); 706 xc(3) -= dx(3); 707 } 708 else { 691 if (_antex->satCoMcorrection(corr->prn, Mjd, xc.Rows(1,3), dx) != success) { 692 dx = 0; 709 693 cout << "antenna not found" << endl; 710 694 } 711 695 } 712 696 713 out << 'P' << corr->prn.toAscii().data() 714 << setw(14) << setprecision(6) << xc(1) / 1000.0 715 << setw(14) << setprecision(6) << xc(2) / 1000.0 716 << setw(14) << setprecision(6) << xc(3) / 1000.0 717 << setw(14) << setprecision(6) << xc(4) * 1e6 718 << setw(14) << setprecision(6) << relCorr * 1e6 719 << setw(8) << setprecision(3) << dx(1) 720 << setw(8) << setprecision(3) << dx(2) 721 << setw(8) << setprecision(3) << dx(3) 722 << setw(8) << setprecision(3) << dcbP1C1 723 << setw(8) << setprecision(3) << dcbP1P2 724 << setw(6) << setprecision(1) << dT 725 << setw(14) << setprecision(6) << (xc(1) + vv(1) * dT) / 1000.0 726 << setw(14) << setprecision(6) << (xc(2) + vv(2) * dT) / 1000.0 727 << setw(14) << setprecision(6) << (xc(3) + vv(3) * dT) / 1000.0 << endl; 728 697 out << corr->prn.toAscii().data() 698 << " APC 3 " 699 << setw(15) << setprecision(4) << xc(1) 700 << setw(15) << setprecision(4) << xc(2) 701 << setw(15) << setprecision(4) << xc(3) 702 << " Clk 1 " 703 << setw(15) << setprecision(4) << xc(4) * t_CST::c 704 << " Vel 3 " 705 << setw(15) << setprecision(4) << vv(1) 706 << setw(15) << setprecision(4) << vv(2) 707 << setw(15) << setprecision(4) << vv(3) 708 << " CoM 3 " 709 << setw(15) << setprecision(4) << xc(1) - dx(1) 710 << setw(15) << setprecision(4) << xc(2) - dx(2) 711 << setw(15) << setprecision(4) << xc(3) - dx(2) 712 << endl; 713 729 714 QString line; 730 715 int messageType = COTYPE_GPSCOMBINED;
Note:
See TracChangeset
for help on using the changeset viewer.