Changeset 7260 in ntrip
- Timestamp:
- Aug 28, 2015, 1:21:05 PM (9 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/BNC/src/bncutils.cpp
r7251 r7260 414 414 415 415 const double bell = t_CST::aell*(1.0-1.0/t_CST::fInv) ; 416 const double e2 416 const double e2 = (t_CST::aell*t_CST::aell-bell*bell)/(t_CST::aell*t_CST::aell) ; 417 417 double Ell[3]; 418 418 if (xyz2ell(XYZ, Ell) != success) { … … 420 420 } 421 421 double rho = sqrt(XYZ[0]*XYZ[0]+XYZ[1]*XYZ[1]+XYZ[2]*XYZ[2]); 422 double Rn = t_CST::aell/sqrt(1- (e2)*pow(sin(Ell[0]),2));422 double Rn = t_CST::aell/sqrt(1-e2*pow(sin(Ell[0]),2)); 423 423 424 424 Geoc[0] = atan((1-e2 * Rn/(Rn + Ell[2])) * tan(Ell[0]));
Note:
See TracChangeset
for help on using the changeset viewer.