Changeset 6109 in ntrip for trunk/BNC/src/rinex/reqcanalyze.cpp


Ignore:
Timestamp:
Sep 8, 2014, 3:49:16 PM (10 years ago)
Author:
mervart
Message:
 
File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/BNC/src/rinex/reqcanalyze.cpp

    r6017 r6109  
    554554         
    555555          if (eph) {
    556             double xSat, ySat, zSat, clkSat;
    557             eph->position(oneObs->_GPSWeek, oneObs->_GPSWeeks,
    558                           xSat, ySat, zSat, clkSat);
     556            ColumnVector xc(4);
     557            ColumnVector vv(3);
     558            eph->getCrd(bncTime(oneObs->_GPSWeek, oneObs->_GPSWeeks), xc, vv, false);
    559559         
    560560            double rho, eleSat, azSat;
    561             topos(xyzSta(1), xyzSta(2), xyzSta(3),
    562                   xSat, ySat, zSat, rho, eleSat, azSat);
     561            topos(xyzSta(1), xyzSta(2), xyzSta(3), xc(1), xc(2), xc(3), rho, eleSat, azSat);
    563562         
    564563            aziDeg = azSat * 180.0/M_PI;
     
    764763    if (eph) {
    765764      ++nSatUsed;
    766       ColumnVector xSat(3);
    767       double clkSat;
    768       eph->position(_currEpo->tt.gpsw(), _currEpo->tt.gpssec(),
    769                     xSat(1), xSat(2), xSat(3), clkSat);
    770       ColumnVector dx = xSat - xyzSta;
     765      ColumnVector xSat(4);
     766      ColumnVector vv(3);
     767      eph->getCrd(_currEpo->tt, xSat, vv, false);
     768      ColumnVector dx = xSat.Rows(1,3) - xyzSta;
    771769      double rho = dx.norm_Frobenius();
    772770      AA(nSatUsed,1) = dx(1) / rho;
Note: See TracChangeset for help on using the changeset viewer.