- Timestamp:
- Nov 6, 2013, 11:23:49 AM (11 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/src/PacpusTools/src/geodesie.cpp
r208 r209 18 18 #endif //_MSC_VER 19 19 20 namespace Geodes ie20 namespace Geodesy 21 21 { 22 22 … … 170 170 } 171 171 172 } // namespace Geodes ie172 } // namespace Geodesy 173 173 174 174 /// //////////////////////////////////////////////////////////////////// … … 177 177 /// //////////////////////////////////////////////////////////////////// 178 178 //ALGO0001 179 double Geodes ie::LatitueIsometrique(double latitude, double e)179 double Geodesy::LatitueIsometrique(double latitude, double e) 180 180 { 181 181 using ::std::log; … … 189 189 /// //////////////////////////////////////////////////////////////////// 190 190 //ALGO0002 191 double Geodes ie::LatitueIsometrique2Lat(double latitude_iso,double e,double epsilon)191 double Geodesy::LatitueIsometrique2Lat(double latitude_iso,double e,double epsilon) 192 192 { 193 193 using ::std::atan; … … 213 213 } 214 214 /// //////////////////////////////////////////////////////////////////// 215 void Geodes ie::Geo2ProjLambert(215 void Geodesy::Geo2ProjLambert( 216 216 double lambda,double phi, 217 217 double n, double c,double e, … … 225 225 /// //////////////////////////////////////////////////////////////////// 226 226 //ALGO0004 227 void Geodes ie::Proj2GeoLambert(227 void Geodesy::Proj2GeoLambert( 228 228 double X,double Y, 229 229 double n, double c,double e, … … 241 241 } 242 242 /// //////////////////////////////////////////////////////////////////// 243 double Geodes ie::ConvMerApp(double longitude) {243 double Geodesy::ConvMerApp(double longitude) { 244 244 double phi0_Lambert93 = Deg2Rad(46.5); 245 245 double lambda0_Lambert93 = Deg2Rad(3.0); … … 249 249 250 250 //////////////////////////////////////////////////////////////////// 251 void Geodes ie::Geographique_2_Lambert93(const Raf98& raf98,double lambda,double phi,double he,Matrice in,double& E,double& N,double& h,Matrice& out)251 void Geodesy::Geographique_2_Lambert93(const Raf98& raf98,double lambda,double phi,double he,Matrice in,double& E,double& N,double& h,Matrice& out) 252 252 { 253 253 Matrice passage; 254 double conv=Geodes ie::ConvMerApp(lambda);254 double conv=Geodesy::ConvMerApp(lambda); 255 255 double c_=cos(conv); 256 256 double s_=sin(conv); … … 273 273 h=he-diff_h; 274 274 275 Geodes ie::Geo2ProjLambert(275 Geodesy::Geo2ProjLambert( 276 276 lambda,phi, 277 277 n_Lambert93,c_Lambert93,e_Lambert93, … … 280 280 } 281 281 //////////////////////////////////////////////////////////////////////// 282 void Geodes ie::Geographique_2_Lambert93(const Raf98& raf98,double lambda,double phi,double he,double& E,double& N,double& h) {283 Geodes ie::Geo2ProjLambert(282 void Geodesy::Geographique_2_Lambert93(const Raf98& raf98,double lambda,double phi,double he,double& E,double& N,double& h) { 283 Geodesy::Geo2ProjLambert( 284 284 lambda,phi, 285 285 n_Lambert93,c_Lambert93,e_Lambert93, … … 293 293 294 294 /// Converts Lambert93 coordinates (East, North, Height) into geographical coordinates in radians (Longitude = Rad2Deg(lambda), Latitude = Rad2Deg(phi), Height) 295 void Geodes ie::Lambert93_2_Geographique(const Raf98& raf98,double E,double N,double h,double& lambda,double& phi,double& he) {296 Geodes ie::Proj2GeoLambert(295 void Geodesy::Lambert93_2_Geographique(const Raf98& raf98,double E,double N,double h,double& lambda,double& phi,double& he) { 296 Geodesy::Proj2GeoLambert( 297 297 E,N, 298 298 n_Lambert93,c_Lambert93,e_Lambert93, … … 307 307 308 308 //////////////////////////////////////////////////////////////////////// 309 void Geodes ie::Lambert93_2_Geographique(const Raf98& raf98,double E,double N,double h,Matrice in,double& lambda,double& phi,double& he,Matrice& out) {310 Geodes ie::Proj2GeoLambert(309 void Geodesy::Lambert93_2_Geographique(const Raf98& raf98,double E,double N,double h,Matrice in,double& lambda,double& phi,double& he,Matrice& out) { 310 Geodesy::Proj2GeoLambert( 311 311 E,N, 312 312 n_Lambert93,c_Lambert93,e_Lambert93, … … 316 316 317 317 Matrice passage; 318 double conv=Geodes ie::ConvMerApp(lambda);318 double conv=Geodesy::ConvMerApp(lambda); 319 319 double c_=cos(conv); 320 320 double s_=sin(conv); … … 340 340 341 341 //////////////////////////////////////////////////////////////////////// 342 void Geodes ie::Geographique_2_ECEF(double longitude,double latitude,double he,double& x,double& y,double& z) {342 void Geodesy::Geographique_2_ECEF(double longitude,double latitude,double he,double& x,double& y,double& z) { 343 343 const double n = GRS_a / sqrt(1.0 - pow(GRS_e,2) * pow(sin(latitude),2)); 344 344 x = (n + he) * cos(latitude) * cos(longitude); … … 348 348 349 349 //////////////////////////////////////////////////////////////////////// 350 void Geodes ie::ECEF_2_ENU(double x,double y,double z,double& e,double& n,double& u,double lon0,double lat0,double he0) {350 void Geodesy::ECEF_2_ENU(double x,double y,double z,double& e,double& n,double& u,double lon0,double lat0,double he0) { 351 351 double slat = std::sin(lat0); 352 352 double clat = std::cos(lat0); … … 354 354 double clon = std::cos(lon0); 355 355 356 Geodes ie::Matrice C;356 Geodesy::Matrice C; 357 357 C.c0_l0 = -slon; 358 358 C.c1_l0 = clon; … … 376 376 } 377 377 378 QMatrix4x4 Geodes ie::yprenuToMatrix(QVector3D angle, QVector3D position)378 QMatrix4x4 Geodesy::yprenuToMatrix(QVector3D angle, QVector3D position) 379 379 { 380 380 float c1 = cos(angle.x());
Note:
See TracChangeset
for help on using the changeset viewer.