Changeset 6213 in ntrip for trunk/BNC/src/rinex
- Timestamp:
- Oct 5, 2014, 5:15:25 PM (10 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/BNC/src/rinex/reqcanalyze.cpp
r6137 r6213 538 538 ColumnVector xc(4); 539 539 ColumnVector vv(3); 540 eph->getCrd(bncTime(oneObs->_GPSWeek, oneObs->_GPSWeeks), xc, vv, false); 541 542 double rho, eleSat, azSat; 543 topos(xyzSta(1), xyzSta(2), xyzSta(3), xc(1), xc(2), xc(3), rho, eleSat, azSat); 544 545 aziDeg = azSat * 180.0/M_PI; 546 zenDeg = 90.0 - eleSat * 180.0/M_PI; 547 zenFlag = true; 540 if (eph->getCrd(bncTime(oneObs->_GPSWeek, oneObs->_GPSWeeks), xc, vv, false) == success) { 541 double rho, eleSat, azSat; 542 topos(xyzSta(1), xyzSta(2), xyzSta(3), xc(1), xc(2), xc(3), rho, eleSat, azSat); 543 aziDeg = azSat * 180.0/M_PI; 544 zenDeg = 90.0 - eleSat * 180.0/M_PI; 545 zenFlag = true; 546 } 548 547 } 549 548 } … … 742 741 } 743 742 if (eph) { 744 ++nSatUsed;745 743 ColumnVector xSat(4); 746 744 ColumnVector vv(3); 747 eph->getCrd(_currEpo->tt, xSat, vv, false); 748 ColumnVector dx = xSat.Rows(1,3) - xyzSta; 749 double rho = dx.norm_Frobenius(); 750 AA(nSatUsed,1) = dx(1) / rho; 751 AA(nSatUsed,2) = dx(2) / rho; 752 AA(nSatUsed,3) = dx(3) / rho; 753 AA(nSatUsed,4) = 1.0; 745 if (eph->getCrd(_currEpo->tt, xSat, vv, false) == success) { 746 ++nSatUsed; 747 ColumnVector dx = xSat.Rows(1,3) - xyzSta; 748 double rho = dx.norm_Frobenius(); 749 AA(nSatUsed,1) = dx(1) / rho; 750 AA(nSatUsed,2) = dx(2) / rho; 751 AA(nSatUsed,3) = dx(3) / rho; 752 AA(nSatUsed,4) = 1.0; 753 } 754 754 } 755 755 }
Note:
See TracChangeset
for help on using the changeset viewer.