Changeset 6213 in ntrip for trunk/BNC/src/rinex


Ignore:
Timestamp:
Oct 5, 2014, 5:15:25 PM (10 years ago)
Author:
mervart
Message:
 
File:
1 edited

Legend:

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

    r6137 r6213  
    538538            ColumnVector xc(4);
    539539            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            }
    548547          }
    549548        }
     
    742741    }
    743742    if (eph) {
    744       ++nSatUsed;
    745743      ColumnVector xSat(4);
    746744      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      }
    754754    }
    755755  }
Note: See TracChangeset for help on using the changeset viewer.