VC高斯投影与坐标转换的源代码.rar

#include "stdafx.h" #include #define Pi 3.1415926535897932384626433832795028841971693993751058209749445923078164 const double awgs = 6378137.0; // WGS84 Semi-Major Axis = Equatorial Radius in meters const double bwgs = 6356752.314; // WGS84 Semi-Minor Axis = Polar Radius in meters const double abes = 6377397.155; // Bessel Semi-Major Axis = Equatorial Radius in meters const double bbes = 6356078.962; // Bessel Semi-Minor Axis = Polar Radius in meters const double cbes = 111120.6196; // Bessel latitude to Gauss-Krueger meters const double dx = -585.7; // Translation Parameter 1 const double dy = -87.0; // Translation Parameter 2 const double dz = -409.2; // Translation Parameter 3 const double rotx = 2.540423689E-6; // Rotation Parameter 1 const double roty = 7.514612057E-7; // Rotation Parameter 2 const double rotz = -1.368144208E-5; // Rotation Parameter 3 // const double sc = 1/0.99999122; // Scaling Factor wrong! // Maik Stoeckmann reported this error on Nov 12th 2002. Thank you, Maik! const double sc = 0.99999122; // Scaling Factor